inst/demo_cjk/ex_3_single_match.r

In olssol/psrwe: PS-Integrated Methods for Incorporating RWE in Clinical Studies

### Example of single arm and matching
suppressMessages(library(psrwe, quietly = TRUE))
options(digits = 3)
data(ex_dta)

### Obtain PSs.
dta_ps_single <- psrwe_est(ex_dta,
                     v_covs = paste("V", 1:7, sep = ""),
                     v_grp = "Group", cur_grp_level = "current",
                     ps_method = "logistic")

### PS matching.
dta_ps_match <- psrwe_match(dta_ps_single, ratio = 2, strata_covs = "V1")

### Obtain discounting parameters.
ps_bor_match <- psrwe_borrow(dta_ps_match, total_borrow = 30)

### PSCL, single arm study, binary outcome.
rst_cl <- psrwe_compl(ps_bor_match,
                      outcome_type = "binary",
                      v_outcome    = "Y_Bin")

### Outcome analysis.
oa_cl <- psrwe_outana(rst_cl, method_ci = "wilson", mu = 0.40)
oa_cl

### Use complex Jackknife stderr. This may take very long to finish.
rst_cl_cjk <- psrwe_compl(ps_bor_match,
                          outcome_type = "binary",
                          v_outcome = "Y_Bin",
                          stderr_method = "cjk")
oa_cl_cjk <- psrwe_outana(rst_cl_cjk, mu = 0.4)
oa_cl_cjk