E.g., y <- rnorm(12); a=gl(4,1,12); b=gl(3,4,12). Can I fit an additive model like y~a+b but with a and b random?
Everything I've tried gives an error:
> lme(y~1,random=~1|(a+b))
Error in switch(mode(object), name = , numeric = , call = object, character = as.name(object), :
[[ cannot be of mode (
> lme(y~1,random=~a+b)
Error in getGroups.data.frame(dataMix, groups) :
Invalid formula for groups
Thanks Gordon
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help