Displaying 1 result from an estimated 1 matches for "mssite".
Did you mean:
mysite
2012 Feb 14
2
how to test the random factor effect in lme
...18
13 B 5 19
14 B 5 20
15 B 5 21
16 B 6 19
17 B 6 23
18 B 6 21
By hand calculation, the result should be:
Source SS DF MS
Areas 18.00 1 18.00
Sites 14.44 4 3.61
Error 20.67 12 1.72
Total 53.11 17 ---
MSareas/MSsites = 4.99 --- matching the R output
MSsites/MSE = 2.10
Conclusion is that Neither of Areas nor Sites make differences.
My R code so far can only test the fixed effect :
anova(lme(PCB~Area,random=~1|Sites, data = PCBdata))
numDF denDF F-value p-value
(Intercept) 1 12 1841.7845...