StMoMo (version 0.4.1)

apc: Create an Age-Period-Cohort mortality model

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.

Value

An object of class "StMoMo".

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 $$\eta_{xt} = \alpha_x + \kappa_t + \gamma_{t-x}.$$

To ensure identifiability we follow Cairns et al. (2009) and impose constraints $$\sum_c \gamma_c = 0$$ and $$\sum_c c\gamma_c = 0$$

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

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

# }