Skip to contents

Sample Population Life Table Data

Source: R/hazard-plot.R
sample_life_table.Rd

Generates age-group-specific population survival curves using Gompertz mortality, matching the US population life table overlays used in tp.hp.dead.age_with_population_life_table.sas and tp.hp.dead.uslife.stratifed.sas.

Usage

sample_life_table(
  age_groups = c("<65", "65-80", "≥80"),
  age_mids = c(55, 72, 85),
  time_max = 10,
  n_points = 100
)

Arguments

age_groups

Character vector of age group labels. Default c("<65", "65-80", "\u226580").

age_mids

Numeric vector of representative ages (years) for each group, same length as age_groups. Default c(55, 72, 85).

time_max

Upper end of the time axis (years). Default 10.

n_points

Number of time points. Default 100.

Value

A data frame with columns time, survival, and group.

Details

SAS column mapping:

  • time ← prediction time grid (years)

  • survival ← SMATCHED (age-group-specific survivorship, 0–100 %)

  • group ← age group label (e.g. "<65", "65-80", "\u226580")

See also

hv_hazard(), sample_hazard_data()

Examples

# Default: three age groups (<65, 65-80, ≥80) using Gompertz mortality
lt <- sample_life_table(time_max = 10)
head(lt)
#>        time  survival group
#> 1 0.0000000 100.00000   <65
#> 2 0.1010101  99.94562   <65
#> 3 0.2020202  99.89081   <65
#> 4 0.3030303  99.83555   <65
#> 5 0.4040404  99.77985   <65
#> 6 0.5050505  99.72370   <65
nlevels(lt$group)    # 3 age groups
#> [1] 3
range(lt$survival)   # 0-100 % survivorship scale
#> [1]  33.88395 100.00000

# Custom strata — two age groups, 15-year follow-up
lt2 <- sample_life_table(
  age_groups = c("Under 70", "70 and over"),
  age_mids   = c(60, 78),
  time_max   = 15
)
levels(lt2$group)
#> [1] "Under 70"    "70 and over"