lifetable: Construct lifetable from mortality rates

Description

Creates lifetable replicating the excel spreadsheet life table. a0 is set to be $0.07 + (1.7*m0)$ as per Keyfitz. $qx = mx/(1 + ((1-ax) * mx))$ as per Chiang (1984). lifetable constructs a lifetable for multiple years.

Usage

lifetable(data, series = names(data$rate)[1], years = data$year,
                 ages = data$age, max.age = min(100, max(ages)),
                 type = c("period", "cohort"))

Arguments

data

Demogdata object such as obtained from read.demogdata, forecast.fdm or forecast.lca.

series

Name of series to use. Default is the first series in data$rate.

years

Vector indicating which years to include in table. Default is all years available.

ages

Vector indicating which ages to include in table. Default is all ages available.

max.age

Age for last row. Ages beyond this are combined.

type

Type of lifetable: period or cohort.

Value

Object of class lifetable containing the following components:
labelName of region from which data are taken.
seriesName of series
ageAges for lifetable
yearPeriod years or cohort years
mxDeath rate at age x.
qxThe probability that an individual of exact age x will die before exact age x+1.
lxNumber of survivors to exact age x. The radix is 1.
dxThe number of deaths between exact ages x and x+1.
LxNumber of years lived between exact age x and exact age x+1.
TxNumber of years lived after exact age x.
exRemaining life expectancy at exact age x.

References

Chiang CL. (1984) The life table and its applications. Robert E Krieger Publishing Company: Malabar. Keyfitz, N, and Caswell, H. (2005) Applied mathematical demography, Springer-Verlag: New York.

Examples

Run this code

france.lt <- lifetable(fr.mort)
plot(france.lt)

france.LC <- lca(fr.mort)
france.fcast <- forecast(france.LC)
france.lt.f <- lifetable(france.fcast)
plot(france.lt.f,years=2010)