
   Hello:

   I have a doubt about the syntax for random effects, and I hope you can help me with that. 

   I am fitting a mixed effects model with 3 random factors, let's call them A, B and C for simplicity. A and B are nested factors, B is nested into A, and C is not nested, it's just another random effect independent from the other 2. I guess the correct way to specify this in "lmer" would be something like:

   model1<-lmer(y~fixed1*fixed2+(1|A/B)+(1|C))

   but I am not sure how I could do this with "lme". Maybe:

   model1<-lme(y~fixed1*fixed2,random=~1|A/B,random=~1/C) ???
 
   or perhaps:

   model1<-lme(y~fixed1*fixed2,random=~C|A/B)    ???

   Are any of these correct? 

   Thank you very much.

   Iker


  
   


   



      
	[[alternative HTML version deleted]]


