anova.mdm: Analysis of Deviance, Entropy and Diversity Tables

Description Usage Arguments Details Value Note References See Also Examples

View source: R/anova.mdm.R

Description

Provides an analysis of deviance, entropy and diversity for a collection of diversity models (outputs from mdm).

Usage

1
2
## S3 method for class 'mdm'
anova(object, ..., topnote = TRUE, cols = c("df","dev","ent","div")[1:4])

Arguments

object, ...

objects of class mdm, usually, a result of a call to mdm.

topnote

If TRUE then model descriptions appear above the anova table, if FALSE they appear as the first column of the table

cols

The list of colums to print out. Defaults to all columns.

Details

Specifying a single object gives a sequential analysis of deviance table for that fit. That is, the reductions in the residual sum of deviances, plus the residual deviance The deviances are converted into entropies and diversities, and differences in deviances are converted into differnces in entropies and diversities.

If more than one object is specified, the table has a row for the residual degrees of freedom and sum of deviances for each model. For all but the first model, the change in degrees of freedom and deviances. This only makes statistical sense if the models are nested. It is conventional to list the models from smallest to largest, but this is up to the user.

Value

The analysis of deviance, entropy and diversity for a collection of diversity models.

Note

mdm is a modifed version of multinom in the nnet package.

References

De'ath, G. (2011) The Multinomial Diversity Model: Linking Shannon Diversity To Multiple Predictors.

Ripley, B. D. (1996) Pattern Recognition and Neural Networks. Cambridge.

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.

See Also

multinom, nnet

Examples

1
2
3
4
5
6
7
library(MDM)
data(spider6)
fit0 <- mdm(y2p(spider6[,1:6])~1,data=spider6)
fit1 <- mdm(y2p(spider6[,1:6])~Water,data=spider6)
fit2 <- mdm(y2p(spider6[,1:6])~Water+Herbs,data=spider6)
fit3 <- mdm(y2p(spider6[,1:6])~Site,data=spider6)
anova(fit0,fit1,fit2,fit3)

MDM documentation built on May 30, 2017, 4:58 a.m.