data(phyllostomid)
plot(phyllostomid$tree)
head(phyllostomid$mandible)
# \donttest{
# Fit a linear model by PL
fit1 <- mvgls(mandible~grp1, data=phyllostomid, phyllostomid$tree, model="lambda", method="LOO")
# regularized MANOVA test
(manova.gls(fit1, test="Wilks", verbose=TRUE))
# }
Run the code above in your browser using DataLab