[R-sig-ME] lme: count the number extra parameters estimated for variance or covariances

Simon Harmel @|m@h@rme| @end|ng |rom gm@||@com
Sun Nov 1 10:02:15 CET 2020


Hello All,

In addition to fixed and random effects, is there a way to extract how many
other parameters (for modeling residual variances or covariances) an
"lme()" object has estimated?

Here is a reproducible example:

library(nlme)

hsb <- read.csv('
https://raw.githubusercontent.com/rnorouzian/e/master/hsb.csv')
hsb$female <- as.factor(hsb$female)

fit <- lme(math ~ female, random = ~ 1|sch.id, data = hsb, weights =
varIdent(form = ~1 |female))

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list