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)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.