getAncStates {geiger} | R Documentation |
Estimates ancestral character states for continuous characters under a Brownian motion model
getAncStates(x, phy)
x |
Data vector |
phy |
Phylogenetic tree in phylo format |
This function uses an algorithm to calculate Maximum-likelihood estimates of ancestral character states for continuous characters.
The ape function ace carries out a similar function but finds estimates by actually maximizing the likelihood across all ancestors simultaneously.
I have implemented an algorithm that is faster and sometimes more reliable especially for large trees.
Vector of ancestral character states at each node. Nodes are numbered following the order in phy edge
Luke J. Harmon
ace
data(geospiza) attach(geospiza) tt<-drop.tip(geospiza.tree, "olivacea") getAncStates(geospiza.data[,1], tt)