lillies (version 0.2.5)

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

Description

lyl_2plot creates a figure of Life Years Lost at one specific age for two different populations.

Usage

lyl_2plot(x, y, 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).

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.

See Also

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

• lyl_diff to compare Life Years Lost for two populations.

Examples

# NOT RUN {
# 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_2plot(lyl_estimation1, lyl_estimation)
lyl_2plot(lyl_estimation1, lyl_estimation,
          labels = c("Population with a disease", "General population"))

# The plot can be modified with a usual ggplot2 format
lyl_2plot(lyl_estimation1, lyl_estimation) +
  ggplot2::xlab("Age [in years]") +
  ggplot2::ggtitle("Differences in Life Years Lost at age 45 years")
# }