demo/phyloPCA_demo.R

#library(ape)

tree<-pbtree(n=100, scale=1)
mm<-cbind(c(1, -0.5), c(-0.5, 1))
z<-sim.char(tree, mm, model="BM")[,,1]
colnames(z)<-c("x", "y")
td<-make.treedata(tree, z)

phylopcaArbor(td)
arborworkflows/aRbor documentation built on June 11, 2020, 7:44 p.m.