packrat/lib-R/x86_64-w64-mingw32/3.6.1/nlme/tests/lmList.R

In jmcascalheira/LGMIberiaCluster: Research compendium for the paper "Human adaptations and palaeoenvironmental record during the Last Glacial Maximum in the Iberian Peninsula"

## PR#13788

library(nlme)
packageDescription("nlme")# gets the .Library one even in R CMD checking < FIXME !

qm <-  lmList(height ~ age | Subject, data = Oxboys)
nd <- with(Oxboys,
           expand.grid(age = seq(min(age),max(age),length=50),
                       Subject = levels(Subject))
           )

## failed in 3.1-92
res <- predict(qm, nd, se=TRUE)
stopifnot(is.data.frame(res), dim(res) == c(1300, 3),
          identical(names(res), c("Subject", "fit", "se.fit")))

req <- ranef(qm)
(p.req <- plot(req, xlab = "R.Eff.", ylab = "Subj"))
# Fails (p.re2 <- plot(req, age ~ fitted(.)))
iqm <- intervals(qm)
stopifnot(is.array(iqm), dim(iqm) == c(26,3,2))
p.iq <- plot(iqm, ylab = "Subject [factor]")
## Error: formal argument "ylab" matched by multiple .. in 3.1.137
stopifnot(inherits(p.iq, "trellis"),
          inherits(p.req, "trellis"),
          identical( unclass(p.req)[c("xlab","ylab")],
                    list(xlab = "R.Eff.", ylab = "Subj")),
          formula(p.iq) == (group ~ intervals | what))
p.iq