demo/old/boot_control.R

In harvard-P01/causalMediation: Causal Mediation Analysis

library(causalMediation)

df <- data.frame('smoking'    = c(0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0),
                 'lbw'        = c(0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0),
                 'death'      = c(0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
                 'drinking'   = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0),
                 'agebelow20' = c(0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
)


result_delta_bin_bin_noint_boot <- causalMediationAll(data = df, 
                                                      outcome = 'death', 
                                                      treatment = 'smoking',
                                                      mediator = 'lbw',
                                                      covariates = 'drinking',
                                                      vecc = 1,
                                                      interaction = FALSE,
                                                      yreg = "logistic", mreg = "logistic",
                                                      boot = TRUE)
result_delta_bin_bin_noint_boot

boot(data = df, statistic = causalMediationAll, R = 10,
     outcome = 'death', 
     treatment = 'smoking',
     mediator = 'lbw',
     covariates = 'drinking',
     vecc = 1,
     interaction = FALSE,
     yreg = "logistic", mreg = "logistic",
     boot = TRUE)

result_delta_bin_bin_int <- causalMediationAll(data = df, 
                                               outcome = 'death', 
                                               treatment = 'smoking',
                                               mediator = 'lbw',
                                               covariates = 'drinking',
                                               vecc = 1,
                                               interaction = TRUE,
                                               yreg = "logistic", mreg = "logistic",
                                               boot = FALSE)
result_delta_bin_bin_int

result_delta_cont_bin_noint <- causalMediationAll(data = df, 
                                                  outcome = 'death', 
                                                  treatment = 'smoking',
                                                  mediator = 'lbw',
                                                  covariates = 'drinking',
                                                  vecc = 1,
                                                  interaction = FALSE,
                                                  yreg = "linear", mreg = "logistic",
                                                  boot = FALSE)
result_delta_cont_bin_noint

result_delta_cont_bin_int <- causalMediationAll(data = df, 
                                                outcome = 'death', 
                                                treatment = 'smoking',
                                                mediator = 'lbw',
                                                covariates = 'drinking',
                                                vecc = 1,
                                                interaction = TRUE,
                                                yreg = "linear", mreg = "logistic",
                                                boot = FALSE)
result_delta_cont_bin_int

result_delta_bin_cont_noint <- causalMediationAll(data = df, 
                                                  outcome = 'death', 
                                                  treatment = 'smoking',
                                                  mediator = 'lbw',
                                                  covariates = 'drinking',
                                                  vecc = 1,
                                                  interaction = FALSE,
                                                  yreg = "logistic", mreg = "linear",
                                                  boot = FALSE)
result_delta_bin_cont_noint

result_delta_bin_cont_int <- causalMediationAll(data = df, 
                                                outcome = 'death', 
                                                treatment = 'smoking',
                                                mediator = 'lbw',
                                                covariates = 'drinking',
                                                vecc = 1,
                                                interaction = TRUE,
                                                yreg = "logistic", mreg = "linear",
                                                boot = FALSE)
result_delta_bin_cont_int

result_delta_cont_cont_noint <- causalMediationAll(data = df, 
                                                   outcome = 'death', 
                                                   treatment = 'smoking',
                                                   mediator = 'lbw',
                                                   covariates = 'drinking',
                                                   vecc = 1,
                                                   interaction = FALSE,
                                                   yreg = "linear", mreg = "linear",
                                                   boot = FALSE)
result_delta_cont_cont_noint

result_delta_cont_cont_int <- causalMediationAll(data = df, 
                                                 outcome = 'death', 
                                                 treatment = 'smoking',
                                                 mediator = 'lbw',
                                                 covariates = 'drinking',
                                                 vecc = 1,
                                                 interaction = TRUE,
                                                 yreg = "linear", mreg = "linear",
                                                 boot = FALSE)
result_delta_cont_cont_int