codhump is used to fit a constrained non-parametric estimation of the cause- and age-specific contributions to the young adult mortality hump.
list produced with
list containing the cause-of-death typology to use (see details)
age range to consider for the analysis
maximum number of iterations of the optimization algorithm
assumed end of the hump used to estimate the starting values (see details)
logical value to automatically correct incoherent starting values
smoothing parameter for the senesence component
The estimation uses simultaneous constrained splines to estimate a Sum of Smooth Exponentials (SSE) model on cause-deleted forces of mortality. Briefly, the SSE model describes the observed force of mortality over age mu as the sum of three vectors mu1, mu2, mu3 over age. In other words, it assumes that deaths are realizations from Poisson distributions with mean composed of three parts: infant, early-adulthood and old-age mortality. For more information on the SSE model, see Camarda et al. (2016) and sse.fit. For the purpose of the study of the young adult mortality hump, the SSE model is here reduced to two components capturing the hump and the senescence parts of the force of mortality.
In order to decompose the hump by cause of death, this model uses a constrained approach on cause-deleted forces of mortality that can be summarized in four steps.
Identify manually those causes of death that contribute to the young adult hump component
Estimate an SSE model on the overall mortality in order to separate the senescent and young adult hump components
Construct cause-deleted datasets by removing separately deaths from each cause that was identified in step 1.
Simultaneously estimate SSE models for each of these cause-deleted datasets, interpreting the diminution of each component as the contribution of this cause to that component, and constraining the sum of all these contributions to be equal to the components estimated in step 2.
The SSE model on which this algorithm is based is more adaptive to specific mortality schedules than parametric models such as the Heligman-Pollard. It is thus designed to converge to meaningful results in the majority of cases. It however sometimes needs some fine tuning in order to reach coherent results. Please pay special attention to the following parameters in order to maximize the chance to get a meaningful result.
typ argument defines the typology of causes of death that are assumed to contribute to the young adult mortality hump.
Each element of the list is a vector containing one or several numerical values corresponding to the columns of the
mxc data frame from the
If an element of
typ contains only one cause, this cause will be considered on its own. If several causes are mentioned in the same element of the
typ object, these causes will be grouped and considered together in the analysis.
The names of the elements of the
typ object are recycled as the names of the causes of death in the typology.
The choice of the causes of death included in the
typ argument has a profound influence on the results, and should therefore be made with caution.
As each case calls for a specific selection of causes of death, there is no general rule as for which causes should be assumed to contribute to the hump. However, the model assumes that the list accounts for all of the hump.
A good way to test this assumption is to plot the force of mortality after removing all of the deaths from the causes included in the list and check that this leftover category does not display a hump.
Failure to include sufficient causes will result in a probable underestimation of the hump, as well as a difficulty for the algorithm to converge since the all-cause hump will not be totally accounted for by the cause-specific contributions to the hump.
On the other hand, specifying too many causes as contributing to the hump may result in difficulties to estimate cause-specific contributions that are not based on sufficient empirical evidence.
typ object will include somewhere between 2 and 6 causes (or groups of causes) of death depending on the context and the number of causes available in the dataset.
x.range argument defines the age range on which the analysis takes place. By default, the bottom boundary corresponds to the age with the lowest observed death rate, which is usually situated around 10 years of age.
The top of the range should include most of the adult years, but should also avoid the so-called mortality plateau often observed after age 100. It is defined by default to 90, but the user may want to change these values if necessary.
maxit argument defines the maximum number of iterations used for the step 4 of the algorithm. This step uses a Penalized Composite Link Model (PCLM) along with an iterative re-weighted least squares (IRWLS) optimization procedure.
maxit argument will therefore determine the maximum number of iterations used in the IRWLS loop, but also the speed of convergence since the re-weighting defines updated solution as
new = old * (1 - it/maxit) + new * (it/maxit)
maxit argument is defined at 200 by default, but it can be increased in case of an absence of convergence, even if the algorithm stopped before reaching
maxit number of iterations.
x.hump argument is used to compute the starting values of the independant SSE models. More specifically, it is used to determined the apriori age range to be considered for the young adult mortality hump.
In most cases, especially when the hump is obvious, the model converges even with relatively bad starting values. However, in case of a very narrow hump, it is advised to use a small value of
x.hump (20 to 25, but not lower), and in case of a wide hump it may be useful to consider larger values for
x.hump (30 and over).
Inadequate values of
x.hump may result in incoherent starting values for the cause-deleted SSE models and a lack of convergence in the step 4 of the algorithm (see
lambda parameter controls the amount of smoothing imposed on the senescence component of the SSE model. Typically, a large value is advisable since this is the part of the force of mortality that the model aims at removing in order to reveal the hump.
However, in some cases, especially when using abridged datasets, it may be useful to consider smaller values of
lambda such as 10 or 100. A bad choice of
lambda may result in poor starting values for the SSE and a lack of convergence in the step 4 of the algorithm.
start.correct argument is conceived as a safeguard against misspecified starting values of the SSE model. Specifically, while estimating the SSE model on the cause-deleted forces of mortality, if the hump component peaks after the middle of the
x.range interval, the starting values are replaced with the all-cause components.
This parameter is designed as an attempt to salvage a bad choice in other arguments, especially
lambda, but remains an emergency safeguard. In case of a lack of convergence, it is advised to change the values of the other parameters instead on relying on the
start.correct argument to guarantee convergence.
Returns an object of class codhump that includes
Data frame containing for each age the overall mortality rates, the cause-deleted rates and the rates of causes that do not contribute to the hump.
Fit of the
sse model on the all-cause mortality.
Fit of the
sse model for each cause-deleted mortality rates before constraining.
Fit of the
sse model for each cause-deleted mortality rates after constraining.
Age- and cause-specific contributions to the hump.
Percents of negative contributions after each iteration.
Maximum relative change among all coefficients of the components of the
sse model after each iteration.
List of parameters provided to fit the model.
Camarda, C. G., Eilers, P. H. C., & Gampe, J. (2016). Sums of smooth exponentials to decompose complex series of counts. Statistical Modelling.
Remund, A., Riffe, T., & Camarda, C. (2018). A cause-of-death decomposition of young adult excess mortality. Demography.
sse.fit, summary.codhump, plot.codhump, HCD2MH
1 2 3 4 5 6 7 8
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.