[R] Rounding of lme coefficients: Mac vs Windows
David Afshartous
dafshartous at med.miami.edu
Sat Oct 6 20:14:45 CEST 2007
All,
I have an lme model estimated in R 2.5.1 on my Mac; when I estimate the same
model on Windows, the parameter coefficients are rounded to integers. Below
is a similar example for the Orthodont data. Is there some option I need to
set in the Windows version to prevent rounding? Didn't see this in the
archives.
Thanks,
David Afshartous
> sessionInfo()
R version 2.5.1 (2007-06-27)
i386-apple-darwin8.9.1
locale:
en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "datasets" "methods"
"base"
other attached packages:
nlme
"3.1-84"
> fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age
>
> fm1
Linear mixed-effects model fit by REML
Data: Orthodont
Log-restricted-likelihood: -221.3183
Fixed: distance ~ age
(Intercept) age
16.7611111 0.6601852
Random effects:
Formula: ~age | Subject
Structure: General positive-definite
StdDev Corr
(Intercept) 2.3270338 (Intr)
age 0.2264276 -0.609
Residual 1.3100399
Number of Observations: 108
Number of Groups: 27
> sessionInfo()
R version 2.6.0 (2007-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] nlme_3.1-85 foreign_0.8-23 arm_1.0-33 R2WinBUGS_2.1-6
[5] coda_0.12-1 lme4_0.99875-8 Matrix_0.999375-2 lattice_0.16-5
[9] MASS_7.2-36
loaded via a namespace (and not attached):
[1] grid_2.6.0 tools_2.6.0
> fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age
> fm1
Linear mixed-effects model fit by REML
Data: Orthodont
Log-restricted-likelihood: -221
Fixed: distance ~ age
(Intercept) age
16.76 0.66
Random effects:
Formula: ~age | Subject
Structure: General positive-definite
StdDev Corr
(Intercept) 2.33 (Intr)
age 0.23 -0.61
Residual 1.31
Number of Observations: 108
Number of Groups: 27
>
More information about the R-help
mailing list