# Small example
data("OnoBurden")
OnoBurden_data_pr_s <- OnoBurden$OnoBurden_data_pr[1:500, ]
# randomization based on marginal population design
target_dist_marginal_s <- OnoBurden$target_dist_marginal[c("gender", "age")]
# design-based estimation
out_design_mar_s <-
design_pAMCE(formula = Y ~ gender + age,
factor_name = "gender",
data = OnoBurden_data_pr_s,
pair_id = OnoBurden_data_pr_s$pair_id,
cluster_id = OnoBurden_data_pr_s$id,
target_dist = target_dist_marginal_s, target_type = "marginal")
summary(out_design_mar_s)
# \donttest{
# Example
data("OnoBurden")
OnoBurden_data_pr <- OnoBurden$OnoBurden_data_pr
# randomization based on marginal population design
target_dist_marginal <- OnoBurden$target_dist_marginal
# design-based estimation
out_design_mar <-
design_pAMCE(formula = Y ~ gender + age + family + race + experience + party + pos_security,
factor_name = c("gender", "age", "experience"),
data = OnoBurden_data_pr,
pair_id = OnoBurden_data_pr$pair_id,
cluster_id = OnoBurden_data_pr$id,
target_dist = target_dist_marginal, target_type = "marginal")
summary(out_design_mar)
# }
Run the code above in your browser using DataLab