A likely story

The foundation for most modern statistical inference is the log-likelihood function.  By maximising the value of this function, we find the maximum-likelihood estimate (MLE) for a given parameter, i.e. the most likely value given the model and data.  For models with more than one parameter, we find the set of values which jointly maximise the log-likelihood.

This much is basic statistics.  However, the log-likelihood function can give you more insight than just yielding MLEs.  In particular the shape and curvature of the log-likelihood tells you how much confidence you can have in a particular MLE.  By way of example, consider fitting a simple Makeham model for the force of mortality, μx:

μx = exp(ε) + exp(α + βx)

for which we can find the triplet of values for α, β and ε which jointly maximises the log-likelihood, and therefore provides the joint MLEs for the three parameters. We can also plot the log-likelihood for values around each parameter, which we do below for ε:

graph: MakehamLogL1

which shows that while the MLE for ε is -5.94, there is support for values in the interval (-6.40, -5.62) given the MLEs for α and β.  This suggests that the Makeham parameter for a constant background level of mortality is a valid model for this data set.  One of the handy things about support intervals is that they don't involve any further assumptions and don't have to be symmetric.  In the above example we can see that support on the right falls away faster than on the left.

The classical approach to statistical inference is to assume that the parameter estimate has a Normal distribution around the true value with a standard error of σ (say).  Thus, we would be 95% confident that the true value of ε in our example would lie within ±2σ of the estimated value.  The classical approach assumes a symmetric distribution around the true parameter value, and the usual means of doing this is to approximate the log-likelihood by a quadratic function.  We can see in the above graph that the log-likelihood is actually mildly skewed, but a quadratic approximation will be reasonable for the purpose of getting a rough standard error.

The problem with the classical approach is when the log-likelihood is neither symmetric nor approximately quadratic.  To illustrate this, consider the same profile log-likelihood for ε in the same Makeham model, but this time for a different data set where the Makeham parameter is not significant:

graph: MakehamLogL2

In this case the classical approach leads to a value for σ of 84,432, which is wildly larger than the likely values ε can take. This is caused by the log-likelihood being not remotely quadratic in shape, and certainly not symmetric.  The support-interval approach is more revealing, since the support interval for ε is actually (-infinity, -7.34).  This reflects the fact that the Makeham parameter is not significant, since exp(-infinity) is zero.

Written by: Stephen Richards
Publication Date:
Last Updated:

Log-likelihoods in Longevitas

Longevitas fits all its models using the log-likelihood function.  Classical standard errors are produced for all estimates, but the expert system will warn you if the log-likelihood has a shape which invalidates them.

Previous posts

Do we need standard tables any more?

Actuaries are long used to using standard tables. In the UK these are created by the Continuous Mortality Investigation Bureau (CMIB), and the use of certain tables is often prescribed in legislation. As actuaries increasingly move to using statistical models for mortality, it is perhaps natural that they should first consider incorporating standard tables into these models. But are standard tables necessary, or even useful, in such a context?

Tags: Filter information matrix by tag: standard table, Filter information matrix by tag: GLM, Filter information matrix by tag: survival models

New SAPS tables

On 31st October 2008 the CMIB published the new SAPS S1 mortality tables based on the mortality of defined-benefit pension schemes. These tables follow the previous release of the '00 Series' mortality tables based on the mortality of life-office pensioners. 
Tags: Filter information matrix by tag: standard table

Add new comment

Restricted HTML

  • Allowed HTML tags: <a href hreflang> <em> <strong> <cite> <blockquote cite> <code> <ul type> <ol start type> <li> <dl> <dt> <dd> <h2 id> <h3 id> <h4 id> <h5 id> <h6 id>
  • Lines and paragraphs break automatically.
  • Web page addresses and email addresses turn into links automatically.