Anova method for grouped objects
Performs a Likelihood Ratio Test between two nested grouped models.
## S3 method for class 'grouped' anova(object, object2, ...)
object |
an object inheriting from class |
object2 |
an object inheriting from class |
... |
additional arguments; currently none is used. |
a list of class aov.grouped with the following components:
name0 |
the name of the null model represented by |
L0 |
the log-likelihood under |
df0 |
the number of parameters in |
AIC0 |
the AIC under |
BIC0 |
the BIC under |
name1 |
the name of the alternative model represented by |
L1 |
the log-likelihood under |
df1 |
the number of parameters in |
AIC1 |
the AIC under |
BIC1 |
the BIC under |
L01 |
the value of the likelihood ratio test statistic. |
p.value |
the p-value of the test. |
The function does only partial checking whether the two models are nested; the user is responsible to supply nested models in order to perform a valid test.
anova.grouped performs a likelihood ratio test between two nested models; for simple Wald tests
for the estimated parameters use summary.grouped.
Dimitris Rizopoulos d.rizopoulos@erasmusmc.nl
m1 <- grouped(cbind(lo, up) ~ treat, link = "logit", data = Sdata) m2 <- grouped(cbind(lo, up) ~ treat * x, link = "logit", data = Sdata) anova(m1, m2) m1 <- grouped(equispaced(r, n) ~ x1, link = "logit", data = Seeds) m2 <- grouped(equispaced(r, n) ~ x1 * x2, link = "logit", data = Seeds) anova(m1, m2)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.