Statistical Modelling 13 (1) (2013), 41–67

Clustering in linear mixed models with approximate Dirichlet process mixtures using EM algorithm

Ian D Currie
Department of Actuarial Mathematics and Statistics,
and the maxwell Institute for Mathematical Sciences,
Herriot-Watt University,
Edimburgh
UK
e-mail: I.D.Currie@hw.ac.uk

Abstract:

We consider a generalized linear model (GLM) with canonical link function in which parameters can be subject to (a) a set of linear constraints and (b) smoothing. We apply Lagrange methods to give a general Newton-Raphson algorithm for such a GLM in which parameters are estimated, constraints are applied and smoothing is performed simultaneously. We express the Lee-Carter model, an important model for the forecasting of human mortality, in terms of GLMs, and use our method to estimate the parameters in the model. The smoothing option allows us to improve the forecasting properties of the model. We compare the performance of (a) the Poisson model with log link for the force of mortality and (b) the binomial model with logit link for the probability of death in a calendar year. Examples using UK Office for National Statistics data are provided.

Keywords:

constraints; generalised linear models; identifiability; Lagrange methods; Lee-Carter; mortality; Newton-Raphson; smoothing

Downloads:

Example data in zipped archive
back