project_mortality_rates: Project mortality rates

In RobbenJ/MultiMoMo: Multi-population mortality models

Project mortality rates

Description

This function project mortality rates using the simulations that were set up for the joint time series \mjeqn(K_t^M, \kappa_t^M, K_t^F, \kappa_t^F)ASCII representation.

Usage

project_mortality_rates(fit_M, fit_F, proj_par)

Arguments

fit_M	The calibrated Li-Lee model for males.
fit_F	The calibrated Li-Lee model for females.
proj_par	The (simulations of) the projected time dependent parameters \mjeqn(K_t^M, \kappa_t^M, K_t^F, \kappa_t^F)ASCII representation.

Details

\loadmathjax

The arguments fit_M and fit_F should be the outputs of the function fit_li_lee.

The argument proj_par should be the output of the function project_parameters.

Value

A list containing the following objects

• The projected/simulated male mortality rates ($Male).

• The projected/simulated female mortality rates ($Female).

Examples

lst   <- MultiMoMo::european_mortality_data
dat_M <- lst$M
dat_F <- lst$F
xv    <- 0:90
yv = yvSPEC <- 1970:2018
Countries   <- names(dat_M$UNI)
CountrySPEC <- "BE"
fit_M <- fit_li_lee(xv, yv, yvSPEC, CountrySPEC, dat_M, "NR", TRUE, FALSE)
fit_F <- fit_li_lee(xv, yv, yvSPEC, CountrySPEC, dat_F, "NR", TRUE, FALSE)

arima_spec <- list(K.t_M = "RWD", k.t_M = "AR3.1", K.t_F = "RWD", k.t_F = "AR5.0")
n_ahead    <- 50
n_sim      <- 10000
est_method <- "PORT"
proj_par   <- project_parameters(fit_M, fit_F, n_ahead, n_sim, arima_spec, est_method)
proj_rates <- project_mortality_rates(fit_M, fit_F, proj_par)

RobbenJ/MultiMoMo documentation built on June 28, 2022, 9:29 p.m.