psrwe_survkm: PS-Integrated Kaplan-Meier Estimation

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

PS-Integrated Kaplan-Meier Estimation

Description

Estimate the mean of a survival outcome at a given time point based on PS-integrated Kaplan-Meier approach. Variance can be estimated by Jackknife methods. Apply to the case when there is only one external data source.

Usage

psrwe_survkm(
  dta_psbor,
  pred_tp,
  v_time = "time",
  v_event = "event",
  stderr_method = c("naive", "jk", "sjk", "cjk", "sbs", "cbs", "none"),
  n_bootstrap = 200,
  ...
)

Arguments

dta_psbor	A class PSRWE_BOR object generated by psrwe_borrow.
pred_tp	A numeric value corresponding to time of interest (e.g., 365 days or 1 year)
v_time	Column name corresponding to event time
v_event	Column name corresponding to event status
stderr_method	Method for computing StdErr, see Details
n_bootstrap	Number of bootstrap samples (for bootstrap stderr)
...	Additional Parameters

Details

stderr_method includes naive as default which mostly follows Greenwood formula, jk using Jackknife method within each stratum, sjk using simple Jackknife method for combined estimates such as point estimates in single arm or treatment effects in RCT, or cjk for complex Jackknife method including refitting PS model, matching, trimming, calculating borrowing parameters, and combining overall estimates. Note that sjk may take a while longer to finish and cjk will take even much longer to finish. The sbs and cbs is for simple and complex Bootstrap methods.

Value

A data frame with class name PSRWE_RST. It contains the composite estimation of the mean for each stratum as well as the Jackknife estimation. The results can be further summarized by its S3 method summary. The results can be also analyzed by psrwe_outana for outcome analysis and inference.

Examples

data(ex_dta)
dta_ps <- psrwe_est(ex_dta,
       v_covs = paste("V", 1:7, sep = ""),
       v_grp = "Group",
       cur_grp_level = "current")
ps_borrow <- psrwe_borrow(total_borrow = 30, dta_ps)
rst       <- psrwe_survkm(ps_borrow,
                          pred_tp = 365,
                          v_time = "Y_Surv",
                          v_event = "Status")
rst

olssol/psrwe documentation built on July 17, 2024, 4:06 p.m.