getAncStates {geiger}R Documentation

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

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.

Value

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

Author(s)

Luke J. Harmon

See Also

ace

Examples


data(geospiza)
attach(geospiza)

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


[Package geiger version 1.3-1 Index]