In jbryer/psa: Applied Propensity Score Analysis with R

Bootstrapping {#chapter-bootstrapping}

library(PSAboot)

boot.matching.1to3 <- function(Tr, Y, X, X.trans, formu, ...) {
    return(boot.matching(Tr=Tr, Y=Y, X=X, X.trans=X.trans, formu=formu, M=3, ...))
}


boot_out <- PSAboot(Tr = lalonde$treat == 1, 
                    Y = lalonde$re78, 
                    X = lalonde[,all.vars(lalonde.formu)[-1]], 
                    seed = 2112,
                    methods=c('Stratification' = boot.strata,
                              'ctree' = boot.ctree,
                              'rpart' = boot.rpart,
                              'Matching' = boot.matching,
                              'Matching-1-to-3' = boot.matching.1to3,
                              'MatchIt' = boot.matchit) )

summary(boot_out)

plot(boot_out)

boxplot(boot_out)

matrixplot(boot_out)

boot_balance <- balance(boot_out)
boot_balance

plot(boot_balance)

boxplot(boot_balance) + geom_hline(yintercept=.1, color='red')

Details are available within the returned object

boot_balance$unadjusted
boot_balance$complete
boot_balance$pooled |> head()

jbryer/psa documentation built on Nov. 17, 2023, 8:21 a.m.