I downloaded the package and got it to work with the coding:
model <- function(df) {aov(values ~ WellID, data = twelldata)}
ANOVA1way <- dlply(twelldata, .(Analyte), model)
print(ANOVA1way)
This gives me degrees of freedom and sum of squares for each anova per
analyte. However, I cant get the summary(ANOVA1way) to work so that I can
get p-values, etc... How can I do this?
Jodi
Tobias Verbeke-2 wrote:
>
> baxterj wrote:
>> I have a simple 1 way anova coded like
>>
>> summary(ANOVA1way <- aov(Value ~ WellID, data = welldata))
>>
>> How can I use the BY function to do this ANOVA for each group using
>> another
>> variable in the dataset?? I tried coding it like this, but it doesn't
>> seem
>> to work.
>>
>> summary(ANOVA1way <- by(welldata, Analyte, function(x) aov(Value ~
>> WellID,
>> data = welldata)))
>>
>> In SAS I would code it like this:
>> Proc sort data=welldata; by analyte; run;
>> Proc glm data=welldata;
>> by analyte;
>> class wellid;
>> model value = wellid;
>> run;
>
> Look at the plyr package for a general solution to
> this type of problems:
>
> http://cran.r-project.org/web/packages/plyr/index.html
>
> and its introductory guide on the package home page:
>
> http://had.co.nz/plyr/
>
> HTH,
> Tobias
>
> ______________________________________________
> 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
> and provide commented, minimal, self-contained, reproducible code.
>
>
--
View this message in context:
http://www.nabble.com/Running-an-ANOVA-with-a-BY-tp25645413p25650273.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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
and provide commented, minimal, self-contained, reproducible code.
