Abstract. Recent literature on mortality modeling suggests to include in the dynamics of mortality rates the effects of time, age, the interaction of these two and a term for possible shocks. In this paper we investigate models that use Legendre polynomials for the inclusion of age and cohort effects. In order to capture the dynamics of the shock term it is suggested to consider continuous autoregressive moving average (CARMA) models due to their flexibility in reproducing different autoregressive profiles of the shock term. In order to validate the proposed model, different life tables are considered. In particular the male life tables for New Zealand, Taiwan and Japan are used for the presentation of in-sample fitting. Empirical analysis suggests that the inclusion of more flexible models such as higher-order CARMA(p,q) models leads to better in-sample fitting.
With Asmerilda Hitaj & Edit Rroji
You may also like
YSS2019: Computational and Statistical Methods for Stochastic Process
The first YUIMA Summer School on Computational and Statistical Methods for Stochastic Process This 4 days course aims at introducing researchers, PhD students and practitioners to several aspects of numerical and statistical analysis of time series…
Simulation and Inference for Stochastic Processes with YUIMA
Contains both theory and code with step-by-step examples and figures. Uses YUIMA package to implement the latest techniques available in the literature of inference for stochastic processes. Shows how to create the description of very abstract models in…
The YUIMA Project: A Computational Framework for Simulation and Inference of Stochastic Differential Equations
Abstract. The YUIMA Project is an open source and collaborative effort aimed at developing the R package yuima for simulation and inference of stochastic differential equations. In the yuima package stochastic differential equations can be of very…
On CARMA estimation in YUIMA
Q: Can YUIMA manage with statistical inference the following model? # yuima.model object model <- setModel(drift = c("-3*x1","x1*x2"), diffusion = matrix(c("1","0"),2,1), state.variable = c("x1","x2"),…
R Code for calculating standard errors(p-values) in COGARCH
To determine p-values and or standard errors for estimates, we combine the estimation method for COGARCH in yuima with the bootstrap methodology. Here you can find an example of the code tha we used in our work "Discrete‐Time Approximation of a…