# 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