data(api, package = "survey")
dstrata <- apistrat %>%
as_survey_design(strata = stype, weights = pw)
results <- dstrata %>%
group_by(both) %>%
group_map(~survey::svyglm(api00~api99 + stype, .))
# group_map_dfr calls `bind_rows` on the list returned and includes
# grouping variables. This is most useful with a package like `broom`
# but could also be used with survey package functions.
result_coef <- dstrata %>%
group_by(both) %>%
group_map_dfr(
~data.frame(
api99_coef = coef(survey::svyglm(api00~api99 + stype, .))[["api99"]]
)
)
Run the code above in your browser using DataLab