# NOT RUN {
data("phyllostomid")
# model fit with mandible~"grp2"
fit <- mvgls(mandible~grp2, data=phyllostomid, phyllostomid$tree, model="lambda", method="PL")
# pairwise tests
pairwise.contrasts(fit, term="grp2")
# }
Run the code above in your browser using DataLab