MortCast (version 2.7-0)

cokannisto.estimate: Coherent Kannisto Estimation

Description

Estimate the coherent Kannisto parameters as described in Sevcikova et al. (2016).

Usage

cokannisto.estimate(mxM, mxF, ages, fitted = TRUE)

Arguments

mxM

A vector of male mortality rates.

mxF

A vector of female mortality rates.

ages

A vector of ages corresponding to mxM and mxF.

fitted

Logical. If TRUE the fitted values and residuals are returned.

Value

List of two lists, one for male and one for female. Each of the two lists contains the following components:

coefficients:

named vector with the coherent Kannisto coefficients \(c\) and \(d\). The \(d\) values are the same in both lists.

fitted.values:

the fitted values (not included if fitted is FALSE)

residuals:

input rates minus the fitted values (not included if fitted is FALSE)

Details

Given the Kannisto equation \(logit(m_x) = \log(c) + dx\), the Coherent Kannisto method estimates the \(d\) parameter jointly for male and female data, in order to prevent mortality crossovers in higher ages.

References

Sevcikova H., Li N., Kantorova V., Gerland P., Raftery A.E. (2016). Age-Specific Mortality and Fertility Rates for Probabilistic Population Projections. In: Schoen R. (eds) Dynamic Demographic Analysis. The Springer Series on Demographic Methods and Population Analysis, vol 39. Springer, Cham

See Also

cokannisto, kannisto.predict, kannisto

Examples

# NOT RUN {
data(mxM, mxF, package = "wpp2017")
country <- "Brazil"
mxm <- subset(mxM, name == country)[,"2010-2015"]
mxf <- subset(mxF, name == country)[,"2010-2015"]
cokannisto.estimate(mxm[18:21], mxf[18:21], ages = 18:21)

# }