R: Compare two nested models using D3-statistic

D3 {mice} | R Documentation

## 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–11.

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

### Examples

    # 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)

[Package mice version 3.8.0 Index]