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 |
sorted |
should the models be sorted according to their numbers or parameters? |
Details
The test statistic is computed where
is the likelihood of model
with
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-1 Index]