[R] Trouble phrasing an R command that will run the model i need (ANOVA, nes

David Elliott dalelemu at hotmail.com
Sat Mar 25 03:41:46 CET 2006


Hi John,

Thank-you very much for replying, - in answer to your questions:

*"As I understand what you wrote, all factors but  linelabel are crossed, 
with linelabel nested within all the other  factors?"

That's exactly it,

*" What is in the highest order cells?  Single observations,  or multiple 
scores?"

I would say "single observations" - but i have to confess I'm not entirely 
sure what you mean by "multiple scores".

Any one data point refers to a single animal (an individual fly, to be 
exact) giving a single value for a response variable, and having a single 
state with respect to each factor. One of my datasets involves crosses 
rather than flies, but it's the same set-up for them: each individual cross 
only gives one value for a response variable, and has one state with respect 
to each factor.

Just in the interests of clarifying things as much as possible: one of my 
datasets looks like this(apologies if the formatting gets weird):

foresttype	region	crosstype	experiment	linelabel	yproductivity
l               g               wx            1              al1            
63
l               g               wx            1              al2            
26
... and so on.

Each horizontal line refers to a single observation. The first horizontal 
line, for example, refers to a single cross that was found to have a 
yproductivity value of "63", and that cross happened to be in the "l" 
foresttype, the "g" region, it was crosstype "wx", in experiment "1", and of 
linelabel "al1". And any one cross only appears once in the dataset, so each 
horizontal line refers to a distinct cross, about which a single observation 
has been made.

Is that the kind of information you were looking for?

- David

>From: John Vokey <vokey at uleth.ca>
>To: dalelemu at hotmail.com
>Subject: RE: [R] Trouble phrasing an R command that will run the model i 
>	need	(ANOVA, nested)
>Date: Fri, 24 Mar 2006 09:41:49 -0700
>
>David,
>   I could probably help you with aov(), but I need more information  about 
>your design.  As I understand what you wrote, all factors but  linelabel 
>are crossed, with linelabel nested within all the other  factors?  What is 
>in the highest order cells?  Single observations,  or multiple scores?
>
>   BTW, Anova() has no methods for the output from aov(), you need to  use 
>summary(aov()); hence, the type III SS of Anova() aren't directly  
>available for aov().  But, you wouldn't want them anyway, because  they 
>aren't the type III SS you probably were expecting (e.g., such  as those 
>found in SAS or SuperAnova).
>




More information about the R-help mailing list