waldtest.BinaryEPPM {BinaryEPPM} | R Documentation |
waldtest
is a generic function for comparisons of nested (generalized) linear models via Wald tests.
## S3 method for class 'BinaryEPPM' waldtest(object, ..., vcov = NULL, test = c("Chisq", "F"))
object |
an object of class "BinaryEPPM". |
... |
further object specifications passed to methods. See below for details. |
vcov |
a function for estimating the covariance matrix of the regression coefficients. If only two models are compared it can also be the covariance matrix of the more general model. |
test |
character specifying whether to compute the large sample Chi-squared statistic (with asymptotic Chi-squared distribution) or the finite sample F statistic (with approximate F distribution). |
waldtest
is a generic function for comparisons of nested (generalized)linear models via Wald tests. It does not have the same functionality as the versions of betareg and lmtest with a reduced list of arguments. With these caveats, more details can be obtained from the Details pages of those packages.
An object of class "anova"
which contains the residual degrees of freedom, the difference in degrees of freedom, Wald statistic (either "Chisq"
or "F"
) and corresponding p value.
David M. Smith <smithdm1@us.ibm.com>
Cribari-Neto F, Zeileis A. (2010). Beta Regression in R. Journal of Statistical Software, 34(2), 1-24. doi: 10.18637/jss.v034.i02.
Zeileis A, Hothorn T. (2002). Diagnostic Checking in Regression Relationships. R News, 2(3), 7-10. https://CRAN.R-project.org/doc/Rnews/.
## Not run: data("ropespores.case") output.fn <- BinaryEPPM(data = ropespores.case, number.spores / number.tested ~ 1 + offset(logdilution), model.type = 'p only', model.name = 'binomial') output.fn.one <- BinaryEPPM(data = ropespores.case, number.spores / number.tested ~ 1 + offset(logdilution), model.type = 'p only', model.name = 'beta binomial') waldtest(output.fn, output.fn.one, test = c("Chisq", "F"), vcov = vcov) ## End(Not run)