Dear all,

 

I have an experimental setup in which a random variable is nested within a 
fixed variable; however I have troubles specifying the correct LMM with lme. I 
have searched the lists but haven't been
able to find an example like my setup, which I unfortunately need to get this 
stuff right. Pinheiro & Bates is great but I still can't figure out how to do 
it. 

 

My experimental setup was as follows:

100 measurements per treatment plot

2 treatment plots per site

4 sites: 2 in one area and 2 in another area

 

Both treatment and area are fixed factors,while site is random. I am interested 
in the significance of the fixed effects,less in the magnitude of the random 
effect. 

 

I have tried:

 

mod1 <- lme(response ~ area*treatment, data=data,random= ~1|site)

but now site is not nested in area…

 

mod2 <- lme(response ~ area*treatment, data=data,random= ~1|area/site)

but now area is both a fixed and a random variable, which doesn't seem to make 
sense, plus I run out of df for treatment

 

mod3 <- lme(response ~ area*treatment, data=data,random= ~1| plot)

but here plots are not grouped according to site

 

I hope someone would be willing to help me,
thank you in advance!

 

Robert Buitenwerf

Ecologist

South African Environmental Observation Network
_________________________________________________________________
[[elided Hotmail spam]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to