MortCast (version 2.7-0)

mlt: Model Life Tables Mortality Patterns

Description

Predict age-specific mortality rates using Coale-Demeny and UN model life tables.

Usage

mlt(e0, sex = c("male", "female"), type = "CD_West", nx = 5, ...)
mltj(e0m, e0f, ..., nx = 5)

Arguments

e0

A time series of target life expectancy.

sex

Either "male" or "female".

type

Type of the model life table. Available options are “CD_East”, “CD_North”, “CD_South”, “CD_West”, “UN_Chilean”, “UN_Far_Eastern”, “UN_General”, “UN_Latin_American”, “UN_South_Asian”.

nx

Size of age groups. Should be either 5 or 1.

…

Additional arguments passed to the underlying function.

e0m

A time series of target male life expectancy.

e0f

A time series of target female life expectancy.

Value

Function mlt returns a matrix with the predicted mortality rates. Columns correspond to the values in the e0 vector and rows correspond to age groups. Function mltj returns a list of such matrices, one for each sex.

Details

Given a level of life expectancy (e0), sex and a type of model life table, the function extracts the corresponding mortality pattern from MLTlookup (for abridged LT) or MLT1Ylookup (for 1-year LT), while interpolating between neighboring e0 groups. Function mlt is for one sex, while mltj can be used for both sexes.

References

https://www.un.org/development/desa/pd/data/extended-model-life-tables

Coale, A., P. Demeny, and B. Vaughn. 1983. Regional model life tables and stable populations. 2nd ed. New York: Academic Press.

See Also

mortcast, mortcast.blend, pmd, MLTlookup

Examples

# NOT RUN {
data(e0Fproj, package = "wpp2017")
country <- "Uganda"
# get target e0
e0f <- subset(e0Fproj, name == country)[-(1:2)]
# project into future using life table Cole-Demeny North
mx <- mlt(e0f, sex = "female", type = "CD_North")
# plot first projection in black and the remaining ones in grey 
plot(mx[,1], type = "l", log = "y", ylim = range(mx),
    ylab = "female mx", xlab = "Age", 
    main = paste(country, "5-year age groups"))
for(i in 2:ncol(mx)) lines(mx[,i], col = "grey")

# MLT for 1-year age groups
mx1y <- mlt(e0f, sex = "female", type = "CD_North", nx = 1)
plot(mx1y[,1], type = "l", log = "y", ylim = range(mx1y),
    ylab = "female mx", xlab = "Age", 
    main = paste(country, "1-year age groups"))
for(i in 2:ncol(mx1y)) lines(mx1y[,i], col = "grey")
    
# }