[R] (no subject)

http://www.sangermanomobili.it/mfoto.php [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] how to view the top 20 lines in a long dataset

is this what you want? head(yourdataset,20)   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: Louis Plough lplo...@usc.edu To: r-help@r-project.org Sent: Mon, November 1, 2010 12:40:19 PM

[R] Systat NRM function

Hi: I have some old scripts from when I used to use Systat and have the NRM(value1,value2) that I need to use with R. Does anyone know the R equivalency to this function? Thanks   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California,

[R] drop levels problem

Hi all: I am having trouble dropping levels, got a few hints online without success. Please consider the dataset below:  I was under the inpression that subset(..drop=TRUE) would work but it doesn't library(ggplot2)     library(hmisc) x - structure(list(first = c(38.2086, 43.1768, 43.146,

Re: [R] Filling in missing time samples with na.approx

strange,,I don't see any change either, could it be that we have an older version of zoo?   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: Gabor Grothendieck ggrothendi...@gmail.com To:

Re: [R] drop levels problem

Thanks Joshua, I get it now, levels sometimes drive me loco   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: Joshua Wiley jwiley.ps...@gmail.com To: Felipe Carrillo mazatlanmex

[R] subset

 Hi: I always use subset the same way but now is returning 0 rows. What's wrong with the way I am subsetting? library(ggplot2) structure(list(first = c(38.2086, 43.1768, 43.146, 41.8044, 42.4232, 46.3646, 38.0813, 40.0745, 40.4889, 38.6246, 40.2826, 41.6056, 34.5353, 40.0768), second =

Re: [R] subset

...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Mon, November 29, 2010 1:11:23 PM Subject: Re: [R] subset Hi Felipe, On Mon, Nov 29, 2010 at 1:02 PM, Felipe Carrillo mazatlanmex...@yahoo.com wrote:  Hi: I always use subset the same way but now

[R] How to select the column header with \Sexpr{}

Hi: Since I work with a few different fish runs my column headers change everytime I start a new Year. I have been using \Sexpr{} for my row and columns and now I am trying to use with my report column headers. \Sexpr{1,1} is row 1 column 1, what can I use for headers? I tried \Sexpr{0,1} but

Re: [R] How to select the column header with \Sexpr{}

column headers with \Sexpr{} because I can't find the way to reference run1,run2,run3 and run4. Sorry if I am not explain myself really well.   - Original Message From: Duncan Murdoch murdoch.dun...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch

Re: [R] How to select the column header with \Sexpr{}

Service California, USA - Original Message From: David Winsemius dwinsem...@comcast.net To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: Duncan Murdoch murdoch.dun...@gmail.com; r-h...@stat.math.ethz.ch Sent: Mon, July 12, 2010 3:14:49 PM Subject: Re: [R] How to select the column header

Re: [R] How to select the column header with \Sexpr{}

anymore and returns a NULL when trying to extract the column names. I was not aware that \Sexpr{} only work on dataframes, thanks for your help. - Original Message From: Ista Zahn iz...@psych.rochester.edu To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: David Winsemius dwinsem

[R] Wrap column headers caption

Hi: Using this dataframe with quite long column headers, how can I wrap the text so that the columns are narrower. I was trying to use strwrap without success. Thanks reportDF - structure(list(IDDate = c(3/12/2010, 3/13/2010, 3/14/2010, 3/15/2010), FirstRunoftheYear = c(33 (119 ? 119), n (0 ?

[R] latex table question

Hi: My head is spinning with this latex doc so hopefully after I align my tables to the left of the page my headache are going to be over. I always use: \hspace*{-0.1in} to move my figures horizontally to the left margin of the page but the table below doesn't move at all, but instead it gets

Re: [R] export tables to excel files on multiple sheets with titles for each table

Check the WriteXLS package, I think it does that and also saves each R object on a different excel sheet.   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: eugen pircalabelu

Re: [R] latex table question

To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Sun, July 18, 2010 10:24:39 AM Subject: Re: [R] latex table question Why do you think R-help is a mailing list about LaTeX? Best, Uwe Ligges On 13.07.2010 23:05, Felipe Carrillo wrote: Hi: My head

Re: [R] Figures in Latex

Hannah: I am not sure if this is what you need but you can use an array to do that. Copy and paste the below code to your latex code. \newpage \begin{landscape} \begin{figure}[h] \begin{center}$ \begin{array}{cc} \includegraphics[width=2in]{yourgraphicname}

[R] How to generate a sequence of dates without hardcoding the year

Hi: I have a dataframe named 'spring' and I am trying to add a new variable named 'IdDate' This line of code works fine: spring$idDate - seq(as.Date(2008-07-01),as.Date(2009-06-30),by=week) But I don't want to hardcode the year because it will be used again the following year Is it possible to

Re: [R] Trouble using grid.layout in Sweave

I just run the code below with sweave and works fine It looks like you might be missing the sequence of vplay fig.R,echo=F,fig=T= library(ggplot2) vplay- function(x, y) viewport(layout.pos.row=x, layout.pos.col=y) grid.newpage() p - ggplot(diamonds, aes(x=carat, y=..density..)) +

Re: [R] How to generate a sequence of dates without hardcoding the year

holtman jholt...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Sat, July 24, 2010 4:02:57 PM Subject: Re: [R] How to generate a sequence of dates without hardcoding the year Is this what you want if you want to assume that the date without a year

[R] reordering levels error

Please consider the following dataset: I want to reorder the levels by year but get the following error: Error in tapply(v, x, FUN, ...) : arguments must have same length   I suspect that I need to add the levels before I melt the dataset  but either way I have only use 'reorder' once before and

Re: [R] reordering levels error

-help Help r-h...@stat.math.ethz.ch Sent: Mon, August 30, 2010 2:49:15 PM Subject: Re: [R] reordering levels error On Aug 30, 2010, at 5:25 PM, Felipe Carrillo wrote: Please consider the following dataset: I want to reorder the levels by year but get the following error: Error

Re: [R] reordering levels error

 Thanks Dennis: I always seem to have a hard time defining levels. That's exactly what I needed.    From: Dennis Murphy djmu...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Mon, August 30, 2010 5:41:02 PM Subject: Re: [R] reordering levels error

[R] Controling R from MS Access

HI: I've seen a few threads about this topic but still can't find a straightforward way on this.   Is there a package that can control R within an access form. For example, I want to send a query to R, perform some statistics in R and send the output or summary back to Access and display it on a

[R] gridExtra question

  Hi: I get a couple of warnings  when trying to download gridExtra: install.packages(gridExtra,repos=http://R-Forge.R-project.org)  Warning: unable to access index for repository http://R-Forge.R-project.org/bin/windows/contrib/2.10 Warning message: In getDependencies(pkgs, dependencies,

Re: [R] gridExtra question

Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Fri, October 1, 2010 12:22:22 PM Subject: Re: [R] gridExtra question Is there some reason you don't want the CRAN version? -Ista On Fri, Oct 1, 2010 at 3:08 PM, Felipe Carrillo mazatlanmex...@yahoo.com wrote:   Hi: I getÂ

Re: [R] gridExtra question

desktop and then installed packages from local zip files   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: Erik Iverson er...@ccbr.umn.edu To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h

Re: [R] Count values in a dataframe with respect to groups

Here is an option: grp = c(1,1,1, 1,2, 2,2)  val = c(2,1,5,NA,3,NA,1)  dta = data.frame(grp=grp, val=val) ddply(dta,grp,summarise,count=length(na.omit(val)))   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original

Re: [R] Apply R in Excel through RExcel

ryusuke: It sounds like you need to have (D)COM server to be able to work on the background. As for the foreground, rcom is what you need and it appears to be working OK. I am out of my office right now but I'll back to work next week and I will be able to explain in more detail how the

Re: [R] Placing text in a ggplot

John: The 'year' dataset has 366 rows, I used the 'line.count' without the '1' to come up with 12 rows to match the mlabs and then used line.count to draw the labels. Is this close to what you want? line.count - c(cumsum(as.vector((table(year$monthnum);line.count namposts -

Re: [R] change order of bar plot categories

Is this what you want? temp-c(rep(Low,2),rep(Medium,2),rep(High,2)) light-rep(c(Dark,light),3) avg-dat.avg2[,3] # se-dat.avg2[,4] dat.avg.temp-data.frame(cbind(avg,se)) dat.avg.temp-data.frame(cbind(temp,light,dat.avg.temp)) dat.plot-qplot(light,avg, fill=factor(temp),data=dat.avg.temp, geom=bar,

Re: [R] two graphs 1 x-axis

You can use ggplot2. library(ggplot2) a-c(1,2,3,4,5,6) b-c(3,5,4,6,1,1) c-c(1,1,1,1,1,1) dframe = data.frame(a,b,c);dframe melt.dframe - melt(dframe, id= a);melt.dframe qplot(a,value,data=melt.dframe) + facet_grid(variable~.,scales=free) Felipe D. Carrillo Supervisory Fishery Biologist

Re: [R] RODBC sqlSave does not append the records to a DB2 table

Elaine: Try this. It works for me using my own database. Perhaps you don't have admin rights. Good luck library(RODBC) myDB - odbcConnectAccess(DB2.mdb,uid=admin,pwd=) sqlSave(myDB,se2,rownames=FALSE,append=TRUE) close(myDB) Felipe D. Carrillo Supervisory Fishery Biologist

[R] ggplot2-caption on graphics

Hi: I am trying to add some caption at the bottom of the graphic by using the 'xlab' and just adding a new line to it. Is there another way to do this or this is the correct way? Thanks library(ggplot2) library(effects) df - data.frame( x = c(3, 1, 5), y = c(2, 4, 6), label = c(a,b,c) )

Re: [R] data.frame and ddply

You can do something like this after the output from opfut opfut - data.frame(opfut$CONTRAT,opfut$POSITION,opfut$SETTLEMENT) names(opfut) - c('CONTRAT','POSITION','SETTLEMENT') opfut   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California,

[R] how to lmport this dataframe into R

Hi: I need help with a dataframe(see pic attached). is a mix of dates and text. I want to create a table either using latex function from hmisc or xtable. I already know how to do this but the problem is getting the dataframe into R. I don't have a reproducible example but I am hoping that the pic

[R] how to import a dataframe with mixed text and numbers

Hi: I need help with a dataframe. is a mix of dates and text. I want to create a table either using latex function from hmisc or xtable. I already know how to do this but the problem is getting the dataframe into R. I don't have a reproducible example . If someone is interested in helping with

[R] replace question mark with a dash

Hi: I have the following dataset in R( thanks Gabor for your help) but now the problem is that all the dashes are converted to questions marks (?) I am trying to get the dashes back using: mydf - data.frame(lapply(mydf,function(x) replace(mydf,?(x),-))) but isn't working. I also tried: mydf -

[R] results=hide?

Hi: I am using Sweave and texi2dvi to generate a LaTeX document but can't find the way to hide the graphics while the R chunks are being executed. I thought results=hide would do it but that't not the case. If I do: \begin{figure}[h] figA=true,echo=F,fig=T,results=hide= a rnorm(1000) plot(a) @

[R] Sweave question

Hi: I am using Sweave and texi2dvi to generate a LaTeX document but can't find the way to hide the graphics while the R chunks are being executed. I thought results=hide would do it but that't not the case. If I do: \begin{figure}[h] figA=true,echo=F,fig=T,results=hide= a rnorm(1000) plot(a) @

Re: [R] Sweave question

Message From: Duncan Murdoch murdoch.dun...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Thu, April 29, 2010 4:12:58 AM Subject: Re: [R] Sweave question On 28/04/2010 11:31 PM, Felipe Carrillo wrote: Hi: I am using Sweave and texi2dvi

Re: [R] Question re: interpolation

Are you looking for something like this? data - data.frame(first= c(1,2,3,4,5,6,7,8), zehn = c(15,NA,NA,NA,NA,18,NA,25), second = c(4,NA,7,9,NA,10.2,NA,12),     third= c(6,7,NA,NA,10,12,NA,16))   data     library(zoo)     data2 -na.approx(data,na.rm=F)     data2 ?na.approx   Felipe D.

[R] cumsum output

Hi: Thanks to Dennis and Fernando for your help reordering the levels. Now I have a different issue: I am trying to get the cumulative weekly values using cumsum and it appears to output the wrong values. Here's my dataset: winter -  structure(list(week = c(26L, 27L, 28L, 29L, 30L, 31L, 32L,

Re: [R] cumsum output

I realized that after I hit the send button...I was looking at the wrong dataset,Brr   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA From: jim holtman jholt...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com

[R] predicting without a model

Hello: I have 5 years of weekly passage data and want to predict fish passage for the following year. I don't have a model to use to predict data for the sixth year. Can I somehow still predict based on these five years? I just want to see on the graph what the predicted year would look like and

Re: [R] predicting without a model

: Re: [R] predicting without a model Felipe Carrillo mazatlanmexico at href=http://yahoo.com;yahoo.com writes: ## snip   In the absence of any other information, I would say your best bet would just be to take the weekly average across the previous years.  There are lots of ways to do

Re: [R] a question about latex in Hmisc and .dvi file

HI: If you don't mind me asking this question about latex: Based on Duncan's comment: The print method for latex objects tries to run latex and then display the resulting .dvi file.  (Or maybe it runs pdflatex and displays the .pdf file). Where in MikTex can I change to pdflatex? everytime I use

Re: [R] a question about latex in Hmisc and .dvi file

the latex function?   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA From: RICHARD M. HEIBERGER r...@temple.edu To: Erik Iverson er...@ccbr.umn.edu Cc: Felipe Carrillo mazatlanmex...@yahoo.com; r-help@r-project.org; Charles Dupont

Re: [R] a question about latex in Hmisc and .dvi file

...@temple.edu To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: Erik Iverson er...@ccbr.umn.edu; r-help@r-project.org; Charles Dupont charles.dup...@vanderbilt.edu; Frank E Harrell Jr f.harr...@vanderbilt.edu Sent: Thu, May 13, 2010 1:00:43 PM Subject: Re: [R] a question about latex in Hmisc and .dvi file

Re: [R] Dataframe to word, using R2wd

Hi: Just adding on to this question...Is it a way to add footnotes to tables with R2wd? I started using this nice package about a month ago and is one of the things that I can't figure out yet..   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife

[R] Dataframe to word, using R2wd

  Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: Tal Galili tal.gal...@gmail.com To: Jeremy Miles jeremy.mi...@gmail.com Cc: r-help@r-project.org Sent: Sat, May 15, 2010 1:03:12 AM

[R] Dataframe to word, using R2wd

Hi: Just adding on to this question...Is it a way to add footnotes to tables with R2wd? I started using this nice package about a month ago and is one of the things that I can't figure out yet..   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife

Re: [R] Problem with Sweave not recognising \Sexpr{}

 I had the same problem before and I think it has something to do with the R2HTML package. To take care of that problem simply sweave your rnw file like this: Sweave(yourfile.Rnw,syntax=SweaveSyntaxNoweb) and try Sexpr{} again.   Felipe D. Carrillo Supervisory Fishery Biologist Department of the

Re: [R] yeroon.net/ggplot2 web application v0.11

Hi Jeroen: This is great for someone who is point and click oriented. Can this program be downloaded to be used offline or its just web based? My field crew usually take toughbooks to the field and since they don't know how to program it will be nice to make the graphics just by pointing and

Re: [R] simple data manipulation question

Or like this: d = data.frame(gender=c(m,f,unkown), x=rnorm(300)) ddply(d,gender,summarise,mean=mean(x)) Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA --- On Thu, 12/10/09, smu m...@z107.de wrote: From: smu

Re: [R] Put two plots side by side

Two different ways: library(ggplot2) x=5 size=50 A=data.frame(X=sample(x, size, replace=T), Y=sample(x, size, replace=T),a=rep(1:2,each=25));A # Facetting qplot(X,Y,data=A) + geom_jitter(position=position_jitter(width=.03)) + facet_grid(.~a) # Or with vp p=qplot(X, Y, data=A) +

Re: [R] Can not save plot to png file correctly

You can save as png like this too: library(ggplot2) data=data.frame(   X=sample(10,1000,replace=T)   , Y=letters[1:10]) png(mypng.png) qplot(X, data=data, geom='histogram') + facet_wrap( ~ Y) dev.off()   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife

[R] Can not save plot to png file correctly

With ggsave the graph windows pops up but using: png(mypng.png) qplot(X, data=data, geom='histogram') + facet_wrap( ~ Y) dev.off() The graph is saved in the background Is there a way to hide the graph window when using ggsave? You can save as png like this too: library(ggplot2) data=data.frame(  

Re: [R] cumsum function with data frame

You can also use ddply from the plyr package: library(plyr) Data - read.table(textConnection(variable    Year  value EC01    2005    5 EC01    2006    10 AAO1    2005  2 AAO1  2006  4),header=T) Data

Re: [R] cumsum function with data frame

,CUMSUM=cumsum(value))   - Original Message From: Felipe Carrillo mazatlanmex...@yahoo.com To: Joris Meys jorism...@gmail.com; n.via...@libero.it n.via...@libero.it Cc: r-help@r-project.org Sent: Thu, June 3, 2010 11:28:58 AM Subject: Re: [R] cumsum function with data frame You can

Re: [R] Patch for legend.position={left,top,bottom} in ggplot2

Hi: I downloaded the patch, how do I incorporate it to my current version of ggplot2?   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA - Original Message From: Hadley Wickham had...@rice.edu To: Karsten Loesing

[R] Latex: Date Format conversion

Hi: Can't find a way to convert from shortDate to LongDate format. I got: 3/10/10 that I want to convert to March 10, 2010. I am using: \documentclass[11pt]{article} \usepackage{longtable,verbatim} \usepackage{ctable} \usepackage{datetime} \title{my title} \begin{document}   % Convert date

[R] Date conversion

Hi: Can't find a way to convert from shortDate to LongDate format. I got: 3/10/10 that I want to convert to March 10, 2010. I am using: \documentclass[11pt]{article} \usepackage{longtable,verbatim} \usepackage{ctable} \usepackage{datetime} \title{my title} \begin{document}   % Convert date

Re: [R] Latex: Date Format conversion

this pseudo-code: Report from \longdate\Sexpr{report[1,1]} - \longdate\Sexpr{report[1,15]} Where long date will be the format that converts 6/1/10 to June 01, 2010 Thanks for helping. - Original Message From: Marc Schwartz marc_schwa...@me.com To: Felipe Carrillo mazatlanmex

Re: [R] Latex: Date Format conversion

Service California, USA - Original Message From: Marc Schwartz marc_schwa...@me.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Thu, June 10, 2010 10:19:40 AM Subject: Re: [R] Latex: Date Format conversion Felipe, I would not do the processing

Re: [R] Date conversion

jwiley.ps...@gmail.com To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Thu, June 10, 2010 1:18:27 PM Subject: Re: [R] Date conversion Hello Felipe, Is this what you want? format(as.Date(3/10/10, format=%m/%d/%y), %B %d, %Y) Josh On Thu, Jun 10, 2010 at 8:29

Re: [R] Plotting Data on a Map

, Tom On Mon, Jun 21, 2010 at 11:57 PM, Felipe Carrillo mazatlanmex...@yahoo.com wrote: Hi Tom: I am just starting to use rgdal and maptools but I have a long way to go. I went to a training a couple of weeks ago and the instructor showed us a csv file and a shapefile with wolf data from

Re: [R] Plotting Data on a Map

For some reason the shapefile can't get attached. The shapefile is too large to put it in dput..Is there another way to do this? - Original Message From: Felipe Carrillo mazatlanmex...@yahoo.com To: Tom Hopper tomhop...@gmail.com Cc: r-h...@stat.math.ethz.ch; ggpl

[R] rgdal-maptools

The shapefile data can be downloaded from the link below: download all the six files and save them on your working directory and make sure the dsn path is set to where the files are saved. My shapefiles are saved on C:/Data. https://secure.filesanywhere.com/fs/v.aspx?v=897263875a6472a99baa Hi:

[R] predict newdata question

Hi: I am using a subset of the below dataset to predict PRED_SUIT for the whole dataset but I am having trouble with 'newdata'. The model was created with 153 records and want to predict for 208 records. wolf2 - structure(list(gridcell = c(367L, 444L, 533L, 587L, 598L, 609L, 620L, 629L, 641L,

Re: [R] predict newdata question

-project.org; href=mailto:r-help-boun...@r-project.org;r-help-boun...@r-project.org [mailto: href=mailto:r-help-boun...@r-project.org;r-help-boun...@r-project.org] On Behalf Of Felipe Carrillo Sent: Saturday, 26 June 2010 10:35 AM To: ymailto=mailto:r-h...@stat.math.ethz.ch; href=mailto:r-h

Re: [R] predict newdata question

: href=mailto: ymailto=mailto:r-help-boun...@r-project.org; href=mailto:r-help-boun...@r-project.org;r-help-boun...@r-project.org ymailto=mailto:r-help-boun...@r-project.org; href=mailto:r-help-boun...@r-project.org;r-help-boun...@r-project.org] On Behalf Of Felipe Carrillo Sent: Saturday, 26

[R] R2wd- how to open an existing document

Hi: How can one open an existing word document with wdGet() I am getting an error message when trying to open it like this: wdGet(filename=myDoc.doc,path=c/mydata) Looking at Tal Galili's website example it appears that double backslashes are used but i tried it and didn't work either. Thanks  

Re: [R] R2wd- how to open an existing document

California, USA - Original Message From: Erik Iverson er...@ccbr.umn.edu To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Sent: Tue, July 6, 2010 3:03:31 PM Subject: Re: [R] R2wd- how to open an existing document Having no clue about the specifics here

Re: [R] median of grouped data

Like this? d = read.table(textConnection(ID PM adt_01 25.255 adt_01 225.36 adt_01 14.2325 adt_02 15 adt_02 12.3 adt_03 15.2 adt_03 148.3 adt_03 25.5 adt_03 14.25),head=T) d d$PM - as.numeric(d$PM) ddply(d,.(ID),numcolwise(median)) Felipe D. Carrillo Supervisory Fishery Biologist

Re: [R] how to read this data file into R?

Convert the file into a csv file( myfile.csv), delete some of the unwanted rows before importing into R. Save the file on your working directory and then import into R: x - read.csv(myfile.csv) # Change column names names(x) - LETTERS[1:ncol(x)];x Change your column names to whatever you want

Re: [R] how to read this data file into R?

That's strange because I was able to read the file into R just fine --- On Tue, 2/2/10, David Winsemius dwinsem...@comcast.net wrote: From: David Winsemius dwinsem...@comcast.net Subject: Re: [R] how to read this data file into R? To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-help@r

Re: [R] diagnostic plots

Or just click on the Enter key to see the following graphic Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA --- On Wed, 2/3/10, jim holtman jholt...@gmail.com wrote: From: jim holtman jholt...@gmail.com Subject:

[R] ggplot2 X axis levels

Hi all: I've done this before with factors but can't figure how to do it with a continuous variable. I am trying to reorder the sequence of my weeks along the X axis. I want to start with week 27 to 52 and then 1 to 26. I guess I could use levels along with seq() but doesn't seem to work for me.

Re: [R] ggplot2 X axis levels

...@menne-biomed.de To: r-help@r-project.org Sent: Sat, February 20, 2010 2:56:14 AM Subject: Re: [R] ggplot2 X axis levels Felipe Carrillo wrote: Hi all: I've done this before with factors but can't figure how to do it with a continuous variable. I am trying to reorder the sequence of my

[R] hide or minimize graphics

Hello all: Is there a way to create a graph and then hide it or just make a graph in the background (invisible). Thanks Felipe D. Carrillo Fishery Biologist US Fish Wildlife Service Red Bluff, California 96080 __ R-help@r-project.org mailing

Re: [R] adding the mean and standard deviation to boxplots

Tom: You can do this with ggplot2. The code below puts 95% CI,a smooth line and the mean(blue point)on the same plot. Felipe library(ggplot2) r - ggplot(ToothGrowth, aes(y=len, x=factor(dose))) r$background.fill = cornsilk r + geom_boxplot(aes(colour=supp)) +

[R] How to estimate weekly Variance

Dear All: Using the Excel formula below I am able to estimate daily variance(Daily VarPD)for PD,but I can't figure out how to estimate the Weekly variance for PD. Any hints? Daily VarPD = (PD * (1 - TD))/TD + VarTD * ((PD*(1-TD)+ PD^2*TD))/TD^3 WK TDVARTD PD Daily

[R] Bootstrapping data with a regression model

Hello all: Just wondering if I can get advice on what kind of bootstrapping I should use when using a regression model to estimate juvenile fish passage data. I use rotary screw traps to do fish mark-recapture trials and the efficiency of every trial is added to the graph generating a different

[R] How to hold a value(Mean sq) with a string

Hi all: Can someone advice me on how to hold the residuals Mean sq value on a string so it can be used in other calculations. I was trying something like this: Msquare-dfr$Mean sq but fails..Thanks dfr - read.table(textConnection(percentQ Efficiency 1.5650.0125 1.94 0.0213

[R] How to Estimate Covariance by Week based on a linear regression model

Hi all: I have always used SPSS to estimate weekly covariance based on a linear regression model but have to hard code the model Std. Error and the Mean-Square and then execute one week a the time. I was wondering if someone could give me an idea on how to estimate weekly(WK) covariance

[R] transpose dataset

Hi all: How can I transpose this dataset from column to row idweek value 1 5 51 2 6 73 3 7 41 4 8 22 5 9 83 6 10 55 7 11 42 to something like this... id1 2 3 4 5 6 7 week 5 6 7 8 9 10 11 value 51 73 41

[R] How to loop through all the columns in dataframe

Hi: Can anyone advice me on how to loop and perform a calculation through all the columns. here's my data xd- c(2.2024,2.4216,1.4672,1.4817,1.4957,1.4431,1.5676) pd- c(0.017046,0.018504,0.012157,0.012253,0.012348,0.011997,0.012825) td- c(160524,163565,143973,111956,89677,95269,81558)

Re: [R] How to loop through all the columns in dataframe

--- jim holtman [EMAIL PROTECTED] wrote: Glad I could help. You might want to post it back to R-Help so that others can see what was done. On Sun, Mar 16, 2008 at 6:02 PM, Felipe Carrillo [EMAIL PROTECTED] wrote: Jim: I owe you man, this is great,I never thought that I could

Re: [R] I need help integrating ggplot2 into Excel

Create an R script,name it Myggplot2 and place it on your desktop.In excel open the vb editor and write this code: Sub Myggplot2() 'start R Call rinterface.StartRServer 'Put the dataframe into R,assuming that your 'data is in sheet1 cells A1:D1 Call rinterface.PutDataframe(yourdataframename,

[R] counting values on one colum only

Hi all: Can someone help me count the number of rows with values in colum a only. assume the name of my dataframe is weekly I was trying i- nrows(weekly$a) i but returns 7 when it should be 4. Thanks a b c d 27.000 27.000 1.569 0.013

Re: [R] counting values on one colum only

: I basically need to count the number of rows that have values only. I tried your example and it crashed.I have been reading and looking for a function that does that but no luck... --- [EMAIL PROTECTED] wrote: From: Felipe Carrillo [EMAIL PROTECTED] Date: 2008/03/22 Sat PM 06:16:59

[R] Inverse t-distribution

Hello all: Is there a function in R to estimate the Inverse t-distribution(tif in Systat).If so how can I see an example on how is used? Thanks Felipe D. Carrillo Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA

Re: [R] Inverse t-distribution(TINV function in Excel)

-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Felipe Carrillo Sent: Sunday, March 23, 2008 6:18 PM To: [EMAIL PROTECTED] Subject: [R] Inverse t-distribution Hello all: Is there a function in R to estimate the Inverse t-distribution(tif in Systat).If so

[R] cbind and mean by week

Hi: I have been able to finally crunch my data by importing it by week(thank you all for your help),but here we go again.. Now I'am trying to do it for the the whole year. Since the dataset is huge I'm only making a 3 weeks dataframe. - I want to get the mean of pd by week - I want to count

Re: [R] GGPlot 2 and odfweave

Hi John: Use print before qplot or ggplot. Something like this: TV,fig=TRUE,echo=FALSE= print(qplot(subscriber,data=TV09,geom=bar,fill=subscriber, xlab=NULL,ylab=NULL)+coord_flip()+ opts(legend.position=none)) @ Felipe D. Carrillo  Supervisory Fishery Biologist  Department of the

Re: [R] USGS stream flow data automatic download R

An Excel Web query will also do what you want. From an Excel Workbook create a web query and link it to the USGS website. You can schedule downloads or the query will download data up-to-date everytime the workbook is open. Felipe D. Carrillo Supervisory Fishery Biologist Department of

Re: [R] Scatter Plot

Tri this for both, points and text.. library(ggplot2) qplot(x=1:5,y=1:5) + geom_text(label=1:5) Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA --- On Thu, 7/30/09, amna khan amnakhan...@gmail.com wrote: From:

Re: [R] making scatter plot points fill semi-transparent

You can also use ggplot2: library(ggplot2) x - rnorm(1);y - rnorm(x);myplot - data.frame(x,y) qplot(x,y,data= myplot,colour=I(alpha(blue,1/25))) Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA --- On Wed,

[R] ggplot2-ddply question

Hi all: I am trying to use the ddply function to estimate the mean of 'Total','Fry','Smolt' and 'Fry.Eq' columns without success. I have the dput of my dataset below. I wonder if someone can give me a hand with this function. # dput(winter) winter -structure(list(IDDate = structure(c(37L, 48L,

  1   2   3   >