tests/varIdent.R

## test for PR#9765
library(nlme)
Orth <- subset(Orthodont, Subject %in% c("M01","F01"))
# Argument fixed in varIdent is ignored
vf  <- varIdent(form=~1|Sex,fixed=c(Female=0.5))
vf <- Initialize(vf, data=Orth)
stopifnot(varWeights(vf) == rep(c(1,2), each=4))
bbolker/nlme documentation built on April 11, 2023, 8:29 p.m.