[R] basic question on getGroups for lme analyses

Jose A. Andres
Thu, 22 May 2003 12:02:44 -0000

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