Skip to contents

This function calculates the incidence density of both adverse events and specified competing events observed in [0, tau] and then combines and transforms the incidence densities on a probability scale. Please also refer to formulas (4) and (5) in Stegherr et al. (2021) .

Usage

prop_trans_inc_dens_ce(data, ce, tau)

Arguments

data

(data.frame)
with columns including

  • time_to_event: Time to the first AE, death or soft competing event.

  • type_of_event: 0 for censored, 1 for AE, 2 for death, 3 for soft competing event.

ce

(number)
code for competing event.

tau

(number)
milestone at which Probability Transform Incidence Density is computed.

Value

(vector)
with the following entries:

  • ae_prob: Estimated probability of AE.

  • ae_prob_var: Variance of that estimate.

References

Stegherr R, Schmoor C, Beyersmann J, Rufibach K, Jehl V, Brückner A, Eisele L, Künzel T, Kupas K, Langer F, Leverkus F, Loos A, Norenberg C, Voss F, Friede T (2021). “Survival analysis for AdVerse events with VarYing follow-up times (SAVVY)—estimation of adverse event risks.” Trials, 22(1), 420. ISSN 1745-6215, doi:10.1186/s13063-021-05354-x .

Examples

set.seed(123)
dat <- generate_data(n = 5, cens = c(2, 5), haz_ae = 2, haz_death = 3, haz_soft = 5)
prop_trans_inc_dens_ce(dat, ce = 2, tau = 4)
#>     ae_prob ae_prob_var 
#>       0.200       0.032