Re: [R] How to obtain a mean by more than one other variable

2014-03-09 Thread Andrew Koeser
One function I like is from package psych. library(psych) describe.by(ARK$TotalWater, group=Tray:Pot) Andrew Koeser On 3/9/2014 10:39 PM, Erynn Call wrote: I'm trying to create a new variable (meanABUNDHA) by obtaining a mean value *by* other variables. I surveyed birds (multiple specie

Re: [R] post hoc test to a significant Kruskal-Wallis test package asbio function pairw.kw()

2013-10-21 Thread Andrew Koeser
Luigi, agricolae has a nice function for this with the post hoc test as part of the output. AKK On 10/21/2013 11:12 PM, arun wrote: Hi, Try: pairw.kw(y,factor(x),conf=.95) A.K. On Monday, October 21, 2013 7:10 PM, luigi marongiu wrote: Dear all, I am trying to apply a post hoc test t

[R] Fwd: Re: How to reference to the `stats` package in academical paper

2013-03-05 Thread Andrew Koeser
Julien, I would just try your best given the journal's style guide and wait for them to change it. For what it is worth, my last paper was corrected as follows. (In text)"...one-way analysis of variance (ANOVA) in R [version 2.14.2 (R Core Team, 2012)]. (reference section) R Core Team. 201

Re: [R] Learning the R way – A Wish

2013-03-04 Thread Andrew Koeser
The book that helped me break into R and more advanced texts was Crawley's "Statistics: An Introduction with R." Very light read that assumes no prior knowledge with stats or R. I am using it to teach my fellow grad students R and all agree it was worth scrimping pennies to get. He also has a

Re: [R] Getting the correct factor level as Dunnett control in glht()

2013-02-27 Thread Andrew Koeser
Fantastic! Thanks, Don. I saw contrMat() and base option, but for whatever reason I didn't see how/if they fit together. Your help is greatly appreciated. Cheers, Andrew On 02/27/2013 03:07 PM, MacQueen, Don wrote: I think you can probably specify the base level using the contrMat() functio

[R] Getting the correct factor level as Dunnett control in glht()

2013-02-26 Thread Andrew Koeser
Hello all, I would like to do a Dunnett test in glht(). However, the factor level I want to use as the control is not the first. dunn1<-glht(model3, linfct = mcp(Container = "Dunnett"), alternative = "less") The factor container has 8 levels, so it would be nice not to manually enter in all o