[R] nested ANOVA with fixed and random variables missing data

[Adam D. I. Kramer]

On Mon, 8 Sep 2008, Lynise.Wearne at csiro.au wrote:

> I have a nested ANOVA, with a fixed factor "tmt" nested within "site"
> (random). There are missing values in the data set.
>
> aeucs, tmt and site have been defined as objects
>
> I have tried:
>
> model1=lme(aeucs~tmt,random=~1|tmt/site)

I think you want lme(aeucs~tmt,random=~tmt|site)

...as tmt is a fixed factor, you don't want to use it for grouping. Instead,
you want to estimate a separate tmt effect for each site...so estimate (~) a
fixed tmt effect (tmt) within (|) sites (site) = ~tmt|site.

--Adam

>
> I get the following error message
>
> Error in na.fail.default(list(aeucs = c(0.833333333, 1, 1, 1, 1, 1, 0.666666667,  :
>  missing values in object
>
> I need to know where I am going wrong.
>
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.
>