[R] lme in R (WinXP) vs. Splus (HP UNIX)

[Peter Dalgaard]

Greg Tarpinian <sasprog474474 at yahoo.com> 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
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 35327907