0byt3m1n1
Path:
C:
/
Program Files
/
R
/
R-4.4.1
/
library
/
nlme
/
tests
/
[
Home
]
File: corFactor.R
## test for PR#16110 library(nlme) Orth <- subset(Orthodont, Subject %in% c("M01","F01")) cs1CompSymm <- corCompSymm(value = 0.3, form = ~ 1 | Subject) cs1CompSymm <- Initialize(cs1CompSymm, data = Orth) ## corFactor should return corMatrix(, corr = FALSE) as a vector Linvt <- corMatrix(cs1CompSymm, corr = FALSE) stopifnot(all.equal(unlist(Linvt, use.names = FALSE), as.vector(corFactor(cs1CompSymm)))) ## failed in 3.1-145 because the corFactor.corCompSymm method was ## misnamed corFactor.compSymm (a non-existent class), such that the ## general corFactor.corStruct method was called instead, which returned ## a different solution for the (transpose inverse) square-root factor