Displaying 1 result from an estimated 1 matches for "216x5".
Did you mean:
2165
2005 Dec 09
1
lmer for 3-way random anova
...this using lmer from lme4.
I have 3 main effects, tim, trt, ctr, and all the interaction effects tim*trt*ctr. The response variable is ge.
Here is my lme code:
dat <- data.frame(ge=ge,trt=factor(trt),tim=factor(tim),ctr=factor(ctr))
dat$grp = as.factor(rep(1, nrow(dat)))
# dim(dat) = 216x5
w <- lme(ge ~ 1,data=dat,
random = list(
grp = pdBlocked(list( pdIdent(~ trt - 1)
,pdIdent(~ tim - 1)
,pdIdent(~ ctr - 1)
,pdIdent(~ trt:tim - 1)
,pdIdent(~ trt:ctr - 1)...