[R-sig-ME] Random effects variances in R and SPSS not matching
Simon Harmel
@|m@h@rme| @end|ng |rom gm@||@com
Wed Mar 31 22:15:24 CEST 2021
Dear All,
For my reproducible model below, SPSS gives the variance component of
119.95 for Y1, and 127.90 for Y2.
But in `nlme::lme()` my variance components are 105.78 for Y1 and 113.73
for Y2.
Can we make the `lme()` reproduce the SPSS's variance components?
#======= Data and R code:
dat <- read.csv('https://raw.githubusercontent.com/hkil/m/master/mv.l.csv')
library(nlme)
m2 <- lme(value ~0 + name, random = ~0 + name| Student, data = dat, method
= "ML")
Random effects variance covariance matrix
nameY1 nameY2
nameY1 105.780 60.869
nameY2 60.869 113.730
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list