[R-sig-ME] Peijian Shi
shi_peijian at 163.com
Mon Mar 23 04:37:03 CET 2015
Dear Sir / Ms.,
I am using lme4 to analyze the effects of three treatment levels, three seed types, six sites on the removal percentage of seeds over 20 days. Each day, six sites were randomly chose, so the six sites' locations are different from another day's. There are total 10 seed for each combination (Treatment * Type) for oberving the removal number.
I think, Day should be a random effect and Site also a random effect. Then I use the following code:
glmer(cbind(Removal, 10-Removal) ~ Type + Treatment + (1 | Day / Site), family=binomial)
My quesiton is whether (1|Day/Site) correct for my experiment? Each day, there were six random sites, and we observed the removal number per 10 seeds for the combinations (Treatment * Type).
Thank you very much!
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models