R: Compare Dirichlet Regression Models using an LRT

anova.DirichletRegModel {DirichletReg} | R Documentation

## Compare Dirichlet Regression Models using an LRT

### Description

This function allows for pairwise tests of Dirichlet regression models using a likelihood ratio test (LRT).

### Usage

    ## S3 method for class 'DirichletRegModel'
    anova(object, ..., sorted = FALSE)

### Arguments

`object` | the model to be compared against those listed in ...
`...` | models to be tested against the one specified as `object`
`sorted` | should the models be sorted according to their numbers or parameters?

### Details

The test statistic is computed LR=-2[log(La)-log(Lb)] where Li is the likelihood of model i with df equal to the difference of the number of parameters in the models.

### Author(s)

Marco J. Maier

### Examples

    ALake <- ArcticLake
    ALake$AL <- DR_data(ArcticLake[,1:3])
    mod0 <- DirichReg(AL ~ 1, ALake)
    mod1 <- DirichReg(AL ~ depth, ALake)
    mod2 <- DirichReg(AL ~ depth + I(depth^2), ALake)
    anova(mod1, mod0, mod2, sorted = TRUE)

[Package DirichletReg version 0.7-0 Index]