# NOT RUN {
data(shorebird)
shorebird <- comparative.data(shorebird.tree, shorebird.data, Species, vcv=TRUE, vcv.dim=3)
mod1 <- pgls(log(Egg.Mass) ~ log(M.Mass) * log(F.Mass), shorebird)
print(mod1)
mod1.sum <- summary(mod1)
print(mod1.sum)
# }
Run the code above in your browser using DataLab