Displaying 1 result from an estimated 1 matches for "pglm_brain".
2011 Aug 04
0
phyres function in caper package
...unning following phylogenetic analyses with the caper package:
data=read.table(file="data.txt",header=T,sep="\t")
tree = read.nexus("Tree.nex")
primate = comparative.data(phy=tree, data=data,
names.col=Species, vcv=TRUE, na.omit=FALSE, warn.dropped=TRUE)
PGLM_brain=pgls(ln_brain~ln_body,data=primate, lambda="ML")
summary(PGLM_brain)
Now I would like the get the phylogenetic residuals of this model, which I should be able to do with the phyres function:
resbrain=phyres(PGLM_brain)
But then I get the following error message:
Error: could not f...