## Not run:
#
# require(lme4)
#
# names <- colnames(corn)
# X <- t(matrix(unlist(strsplit(names, ".", fixed=T)), 3, 36))
# X <- data.frame(X)
# colnames(X) <- c("spec", "temp", "plant")
#
# summary(X)
#
# y <- corn[4662,]
# lmer(y~spec*temp + (1|plant:spec:temp), data=X)
#
# ## End(Not run)
Run the code above in your browser using DataLab