context("blocks")
imp <- mice(nhanes, blocks = make.blocks(list(c("bmi", "chl"), "bmi", "age")), m = 10, print = FALSE)
# plot(imp)
test_that("removes variables from 'where'", {
expect_identical(sum(imp$where[, "hyp"]), 0L)
})
# reprex https://github.com/amices/mice/issues/326
imp1 <- mice(nhanes, seed = 1, m = 1, maxit = 2, print = FALSE)
imp2 <- mice(nhanes, blocks = list(c("bmi", "hyp"), "chl"), modeltype = c("pred", "pred"), m = 1, maxit = 2, seed = 1, print = FALSE)
test_that("expands a univariate method to all variables in the block", {
expect_identical(complete(imp1, 1), complete(imp2, 1))
})
imp3 <- mice(nhanes, blocks = list(c("hyp", "bmi"), "chl"), m = 1, maxit = 2, seed = 1, print = FALSE)
imp4 <- mice(nhanes, visitSequence = c("hyp", "bmi", "chl"), m = 1, maxit = 2, seed = 1, print = FALSE)
test_that("blocks alter the visit sequence", {
expect_identical(complete(imp3, 1), complete(imp3, 1))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.