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


