skip_on_cran()
skip_if_not_installed("mice")
skip_if_not_installed("nnet")
test_that("param", {
set.seed(1234)
d <- suppressWarnings(mice::ampute(mtcars)) ## Ampute mtcars and impute two data sets
imp <- suppressWarnings(mice::mice(d$amp, m = 2, printFlag = FALSE))
imp.l <- mice::complete(imp, action = "long")
model <- list() ## Fit and pool models
for (i in 1:2) {
capture.output({
model[[i]] <- nnet::multinom(cyl ~ disp + hp, data = imp.l, subset = .imp == i)
})
}
pooled <- mice::pool(model)
out1 <- get_parameters(pooled)
out2 <- get_statistic(pooled)
out1$Estimate <- round(out1$Estimate, 4)
out2$Statistic <- round(out2$Statistic, 4)
expect_snapshot(out1)
expect_snapshot(out2)
expect_identical(find_parameters(pooled), list(conditional = c("(Intercept)", "disp", "hp")))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.