[R-sig-ME] Peijian Shi

[shi_peijian]

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!
  
Best wishes,
 
Peijian Shi

	[[alternative HTML version deleted]]