data(admin)
data(jvs)
jvs_svy <- svydesign(ids = ~ 1, weights = ~ weight,
strata = ~ size + nace + region, data = jvs)
ipw_est1 <- nonprob(selection = ~ region + private + nace + size,
target = ~ single_shift,
svydesign = jvs_svy,
data = admin, method_selection = "logit"
)
ipw_est2 <- nonprob(
selection = ~ region + private + nace + size,
target = ~ single_shift,
svydesign = jvs_svy,
data = admin, method_selection = "logit",
control_selection = control_sel(est_method = "gee", gee_h_fun = 1))
## check the balance for the standard IPW
check_balance(~size+private, ipw_est1)
## check the balance for the calibrated IPW
check_balance(~size+private, ipw_est2)
## check balance for a more complicated example
check_balance(~ I(size=="M") + I(nace == "C"), ipw_est1)
Run the code above in your browser using DataLab