lyl_2plot: Plot Life Years Lost at one specific age for two different...

View source: R/lyl_plot.R

lyl_2plotR Documentation

Plot Life Years Lost at one specific age for two different populations

Description

lyl_2plot was used to create a figure of Life Years Lost at one specific age for two different populations. Please use lyl_compare_plot instead.

Usage

lyl_2plot(
  x,
  y,
  color_alive = NA,
  colors = NA,
  labels = c("Population of interest", "Reference population"),
  ...
)

Arguments

x

An object of class lyl (obtained with function lyl).

y

An object of class lyl (obtained with function lyl).

color_alive

Color to be used for the censoring category. Default is NA, and default color is "white".

colors

Vector with one color for each cause of death. Default is NA, and default colors are used.

labels

Vector with labels for the two populations (default are "Population of interest" for x, and "Reference population" for y)

...

Additional arguments affecting the plot produced.

Value

A plot with survival function and stacked cause-specific cumulative incidences for two populations side by side.

References

  • Plana-Ripoll et al. lillies – An R package for the estimation of excess Life Years Lost among patients with a given disease or condition. PLoS ONE. 2020;15(3):e0228073.

See Also

  • lyl for estimation of Life Years Lost at one specific age.

  • lyl_diff to compare Life Years Lost for two populations.

Examples

# Load simulated data as example
data(simu_data)

# Estimate remaining life expectancy and Life Years
# Lost after age 45 years and before age 95 years
lyl_estimation <- lyl(data = simu_data, t = age_death, status = cause_death,
                      age_specific = 45, tau = 95)

# Same estimate for those with a specific disease
diseased <- simu_data[!is.na(simu_data$age_disease), ]

lyl_estimation1 <- lyl(data = diseased, t0 = age_disease,
                       t = age_death, status = cause_death,
                      age_specific = 45, tau = 95)

# Plot the data
lyl_compare_plot(list(lyl_estimation1, lyl_estimation))

lillies documentation built on Aug. 7, 2023, 5:09 p.m.

Related to lyl_2plot in lillies...