stat.anova.jmdem: JMDEM Anova Statistics

In jmdem: Fitting Joint Mean and Dispersion Effects Models

Description

This is a utility function, used in jmdem method for anova(..., test != NULL) and should not be used by the average user.

Usage

stat.anova.jmdem(table, test = c("Rao", "Wald")) 

Arguments

table	numeric matrix as results from anova.jmdem(..., test = NULL, print.results = FALSE) saved as the attributes table.x or table.z.
test	a character string, partially matching one of "Rao" or "Wald".

Value

A matrix which is the original table, augmented by a column of test statistics, depending on the test argument.

Author(s)

Karl Wu Ka Yui ([e-mail])

References

Hastie, T. J. and Pregibon, D. (1992). Generalized linear models. Chapter 6 of Statistical Models in S eds J. M. Chambers and T. J. Hastie, Wadsworth & Brooks/Cole.

See Also

anova.jmdem

Examples

## Example in jmdem(...)
MyData <- simdata.jmdem.sim(mformula = y ~ x, dformula = ~ z, 
                            mfamily = poisson(), 
                            dfamily = Gamma(link = "log"), 
                            beta.true = c(0.5, 4), 
                            lambda.true = c(2.5, 3), n = 100)
                            
fit <- jmdem(mformula = y ~ x, dformula = ~ z, data = MyData, 
             mfamily = poisson, dfamily = Gamma(link = "log"), 
             dev.type = "deviance", method = "CG")
              
stat.anova.jmdem(anova(fit, test = "Rao", print.results = FALSE)$table.x)

jmdem documentation built on March 13, 2020, 2:20 a.m.