# anova.mdm: Analysis of Deviance, Entropy and Diversity Tables In MDM: Multinomial Diversity Model

## 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.

`multinom`, `nnet`
 ```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) ```