Displaying 1 result from an estimated 1 matches for "mdl1g".
Did you mean:
mdl1
2007 May 21
1
can I get same results using lme and gls?
..."V20","V21","V22","V23","V24","V25"))
colnames(dta1)[1] = "schoolNR"
dta2 = dta1[order(dta1$id),]
head(dta2)
timef = factor(dta2$time)
summary(mdl1l <- lme(score~timef-1, dta2, ~timef-1|schoolNR/id,,,,"ML"))
summary(mdl1g <- gls(score~timef-1, dta2, corCompSymm(, ~timef|schoolNR/id),
varIdent(, ~1|id*timef),,"ML"))