observed_survival: Calculate observed survival informations

In ampel-leipzig/ameld: Data and Model of End-Stage Liver Disease used in the AMPEL Project

Calculate observed survival informations

Description

This functions calculate some survival information (e.g. survival rate, mortality rate).

Usage

observed_survival(x, f, times)

observed_mortality(x, f, times)

observed_events(x, f, times, cumulative = FALSE)

Arguments

x	Surv, object
f	factor, strata/splitting factor
times	numeric, information at times
cumulative	logical, if TRUE (default FALSE) and multiple times are given the cumulative events are reported instead of the events in the time periods.

Details

In contrast to survival::summary.survfit() these functions do not remove empty strata (e.g. important for bootstrapping).

Value

double, with an element for each level in f, if times is a numeric of length 1; or a matrix with as much rows as elements in times and a column for each level in f.

Examples

sv <- Surv(1:4, c(0, 1, 1, 0))
f <- factor(c("a", "c", "a", "c"), levels = c("a", "b", "c"))
observed_survival(sv, times = 3)
observed_survival(sv, f = f, times = 3)
observed_survival(sv, f = f, times = 1:4)

observed_mortality(sv, f = f, times = 3)

observed_events(sv, f = f, times = 3)
observed_events(sv, f = f, times = 1:4)
observed_events(sv, f = f, times = 1:4, cumulative = TRUE)

ampel-leipzig/ameld documentation built on Aug. 23, 2024, 7:31 p.m.