Description Usage Arguments Details Value References See Also Examples
View source: R/standardModels.R
Utility function to initialise a StMoMo
object representing the
M7 extension of the Cairns-Blake-Dowd mortality model introduced
in Cairns et al (2009).
1 |
link |
defines the link function and random component associated with
the mortality model. |
The created model is either a logit-Binomial or a log-Poisson version of the M7 model which has predictor structure
η_{xt} = κ_t^{(1)} + (x-\bar{x})κ_t^{(2)} + ((x-\bar{x})^2 - \hat{σ}^2_x)κ_t^{(2)} + γ_{t-x},
where \bar{x} is the average age in the data and \hat{σ}^2_x is the average value of (x-\bar{x})^2.
Identifiability of the model is accomplished by applying parameters constraints
∑_cγ_c = 0, ∑_c cγ_c = 0, ∑_c c^2γ_c = 0
which ensure that the cohort effect fluctuates around zero and has no linear or quadratic trend. These constraints are applied using the strategy discussed in Appendix A of Haberman and Renshaw (2011).
An object of class "StMoMo"
.
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.
Haberman, S., & Renshaw, A. (2011). A comparative study of parametric mortality projection models. Insurance: Mathematics and Economics, 48(1), 35-55.
StMoMo
, central2initial
,
cbd
, m6
, m8
1 2 3 4 | M7 <- m7()
wxt <- genWeightMat(55:89, EWMaleData$years, clip = 3)
M7fit <- fit(M7, data = central2initial(EWMaleData), ages.fit = 55:89)
plot(M7fit, parametricbx = FALSE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.