R help - Sep 2008 - nested ANOVA with fixed and random variables missing data

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 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]]

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.
>