varcomp: Compute Variance Component Estimates

In gjuggler/ape: Analyses of Phylogenetics and Evolution

Description

Get variance component estimates from a fitted lme object.

Usage

varcomp(x, scale = FALSE, cum = FALSE)

Arguments

x	A fitted lme object
scale	Scale all variance so that they sum to 1
cum	Send cumulative variance components.

Details

Variance computations is done as in Venables and Ripley's book.

Value

A named vector of class varcomp with estimated variance components.

Author(s)

Julien Dutheil [e-mail]

References

Venables, W. N. and Ripley, B. D. (2002) Modern applied statistics with S (fourth edition). New York: Springer-Verlag.

See Also

lme

Examples

data(carnivora)
library(nlme)
m <- lme(log10(SW) ~ 1, random = ~ 1|Order/SuperFamily/Family/Genus, data=carnivora)
v <- varcomp(m, TRUE, TRUE)
plot(v)

gjuggler/ape documentation built on May 17, 2019, 6:03 a.m.