Lucy D’Agostino McGowan
@lucystats.bsky.social
5.6K followers 178 following 347 posts
Biostatistician • Associate Prof @ Wake Forest University • former postdoc @ Hopkins Biostat • PhD @ Vandy Biostat • 🎙 Casual Inference • lucymcgowan.com
Biostatistician • Associate Prof @ Wake Forest University • former postdoc @ Hopkins Biostat • PhD @ Vandy Biostat • 🎙 Casual Inference • lucymcgowan.com
</sup>](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:xzbbgzmcitcbi6nmi3w2poix/bafkreibh4cuhy6nvnnkn5wtkb7ca43afbtsqlpsdvywxhidbi4qo6jaoku)
![sim <- function() {
n <- 100
x1 <- rnorm(n)
x2 <- rnorm(n)
y <- x1 + x2 + rnorm(n)
x1_obs <- ifelse(x1 > 0, NA, x1)
x2_obs <- ifelse(x2 > 0, NA, x2)
d <- data.frame(y, x1_obs, x2_obs)
cc <- broom::tidy(lm(y ~ x1_obs + x2_obs, data = d))
m <- mice::mice(d, print = F, method = "norm")
imp <- broom::tidy(mice::pool(with(m, lm(y ~ x1_obs + x2_obs))), conf.int = T)
data.frame(
cc = cc[2, "std.error"],
imp = imp[2, "std.error"]
)
}
o <- purrr::map_df(1:100, ~sim())
set.seed(1)
mean(o$std.err)
#> [1] 0.3642946
mean(o$imp)
#> [1] 0.3741719](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:xzbbgzmcitcbi6nmi3w2poix/bafkreicvkpq2y4qchuuh2szuxmaebq5itynkqxi7fi3t7tks5fbexvnbra)
