anova.grouped {grouped} | R Documentation |
Performs a Likelihood Ratio Test between two nested grouped
models.
## S3 method for class 'grouped': anova(object, object2, ...)
object |
an object inheriting from class grouped , nested to object2 . |
object2 |
an object inheriting from class grouped . |
... |
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 object . |
L0 |
the log-likelihood under object . |
df0 |
the number of parameters in object . |
AIC0 |
the AIC under object . |
BIC0 |
the BIC under object . |
name1 |
the name of the alternative model represented by object2 . |
L1 |
the log-likelihood under object2 . |
df1 |
the number of parameters in object2 . |
AIC1 |
the AIC under object2 . |
BIC1 |
the BIC under object2 . |
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)