`codhump`

is used to fit a constrained non-parametric estimation of the cause- and age-specific contributions to the young adult mortality hump.

1 2 |

`data` |
list produced with |

`typ` |
list containing the cause-of-death typology to use (see details) |

`x.range` |
age range to consider for the analysis |

`maxit` |
maximum number of iterations of the optimization algorithm |

`x.hump` |
assumed end of the hump used to estimate the starting values (see details) |

`start.correct` |
logical value to automatically correct incoherent starting values |

`lambda` |
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.

The `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 `data`

object.
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.
Typically, the `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.

The `x.range`

argument defines the age range on which the analysis takes place. The bottom boundary should correspond approximately to the lowest observed age-specific 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. Therefore, the `x.range`

argument is defined by default as `10:90`

, but the user may want to change these values if necessary.

The `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.
The `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)*

The `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.

The `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 case of a very narrow hump, it is advised to reduce the value of `x.hump`

to 25 or even 20 (but not lower), wherease in case of an especially wide hump it may be useful to consider larger values for `x.hump`

.
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 `start.correct`

).

The `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 on which the model is not focused.
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.

The `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 `typ`

, `maxit`

, `x.hump`

and `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

- mxc
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.

- all
Fit of the

`sse`

model on the all-cause mortality.- start
Fit of the

`sse`

model for each cause-deleted mortality rates before constraining.- fits
Fit of the

`sse`

model for each cause-deleted mortality rates after constraining.- decomp
Age- and cause-specific contributions to the hump.

- neg
Percents of negative contributions after each iteration.

- dif
Maximum relative change among all coefficients of the components of the

`sse`

model after each iteration.- par
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. (2016). A Non-Parametric Approach to Decompose the Young Adult Mortality Hump by Causes of Death. Poster presented at the 2016 annual meeting of the Population Association of America. March 2016, Washington DC.

sse.fit, summary.codhump, plot.codhump, HCD2MH

1 2 3 4 5 6 7 8 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.