Compare several nested models
Performs model comparisons for a series of nested statistical models fitted using with.mitml.list
.
## S3 method for class 'mitml.result' anova(object, ..., method = c("D3", "D4", "D2"), ariv = c("default", "positive", "robust"), data = NULL)
object |
An object of class |
... |
Additional objects of class |
method |
A character string denoting the method used for the model comparison. Can be |
ariv |
A character string denoting how the ARIV is calculated. Can be |
data |
(optional) A list of imputed data sets (see 'Details'). |
This function performs likelihood-based comparisons between multiple statistical models fitted with with.mitml.list
.
If possible, the models are compared using the D_3 statistic (Meng & Rubin, 1992).
If this method is unavailable, the D_4 or D_2 statistic is used instead (Chan & Meng, 2019; Li, Meng, Raghunathan, & Rubin, 1991).
This function is essentially a wrapper for testModels
with the advantage that several models can be compared simultaneously.
For a list of supported models and further options for more specific model comparisons, see testModels
.
The ariv
argument affects how the average relative increase in variance is calculated (see also testModels
).
Note that the D_4 method can fail if the data to which the model was fitted cannot be found.
In such a case, the data
argument can be used to specify the list of imputed data sets directly (see also testModels
).
A list containing the results of each model comparison.
A print
method is used for more readable output.
Simon Grund
Meng, X.-L., & Rubin, D. B. (1992). Performing likelihood ratio tests with multiply-imputed data sets. Biometrika, 79, 103-111.
Laird, N., Lange, N., & Stram, D. (1987). Maximum likelihood computations with repeated measures: Application of the em algorithm. Journal of the American Statistical Association, 82, 97-105.
Li, K. H., Raghunathan, T. E., & Rubin, D. B. (1991). Large-sample significance levels from multiply imputed data using moment-based statistics and an F reference distribution. Journal of the American Statistical Association, 86, 1065-1073.
require(lme4) data(studentratings) fml <- ReadDis + SES ~ ReadAchiev + (1|ID) imp <- panImpute(studentratings, formula = fml, n.burn = 1000, n.iter = 100, m = 5) implist <- mitmlComplete(imp) # simple comparison (same as testModels) fit0 <- with(implist, lmer(ReadAchiev ~ (1|ID), REML = FALSE)) fit1 <- with(implist, lmer(ReadAchiev ~ ReadDis + (1|ID), REML = FALSE)) anova(fit1, fit0) ## Not run: # multiple comparisons fit2 <- with(implist, lmer(ReadAchiev ~ ReadDis + (1 + ReadDis|ID), REML = FALSE)) anova(fit2, fit1, fit0) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.