blendR (version 1.0.0)

ext_surv_sim: Create an external survival data based on expert opinion

Description

Generally, the sampling is done is two steps $$ p(T) = p(T | interval i) p(interval i) $$

Usage

ext_surv_sim(t_info, S_info, T_max, n = 100)

Value

Dataframe of times and censoring status

Arguments

t_info

A vector of times for which expert opinion is elicited

S_info

A vector of mean survival probabilities estimated by experts corresponding to time points in t_info

T_max

The maximum survival time to be used

n

The number of patients to construct the artificial external data set; default 100

Details

In particular \(T ~ U(x_{i}, x_{i+1})\) \(i ~ multinomial(\hat{\pi})\)

Examples

dat <- ext_surv_sim(t_info = c(10,20,50),
                    S_info = c(0.9, 0.8, 0.2),
                    T_max = 100, n = 100)
if (require(survival)) {
    km_fit <- survfit(Surv(time, event) ~ 1, data = dat)
    plot(km_fit)
}