geiger (version 1.3-1)

getAncStates: Estimate ancestral character states for continuous characters

Description

Estimates ancestral character states for continuous characters under a Brownian motion model

Usage

getAncStates(x, phy)

Arguments

x

Data vector

phy

Phylogenetic tree in phylo format

Value

• Vector of ancestral character states at each node. Nodes are numbered following the order in phy edge

Details

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.

See Also

ace

Examples

data(geospiza)
attach(geospiza)

tt<-drop.tip(geospiza.tree, "olivacea")
getAncStates(geospiza.data[,1], tt)