[R-sig-ME] Random effects variances in R and SPSS not matching

[Simon Harmel]

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]]