Part I – doi:10.1017/S1748499500000713, Published by Cambridge University Press, 10 May 2010
Part II – doi:10.1017/S1748499500000725, Published by Cambridge University Press, 10 May 2010
Part III – doi:10.1017/S1748499500000737, Published by Cambridge University Press, 10 May 2010
In a series of three papers, Ling et al. (Reference Ling, Waters and Wilkie2009a, Reference Ling, Waters and Wilkie2009b & Reference Ling, Waters and Wilkie2009c), the authors set out formulae for the calculation of the recovery and mortality intensities for Income Protection Insurance (IPI) claimants. These formulae, which allow for duration of sickness, age, cause of claim and other factors, are, unsurprisingly, complex. Regrettably, there are some typographical errors in the published formulae. This short note corrects these errors and provides some clarification regarding the calculation of the intensities. The errors in the original papers were all typographical and did not affect the numerical results or the conclusions presented in those papers.
The corrections/clarifications all relate to Papers I and II and are set out in Paper/Section order.
1. Paper I
1.1. Paragraph 4.3
The definition of xyear should be:
1.2. Table 4
There is no need to include αDP 1 as a retained covariate for the three models since DP1 is the baseline category for deferred period and hence is included by default in the three models.
αage 2 is not a retained covariate in the model for cs36. See Table 6.
1.3. Paragraph 4.5
The duration of sickness, z, is measured in years but the break points, for example , are expressed in days, with 365 days in each year. Also, in formulae (iii) and (iv), the i−th power is taken after taking the positive part. Hence, the correct illustrative formulae for (i), (iii) and (iv) should be:
(i)
(iii)
(iv)
Formula (ii) is unchanged.
For example, the terms to be included in the model for cs32, cs33, cs34, cs37 and cs38 should be:
where and .
As a further example, the term to be included in the model for cause cs27 is:
See Table 11 for the parameter values.
1.4. Paragraph 4.6
For the calculation of the baseline hazard, ρ 0(z), the duration of sickness, z, is measured in years.
1.5. Paragraph 4.7
Formula (3) should include within the second exponential a term:
since, where several causes are being modelled within the same sickness category, each cause may have its own cause parameter. This is the case for causes cs32, cs33, cs34 and cs37 in Category G6, with cs38 being the baseline cause for this group.
Formula (3) contains the terms and . Note that, as explained in Section 4.3:
is , where is the Chebycheff polynomial of order and , and,
is , where Ci() is the Chebycheff polynomial of order and .
1.6. Figure 6
The caption and legend for Figure 6 should specify that the recovery intensities for cs36 are for males. This is not required for cs38 since for this cause the recovery intensities do not depend on sex.
1.7. Paragraph 5.1
The coefficient −0.2885 in the second formula in point (iv) should be −0.2881. The final result, 0.7646, is unchanged.
The final two paragraphs of point (v) has some numerical errors; in particular, the denominator in the factor for Year should be 16 rather than 13. Hence these paragraphs should be:
For cs36 we can adjust the intensities so that they relate to other years by multiplying by a factor . For example, intensities for the calendar year 1975 would be a multiple 1.4513 and intensities for 2002 would be a multipe 0.6696 of those shown in Figure 6. This deterioration in recovery rates over the range of our data is consistent with the values in CMIR18 (2000, Table 4A) and CMIR22 (2005, Table 4A) for all causes combined and for deferred periods greater than one week. The recovery intensities for 2011, which is well beyond the range of our data, would be a multiple 0.5174 of those shown in Figure 6. Whether this continued deterioration is realistic is beyond the scope of this paper. A pragmatic approach for anyone requiring recovery intensities applicable to years later than 2002 might be to use the 2002 rates as produced by our models, but, given the deterioration in these rates over the period 1975 to 2002, caution would need to be exercised.
For cs38 the Year adjustment is a function of the duration of sickness z (in years), given by
Hence, for z = 26 weeks and for calendar year 2002, the recovery intensities for both DP1 and DP4 would be a multiple 0.6407 of those shown in Figure 6.
1.8. Table 11
The value of the parameter is zero so that the formula for is:
1.9. Table 14
For cause 50, the parameter has the value 43.5 days.
2. Paper II
2.1. Paragraph 4.3
The final model for the mortality of male assured lives presented in Paper II (formula (5)) was:
This is incorrect. The correct formula is:
The difference is that there are separate spline functions for the age-dependent coefficients of the year adjustment, labelled i = 1, 2, 3, 4. Each of these splines has six knots and the knots are placed as shown in Table 1.
2.2. Paragraph 5
The final model for the mortality of female assured lives presented in Paper II (formula (7)) was:
This is incorrect. The correct formula is:
The authors regret any inconvenience caused by the errors in the original papers and are grateful to Andrew Stott for identifying many of the errors.