[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')


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