Hi, I'm a SAS user trying to convert myself to R but I still have problems with some pretty simple commands. First I wanted to add up a number of red-tailed hawks seen per day (julian day) per year. So I tried: RTyrday <- tapply(RThr,list(year,julian),sum) And then I tried the following regression: mod1 <- glm(RTyrday~julian+year, family=gaussian (link=identity),data=RT) Wich didn't work since my vector RTyrday and julian don't have the same length. My question is: How can I create a new data sheet with the output of my function "tapply" ? Something I could have done in SAS by giving an "output out" to my "proc summary" Thank you, Emilie Emilie Berthiaume Graduate Student Biology Department Sherbooke University Sherbrooke, Québec CANADA emilie.berthiaume@USherbrooke.ca [[alternative HTML version deleted]]
Emilie Berthiaume wrote:> Hi, > > I'm a SAS user trying to convert myself to R but I still have problems with somepretty simple commands.> > First I wanted to add up a number of red-tailed hawks seen per day (julian day) per year. So I tried:> > RTyrday <- tapply(RThr,list(year,julian),sum)Try: RT <- as.data.frame.table(RTyrday)> > And then I tried the following regression: > > mod1 <- glm(RTyrday~julian+year, family=gaussian (link=identity),data=RT)here it is simpler with lm() Kjetil> > Wich didn't work since my vector RTyrday and julian don't have the same length. My question is: How can I create a new data sheet with the output of my function "tapply" ? Something I could have done in SAS by giving an "output out" to my "proc summary" > > Thank you, > > Emilie > > > > Emilie Berthiaume > Graduate Student > Biology Department > Sherbooke University > Sherbrooke, Qu?bec > CANADA > > emilie.berthiaume at USherbrooke.ca > [[alternative HTML version deleted]] > > > > ------------------------------------------------------------------------ > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Emilie Berthiaume wrote:> Hi, > > I'm a SAS user trying to convert myself to R but I still have problems with some pretty simple commands. > > First I wanted to add up a number of red-tailed hawks seen per day (julian day) per year. So I tried: > > RTyrday <- tapply(RThr,list(year,julian),sum) > > And then I tried the following regression: > > mod1 <- glm(RTyrday~julian+year, family=gaussian (link=identity),data=RT) > > Wich didn't work since my vector RTyrday and julian don't have the same length. My question is: How can I create a new data sheet with the output of my function "tapply" ? Something I could have done in SAS by giving an "output out" to my "proc summary" > > Thank you, > > EmilieOne of many ways is at http://biostat.mc.vanderbilt.edu/SasByMeansExample (especially for SAS users). The example uses the R Hmisc package. Frank> > > > Emilie Berthiaume > Graduate Student > Biology Department > Sherbooke University > Sherbrooke, Qu?bec > CANADA > > emilie.berthiaume at USherbrooke.ca > [[alternative HTML version deleted]] > > > > ------------------------------------------------------------------------ > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html-- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University
Emilie Berthiaume wrote:> > I found out that for want I want to do I must use the function > "summarize". My only problem now is that I don't understand how to > summarize two vectors at a time. For exemple, I wanted to get the sum > of RThr and effortwt per day per year so I've tried this: > > RTyrday <- with(RT.sub,summarize(cbind(RThr,effortwt), > llist(year,jjulian), FUN=sum)) > > But the output combine RThr and effortwt instead of doing two separate > column with those vectors. > What am I missing? > > Thanks for your help, > > EmilieSee the example I posted 2 hours ago, especially how you need to specify FUN. FH> > ----- Original Message ----- From: "Frank E Harrell Jr" > <f.harrell at vanderbilt.edu> > To: "Emilie Berthiaume" <emilie.berthiaume at usherbrooke.ca> > Cc: <R-help at stat.math.ethz.ch> > Sent: Wednesday, March 01, 2006 12:16 PM > Subject: Re: [R] How to do a "proc summary" in R? > > > Emilie Berthiaume wrote: > >> Hi, >> >> I'm a SAS user trying to convert myself to R but I still have problems >> with some pretty simple commands. >> >> First I wanted to add up a number of red-tailed hawks seen per day >> (julian day) per year. So I tried: >> >> RTyrday <- tapply(RThr,list(year,julian),sum) >> >> And then I tried the following regression: >> >> mod1 <- glm(RTyrday~julian+year, family=gaussian >> (link=identity),data=RT) >> >> Wich didn't work since my vector RTyrday and julian don't have the >> same length. My question is: How can I create a new data sheet with >> the output of my function "tapply" ? Something I could have done in >> SAS by giving an "output out" to my "proc summary" >> >> Thank you, >> >> Emilie > > > One of many ways is at > http://biostat.mc.vanderbilt.edu/SasByMeansExample (especially for SAS > users). The example uses the R Hmisc package. > > Frank > >> >> >> >> Emilie Berthiaume >> Graduate Student >> Biology Department >> Sherbooke University >> Sherbrooke, Qu?bec >> CANADA >> >> emilie.berthiaume at USherbrooke.ca >> [[alternative HTML version deleted]] >> >> >> >> ------------------------------------------------------------------------ >> >> ______________________________________________ >> R-help at stat.math.ethz.ch mailing list >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide! >> http://www.R-project.org/posting-guide.html > > >-- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University