Dear R-help group,
How can I define a lme with 3 factors(a,b,c), where c is nested in b,
and a is crossed with b/c?
I think that:
lme(response ~ ..., data = Data,
random = pdBlocked(list(pdIdent(~ a - 1), pdIdent(~ b - 1))))
is one part of the answer and:
lme(response~..., data=Data, random=~1|b/c)
is the other part of the answer but how can I combine them??
Could anybody please help me?
Many greetings,
Thomas
PS: Is there good literatur available which could help working with
mixed random effects (I just know Pinheiro and Bates 2000), perhaps
even on the internet?
[[alternative HTML version deleted]]