caper (version 1.0.3)

anova.pgls: Anova and AIC tables for 'pgls' models.

Description

The 'anova' function creates ANOVA tables for a 'pgls' models using sequential sums of squares.

Usage

# S3 method for pgls
anova(object, ...)
# S3 method for pglslist
anova(object, ..., scale = 0, test = "F")

Value

A table of class 'anova' and 'data.frame' that employs the generic plot methods for 'anova' tables.

Arguments

object

A 'pgls' model object.

...

Additional 'pgls' models.

scale

A character string specifying the test statistic to be used. Can be one of "F", "Chisq" or "Cp", with partial matching allowed, or NULL for no test.

test

numeric. An estimate of the noise variance sigma^2. If zero this will be estimated from the largest model considered.

Author

Rob Freckleton, David Orme

Details

The sequential sums of squares are calculated by refitting the model in the order of the terms of the formula and so can take a little time to calculate. Branch length transformations are held at the values of the initial object. The 'logLik.pgls' provides a simple accessor function that allows the use of AIC model comparisons. Note that the generic AIC methods do no checking to ensure that sensible models are being compared.

See Also

pgls

Examples

data(shorebird)
shorebird <- comparative.data(shorebird.tree, shorebird.data, Species, vcv=TRUE, vcv.dim=3)

mod1 <- pgls(log(Egg.Mass) ~ log(M.Mass) * log(F.Mass), shorebird) 
anova(mod1)

mod2 <- pgls(log(Egg.Mass) ~ log(M.Mass) + log(F.Mass), shorebird)  
mod3 <- pgls(log(Egg.Mass) ~ log(M.Mass) , shorebird)
mod4 <- pgls(log(Egg.Mass) ~ 1, shorebird)

anova(mod1, mod2, mod3, mod4)
AIC(mod1, mod2, mod3, mod4)