OUprior: Fit OU model for univariate data

In ouxy: Model of Adaptive Trait Evolution

Description

Fit OU model given tree and trait

Usage

OUprior(tree = tree, trait = trait, model = model)

Arguments

tree	An ape: tree object stored in phylo format
trait	a univaraite trait
model	specified model preset "OU".

Details

Parameter estimates α, θ, σ^2 are estimated by BM (when α = 0) or OU model from geiger for the next step analysis with function HyperParam to get the reasonable range of the hyper parameter as well as the ancestral value.

Value

MLE parameter estimates α, θ, σ^2.

References

1. Jhwueng, D-C. (2019) Statistical modeling for adaptive trait evolution. Under review.

2. Jhwueng, D-C., and Vasileios Maroulas. "Adaptive trait evolution in random environment." Journal of Applied Statistics 43.12 (2016): 2310-2324.

3. Hansen, Thomas F., Jason Pienaar, and Steven Hecht Orzack. "A comparative method for studying adaptation to a randomly evolving environment." Evolution: International Journal of Organic Evolution 62.8 (2008): 1965-1977.

Examples

library(ape)
tree<-rcoal(3)
trait<-rnorm(3)
names(trait)<-tree$tip.label
model <- "OU"
OUprior(tree=tree,trait=trait,model=model)

ouxy documentation built on July 2, 2020, 4:05 a.m.