Greg Tarpinian <[EMAIL PROTECTED]> writes: > R2.2 for WinXP, Splus 6.2.1 for HP 9000 Series, HP-UX 11.0. > > > I am trying to get a handle on why the same lme( ) code gives > such different answers. My output makes me wonder if the > fact that the UNIX box is 64 bits is the reason. The estimated > random effects are identical, but the fixed effects are very > different. Here is my R code and output, with some columns > and rows deleted for space considerations: > > FOO.lme1 <- lme(fixed = Y ~ X1*X2, > random = ~ X2 | X3, data = FOO, > method = "REML", > control = list(maxIter = 200, msMaxIter = 200, > tolerance = 1e-8, niterEM = 200, > msTol = 1e-8, msVerbose = TRUE, > optimMethod = "Nelder-Mead"))
Try FOO.lme2 <- update(FOO.lme1,.~C(X1,"contr.helmert")*X2) summary(FOO.lme2) and see if enlightenment does not ensue. -- O__ ---- Peter Dalgaard Øster Farimagsgade 5, Entr.B c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
