Hi all!

I am working on a nested lme model with one fixed effect ("treatment", which 3 levels) and two random effects for "Individuals" (four of them) within "treatment" and "replicate -2 levels-" within "individual" within "treatment". For doing so, I´ve been trying to create a factor for Individual%in%Treatment, say IT 
by

Hongos$IT=getGroups(~1|Treatment/Individual, data=hongos1, level=4) 

but I got the following error:

Error in getGroups(~1 | Treatment/Individual, data = hongos1, level = 4) : 
        no applicable method for "getGroups"

I've been trying some other things but nothing seem to work, probably because I cannot fully understand how o use the getGroups command. Any suggestions?

Thanks a lot!

/Jose



	[[alternate HTML version deleted]]


