brown.fit {ouch} | R Documentation |
These functions relate to the Brownian motion model for phylogenetic evolution.
brown.fit(data, node, ancestor, times) brown.dev(n = 1, node, ancestor, times, sigma, theta)
data |
Phenotypic data for extant species, i.e., at the terminal ends of the phylogenetic tree. |
node |
Specification of the names of the nodes. |
ancestor |
Specification of the topology of the phylogenetic tree. This is in the form of a character vector of node names, one for each node in the tree. The i-th name is that of the ancestor of the i-th node. The root node is distinguished by having no ancestor (i.e., NA). |
times |
A vector of nonnegative numbers, one per node in the tree, specifying the time at which each node is located. The root node should be assigned time 0. |
n |
the number of simulated data sets to generate. |
sigma |
the value of sigma to be used in the simulations. |
theta |
the value of theta to be used in the simulations. |
brown returns a list of the following elements:
sigma |
Maximum likelihood estimate of sigma. |
theta |
Maximum likelihood estimate of theta. |
u |
-2 log likelihood. |
aic |
Akaike information criterion. |
sic |
Schwartz information criterion (=BIC) |
df |
Number of parameters estimated (= 2). |
brown.dev returns a data frame containing simulated data sets. Each
realization is a row.
Aaron A. King <kingaa at umich dot edu>
Butler, M.A. and A.A. King (2004) Phylogenetic comparative analysis: a modeling approach for adaptive evolution. American Naturalist, in press.