mg <- as.MixedGaussian(PCMBaseTestObjects$model.ab.123.bSigmae_x)
stopifnot(
PCMLik(
X = PCMBaseTestObjects$traits.ab.123,
PCMBaseTestObjects$tree.ab,
PCMBaseTestObjects$model.ab.123.bSigmae_x) ==
PCMLik(
X = PCMBaseTestObjects$traits.ab.123,
PCMBaseTestObjects$tree.ab,
mg))
Run the code above in your browser using DataLab