Title: | Objects to Simulate Survival Times |
---|---|
Description: | Generate objects that simulate survival times. Random values for the distributions are generated using the method described by Bender (2003) <https://epub.ub.uni-muenchen.de/id/eprint/1716> and Leemis (1987) in Operations Research, 35(6), 892–894. |
Authors: | Aponte John [aut, cre, cph] |
Maintainer: | Aponte John <[email protected]> |
License: | GPL (>= 3) |
Version: | 3.1.1 |
Built: | 2024-10-13 06:51:40 UTC |
Source: | https://github.com/johnaponte/survobj |
if censor_time < time, event is change to 0, otherwise not changed
censor_event(censor_time, time, event) censor_time(censor_time, time)
censor_event(censor_time, time, event) censor_time(censor_time, time)
censor_time |
the time to censor |
time |
the time variable where the censor_time is applied |
event |
the variable with the event. It can be initialized in the call with a value for all times. |
if censor_time < time, time is changed to censor_time, otherwise no change
Be careful and do not overwrite the time variable with the censor time variable to not loose track of the events
censored time or event
censor_time()
: Censor time
# Typical workflow in a simulation of survival time. # Simulate time to event (sim_t_event) # and simulates the time to lost to follow up (tim_t_ltof) # the simulation time frame is 1, so everything after 1 is censored require(dplyr) data.frame( sim_t_event = c(0.5,0.6,1,10,20), sim_t_ltof = c(2,0.5,2,2,0.8) ) |> mutate(sevent = censor_event(1,sim_t_event,sim_event=1)) |> mutate(stime = censor_time(1,sim_t_event)) |> mutate(event = censor_event(sim_t_ltof, stime, sevent)) |> mutate(timeto = censor_time(sim_t_ltof, stime))
# Typical workflow in a simulation of survival time. # Simulate time to event (sim_t_event) # and simulates the time to lost to follow up (tim_t_ltof) # the simulation time frame is 1, so everything after 1 is censored require(dplyr) data.frame( sim_t_event = c(0.5,0.6,1,10,20), sim_t_ltof = c(2,0.5,2,2,0.8) ) |> mutate(sevent = censor_event(1,sim_t_event,sim_event=1)) |> mutate(stime = censor_time(1,sim_t_event)) |> mutate(event = censor_event(sim_t_ltof, stime, sevent)) |> mutate(timeto = censor_time(sim_t_ltof, stime))
This set of functions help in the time conversion, taking into account generic times and not specific times. The conversions are based on the assumption that 1 year is 365.25 days and is 12 months. There is no adjustment for lap days or ours or difference of days between months
dtom(x) mtod(x) dtoy(x) ytod(x) mtoy(x) ytom(x)
dtom(x) mtod(x) dtoy(x) ytod(x) mtoy(x) ytom(x)
x |
the time to convert |
the converted time
dtom()
: convert days to months
mtod()
: convert months to days
dtoy()
: convert days to years
ytod()
: convert years to days
mtoy()
: convert months to year
ytom()
: convert years to months
dtom(365.25) mtod(12) dtoy(165.25) ytod(1) mtoy(12) ytom(365.25)
dtom(365.25) mtod(12) dtoy(165.25) ytod(1) mtoy(12) ytom(365.25)
Evaluates if the argument is a single number
is_single_number(x)
is_single_number(x)
x |
a variable to evaluate |
TRUE
if is a single number, FALSE
otherwise
is_single_number(3) #TRUE is_single_number(c(3,3,3)) #FALSE is_single_number(list(a=3)) #FALSE is_single_number("3") #FALSE
is_single_number(3) #TRUE is_single_number(c(3,3,3)) #FALSE is_single_number(list(a=3)) #FALSE is_single_number("3") #FALSE
Generate a random draws from the distribution of recurrent set of survival times under a a non homogeneous Poisson process following Leemis (1987)
nhpphr(SURVIVAL, hr, prevtime) nhppaft(SURVIVAL, aft, prevtime)
nhpphr(SURVIVAL, hr, prevtime) nhppaft(SURVIVAL, aft, prevtime)
SURVIVAL |
Object of survival class |
hr |
Vector of hazard ratios |
prevtime |
Vector of previous survival times |
aft |
Vector of accelerated failure ratios |
Vector of survival times
nhpphr()
: Recurrent episodes under a proportional hazard model
nhppaft()
: Recurrent episodes under an accelerated failure time model
s_obj <- s_exponential(fail = 0.4, t = 1) hr <- c(1,1,0.5,0.5) time1 <- rsurvhr(s_obj, hr) time2 <- nhpphr(s_obj,hr, time1) s_obj <- s_exponential(fail = 0.4, t = 1) aft <- c(1,1,0.5,0.5) timea <- rsurvaft(s_obj, aft) timeb <- nhppaft(s_obj, aft, timea)
s_obj <- s_exponential(fail = 0.4, t = 1) hr <- c(1,1,0.5,0.5) time1 <- rsurvhr(s_obj, hr) time2 <- nhpphr(s_obj,hr, time1) s_obj <- s_exponential(fail = 0.4, t = 1) aft <- c(1,1,0.5,0.5) timea <- rsurvaft(s_obj, aft) timeb <- nhppaft(s_obj, aft, timea)
Generate a random draws from the distribution of recurrent set of survival times under a renewal Poisson process following Leemis (1987)
renewhr(SURVIVAL, hr, prevtime) renewaft(SURVIVAL, prevtime, aft)
renewhr(SURVIVAL, hr, prevtime) renewaft(SURVIVAL, prevtime, aft)
SURVIVAL |
Object of survival class |
hr |
Vector of hazard ratios |
prevtime |
Vector of previous survival times |
aft |
Vector of accelerated failure time ratios |
Vector of survival times
Vector of survival times
renewhr()
: Recurrent episodes under a proportional hazard model
renewaft()
: Recurrent episodes under an accelerated failure time model
s_obj <- s_exponential(fail = 0.4, t = 1) hr <- c(1,1,0.5,0.5) time1 <- rsurvhr(s_obj, hr) time2 <- renewhr(s_obj, hr, time1) s_obj2 <- s_exponential(fail = 0.4, t = 1) aft <- c(1,1,0.5,0.5) timea <- rsurvaft(s_obj2, aft) timeb <- renewaft(s_obj2, aft, timea)
s_obj <- s_exponential(fail = 0.4, t = 1) hr <- c(1,1,0.5,0.5) time1 <- rsurvhr(s_obj, hr) time2 <- renewhr(s_obj, hr, time1) s_obj2 <- s_exponential(fail = 0.4, t = 1) aft <- c(1,1,0.5,0.5) timea <- rsurvaft(s_obj2, aft) timeb <- renewaft(s_obj2, aft, timea)
Creates a SURVIVAL object with an Exponential distribution.
s_exponential(...)
s_exponential(...)
... |
Parameters to define the distribution. See the Parameters for details |
a SURVIVAL object of the exponential distribution family. See the
documentation of s_factory
for the methods available for SURVIVAL objects
To create an exponential survival object the following options are available:
lambda
to specify the canonical parameter of the distribution, or
surv
and t
for the proportion surviving (no events) at time t, or
fail
and t
for the proportion failing (events) at time t
lambda = -log(surv)/t
lambda = -log(1-fail)/t
The parameters should be spell correctly as partial matching is not available
s_exponential(lambda = 3) s_exponential(surv = 0.4, t = 2) s_exponential(fail = 0.6, t = 2)
s_exponential(lambda = 3) s_exponential(surv = 0.4, t = 2) s_exponential(fail = 0.6, t = 2)
Create objects of the class SURVIVAL
s_factory(s_family, ...)
s_factory(s_family, ...)
s_family |
a factory for a specific distribution |
... |
parameters to define the survival distribution |
The objects of the class SURVIVAL define different distributions of survival times. Each class has its own set of parameters but once the SURVIVAL object is defined, they have access to the same functions to calculate:
survival time function: sfx()
,
hazard time function: hfx()
,
cumulative hazard function: Cum_Hfx()
the inverse of the cumulative hazard function: invCum_Hfx()
.
generate random survival times: rsurv()
generate random survival times under proportional hazard ratio: rsurvhr()
.
There several functions to plot the distributions
generic S3: plot.SURVIVAL()
plot_survival()
: to plot the functions
ggplot_survival_random()
: to ggplot random draws from the distribution
compare_survival()
: to compare the functions of two SURVIVAL objects
a SURVIVAL object
The current factories are implemented:
s_exponential()
: for Exponential distributions
s_weibull()
: for Weibull distributions
s_gompertz()
: for Gompertz distributions
s_picewise()
: for Piecewise exponential distributions
# Define a SURVIVAL object obj <- s_factory(s_exponential, lambda = 2) # Survival, Hazard and Cumulative hazard at time 0.4 sfx(SURVIVAL = obj, t= 0.4) hfx(SURVIVAL = obj, t = 0.4) Cum_Hfx(SURVIVAL = obj, t = 0.4) # Time when the Cumulative hazard is 0.8 invCum_Hfx(SURVIVAL = obj, H = 0.8) # Draw one random survival time from the distribution rsurv(SURVIVAL = obj, n = 1) # Draw one random survival time from the distribution, with hazard ratio 0.5 rsurvhr(SURVIVAL = obj, hr = 0.5) # Plot the survival functions plot(obj)
# Define a SURVIVAL object obj <- s_factory(s_exponential, lambda = 2) # Survival, Hazard and Cumulative hazard at time 0.4 sfx(SURVIVAL = obj, t= 0.4) hfx(SURVIVAL = obj, t = 0.4) Cum_Hfx(SURVIVAL = obj, t = 0.4) # Time when the Cumulative hazard is 0.8 invCum_Hfx(SURVIVAL = obj, H = 0.8) # Draw one random survival time from the distribution rsurv(SURVIVAL = obj, n = 1) # Draw one random survival time from the distribution, with hazard ratio 0.5 rsurvhr(SURVIVAL = obj, hr = 0.5) # Plot the survival functions plot(obj)
Creates a SURVIVAL object with an Gompertz distribution.
s_gompertz(...)
s_gompertz(...)
... |
Parameters to define the distribution. See the Parameters for details |
a SURVIVAL object of the Gompertz distribution family. See the
documentation of s_factory
for the methods available for SURVIVAL objects
To create an exponential survival object the following options are available:
scale
and shape
to specify the canonical parameter of the distribution, or
surv
, t
and shape
for the proportion surviving (no events) at time t and shape, or
fail
and t
and shape
for the proportion failing (events) at time t and shape.
scale = -log(surv)·shape/(exp(shape·t))
scale = -log(1-fail)·shape/(exp(shape·t))
The parameters should be spell correctly as partial matching is not available
s_gompertz(scale = 1, shape = 1.5) s_gompertz(surv = 0.4, t = 2, shape = 1.5) s_gompertz(fail = 0.6, t = 2, shape = 1.5)
s_gompertz(scale = 1, shape = 1.5) s_gompertz(surv = 0.4, t = 2, shape = 1.5) s_gompertz(fail = 0.6, t = 2, shape = 1.5)
Creates a SURVIVAL object with a Log Logistic distribution.
s_loglogistic(...)
s_loglogistic(...)
... |
Parameters to define the distribution. See the Parameters for details |
a SURVIVAL object of the log-logistic distribution family. See the
documentation of s_factory
for the methods available for SURVIVAL objects
To create an exponential survival object the following options are available:
scale
and shape
to specify the canonical parameters of the distribution, or
surv
, t
and shape
for the proportion surviving (no events) at time t and the shape parameter, or
fail
, t
and shape
for the proportion failing (events) at time t and the shape parameter or
intercept
and scale
for the parameters returned by survreg(.., dist = "loglogistic")
models.
The parameters should be spell correctly as partial matching is not available
s_loglogistic(scale = 2,shape = 2) s_loglogistic(surv = 0.6, t= 12, shape = 0.5) s_loglogistic(fail = 0.4, t = 12, shape =0.5) s_loglogistic(intercept = 0.4, scale = 0.5)
s_loglogistic(scale = 2,shape = 2) s_loglogistic(surv = 0.6, t= 12, shape = 0.5) s_loglogistic(fail = 0.4, t = 12, shape =0.5) s_loglogistic(intercept = 0.4, scale = 0.5)
Creates a SURVIVAL object with a Log Normal distribution.
s_lognormal(...)
s_lognormal(...)
... |
Parameters to define the distribution. See the Parameters for details |
a SURVIVAL object of the log-normal distribution family. See the
documentation of s_factory
for the methods available for SURVIVAL objects
To create an exponential survival object the following options are available:
scale
and shape
to specify the canonical parameters of the distribution, or
surv
, t
and shape
for the proportion surviving (no events) at time t and the shape parameter, or
fail
, t
and shape
for the proportion failing (events) at time t and the shape parameter or
intercept
and shape
for the parameters returned by survreg(.., dist = "lognormal")
models.
The scale parameter is the median value of the distribution, and the shape is the log standard deviation
The parameters should be spell correctly as partial matching is not available
s_lognormal(scale = 2,shape = 2) s_lognormal(surv = 0.6, t= 12, shape = 0.5) s_lognormal(fail = 0.4, t = 12, shape =0.5) s_lognormal(intercept = 0.4, scale = 0.5)
s_lognormal(scale = 2,shape = 2) s_lognormal(surv = 0.6, t= 12, shape = 0.5) s_lognormal(fail = 0.4, t = 12, shape =0.5) s_lognormal(intercept = 0.4, scale = 0.5)
Creates a SURVIVAL object with an Piecewise Exponential distribution.
s_piecewise(...)
s_piecewise(...)
... |
Parameters to define the distribution. See the Parameters for details |
a SURVIVAL object of the piecewise exponential distribution family. See the
documentation of s_factory
for the methods available for SURVIVAL objects
To create an piecewise exponential survival object the following options are available:
breaks
and hazards
to specify the exponential (constant) hazard until each break, or
surv
, breaks
and segments
for the proportion surviving (no events) at the end of last segment or
fail
, breaks
and segments
for the proportion failing (events) at the end of last segment
If surv
or fail
parameters are indicated, the segments
are scaled to hazards in order
to mach the surviving or failing proportion at the end of the last segment.
Define the last break point as Inf
to fully define the distribution, otherwise
an error will be produce if function after the last break is requested
The parameters should be spell correctly as partial matching is not available
s_piecewise(breaks = c(1,2,3,Inf), hazards = c(0.5,0.6,0.5,0.1)) s_piecewise(surv = 0.4, breaks = c(1,2,3,Inf), segments = c(1,2,1,2)) s_piecewise(fail = 0.6, breaks = c(1,2,3,Inf), segments = c(1,2,1,2))
s_piecewise(breaks = c(1,2,3,Inf), hazards = c(0.5,0.6,0.5,0.1)) s_piecewise(surv = 0.4, breaks = c(1,2,3,Inf), segments = c(1,2,1,2)) s_piecewise(fail = 0.6, breaks = c(1,2,3,Inf), segments = c(1,2,1,2))
Creates a SURVIVAL object with a Weibull distribution.
s_weibull(...)
s_weibull(...)
... |
Parameters to define the distribution. See the Parameters for details |
a SURVIVAL object of the Weibull distribution family. See the
documentation of s_factory
for the methods available for SURVIVAL objects
To create an exponential survival object the following options are available:
scale
and shape
to specify the canonical parameters of the distribution, or
surv
, t
and shape
for the proportion surviving (no events) at time t and the shape parameter, or
fail
, t
and shape
for the proportion failing (events) at time t and the shape parameter or
intercept
and scale
for the parameters returned by survreg(.., dist = "weibull")
models.
scale = -log(surv)/(t^shape)
scale = -log(1-fail)/(t^shape)
The parameters should be spell correctly as partial matching is not available
s_weibull(scale = 2,shape = 2) s_weibull(surv = 0.6, t= 12, shape = 0.5) s_weibull(fail = 0.4, t = 12, shape =0.5) s_weibull(intercept = 0.4, scale = 0.5)
s_weibull(scale = 2,shape = 2) s_weibull(surv = 0.6, t= 12, shape = 0.5) s_weibull(fail = 0.4, t = 12, shape =0.5) s_weibull(intercept = 0.4, scale = 0.5)
All the SURVIVAL objects have access to the functions described here
sfx(SURVIVAL, t) hfx(SURVIVAL, t) Cum_Hfx(SURVIVAL, t) invCum_Hfx(SURVIVAL, H) rsurv(SURVIVAL, n) rsurvhr(SURVIVAL, hr) rsurvaft(SURVIVAL, aft) rsurvah(SURVIVAL, aft, hr) plot_survival(SURVIVAL, timeto, main) ggplot_survival_random(SURVIVAL, timeto, subjects, nsim, alpha = 0.1) compare_survival(SURVIVAL1, SURVIVAL2, timeto, main) ggplot_survival_hr(SURVIVAL, hr, timeto, subjects, nsim, alpha = 0.1) ggplot_survival_aft(SURVIVAL, aft, timeto, subjects, nsim, alpha = 0.1) ggplot_survival_ah(SURVIVAL, aft, hr, timeto, subjects, nsim, alpha = 0.1)
sfx(SURVIVAL, t) hfx(SURVIVAL, t) Cum_Hfx(SURVIVAL, t) invCum_Hfx(SURVIVAL, H) rsurv(SURVIVAL, n) rsurvhr(SURVIVAL, hr) rsurvaft(SURVIVAL, aft) rsurvah(SURVIVAL, aft, hr) plot_survival(SURVIVAL, timeto, main) ggplot_survival_random(SURVIVAL, timeto, subjects, nsim, alpha = 0.1) compare_survival(SURVIVAL1, SURVIVAL2, timeto, main) ggplot_survival_hr(SURVIVAL, hr, timeto, subjects, nsim, alpha = 0.1) ggplot_survival_aft(SURVIVAL, aft, timeto, subjects, nsim, alpha = 0.1) ggplot_survival_ah(SURVIVAL, aft, hr, timeto, subjects, nsim, alpha = 0.1)
SURVIVAL |
a SURVIVAL object |
t |
Time |
H |
cumulative hazard |
n |
number of observations |
hr |
hazard ratio |
aft |
accelerated failure time |
timeto |
plot the distribution up to timeto |
main |
title of the graph |
subjects |
number of subjects per group to simulate in each simulation |
nsim |
number of simulations |
alpha |
alpha value for the graph |
SURVIVAL1 |
a SURVIVAL object |
SURVIVAL2 |
a SURVIVAL object |
Depending on the function a proportion surviving, hazard, cumulative hazard, inverse of the cumulative hazard, a random draw or a plot
sfx()
: Survival function
hfx()
: Hazard function
Cum_Hfx()
: Cumulative Hazard function
invCum_Hfx()
: Inverse of the Cumulative Hazard function
rsurv()
: Generate random values from the distribution
rsurvhr()
: Generate random values from the distribution under proportional hazard ratios
rsurvaft()
: Generate random values from the distribution under accelerated failure time ratios
rsurvah()
: Generate random values from the distribution under accelerated hazard ratios
plot_survival()
: Plot of the survival functions
ggplot_survival_random()
: ggplot of the simulation of survival times
compare_survival()
: Compare graphically two survival distributions
ggplot_survival_hr()
: ggplot of the simulation of survival times with hazard ratios
ggplot_survival_aft()
: ggplot of the simulation of survival times with accelerated time failures
ggplot_survival_ah()
: ggplot of the simulation of survival times with accelerated hazard
#' # Define a SURVIVAL object obj <- s_factory(s_weibull, surv = 0.8, t = 2, shape = 1.2) # Survival, Hazard and Cumulative hazard at time 0.4 sfx(SURVIVAL = obj, t= 0.4) hfx(SURVIVAL = obj, t = 0.4) Cum_Hfx(SURVIVAL = obj, t = 0.4) # Time when the Cumulative hazard is 0.8 invCum_Hfx(SURVIVAL = obj, H = 0.8) # Draw one random survival time from the distribution rsurv(SURVIVAL = obj, n = 1) # Draw one random survival time from the distribution under Proportional # hazard, Accelerated time failure or Accelerated hazard. rsurvhr(SURVIVAL = obj, hr = 0.5) rsurvaft(SURVIVAL = obj, aft = 2) rsurvah(SURVIVAL = obj, aft = 2, hr = 0.5) # Plot the survival functions plot_survival(SURVIVAL = obj, timeto = 2, main = "Example of Weibull distribution" )
#' # Define a SURVIVAL object obj <- s_factory(s_weibull, surv = 0.8, t = 2, shape = 1.2) # Survival, Hazard and Cumulative hazard at time 0.4 sfx(SURVIVAL = obj, t= 0.4) hfx(SURVIVAL = obj, t = 0.4) Cum_Hfx(SURVIVAL = obj, t = 0.4) # Time when the Cumulative hazard is 0.8 invCum_Hfx(SURVIVAL = obj, H = 0.8) # Draw one random survival time from the distribution rsurv(SURVIVAL = obj, n = 1) # Draw one random survival time from the distribution under Proportional # hazard, Accelerated time failure or Accelerated hazard. rsurvhr(SURVIVAL = obj, hr = 0.5) rsurvaft(SURVIVAL = obj, aft = 2) rsurvah(SURVIVAL = obj, aft = 2, hr = 0.5) # Plot the survival functions plot_survival(SURVIVAL = obj, timeto = 2, main = "Example of Weibull distribution" )