library(magrittr)
lbpmilr %>%
df2milist(impvar="Impnr") %>%
with(expr=prop_wald(Radiation ~ 1)) %>%
pool_prop_wilson()
# Same as
imp_dat <- df2milist(lbpmilr, impvar="Impnr")
ra <- with(imp_dat, expr=prop_wald(Radiation ~ 1))
res <- pool_prop_wilson(ra)
Run the code above in your browser using DataLab