tests/testthat/test-mira.R

In parameters: Processing of Model Parameters

skip_if_not_installed("mice")
data("nhanes2", package = "mice")
imp <- mice::mice(nhanes2, printFlag = FALSE)
fit <- with(data = imp, exp = lm(bmi ~ age + hyp + chl))

mp1 <- model_parameters(fit)
mp2 <- summary(mice::pool(fit))

test_that("param", {
  expect_equal(mp1$Parameter, as.vector(mp2$term))
})

test_that("coef", {
  expect_equal(mp1$Coefficient, mp2$estimate, tolerance = 1e-3)
})

test_that("se", {
  expect_equal(mp1$SE, mp2$std.error, tolerance = 1e-3)
})