[R] split plot design with missing plots

[Martin Henry H. Stevens]

Dear Peter,
I did double check whether all factors had "factor" as their class. "plot" 
is a factor, with 48 levels. I used Professor Ripley's suggestion to go 
with lme() which can handle unbalanced designs. It had no problem (that I 
was aware of!). Its format and methods, however, are less familiar to me, 
and so I was resisting it. Now, as a consequence, however, I am relearning 
(again!) the relations among variance components, Mean squares, etc. I am 
still not clear, though, how lme() tests fixed effects. I'll get it yet.
Thank you for your feed back.

At 04:45 PM 6/19/2002, Peter Dalgaard BSA wrote:
>"Martin Henry H. Stevens" <HStevens at muohio.edu> writes:
>
> > The output, below, tests All factors over Error: Within. I don't
> > understand. I expected Fert to be tested over plot residuals, Seed and
> > Seed:Fert tested over plot:Seed residuals, and all else tested over
> > Within residuals.
> >
> > Can any body help?
> >
> > summary(fmc)
> >
> > Error: plot
> >       Df  Sum Sq Mean Sq
> > Fert  1 0.98913 0.98913
> >
> > Error: plot:Seed
> >       Df  Sum Sq Mean Sq
> > Fert  1 0.02016 0.02016
>
>1DF total for these two strata makes me suspect that you put in "plot"
>as a numerical variable rather than as a factor.
>
>--
>    O__  ---- Peter Dalgaard             Blegdamsvej 3
>   c/ /'_ --- Dept. of Biostatistics     2200 Cph. N
>  (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
>~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907

Martin Henry H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056

Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.muohio.edu/~botcwis/bot/henry.html

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._