Honorable newslisters:
S-Plus 6.0 for Windows.
On p. 165 of Bates and Pinheiro's Mixed Effect Modeling (2000) book it
explains how to specify a model to lme() in which (some of) random effects
are crossed. The random=specification uses a class pdBlocked matrix and
works as given only if the data are a groupedData object. This worked fine
for me, but I wonder how one would specify the random object if one uses a
plain data frame that is not of class groupedData. Maybe it's obvious, but a
couple of things I tried didn't work, so I obviously don't get it. Any help
would be appreciated.
Bert Gunter
Biometrics Research RY 70-38
Merck & Company
P.O. Box 2000
Rahway, NJ 07065-0900
Phone: (732) 594-7765
mailto: bert_gunter@merck.com
"The business of the statistician is to catalyze the scientific learning
process." -- George E.P. Box
|