tests/glssoln.R

In kingaa/ouch: Ornstein-Uhlenbeck Models for Phylogenetic Comparative Hypotheses

library(ouch)

try(ouch:::glssoln(0,1,0))
try(ouch:::glssoln(matrix(c(1,1,0,1),2,2),c(1,1),matrix(c(1,1,1,1),2,2)))
ouch:::glssoln(matrix(c(1,1,0,1),2,2),c(1,1),matrix(c(2,1,1,2),2,2)) -> x
stopifnot(
  coeff=all.equal(as.numeric(x$coeff),c(1,0)),
  residuals=all.equal(as.numeric(x$residuals),c(0,0))
)