mT.fitExtrapolationLaw {MortalityTables}

Fit interpolation law to a mortality table and extrapolate

Description

Fit an extrapolation law (from the MortalityLaws Package to the base table of the mortality table and use it for extrapolation.

Usage

mT.fitExtrapolationLaw(
  table,
  method = "LF2",
  law = "HP",
  fit = 75:99,
  extrapolate = 80:120,
  fadeIn = 80:90,
  fadeOut = NULL,
  raw = NULL
)

Arguments

table

A life table object (instance of a mortalityTable class) or a list, table or array of mortalityTable objects

method

The fitting method (passed on to [MortalityLaw])

law

The mortality law fitted to the data(passed on to [MortalityLaw])

fit

Age range to use for the fit

extrapolate

Desired age range of the extrapolation (i.e. only those ages will be extrapolated and added to the base table)

fadeIn

age range to linearly fade in from the existing base table's values to the extrapolated

fadeOut

age range to linearly fade out from the extrapolated base table's values to the existing

raw

(optional) raw data to use for fitting. If not given, the raw probabilities of the table (stored in table@data$rawProbs) or the table's base table (table@deathProbs) is used by default.

Details

The fit is done using the MortalityLaws::MortalityLaw function, with the ages, death counts, exposures and death rates taken from the table mortality table object. The law and the fitting method can be given in the mT.fitExtrapolationLaw with the law and the fitting method

The age range fit is used to fit the law, while extrapolation is applied only to ages given in parameter extrapolate. As fitting does usually not result a smooth transition, a linear fade in or fade out range can also be provided.

Examples

mortalityTables.load("Austria_Census")
# use Austrian population mortalities for ages 18-95 and exponentially
# extrapolate them to lower ages
AT11.lowAges = mT.fitExtrapolationLaw(mort.AT.census.2011.male, law = "opperman",
                                      fit = 5:15, extrapolate = 0:15,
                                      fadeIn = NULL, fadeOut = 5:15)
AT11.oldAges = mT.fitExtrapolationLaw(mort.AT.census.2011.male, law = "HP",
                                      fit = 75:90, extrapolate = 75:120)
plotMortalityTables(mT.setName(AT11.lowAges, "Low ages fitted (ages 5-15 used)"),
                    mT.setName(AT11.oldAges, "old ages fitted (ages 75-90 used)"),
                    mort.AT.census.2011.male)