set.seed(100)
logit_fit = specify(a = ~ sample(0:1, size = n, replace = TRUE),
b = ~ a + rlnorm(n)) %>%
define(n = c(40, 50)) %>%
generate(1) %>%
fit(logit = ~ glm(a ~ b, family = "binomial"))
logit_fit %>%
apply_fits(broom::augment)
## Arguments to the function can be passed in ...
logit_fit %>%
apply_fits(broom::augment, se_fit = TRUE)
## Using a purrr-style lambda function
logit_fit %>%
apply_fits(~ summary(.)$cov.scaled)
Run the code above in your browser using DataLab