mice (version 3.8.0)

D3: Compare two nested models using D3-statistic

Description

The D3-statistics is a likelihood-ratio test statistic.

Usage

D3(fit1, fit0 = NULL, df.com = Inf, ...)

Arguments

fit1

An object of class mira, produced by with().

fit0

An object of class mira, produced by with(). The model in fit0 is a nested within fit1. The default null model fit0 = NULL compares fit1 to the intercept-only model.

df.com

A single number or a numeric vector denoting the complete-data degrees of freedom for the hypothesis test. If not specified, it is set equal to df.residual of model fit1.

...

Not used.

References

Meng, X. L., and D. B. Rubin. 1992. Performing Likelihood Ratio Tests with Multiply-Imputed Data Sets. Biometrika, 79 (1): 103<U+2013>11.

https://stefvanbuuren.name/fimd/sec-multiparameter.html#sec:likelihoodratio

Examples

# NOT RUN {
# Compare two linear models:
imp <- mice(nhanes2, seed = 51009, print = FALSE)
mi1 <- with(data = imp, expr = lm(bmi ~ age + hyp + chl))
mi0 <- with(data = imp, expr = lm(bmi ~ age + hyp))
D3(mi1, mi0)

# Compare two logistic regression models
imp  <- mice(boys, maxit = 2, print = FALSE)
fit1 <- with(imp, glm(gen > levels(gen)[1] ~ hgt + hc + reg, family = binomial))
fit0 <- with(imp, glm(gen > levels(gen)[1] ~ hgt + hc, family = binomial))
D3(fit1, fit0)
# }