apc: Create an Age-Period-Cohort mortality model

In amvillegas/StMoMo: Stochastic Mortality Modelling

Description

Utility function to initialise a StMoMo object representing an Age-Period-Cohort mortality model.

Usage

apc(link = c("log", "logit"))

Arguments

link	defines the link function and random component associated with the mortality model. "log" would assume that deaths follow a Poisson distribution and use a log link while "logit" would assume that deaths follow a Binomial distribution and a logit link.

Details

The created model is either a log-Poisson or a logit-Binomial version of the classical age-period-cohort mortality model which has predictor structure

η_{xt} = α_x + κ_t + γ_{t-x}.

To ensure identifiability we follow Cairns et al. (2009) and impose constraints

∑_c γ_c = 0

and

∑_c cγ_c = 0

Value

An object of class "StMoMo".

References

Cairns, A. J. G., Blake, D., Dowd, K., Coughlan, G. D., Epstein, D., Ong, A., & Balevich, I. (2009). A quantitative comparison of stochastic mortality models using data from England and Wales and the United States. North American Actuarial Journal, 13(1), 1-35.

See Also

StMoMo, rh

Examples

APC <- apc()
wxt <- genWeightMat(EWMaleData$ages,  EWMaleData$years, clip = 3)
APCfit <- fit(APC, data = EWMaleData, wxt = wxt)
plot(APCfit, parametricbx = FALSE, nCol = 3)

amvillegas/StMoMo documentation built on Nov. 7, 2019, 5:39 a.m.