ct42005 2009
ct42005 2009
ct42005 2009
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 6 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (103) A2005 Institute of Actuaries
1 (i) Define each of the following examples of a stochastic process
(a) a symmetric simple random walk
(b) a compound Poisson process
[2]
(ii) For each of the processes in (i), classify it as a stochastic process according to
its state space and the time that it operates on. [2]
[Total 4]
2 You have been commissioned to develop a model to project the assets and liabilities
of an insurer after one year. This has been requested following a change in the
regulatory capital requirement. Sufficient capital must now be held such that there is
less than a 0.5% chance of liabilities exceeding assets after one year.
The company does not have any existing stochastic models, but estimates have been
made in the planning process of worst case scenarios.
Set out the steps you would take in the development of the model. [6]
3 Let Y1, Y3, Y5, , be a sequence of independent and identically distributed random
variables with
1
P Y2k 1 = 1 = P Y2k 1 = 1 = , k = 0, 1, 2,...
2
Hint: You may use the fact that, if X, Y are two variables that take only two
values and E XY E X E (Y ), then X, Y are independent. [4]
(b) Hence show that these probabilities do not depend on the current state
and that they satisfy the Chapman-Kolmogorov equations.
[3]
[Total 8]
the mortality rate for any individual is equivalent to a transition rate of 0.025 per
annum, independent of marital status
The state space of the process consists of five states: Never Married (NM),
Married (M), Widowed (W), Divorced (DIV) and Dead (D).
Px is the probability that a person currently in state x, and who has never previously
been widowed, will die without ever being widowed.
(ii) Show, by general reasoning or otherwise, that PNM equals PDIV . [1]
1 4
PNM PM
5 5
1 1
PM PDIV
4 2
[2]
(iv) Calculate the probability of never being widowed if currently in state NM. [2]
(v) Suggest two ways in which the model could be made more realistic. [1]
[Total 8]
Level 1: 0% discount
Level 2: 25% discount
Level 3: 40% discount
Level 4: 60% discount
Following a year with no claims, move to the next higher level, or remain at
level 4.
Following a year with one claim, move to the next lower level, or remain at
level 1.
Following a year with two or more claims, move back two levels, or move to
level 1 (from level 2) or remain at level 1.
For a given policyholder the probability of no claims in a given year is 0.85 and the
probability of making one claim is 0.12.
(ii) Calculate the probability that a policyholder who is currently at level 2 will be
at level 2 after:
The probability of the machine breaking down has been estimated to follow an
exponential distribution with the following annualised frequencies, :
The washing machine has never broken down at the start of the year (time t = 0).
Pi(t) is the probability that the machine has suffered i breakdowns by time t.
(i) Draw a transition diagram for the process defined by the number of
breakdowns occurring up to time t. [1]
t t
(b) demonstrate that P1 (t ) = e 10 e 5. [3]
(v) Calculate the expected number of claims under the policy. [3]
[Total 12]
END OF PAPER
CT4(103) A2005 5
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 7 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (104) A2005 Institute of Actuaries
1 (i) Write down the equation of the Cox proportional hazards model in which the
hazard function depends on duration t and a vector of covariates z. You
should define all the other terms that you use. [2]
(ii) Explain why the Cox model is sometimes described as semi-parametric . [1]
[Total 3]
qx
x t = ,
1 tq x
this implies that deaths between exact ages x and x + 1 are uniformly distributed. [4]
3 An investigation of mortality over the whole age range produced crude estimates of qx
for exact ages x from 2 years to 93 years inclusive. The actual deaths at each age
were compared with the number of deaths which would have been expected had the
mortality of the lives in the investigation been the same as English Life Table 15
(ELT15). 53 of the deviations were positive and 39 were negative.
Test whether the underlying mortality of the lives in the investigation is represented
by ELT15. [5]
4 A life insurance company has investigated the recent mortality experience of its male
term assurance policy holders by estimating the mortality rate at each age, qx . It is
proposed that the crude rates might be graduated by reference to a standard mortality
table for male permanent assurance policy holders with forces of mortality s 1 , so
x
2
that the forces of mortality 1 implied by the graduated rates q x are given by the
x
2
function:
s
1 = 1 k,
x x
2 2
where k is a constant.
(i) Describe how the suitability of the above function for graduating the crude
rates could be investigated. [2]
(ii) (a) Explain how the constant k can be estimated by weighted least squares.
(iii) Explain how the smoothness of the graduated rates is achieved. [1]
[Total 7]
The following table shows the Kaplan-Meier estimate of the survival function, based
on data from the 12 insects.
t (weeks) S(t)
0 t< 1 1.0000
1 t< 3 0.9167
3 t< 6 0.7130
6 t 0.4278
(i) Calculate the number of insects dying at durations 3 and 6 weeks. [6]
(ii) Calculate the number of insects whose history was censored. [1]
[Total 7]
x = total number of policies under which death claims are made when the
policyholder is aged x last birthday in each calendar year
Px(t) = number of in-force policies where the policyholder was aged x nearest
birthday on 1 January in year t
(ii) Obtain an expression, in terms of the Px(t), for the central exposed to risk, Exc ,
which corresponds to the claims data and which may be used to estimate the
force of mortality in year t at each age x, x . State any assumptions you
make. [4]
(iii) Comment on the effect on the estimation of the fact that the x relate to claims,
rather than deaths, and the Px (t ) relate to policies, not lives. [4]
[Total 9]
(i) (a) Estimate the constant force of mortality, 70, using a two-state model
and the data for the 8 lives in the table.
(ii) Show that the maximum likelihood estimate of the constant force, 70, using a
Poisson model of mortality is the same as the estimate using the two-state
model. [5]
(iii) Outline the differences between the two-state model and the Poisson model
when used to estimate transition rates. [3]
[Total 15]
END OF PAPER
EXAMINATION
April 2005
EXAMINERS REPORT
Faculty of Actuaries
Institute of Actuaries
Subject CT4 Models April 2005 Examiners report
EXAMINERS COMMENTS
Comments on solutions presented to individual questions for this April 2005 paper are given
below:
103 Part
Question A5 This was very well answered, with the majority of candidates scoring highly.
Question A6 Overall this was not well answered, but the better candidates did score well.
Many candidates produced good answers to part (i) to (iv). In part (iii), a
number of candidates did not verify that the boundary conditions were
satisfied.
Some candidates struggled with part (v) and a significant number did not
attempt this part of the question.
Page 2
Subject CT4 Models April 2005 Examiners report
104 Part
Question B3 This was well answered overall. Many candidates included a continuity
correction. This was not necessary, as there were 92 ages, but candidates who
did so received full credit if they used it correctly.
Page 3
Subject CT4 Models April 2005 Examiners report
103 Part
A1 (i) (a) Let Y1, Y2, , Yj, , be a sequence of independent and identically
distributed random variables with
1
P Yj 1 P Yj 1
2
and define
n
Xn Yj
j 1
Then X n n 1
constitutes a symmetric simple random walk.
Nt
Xt Yj , t 0.
j 1
(ii) (a) A simple random walk operates on discrete time and has a discrete
state space (the set of all integers, Z).
A2
Review the regulatory guidance.
Define the scope of the model, for example which factors need to be modelled
stochastically.
Plan the development of the model, including how the model will be tested and
validated.
Consider alternative forms of model, and decide and document the chosen
approach. Where appropriate, this may involve discussion with experts on the
underlying stochastic processes.
Page 4
Subject CT4 Models April 2005 Examiners report
Collect any data required, for example historic losses or policy data.
Existing worst case scenarios. Discuss with staff who made the estimates,
especially to gauge views on the probability of events occurring.
Debug the program, for example by checking the model behaves as expected for
simple, defined scenarios.
Validate data.
Run model on historic data to compare model s predictions with previous
observations.
Review parameters that have greatest effect on outputs.
Present range of capital requirements for differing parameter inputs.
A3 (i) It is clear that Y2k can only take two values, ±1, with probabilities
1
P Y2k 1 P Y2 k 1 Y2k 1 1 P Y2k 1 Y2k 1 1
2
and
P Y2 k 1
1
P Y2 k 1 1, Y2k 1 1 P Y2 k 1 1, Y2k 1 1
2
Page 5
Subject CT4 Models April 2005 Examiners report
E Y2 k E Y2 k 1 0.
Next,
E Y2 k Y2 k 1
1 1
E Y2 k Y2k 1 | Y2 k 1 1 E Y2k Y2k 1 | Y2 k 1 1
2 2
But
E Y2 k Y2 k 1 | Y2 k 1 1 1 1 0 ( 1) 1,
1 1
E Y2k Y2 k 1 1 1 0.
2 2
Since
E Y2 k E Y2 k 1 E (Y2 k Y2 k 1 )
1
P Y2 k 1 1| Y2 k 1
2
but
P Y2 k 1 1| Y2 k 1, Y2 k 1 1 0.
(iii) (a) Since the Yk are pairwise independent, we see that for all i, j, m, n,
1
pij (n) P Ym n j | Ym i .
2
Page 6
Subject CT4 Models April 2005 Examiners report
(b) The probabilities do not depend on the current state as they are all ½
A4 (i)
0.025
M W
0.1 0.1
0.1
NM 0.05
0.025 DIV
0.025
0.025
0.025
D
(ii) The transitions out of the divorced state are to the same states, and with the
same transition probabilities, as the transitions out of state NM.
Therefore the probability of ever reaching state W is the same from both
states.
0.025 0.1
PNM PD PM
0.125 0.125
0.025 0.05 0.025
PM PD PDIV PW
0.1 0.1 0.1
Page 7
Subject CT4 Models April 2005 Examiners report
0.025 0.1 1 4
PNM PM PM
0.125 0.125 5 5
0.025 0.05 1 1
PM PDIV PDIV
0.1 0.1 4 2
as required.
1 4 1 1
PNM PNM
5 5 4 2
2 2
1 PNM
5 5
2
PNM
3
(v)
Make mortality and marriage rates age dependent.
Divorce rate dependent on duration of marriage.
Divorce rate dependent on whether previously divorced.
Make mortality rate marital status-dependent.
A5 (i)(a) It is clear that X(t) is a Markov chain; knowing the present state, any
additional information about the past is irrelevant for predicting the next
transition.
0.15 0.85 0 0
0.15 0 0.85 0
P=
0.03 0.12 0 0.85
0 0.03 0.12 0.85
Page 8
Subject CT4 Models April 2005 Examiners report
Alternatively
Alternatively
It is also aperiodic;
If currently at either state 1 or 4, it can remain there. This is not true for states
2 and 3, however these are also aperiodic states since the chain may return e.g.
to state 2 after 2 or 3 transitions.
Page 9
Subject CT4 Models April 2005 Examiners report
This reads
4
Discard the first of these equations and use also that i 1 i
1 . Then, we
obtain first from (4) that 0.85 3 0.15 4 or, that 4 17 3 /3
Page 10
Subject CT4 Models April 2005 Examiners report
A6 (i)
1/10 1/5 1/4
No One Two Three
Breakdowns Breakdown Breakdowns Breakdowns
1
(ii) P0 (t ) P0 (t )
10
1 1
P1 (t ) P0 (t ) P1 (t )
10 5
1 1
P2 (t ) P1 (t ) P2 (t )
5 4
t
Using an integrating factor e 5 , we get
t t t
1 1
e5 P1' t P1 (t ) e 10 5
5 10
t t
d 5 1
e P1 (t ) e10
dt 10
t t
e 5 P1 (t ) e 10 const
t t
P1 (t ) e 10 const e 5
Page 11
Subject CT4 Models April 2005 Examiners report
t t
P1 (t ) exp 10 exp 5
Alternatively
Differentiate the suggested solution and verify it obeys the second equation.
t t
1 10 1 5 1
P2 (t ) exp exp * P2 (t )
5 5 4
t 3 1
d 1 t t
exp 4 P2 (t ) (exp 20 exp 20 )
dt 5
t 3 1
4 t t 8
exp 4 P2 (t ) exp 20 4 exp 20
3 3
t t t
4 10 5 4]
P2 (t ) [exp 3 exp 2 exp
3
1/10
P0 (1) exp 0.905
1/10 1/ 5
P1 (1) exp exp 0.0861
1 1 1
4 10 5 4]
P2 (1) [exp 3 exp 2 exp 0.00832896
3
Page 12
Subject CT4 Models April 2005 Examiners report
104 Part
B1 (i) If the hazard for life i is (t ; zi ) , then
l (t ; zi ) l 0 (t ) exp(b ziT ) ,
B2 Since
t
t px exp x s ds ,
0
t
t qx 1 t px 1 exp x s ds .
0
t
q x ds
t qx 1 exp
1 sq x
0
t
t qx 1 exp log(1 sqx ) 0
This is the assumption of a uniform distribution of deaths and implies that deaths
between exact ages x and x + 1 are uniformly distributed.
Page 13
Subject CT4 Models April 2005 Examiners report
B3 The null hypothesis is that the observed rates are a sample from a population in which
English Life Table 15 represents the true rates.
If the null hypothesis is true, then the observed number of positive deviations, P,
will be such that P ~ Binomial (92, ½).
We use the normal approximation to the Binomial distribution because we have > 20
ages
The z-score associated with the probability of getting 53 positive deviations if the null
hypothesis is true is, therefore
53 46 7
1.46 .
23 4.79
We use a two-tailed test, since both an excess of positive and an excess of negative
deviations are of interest.
s
B4 (i) The suitability of a linear relationship between 1 and 1 could be
x x
2 2
investigated by plotting log(1 q x ) against log(1 q xs ) or by plotting
s
1 against 1 and
x x
2 2
looking for a linear relationship.
If data are scarce, too close a fit is not to be expected, especially at extreme
ages.
Page 14
Subject CT4 Models April 2005 Examiners report
wx (qx qx ) 2
x
or
2
wx 1 1
x x
x 2 2
should be found (note that the summations are over all relevant ages x)
At each age there will be a different sample size or exposed to risk, Ex.
This will usually be largest at ages where many term assurances are
sold (e.g. ages 25 to 50 years) and smaller at other ages.
(b) The estimation procedure should pay more attention to ages where
there are lots of data. These ages should have a greater influence on
the choice of k than other ages.
1 1
wx or wx or wx = Ex
var qx var 1
x
2
(iii) The graduated forces of mortality are a linear function of the forces in the
standard table.
Since the forces in the standard table should already be smooth, a linear
function of them will also be smooth.
S (t ) (1 (t j )) ,
tj t
Page 15
Subject CT4 Models April 2005 Examiners report
d1 d
Since (1) , then we have 1 0.0833 ,
n1 n1
and, since all 12 insects are at risk of dying at tj = 1, we must therefore have
d1 = 1 and n1 = 12.
0.9167 0.7130 d3
so that (3) 0.222 .
0.9167 n3
Since we can have at most 11 insects in the risk set at tj = 3, we must have
d3 = 2 and n3 = 9.
Similarly, at tj = 6, we must have 0.4278 0.7130(1 (6)) ,
0.7130 0.4278 d6
so that (6) 0.400 .
0.7130 n6
Since we can have at most 7 insects in the risk set at tj = 6, we must have
d6 = 2 and n6 = 5.
Alternatively
t S(t) = (1- t) t nt dt ct
0 1.0000 0 12 0
1 0.9167 0.0833 12 1 2
3 0.7130 0.22 9 2 2
6 0.4278 0.4 5 2 3
5 7
Since we started with 12 insects, the remaining 7 insects histories were right-
censored.
Page 16
Subject CT4 Models April 2005 Examiners report
B6 (i) The principle of correspondence states that a life alive at time t should be
included in the exposure at age x at time t if and only if were that life to die
immediately, he or she would be counted in the deaths data x at age x.
(ii) Px(t) is the number of policies under observation aged x nearest birthday on
1 January in year t.
To correspond with the claims data, we wish to have policies classified by age
last birthday.
1
Px (t ) Px (t ) Px 1 (t ) .
2
1
E xc Px (t ) Px (t 1) ,
2
and, substituting for the Px (t ) in terms of Px(t) from the equation above
produces
1 1 1
E xc Px (t ) Px 1 (t ) Px (t 1) Px 1 (t 1) .
2 2 2
(iii) The principle of correspondence still holds, because we are dealing with
claims and policies: one policy can only lead to one claim.
However, because one life may have more than one policy it is possible that
two distinct death claims are the result of the death of the same life.
Page 17
Subject CT4 Models April 2005 Examiners report
i2 i
i
,
i i
i
d70
B7 (i)(a) The two-state estimate of is
, where v70 is the total time the members
70
v70
of the sample are under observation between exact ages 70 and 71 years.
v70 v70,i ,
i
where v70,i is the duration that sample member i is under observation between
exact ages 70 and 71 years.
where ENDDATE is the earliest of the date at which the observation of that
member ceases and the date of the member s 71st birthday,
and STARTDATE is the latest of the date at which observation of that
member begins and the date of the member s 70th birthday.
Page 18
Subject CT4 Models April 2005 Examiners report
(ii) The contributions to the Poisson likelihood made by each member are
proportional to the following
Member
1 exp(-0.25 70 )
2 exp(-0.75 70 )
3 70 exp(-0.5 70 )
4 70 exp(-0.25 70 )
5 exp(-0.25 70 )
6 exp(-0.3333 70 )
7 exp(-0.5833 70 )
8 exp(-0.25 70 )
2
L [exp( 3.167 70 )]( 70 ) .
Then
so that
d log L 2
3.167 .
d 70 70
Setting this equal to zero and solving for 70 produces the maximum
likelihood estimate,
which is 2/3.167 = 0.6316
d 2 log L 2
Since 2 2
, which is always negative, we definitely have a
d 70 70
maximum.
Page 19
Subject CT4 Models April 2005 Examiners report
(iii) The Poisson model is not an exact model, since it allows for a non-zero
probability of more than n deaths in a sample of size n.
The variance of the maximum likelihood estimator for the two-state model is
only available asymptotically, whereas that for the Poisson model is available
exactly in terms of the true .
Page 20
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 7 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (103) S2005 Institute of Actuaries
1 An insurance company has a block of in-force business under which policyholders
have been given options and investment-related guarantees. A stochastic model has
been developed which projects option and guarantee costs. You have used the model
to estimate, for the Company Board, the probability of the insurance company having
insufficient assets to honour the payouts under the policies. A Board member has
asked whether there are any factors which could cause this probability to be
inaccurate.
(ii) Stochastic process models can be placed in one of four categories according to
whether the state space is continuous or discrete, and whether the time set is
continuous or discrete. For each of the four categories:
The total time spent in State H is 652 years and the total time spent in State I is 44
years.
(i) Write down the likelihood function for these data. [3]
5 Claims arrive at an insurance company according to a Poisson process with rate per
week.
(i) Show that, given that there is exactly one claim in the time interval [t, t + s],
the time of the claim arrival is uniformly distributed on [t, t + s]. [3]
(ii) State the joint density of the holding times T0, T1, , Tn between successive
claims. [1]
(iii) Show that, given that there are n claims in the time interval [0, t], the number
of claims in the interval [0, s] for s < t is binomial with parameters n and s/t.
[3]
[Total 7]
ii = for 0 i N 1
i,i+1 = for 0 i N 1
ij =0 otherwise
(i) Write down the generator matrix and the Kolmogorov forward equations (in
component form) associated with this process. [3]
t ( t) j i
pij (t ) = e
( j i )!
(iii) Identify the distribution of the holding times associated with the jump process.
[2]
[Total 7]
7 A time-inhomogeneous Markov jump process has state space {A, B} and the
transition rate for switching between states equals 2t, regardless of the state currently
occupied, where t is time.
(i) Calculate the probability that the process remains in state A until at least
time s. [2]
(ii) Show that the probability that the process is in state B at time T, and that it is
T2
in the first visit to state B, is given by T 2 exp . [3]
(c) Calculate the time at which it is most likely that the process is in its
first visit to state B.
[6]
[Total 11]
END OF PAPER
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 6 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (104) S2005 Institute of Actuaries
1 Describe the advantages and disadvantages of graduating a set of observed mortality
rates using a parametric formula. [4]
Some students switch to another course. Others intend to sit the Survival Models
examination but simply stop attending lectures because they are so boring. In this
university, students who decide not to attend a lecture are not permitted to attend any
subsequent lectures.
The table below gives the number of students switching courses and stopping
attending lectures after each of the first 7 lectures of the course.
1 5 1
2 3 0
3 2 3
4 0 1
5 0 2
6 0 1
7 0 0
(i) Calculate the Index of Lecture Boringness for the Survival Models course. [4]
The deaths during the period of investigation, x, have been classified by age x at the
date of death, where
Censuses of the numbers alive on 1 January in each of the years 2002, 2003, 2004 and
2005 have been tabulated and denoted by
respectively, where x is the age last birthday at the date of each census.
(i) State the rate year implied by the classification of deaths, and give the ages of
the lives at the beginning of the rate year. [2]
(ii) Derive an expression for the exposed to risk in terms of the Px(t) (t = 2002,
2003, 2004, 2005) which corresponds to the deaths data and which may be
used to estimate the force of mortality, x+f at age x + f. [4]
(iii) Determine the value of f, stating any assumptions you make. [3]
[Total 9]
Age x dx qx Exposed-to-risk
18 6 0.0012 5,200
19 8 0.0013 5,000
20 12 0.0015 4,800
21 8 0.0017 5,000
22 9 0.0019 3,800
23 6 0.0020 3,600
24 8 0.0021 3,200
(i) Test whether the overall fit of the graduated rates to the crude data is
satisfactory using a chi-squared test. [5]
(iii) (a) Describe three possible shortcomings in a graduation which the chi-
squared test cannot detect, and
(b) State a test which can be used to detect each one. [3]
[Total 9]
5 An investigation was carried out into the effects of lifestyle factors on the mortality of
people aged between 50 and 65 years. The investigation took the form of a
prospective study following a sample of several hundred individuals from their 50th
birthdays until their 65th birthdays and collecting data on the following covariates for
each person:
In addition, data were collected on the age at death for persons who died during the
period of investigation.
(i) Explain why the Gompertz hazard might be appropriate for analysing the
mortality of persons aged between 50 and 65 years. [2]
B = exp( 0 + 1 X1 + 2 X2 + 3 X3),
(iii) Using the Gompertz hazard, the parameter estimates in the proportional
hazards model were as follows:
Sex 1 +0.40
Cigarette smoking 2 +0.75
Alcohol consumption 3 0.20
0 5.00
c +1.10
(a) Describe the characteristics of the person to whom the baseline hazard
applies in this model.
(b) Calculate the estimated hazard for a female cigarette smoker aged 55
years who does not consume alcohol.
(c) Show that, according to this model, a cigarette smoker at any age has a
risk of death roughly equal to that of a non-smoker aged eight years
older. [6]
[Total 11]
(i) Calculate the probability that a light bulb will fail within the first 20 days. [2]
END OF PAPER
EXAMINATION
September 2005
EXAMINERS REPORT
Faculty of Actuaries
Institute of Actuaries
Subject CT4 Models September 2005 Examiners Report
EXAMINERS COMMENTS
Comments on solutions presented to individual questions for this September 2005 paper are
given below:
103 Part
Question A2 This was well answered overall, even by the weaker candidates.
Credit was not given in part (ii)(b) if the examples cited were not likely to be
encountered by an actuary working in a professional capacity.
Question A4 This was very well answered overall, with the majority of candidates scoring
highly.
One common mistake was the omission of the constant term from the
likelihood function in part (i).
Question A5 This was very poorly answered by all but a few candidates.
Some candidates offered general explanations in parts (i) and (iii), which, if
clear enough, were given some credit.
Page 2
Subject CT4 Models September 2005 Examiners Report
104 Part
Question B4 This was very well answered, even by the weaker candidates.
The main areas where candidates lost marks were: not stating the null
hypothesis, or not stating it clearly enough; failure to identify the correct
degrees of freedom to be used in the test; and insufficient or insufficiently
clear descriptions of the shortcomings.
In part (iii), the majority of candidates seemed confused between two issues in
connection with bias. There are two distinct problems. Firstly, if the
consistent bias is only small, the chi-squared test may fail to detect it because
the resulting number (i.e. the sum of the squared deviations) is not large
enough to exceed the critical value. The signs test, which ignores the
magnitude of the bias and looks only at how consistent it is across the ages,
can be used to identify this. The second problem is that even if the consistent
bias is larger and the chi-squared test leads us to reject the null hypothesis,
the test gives no indication of whether the graduated rates are too high or too
low. This is because the deviations are squared and the test statistic always
positive. The signs test is not a solution to this second problem.
Page 3
Subject CT4 Models September 2005 Examiners Report
103 Part
A1
Items to be mentioned include:
Models will be chosen which it is felt give a reasonable reflection of the underlying
real world processes, but this may not turn out to be the case. (Model error.)
The model may be very sensitive to parameters chosen, and the parameters are
estimates because the true underlying parameters cannot be observed. (Parameter
error.)
Sampling error may result from running insufficient simulations. (It should be
possible to give a confidence interval for the error that could result from this source.)
The management actions assumed may not match what would happen in extreme
circumstances.
Policyholder behaviour, such as take-up rates for options, may differ in practice.
There may be future events, such as legislative changes which affect the
interpretation of the policy conditions, which have not been anticipated in the
modelling.
There may be errors in the coding of the model. The model is likely to be complex
and difficult to verify completely.
The model relies on input data, which may be grouped rather than being able to run
every policy. Any errors in the data could cause the output to be inaccurate.
Page 4
Subject CT4 Models September 2005 Examiners Report
A2
(i) (a) The state space is the set of values which it is possible for each random
variable Xt to take.
(b) The time set is the set J, the times at which the process contains a random
variable Xt.
(c) A sample path is a joint realisation of the variables Xt for all t in J, that is a set
of values for Xt (at each time in the time set) calculated using the previous
values for Xt in the sample path.
(a) Simple random walk, Markov chain, or any other suitable example
(b) Any reasonable example. For example: No Claims Discount systems, Credit
Rating at end of each year
(b) Any reasonable example. For example: Claims received by an insurer, Status
of pension scheme member
(b) Any reasonable example. For example: Share prices at end of each trading
day, Inflation index
(b) Any reasonable example. For example: Share prices during trading period,
Value of claims received by insurer
Page 5
Subject CT4 Models September 2005 Examiners Report
A3 (a) Given the current state (the largest outcome or the number of sixes) up to the
nth roll, no additional information is required to predict the status of the chain
after the next roll. Therefore both Bn and Cn have the Markov property.
(d) The chain Bn is clearly aperiodic; if currently at state i, it can remain there if
the next outcome is at most i.
It is not irreducible, as it cannot be reached from j for i < j.
(e) In the long run, Bn will reach state 6 and will remain there; hence in
equilibrium P(Bn = 6) = 1 for sufficiently large n.
Cn cannot decrease and has an infinite state space; therefore, it is certain that it
will escape to infinity with probability one.
Page 6
Subject CT4 Models September 2005 Examiners Report
23 15 3 5
L K exp( 652( )) exp( 44( ))
l 23
652
23
0 652
23
0.0353 p.a.
652
2
l 23
2 2
0
2 1 2
l
(iii) The variance of is 2
,
23
2
which we can estimate by .
23
Page 7
Subject CT4 Models September 2005 Examiners Report
A5 (i) Let Nt denote the number of claims up to time t. Since the Poisson process has
stationary increments, we may take t = 0, so that the required conditional
distribution is
P T0 y, N s 1
P T0 y | Ns 1
P Ns 1
P Ny 1, N s Ny 0
P Ns 1
But Ns Ny is independent of Ny
and has the same distribution as Ns y.
y (s y)
( ye )e y
s
,
se s
(ii) Since holding times are independent, each having an exponential distribution,
their joint density is
n t1 t2 ... tn
e 1 t ,t
1 2 ,...,tn 0.
P Ns k , Nt n
P Ns k | Nt n
P Nt n
P Ns k , Nt Ns n k
P Nt n
Using again that the Poisson process has stationary and independent
increments, and that the number of claims in an interval [0, t] is Poisson ( t),
we derive from above that
Page 8
Subject CT4 Models September 2005 Examiners Report
s
e ( s)k e (t s ) n k
(t s ) n k
k! (n k )!
P Ns k | Nt n t
e ( t )n
n!
t n k
e s (t s ) n k
n!
k !(n k )! t n n
e t
n! s k (t s ) n k
k !(n k )! tktn k
k n k
n s s
1
k t t
Page 9
Subject CT4 Models September 2005 Examiners Report
. .
. .
A ,
. .
0 0
piN (t ) pi , N 1 (t ).
(ii) Differentiating the function given in the question, we get first for i = j,
pii (t ) e t,
We can then check that the above satisfy the forward equations.
t
(iii) For i = j(<N), the solution in (ii) implies that pii (t ) e , so that the
distribution of the holding times T0 , T1 ,..., TN 1 is exponential with parameter
.
For i = N, this is obviously not true; once the chain reaches state N, it stays
there forever.
Page 10
Subject CT4 Models September 2005 Examiners Report
d
A7 (i) P (t ) 2t PAA (t )
dt AA
d
ln PAA (t ) 2t
dt
ln PAA ( s ) s2 constant
s2
Hence, PAA ( s ) exp
T
P(remains in A to time s )
0
P(transition to B in time s, s + ds)
P(remains in B to time T) ds
T
PAA ( s ) 2 s PBB ( s, T )ds
s 0
Using the result from part (i) and the similar result for PBB with boundary
condition PBB(s, s) = 1, this gives us:
T
s2 T 2 s2
e 2s e ds
s 0
T
T2
2s e ds
s 0
T2
e T2
Page 11
Subject CT4 Models September 2005 Examiners Report
Probability
Time
(b) Commentary:
d t2 t2 t2
e t2 2t e 2t 3 e
dt
Page 12
Subject CT4 Models September 2005 Examiners Report
104 Part
B1 Advantages:
The graduated rates will progress smoothly provided the number of parameters is
small.
Can fit the same formula to different experiences and compare parameter values to
highlight differences between them.
Disadvantages:
It can be hard to find a formula to fit well at all ages without having lots of
parameters.
Care is required when extrapolating: the fit is bound to be best at ages where we have
lots of data, and can often be poor at extreme ages.
Page 13
Subject CT4 Models September 2005 Examiners Report
Lecture nj dj cj j 1 j S(j)
j
This is because the students who switched courses were probably less
interested in the subject matter of Survival Models than those who remained
registered.
Therefore they would have been more likely, had they not switched courses,
to cease attending lectures than those who did not switch.
Page 14
Subject CT4 Models September 2005 Examiners Report
A person who dies will be aged x on the birthday in the calendar year of death,
which implies that he or she will be aged x next birthday on 1 January in the
calendar year of death.
Since 1 January is the start of the rate interval, the age range at the start is x
1 to x.
(ii) A census of those aged x next birthday on 1 January in each year would
correspond to the classification of deaths.
However, the number alive aged x next birthday on any date is equal to the
number alive aged x 1 last birthday.
Thus, using the trapezium rule, the correct exposed to risk at age x in year t is
given by
1
Px 1 (t ) Px (t 1) .
2
Over the three calendar years 2002, 2003 and 2004, we have, therefore,
exposed to risk =
1
Px 1 (2002) Px (2003)
2
1
Px 1 (2003) Px (2004)
2
1
Px 1 (2004) Px (2005) .
2
(iii) Assuming birthdays are uniformly distributed over the calendar year, the
average age at the start of the rate interval will be x ½.
Page 15
Subject CT4 Models September 2005 Examiners Report
B4 (i) The null hypothesis is that the observed data come from a population in which
the graduated rates are the true rates.
(d x Ex qx ) 2
.
x Ex qx
( Ex q x Ex qx ) 2
Age Exqx Ex qx ( Ex qx Ex q x ) 2
Ex qx
Since we have 7 ages, we compare this with the tabulated value at the 5%
level at, say, 4 degrees of freedom (since we lose 2 3 degrees for every
10 ages graduated graphically).
Since 4.4673 < 9.488 we have no evidence to reject the null hypothesis.
(ii) On the basis of the chi-squared test, the graphical graduation adheres to the
data satisfactorily.
Page 16
Subject CT4 Models September 2005 Examiners Report
There may be one or two large deviations at particular ages, balanced by lots
of small deviations (as in the example in part (i))
These can be detected by the individual standardised deviations test.
The graduated rates may be too high or too low over the whole of the age
range, but by an amount too small for the chi-squared test to detect.
The signs test or the cumulative deviations test will detect this.
which indicates that the rate of increase of the hazard with age is constant.
Empirically, this is often a reasonable assumption for middle ages and older
ages, which include the age range 50 65 years.
Therefore the ratio between the hazards for any two persons with different
characteristics does not depend on duration, and so the model is a proportional
hazards model.
a female,
non-smoker,
who drinks less than 21 units of alcohol per week.
Page 17
Subject CT4 Models September 2005 Examiners Report
(b) For a female cigarette smoker who does not consume alcohol we have
X1 = 0, X2 = 1, X3 = 0 and x = 5.
5 = exp( 0 + 1 .0 + 2 .1 + 3 .0) . c5
= exp( 5 + 0.75) 1.105
= 0.0230.
u = exp( 0 + 1 X1 + 3 X 3) . c u,
*v = exp( 0 + 1 X1 + 0.75 + 3 X 3) . c v .
which simplifies to
cu = exp(0.75) . cv,
so that
cu/cv = cu v = exp(0.75) = 2.117.
Therefore
u v = log(2.117)/log(1.1)
= 0.75/0.0953 = 7.87.
Page 18
Subject CT4 Models September 2005 Examiners Report
We have:
20 q0 1 20 p0 1 1 p0 .19 p1
1 (1 1 q0 ) exp( 19 )
1 0.95exp( 19 0.01)
1 0.95exp( 0.19)
1 0.95(0.82696)
which is 0.21439.
(ii) (a) The complete expectation of life of a one-day old light bulb, e1 is
given by
e1 t p1dt
0
0.01t
e dt
0
1 0.01t 1
e1 e 0 1
0.01 0 0.01
= 100 days.
1
e0 t p0 dt t p0 dt t p0 dt . (*)
0 0 1
Alternative 1
Page 19
Subject CT4 Models September 2005 Examiners Report
1
1 1 p0
2
1
1 (1 1 q0 )
2
1
(1 0.95) 0.975
2
1 p0 t p1dt 0.95(100)
0
(using the result from part (i) above).
Therefore:
e0 0.975 100 0.95 95.975 days.
Alternative 2
Then
1
1 p0 exp ds exp( )
0
1 1 q0 0.95.
So that
exp( ) 0.95
and
log(0.95) 0.0513.
Page 20
Subject CT4 Models September 2005 Examiners Report
1 1
t p0 dt exp( 0.0513t )dt
0 0
1 1
exp( 0.0513t ) 0
0.0513
1
exp( 0.0513) 1
0.0513
0.97478,
1 p0 t p1dt 0.95(100)
0
(using the result from part (i) above).
So that
(iii) The complete expectation of life of a light bulb at any age is an average of the
future lifetimes of all bulbs which have not failed before that age.
The value of e0 is lower than e1 because the average e0 includes the very
short lifetimes of the relatively large proportion of bulbs which fail in the first
day, which deflate the average, whereas e1 excludes these.
Page 21
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 6 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (103) A2006 Institute of Actuaries
A1 In the context of a stochastic process {Xt : t J}, explain the meaning of the
following conditions:
A2 A savings provider offers a regular premium pension contract, under which the
customer is able to cease paying in premiums and restart them at a later date. In order
to profit test the product, the provider set up the four-state Markov model shown in
the following diagram:
Policy
matured
AD (D) BD
t t
Premium t
AB Premiums
paying ceased/paid up
(A) (B)
BA
t
BC
AC t
t
Policy
surrendered
(C)
t p0AB t p0AA . AB
t t p0AB .( BA
t
BC
t
BD
t ) [5]
t
The insurer estimates that the probability of making at least one claim in a year is 0.1
if the policyholder made no claims the previous year, and 0.25 if they made a claim
the previous year.
(i) Explain why the system with state space {0%, 25%, 40%, 50%} does not form
a Markov chain. [2]
(ii) (a) Show how a Markov chain can be constructed by the introduction of
additional states.
(b) Write down the transition matrix for this expanded system, or draw its
transition diagram.
[4]
(ii) Describe the difference between a stochastic and a deterministic model. [2]
(iii) Outline the factors you would consider in deciding whether to use a stochastic
or deterministic model to study a problem. [3]
(iv) Explain how a deterministic model might be used to validate model outcomes
where a stochastic approach has been selected. [2]
[Total 9]
H M L
2 2
H 1
P M 1 2
L 2 2
1
(ii) Determine the range of values for for which the matrix P is a valid
transition matrix. [2]
(iv) For = 0.2, calculate the proportion of employees who, in the long run, are in
state L. [3]
(v) Given that = 0.2, calculate the probability that an employee s rating in the
third year, X3, is L:
(a) in the case that the employee s rating in the first year, X1, is H
(b) in the case X1 = M
(c) in the case X1 = L
[2]
[Total 11]
(ii) Outline the principal difficulties in fitting a Markov jump process model with
time-inhomogeneous rates. [2]
A company provides sick pay for a maximum period of six months to its employees
who are unable to work. The following three-state, time-inhomogeneous Markov
jump process has been chosen to model future sick pay costs for an individual:
(t)
Healthy Sick
(H) (S)
(t)
(t) (t)
Dead
(D)
Where Sick means unable to work and Healthy means fit to work.
The time dependence of the transition rates is to reflect increased mortality and
morbidity rates as an employee gets older. Time is expressed in years.
(iii) Write down Kolmorgorov s forward equations for this process, specifying the
appropriate transition matrix. [1]
(iv) (a) Given an employee is sick at time w < T, write down an expression for
the probability that he or she is sick throughout the period w < t < T.
(b) Given that a transition out of state H occurred at time w, state the
probability that the transition was into state S.
(v) Using the results of part (iv) or otherwise, derive an expression for the
probability that an employee is sick at time T and has been sick for less than 6
months, given that they were healthy at time < T - ½. Your expression
should be in terms of the transition rates and PHH ( , w) only. [3]
(a) (t) should also depend on the holding time in state S, and
(b) mortality rates can be ignored.
[3]
[Total 14]
END OF PAPER
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 6 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (104) A2006 Institute of Actuaries
B1 A Cox proportional hazards model was estimated to assess the effect on survival of a
person s sex and his or her self-esteem (measured on a three-point scale as low ,
medium or high ). The baseline category was males with low self-esteem.
Write down the equation of the model, using algebraic symbols to represent variables
and parameters and defining all the symbols that you use. [4]
B2 (i) (a) Explain why it is important to sub-divide data when carrying out
mortality investigations.
(b) Describe the problems that can arise with sub-dividing data.
[4]
(ii) List four factors which are often used to sub-divide life assurance data. [2]
[Total 6]
B3 (i) Assume that the force of mortality between consecutive integer ages, y and
y + 1, is constant and takes the value µy.
Let Tx be the future lifetime after age x ( x y ) and Sx(t) be the survival
function of Tx.
Show that:
(ii) An investigation was carried out into the mortality of male life office
policyholders. Each life was observed from his 50th birthday until the first of
three possible events occurred: his 55th birthday, his death, or the lapsing of
his policy. For those policyholders who died or allowed their policies to lapse,
the exact age at exit was recorded.
Using the result from part (i) or otherwise, describe how the data arising from
this investigation could be used to estimate:
(a) 50
(b) 5 q50
[3]
[Total 7]
the age last birthday of the policyholder when the policy was taken out;
the number of premiums the policyholder paid before the policy lapsed.
(ii) Derive an expression for the exposed-to-risk in terms of Px,t * , stating any
assumptions you make. [7]
The ages at which policyholders died or cancelled their policies were as follows:
Age in Age in
years and months years and months
60y 5m 60y 2m
61y 1m 60y 3m
62y 6m 60y 8m
63y 0m 61y 0m
63y 0m 61y 0m
63y 8m 61y 0m
64y 3m 61y 5m
62y 2m
62y 9m
63y 9m
64y 5m
(i) Explain which types of censoring are present in the investigation. [2]
(ii) Calculate the Nelson-Aalen estimate of the integrated hazard for these
policyholders. [5]
(iv) Estimate the probability that a policyholder will survive to age 65. [2]
[Total 11]
x
qx e
Ex qx qx
x Ex qx qx zx
Ex q x 1 q x
(i) Test the graduation for goodness of fit using the chi-squared test. [5]
(ii) (a) By inspection of the data, suggest one aspect of the graduated rates
where adherence to data seems inadequate.
(b) Explain why this may not be detected by the chi-squared test.
(c) Carry out one other test that may detect this deficiency.
[5]
(iii) Suggest how the graduation could be adjusted to correct the deficiency
identified. [2]
[Total 12]
END OF PAPER
EXAMINATION
April 2006
EXAMINERS REPORT
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
M Flaherty
Chairman of the Board of Examiners
June 2006
Faculty of Actuaries
Institute of Actuaries
Subject CT4 Models Core Technical April 2006 Examiners Report
Comments
Comments on solutions presented to individual questions for this April 2006 paper are given
below.
103 Part
104 Part|
Page 2
Subject CT4 Models Core Technical April 2006 Examiners Report
103 Solutions
This means that the statistical properties of the process remain unchanged over
time.
(b) Because strict stationarity is difficult to test fully in real life, we also use the
less stringent condition of weak stationarity.
Weak stationarity requires that the mean of the process, E[Xt] = m(t), is
constant and the covariance, E[(Xs - m(s)) (Xt - m(t))], depends only on the
time difference t s.
dt ptBB 1 dt ptBA dt pt
BC
dt ptBD
BB
Substituting for dt pt
dt ptBA BA
t .dt o(dt )
dt ptBC BC
t .dt o(dt )
dt ptBD BD
t .dt o(dt )
dt ptAB AB
t .dt o(dt )
Page 3
Subject CT4 Models Core Technical April 2006 Examiners Report
Where o(dt) covers the possibility of more than one transition in time dt and
lim o(dt )
0
dt 0 dt
Substituting in:
AB AA AB AB BA BC BD
t dt p0 t p0 . t .dt t p0 (1 t .dt t .dt t .dt ) o(dt )
lim AB AB
AB t dt p0 t p0 AA AB AB BA BC BD
t p0 t p0 . t t p0 ( t t t )
t dt 0 dt
A3 (i) This is not a Markov chain because it does not possess the Markov property,
that is transition probabilities do not depend only on the current state.
Specifically, if you are in the 25% discount level, the transition probability to
state 0% is 0.25 if a claim was made last year and 0.1 if the previous year was
claim free.
(ii) (a) Split the 25% and 40% discount states to include whether the previous
year was claim free.
0% discount
25%NC (no claim last year)
25%C (at least one claim last year)
40%NC (no claim last year)
40%C (at least one claim last year)
50%
Page 4
Subject CT4 Models Core Technical April 2006 Examiners Report
(b)
0.9
0.75
25% 40%
0.9
NC NC
0% 0.1
0.25 0.75 50%
0.75 0.9
0.1
New state
0% 25%C 25%NC 40%C 40%NC 50%
0% 0.25 0 0.75 0 0 0
25%C 0.25 0 0 0 0.75 0
25%NC 0.1 0 0 0 0.9 0
Old State
(iii) In theory, the insurer should just use 2 NCD states according to whether the
policyholder made a claim in the previous year. This is because the company
believes the claims frequency is the same for drivers who have not made a
claim for 1, 2, 3 years (i.e. it remains at 0.1 whether the driver has been
claims-free for 1 or 10 years).
General comments:
More factors should be taken into account (with a suitable example such
as how long the policyholder has been driving).
Page 5
Subject CT4 Models Core Technical April 2006 Examiners Report
A4 (i) Systems with long time frames such as the operation of a pension fund can be
studied in compressed time.
Different future policies or possible actions can be compared to see which best
suits the requirements or constraints of a user.
Often the output from a stochastic model is in the form of many simulated
possible outcomes of a process, so distributions can be studied.
(iii)
If the distribution of possible outcomes is required then stochastic
modelling would be needed, or if only interested in a single scenario then
deterministic.
Budget and time available stochastic modelling can be considerably
more expensive and time consuming.
Nature of existing models.
Audience for the results and the way they will be communicated.
Page 6
Subject CT4 Models Core Technical April 2006 Examiners Report
2
1 1 2
State H State M
2
1
2
State L
1 5 1 5
The solution of the quadratic is the interval , , so all
2 2 2 2
1
conditions are satisfied simultaneously for [0, ].
2
(iii) The chain is both irreducible, as every state can be reached from every other
state, and aperiodic, as the chain may remain at its current state for all H, M,
L.
Page 7
Subject CT4 Models Core Technical April 2006 Examiners Report
(iv) From the result in (iii), a stationary probability distribution exists and it is
unique. Let = ( H, M, L) denote the stationary distribution. Then, can be
determined by solving P = .
Discard the second of these equations and use also that the stationary
probabilities must also satisfy
H + M + L =1 (3)
The relevant entries are those in the last column, so that the answers are:
(a) 0.1008
(b) 0.28
(c) 0.6192.
Page 8
Subject CT4 Models Core Technical April 2006 Examiners Report
(ii) A model with time-inhomogeneous rates has more parameters, and there may
not be sufficient data available to estimate these parameters.
Also, the solution to Kolmogorov s equations may not be easy (or even
possible) to find analytically.
where
(t ) (t ) (t ) (t )
A(t ) (t ) (t ) (t ) (t )
0 0 0
T
(iv) (a) Pr(Waiting time T w Xw S ) exp ( (t ) (t ))dt
w
(b) Given there is a transition from state H at time w, the probabilities that
this is into state S or D are given by the relative transition rates at time
w.
( w)
So Probability into state S =
( w) ( w)
Page 9
Subject CT4 Models Core Technical April 2006 Examiners Report
T
Pr(Transition from state H at w) Pr(Transition toS) Pr(stays in S to time T) dW
T 1/ 2
T T
( w)
= PHH ( , w).( ( w) ( w)). .exp ( (t ) (t ))dt .dw
( w) ( w)
T 1/ 2 w
T T
= PHH ( , w). ( w).exp ( (t ) (t ))dt .dw
T 1/ 2 w
(vi) (a) This is likely to improve the predictive power of the model because:
However this would make the model more complicated to analyse, and
increase the volume of data required to fit parameters reliably.
(b) For individuals in employment mortality rates are likely to be low, and
may be ignorable. It is less likely that mortality out of state S could be
excluded.
Page 10
Subject CT4 Models Core Technical April 2006 Examiners Report
104 Solutions
B1 h(t ) h0 (t ) exp[ 1F 2M 3H ]
where
H is a variable taking the value 1 if the life has high self-esteem and 0 otherwise,
and
B2 (i) (a) The models of mortality we use assume that we can observe a group of
lives with the same mortality characteristics. This is not possible in
practice.
(b) Sub-dividing data using many factors can result in the numbers in each
class being too low.
Page 11
Subject CT4 Models Core Technical April 2006 Examiners Report
Sex
Age
Type of policy
Smoker/Non-smoker status
Level of underwriting
Duration in force
Sales channel
Policy size
Occupation (or social class) of policyholder
Known impairments
Geographical region
t
t py exp y s ds .
0
py exp y .
y 1 x px py . y x px .
Thus
y 1 x px
py ,
y x px
S x ( y 1 x)
py .
S x ( y x)
Page 12
Subject CT4 Models Core Technical April 2006 Examiners Report
Therefore
S x ( y 1 x)
exp( y) ,
S x ( y x)
so that
S x ( y x)
y log log[ S x ( y x)] log[ S x ( y 1 x)] .
S x ( y 1 x)
(ii) (a) Using the result from part (i) and putting x = 50, y = 50 gives
S50 (0)
50 log log[ S50 (1)]
S50 (1)
and, since
(ii) The age classification of the lapsing data is age last birthday on the policy
anniversary prior to lapsing .
This can be calculated by adding the policyholder s age last birthday when the
policy was taken to out to the number of annual premiums paid minus 1
(assuming that the first premium was paid at policy inception).
Define Px,t as the number of policies in force aged x last birthday at the
preceding policy anniversary at time t. This corresponds with the lapsing
data.
Page 13
Subject CT4 Models Core Technical April 2006 Examiners Report
1
E xc Px,t dt ,
0
which, assuming that policy anniversaries are uniformly distributed across the
calendar year,
may be approximated as
1
Exc [ Px,0 Px,1 ] .
2
But we do not observe Px,t directly. Instead we observe Px,t * the number of
policies in force at time t, classified by age last birthday at time t.
But the range of exact ages that could apply to a life aged x last birthday on
the policy anniversary prior to lapsing is (x, x + 2).
Assuming that birthdays are uniformly distributed across the policy year, half
of these lives will be aged x last birthday and half will be aged x+ 1 last
birthday.
Hence,
1
Px,t [ Px,t * Px 1,t
*
].
2
1 1 1
E xc [ Px,0* Px *
1,0 ] [ Px,1* Px *
1,1 ] .
2 2 2
policies might be taken out in large numbers just before the end of the tax
year,
Page 14
Subject CT4 Models Core Technical April 2006 Examiners Report
dj dj
tj nj dj cj nj
j
nj
(years)
0 t 5 100 0 2 0 0
12
5 t 1 112 98 1 4 1/98 0.0102
12
1 112 t 2 612 93 1 2 1/93 0.0210
2 612 t 3 90 1 1 1/90 0.0321
3 t 3 812 88 2 0 2/88 0.0548
3 812 t 4 312 86 1 1 1/86 0.0664
4 312 t 84 1 1 1/84 0.0783
(b)
0.09
Estimated Integrated Hazard
0.08
0.07
0.06
0.05
0.04
0.03
0.02
0.01
0
0 1 2 3 4 5
Duration since 60th birthday
Page 15
Subject CT4 Models Core Technical April 2006 Examiners Report
(iii) Either
Using the results of the calculation in (ii), the survival function can be
estimated by S t exp t .
Or
dj
Using the Kaplan-Meier estimate of S t 1 ,
tj t
nj
we get, for t 4 3/12:
1 1 1 2 1 1
S t 1 1 1 1 1 1
98 93 90 88 86 84
= 0.9243
B6 (i) The null hypothesis is that the crude rates come from a population in which
true underlying rates are the graduated rates.
2
The critical value of the 12 distribution at the 5% level is 21.03
Page 16
Subject CT4 Models Core Technical April 2006 Examiners Report
(ii) (a) The obvious problem with the graduation is one of overall bias. The
graduated rates are consistently too high, resulting in too many
negative deviations.
(b) This is not detected by the 2 test because the test statistic is the sum
of the squared deviations and so information on the sign and some
information on the size of the individual deviations is lost. The 2 test
would detect large bias, but in this case the graduated and crude rates
are close enough that the statistic is below the critical value.
15
1 15 15 15 15
2 0 1 2 3
15
1
1 15 105 455
2
= 0.0176
Ex qx Ex q x
x
Our test statistic is
Ex q x 1 q x
x
Page 17
Subject CT4 Models Core Technical April 2006 Examiners Report
Age
x
Ex qx qx Ex q x 1 q x
40 -6.40146 38.2751
41 -3.0025 42.84188
42 -7.92472 42.7289
43 -7.62982 36.46509
44 -6.08904 32.93758
45 2.63525 42.20447
46 4.2237 30.62388
47 -3.49218 51.31917
48 -4.9133 61.63457
49 -9.1832 51.99181
50 -7.488 58.25669
51 -8.24226 51.00139
52 1.10244 45.70533
53 -8.55647 72.14466
54 -7.87508 61.5123
Total -72.837 719.643
Ex qx Ex q x
x 72.837
2.715
719.643
Ex q x 1 q x
x
Comments:
(iii) The problem is that the graduated rates are too high. There doesn t appear to
be a problem with the overall shape.
Page 18
Subject CT4 Models Core Technical April 2006 Examiners Report
The problem persists across the whole age range, so the first adjustment to try
would be to decrease the value of .
Page 19
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 6 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (103) S2006 Institute of Actuaries
A1 A manufacturer uses a test rig to estimate the failure rate in a batch of electronic
components. The rig holds 100 components and is designed to detect when a
component fails, at which point it immediately replaces the component with another
from the same batch. The following are recorded for each of the n components used
in the test (i = 1,2, ,n):
The test rig was fully loaded and was run for two years continuously.
You should assume that the force of failure, , of a component is constant and
component failures are independent.
(i) Show that the contribution to the likelihood from component i is:
fi
exp ti si [2]
A2 The price of a stock can either take a value above a certain point (state A), or take a
value below that point (state B). Assume that the evolution of the stock price in time
can be modelled by a two-state Markov jump process with homogeneous transition
rates AB , BA .
(i) Write down the generator matrix of the Markov jump process. [1]
(ii) State the distribution of the holding time in each of states A and B. [1]
(iii) If 3, find the value of t such that the probability that no transition to state
B has occurred until time t is 0.2. [2]
(iv) Assuming all the information about the price of the stock is available for a
time interval [0,T], explain how the model parameters and can be
estimated from the available data. [2]
(v) State what you would test to determine whether the data support the
assumption of a two-state Markov jump process model for the stock price. [1]
[Total 7]
(ii) For each of the processes in (i), state whether it operates in continuous or
discrete time and whether it has a continuous or discrete state space. [2]
(iii) For each of the processes in (i), describe one practical situation in which an
actuary could use such a process to model a real world phenomenon. [3]
[Total 8]
A4 The credit-worthiness of debt issued by companies is assessed at the end of each year
by a credit rating agency. The ratings are A (the most credit-worthy), B and D (debt
defaulted). Historic evidence supports the view that the credit rating of a debt can be
modelled as a Markov chain with one-year transition matrix
(i) Determine the probability that a company rated A will never be rated B in the
future. [2]
(ii) (a) Calculate the second order transition probabilities of the Markov chain.
(b) Hence calculate the expected number of defaults within the next two
years from a group of 100 companies, all initially rated A. [2]
(iii) Calculate the expected number of defaults for this investment manager over
the next two years, given that the portfolio initially consists of 100 A-rated
bonds. [2]
(iv) Comment on the suggestion that the downgrade trigger strategy will improve
the return on the portfolio. [2]
[Total 8]
the probability that a policyholder who had claims in both previous years will
make a claim in the current year is 0.25
the probability that a policyholder who had claims in one of the previous two
years will make a claim in the current year is 0.15; and
the probability that a policyholder who had no claims in the previous two years
will make a claim in the current year is 0.1
(i) Construct this as a Markov chain model, identifying clearly the states of the
chain. [2]
(iii) Explain why this Markov chain will converge to a stationary distribution. [2]
(iv) Calculate the proportion of policyholders who, in the long run, make at least
one claim at a given year. [4]
[Total 9]
An insurance company assumes that the arrival of motor insurance claims follows an
inhomogeneous Poisson process.
Data on claim arrival times are available for several consecutive years.
(ii) (a) Describe the main steps in the verification of the company s
assumption.
(b) State one statistical test that can be used to test the validity of the
assumption.
[3]
(iii) The company concludes that an inhomogeneous Poisson process with rate
t 3 cos 2 t is a suitable fit to the claim data (where t is measured in
years).
(a) Comment on the suitability of this transition rate for motor insurance
claims.
( f ( s, t )) j .exp( f ( s, t ))
P0 j ( s, t )
j!
(d) Comment on the form of the solution compared with the case where
is constant.
[8]
[Total 12]
END OF PAPER
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 6 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
Faculty of Actuaries
CT4 (104) S2006 Institute of Actuaries
B1 Calculate 0.25 p80 and 0.25 p80.5 , using the ELT15 (Females) mortality table and
assuming a uniform distribution of deaths. [4]
B2 A national mortality investigation is carried out over the calendar years 2002, 2003
and 2004. Data are collected from a number of insurance companies.
Deaths during the period of the investigation, x, are classified by age nearest at death.
(i) (a) State the rate year implied by the classification of deaths.
(b) State the ages of the lives at the start of the rate interval.
[1]
(ii) Derive an expression for the exposed to risk, in terms of Px(t), which may be
used to estimate the force of mortality in year t at each age. State any
assumptions you make. [3]
(iii) Describe how your answer to (ii) would change if the census information
provided by some companies was Px* t , the number of in-force policies on
1 January each year, where policyholders are classified by age last birthday.
[3]
[Total 7]
B3 An investigation was undertaken into the effect of a new treatment on the survival
times of cancer patients. Two groups of patients were identified. One group was
given the new treatment and the other an existing treatment.
T
hi t h0 t exp z
where: hi t is the hazard at time t, where t is the time since the start of treatment
h0 t is the baseline hazard at time t
A the risk of death for a male patient is 1.02 times that of a female
patient; and
B the risk of death for a patient given the existing treatment is 1.05 times
that for a patient given the new treatment
(ii) Estimate the ratio by which the risk of death for a male patient who has been
given the new treatment is greater or less than that for a female patient given
the existing treatment. [2]
(iii) Determine, in terms of the baseline hazard only, the probability that a male
patient will die within 3 years of receiving the new treatment. [2]
[Total 7]
B4 An investigation took place into the mortality of persons between exact ages 60 and
61 years. The table below gives an extract from the results. For each person it gives
the age at which they were first observed, the age at which they ceased to be observed
and the reason for their departure from observation.
1 60 0 60 6 withdrew
2 60 1 61 0 survived to 61
3 60 1 60 3 died
4 60 2 61 0 survived to 61
5 60 3 60 9 died
6 60 4 61 0 survived to 61
7 60 5 60 11 died
8 60 7 61 0 survived to 61
9 60 8 60 10 died
10 60 9 61 0 survived to 61
(ii) List the strengths and weaknesses of the Binomial model for the estimation of
empirical mortality rates, compared with the Poisson and two-state models.
[3]
[Total 8]
An extract from the data for 12 policyholders is shown in the table below.
(iii) Sketch on a suitably labelled graph the Nelson-Aalen estimate of the survival
function. [2]
[Total 9]
B6 (i) (a) Describe the general form of the polynomial formula used to graduate
the most recent standard tables produced for use by UK life insurance
companies.
(b) Show how the Gompertz and Makeham formulae arise as special cases
of this formula.
[3]
2
o 1 1
x
1 exp b0 b1 x b2 x .
2 2 2
(a) Explain why this might be a sensible formula to choose for this class of
lives.
(b) Suggest two techniques which can be used to perform the graduation.
[3]
(iii) The table below shows the crude and graduated mortality rates for part of the
relevant age range, together with the exposed to risk at each age and the
standardised deviation at each age.
(c) the existence of individual ages at which the graduated rates depart to a
substantial degree from the observed rates
[9]
[Total 15]
END OF PAPER
EXAMINATION
September 2006
EXAMINERS’ REPORT
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
M A Stocker
Chairman of the Board of Examiners
November 2006
© Faculty of Actuaries
© Institute of Actuaries
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
Comments
Comments on solutions presented to individual questions for this September 2006 paper are
given below.
103 Part
Question A1 This was reasonably well answered, even by the weaker candidates.
In part (ii), very few candidates used the information in the question and
n
calculated ∑ (ti − si ) .
i =1
Question A4 This was not well answered overall, but many of the stronger candidates did
score highly.
In part (i), some candidates incorrectly attempted to calculate the long-run
probability of being in state B.
Part (ii) was generally well answered.
In part (iv), the stronger candidates provided good answers, but overall
candidates did not score well here.
Question A5 Overall this was poorly answered, although the stronger candidates did well.
Many candidates failed to split the two states labelled B and C in the solution,
giving instead a 3-state chain. Some marks were still awarded for the long-
run probability calculations in part (iv), but such candidates were not able to
calculate the required final answer.
Question A6 This was poorly answered by most candidates, even though some parts of the
question had been asked in previous (103) exams.
Marks were lost in all parts of the question. Many candidates did not make a
serious attempt at part (iii)(c).
Page 2
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
104 Part
Question B2 This was poorly answered overall, although some of the stronger candidates
did manage to score highly.
In part (ii), the question asked candidates to “derive an expression” and
therefore we were looking for clearly set out steps here. Many candidates lost
marks by not providing sufficient explanation of their working.
Page 3
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
103 Solutions
A1 (i) If the ith component is still working at the end of the test period its
contribution to the likelihood is:
ti − si psi = exp(−μ(ti − si ))
If the ith component fails at time ti its contribution to the likelihood is:
exp(−μ(ti − si )).μ fi
(ii) Denote the total number of components used in the test by n. The likelihood
for n independent components is:
n
L = ∏ exp(−μ(ti − si )).μ fi
i =1
n n
i =1
Now the rig contains 100 components at all times because it is fully loaded
n
and failed components are immediately replaced, so ∑ (ti − si ) = 200(years) .
i =1
n
∑ fi
So L = exp ( −200μ ) ⋅μ i =1
n
ln L = −200μ + ln μ.∑ fi
i =1
n
∂ ln L
∑ fi
= −200 + i =1
∂μ μ
Page 4
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
n
∑ fi
i =1
μˆ =
200
∂ 2 ln L
∑ fi
i =1
2
=− <0
∂μ μ2
⎛ −σ σ⎞
A=⎜
⎝ ρ − ρ ⎟⎠
∞
−σs
∫ σe ds = e−σt .
t
(iv) The time spent in state A before the next visit to B has mean 1/σ.
(v) Testing whether the successive holding times are exponential variables and
independent would be best. Any procedure which does this test is acceptable.
Page 5
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
N0 = 0;
Nt has independent increments;
Nt has stationary increments, each having a Poisson distribution, i.e.
P [ Nt − N s = n ] =
[ λ (t − s ) ] e −λ (t − s )
n
, s < t , n = 0,1, 2,...
n!
Nt
Xt = ∑Yj , t ≥ 0.
j =1
n
X n = ∑Yj
j =1
∞
with initial condition X0 = 0. Then { X n }n=0 constitutes a general
random walk.
(ii) (a) A Poisson process operates in continuous time and has a discrete state
space, the set of nonnegative integers.
(c) A general random walk operates in discrete time. Again, this has a
discrete or continuous state space according to whether the variables Yj
have a discrete or continuous distribution.
Page 6
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
∞
0.03
= 0.03 × ∑ 0.92i = = 0.375
i =0 1 − 0.92
Page 7
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
(iii) Either
⎛ 0.97 0 0.03 ⎞
⎜ ⎟
A′ = ⎜ 0 0 0 ⎟
⎜ 0 0 1 ⎟⎠
⎝
⎛ 0.9409 0 0.0591⎞
2 ⎜ ⎟
A′ = ⎜ 0 0 0 ⎟
⎜ 0 0 1 ⎟⎠
⎝
Or
Required probability is
(iv) The expected number of defaults has been reduced by this strategy. (The
variance of the number of defaults would also reduce.)
However it is not possible to tell whether the overall return is improved as this
depends on the price at which bonds were bought and sold at the end of year 1.
The price of the debt sold may have been depressed by the companies having
been downgraded to rating B, and the manager loses out on any increase in
price if they recover.
The “downgrade trigger” strategy will incur dealing costs, which should be
considered when comparing the returns.
Page 8
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
A5 (i) Consider the following four states that the policyholder might be at the end of
a year:
• the policyholder has made at least one claim both in the year just ended
and the previous one (state A)
• the policyholder has made no claims in the year just ended but s/he made
at least one claim during the previous year (state B)
• the policyholder has made at least one claim in the year just ended but not
in the previous one (state C)
• the policyholder has made no claim during either the year ended or the
previous one (state D)
If the year ended is year n, and Xn denotes the current state of the policyholder,
then Xn constitutes a Markov chain.
⎛ 0.25 0.75 0 0 ⎞
⎜ ⎟
0 0 0.15 0.85 ⎟
P=⎜
⎜ 0.15 0.85 0 0 ⎟
⎜ ⎟
⎝0 0 0.10 0.90 ⎠
(iii) The chain has a finite number of states (A,B,C,D). In order to show that it has
a stationary distribution, it suffices to show that it is irreducible and aperiodic.
It is apparent from the transition matrix above that any state can be reached
from any other; hence the chain is irreducible.
The chain is also aperiodic since for states A, D the state can remain at the
same state after one step, while for states B, C the state may return to its
current state after 2 or 3 steps.
Page 9
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
Using the transition matrix from (ii) above we obtain the equations
Discard the last of these equations and use also that the stationary probabilities
must also satisfy
πA + πB + πC + πD = 1 (4)
Or 5 πA = πC
πB = πC
17
πD = πB .
2
⎛1 17 ⎞ 10
πB ⎜ + 1 + 1 + ⎟ = 1 ⇒ πB = .
⎝5 2⎠ 107
2 10 85
πA = , πC = , πD = .
107 107 107
The proportion of policyholders who, in the long run, make at least one claim
in a given year is
12
π A + πB = .
107
Page 10
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
A6 (i) The probability that an event occurs during the short time interval between t
and t + h is approximately equal to λ(t) h for small h where λ(t) is called the
rate of the process. For a time-inhomogeneous process, λ(t) depends on the
current time t; for a time-homogeneous process it is independent of time.
(ii) (a) Divide the time period into intervals of a suitable size, say one month.
Estimate the arrival rate separately for each time period.
See if the observed data match the pattern which would be expected if
the model were accurate and if the parameters had their values given
by their estimates.
(b) A goodness of fit test, such as the chi-squared test, should be carried
out for each time period chosen.
Tests for serial correlation [e.g. portmanteau test] should use the whole
data set at once.
(iii) (a) This implies that claims are seasonal with period 12 months, and that
claims in the peak (presumably winter) are double those at the low
point of the year.
∂
P ( s, t ) = P( s, t ). A(t ) t≥s
∂t
Where:
⎛ −λ (t ) λ (t ) ⎞
⎜ ⎟
−λ(t ) λ (t )
A(t ) = ⎜ ⎟
⎜ −λ (t ) ⎟
⎜ ⎟
⎝ ⎠
∂
P0 j ( s, t ) = λ (t ).P0, j −1 ( s, t ) − λ (t ).P0 j ( s, t ) (I)
∂t
with P0 j ( s, s ) = 0
Page 11
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
( f ( s, t )) j .exp(− f ( s, t ))
P0 j ( s, t ) =
j!
LHS of I
exp(− f ( s, t )) d
( j.( f ( s, t )) j −1 − f ( s, t ) j ). . f ( s, t )
j! dt
RHS of I
f ( s, t ) j −1 f ( s, t ) j .exp(− f ( s, t ))
λ (t ). .exp(− f ( s, t )) − λ(t ).
( j − 1)! j!
∂
f ( s, t ) = λ (t )
∂t
Now
t t
∫ λ(v)dv = ∫ (3 + cos(2πv))dv
s s
t
⎡ 1 ⎤
= ⎢3v + sin(2πv) ⎥
⎣ 2π ⎦s
1
= 3(t − s ) + [sin(2πt ) − sin(2πs)] ≡ f ( s, t )
2π
∂
P00 ( s, t ) = −λ (t ).P00 ( s, t ) (II)
∂t
Page 12
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
LHS of II
∂
− exp(− f ( s, t )). ( f ( s, t )) = − P00 ( s, t ).λ (t )
∂t
and P00 ( s, s ) = 1
(d) Solution is of the same form, except that for the homogeneous case
f(s,t) = λ(t-s).
Page 13
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
104 Solutions
ALTERNATIVE 1
(t − s )qx
t − s qx+ s = .
1 − sqx
0.25q80
0.75−0.5 q80+ 0.5 = , and so
1 − 0.5q80
0.25q80
0.25 p80.5 = 1− .
1 − 0.5q80
ALTERNATIVE 2
Using t px = s px ⋅ t − s px+ s,
Page 14
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
B2 (i) (a) The age definition changes 6 months before/after each birthday, so this
is a life year rate interval.
(ii) Under the principle of correspondence the age definition of deaths and census
should correspond, which they do here. So we do not need to adjust the
census information.
3
The exposed to risk is given by E xc = ∫ Px ( t ) dt .
0
2
1
E xc = ∑
2
( Px ( t ) + Px ( t + 1) ) , where t is measured from 1 January 2002
0
⎛1 1 ⎞
= ⎜ Px ( 0 ) + Px (1) + Px ( 2 ) + Px ( 3) ⎟
⎝2 2 ⎠
(iii) The age definitions for deaths and census no longer correspond. So, we need
to adjust the census information for those companies who supply details of
Px* ( t ) .
2
1
E xc = ∑
2
( Px ( t ) + Px ( t + 1) )
0
2
1⎛1
( 1
) (
= ∑ ⎜ Px*−1 ( t ) + Px* ( t ) + Px*−1 ( t + 1) + Px* ( t + 1) ⎟
⎞
)
0 2⎝2 2 ⎠
1
4
( 1
2
) ( 1
) (
= Px*−1 ( 0 ) + Px* ( 0 ) + Px*−1 (1) + Px* (1) + Px*−1 ( 2 ) + Px* ( 2 ) + Px*−1 ( 3) + Px* ( 3)
4
)
Page 15
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
h f ( t ) = h0 ( t ) × exp ( 0 + β2 z2 )
hm ( t ) = h0 ( t ) × exp ( β1 × 1 + β2 z2 )
hm ( t ) = 1.02 × h f ( t ) , so that:
( ) (
h0 ( t ) × exp βˆ 1 + βˆ 2 z2 = 1.02 × h0 ( t ) × exp βˆ 2 z2 )
⇒ exp(βˆ 1 ) = 1.02
( ) (
h0 ( t ) × exp βˆ 1 z1 + 0 = 1.05 × h0 ( t ) × exp βˆ 1z1 + βˆ 2 z2 )
⇒ 1 =1.05 × exp βˆ 2 ( )
⇒ βˆ 2 = ln 1 (
1.05 )
= − 0.0488
(ii) The hazard for a male patient who has been given the new treatment is:
hm,n ( t ) = h0 ( t ) × exp ( β1 ×1 + β2 × 1)
= h0 ( t ) × exp ( −0.029 )
= 0.9714 × h0 ( t )
The hazard for a female patient given the existing treatment is the baseline
hazard.
Page 16
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
Hence, the ratio of the hazard for a male patient who has been given the new
treatment to that for a female patient given the existing treatment is:
hm,n ( t )
= 0.9714
h0 ( t )
ALTERNATIVELY
hm,n ( t ) ⎡ h (t ) ⎤ ⎡ h (t ) ⎤ 1
= ⎢ m, x ⎥ × ⎢ x ,n ⎥ = 1.02 × = 0.9714
h f ,e ( t ) ⎣⎢ h f , x ( t ) ⎦⎥ ⎣⎢ hx ,e ( t ) ⎦⎥ 1.05
{ 3
1 − Sm,n ( 3) = 1 − exp − ∫ hm,n ( s )ds
0 }
= 1 − exp − {∫ 3
0 }
0.9714 × h0 ( s )ds
⎧ ⎫
= 1 − exp ⎨0.9714 × ⎛⎜ − ∫ h0 ( s )ds ⎞⎟ ⎬
3
⎩ ⎝ 0 ⎠⎭
0.9714
⎛ − ∫03 h0 ( s )ds ⎞
=1 − ⎜ e ⎟
⎝ ⎠
B4 (i) Let the age individual i enters observation be ai and the age that individual i
leaves observation be bi. Define an indicator variable di such that di = 0 if
individual i is not observed to die and di = 1 if individual i dies.
10
∑ di
i =1
qˆ60 = 10
.
∑ (1 − ai − ⎡⎣(1 − di )(1 − bi )⎤⎦ )
i =1
Page 17
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
The denominator in this formula shows that for persons who do not die
(di = 0) the exposed to risk is bi – ai and for persons who die (di = 1) the
exposed to risk is 1 – ai.
Thus the relevant calculations are shown in the table below (all durations are
in years).
Totals 4 74/12
4
Therefore qˆ60 = = 0.6486 .
74 /12
ALTERNATIVELY
10
Take the central exposed to risk, ∑ (bi − ai ) (in years) and add
1
½d60 to give the initial exposed to risk.
d60 4 4
qˆ60 = = = = 0.5783.
c
E60 + 0.5d 60 (59 /12) + 2 83 /12
[This approach is inferior to the first, as it does not use all the information
available in the data, and involves the assumption that the deaths take place,
on average, half way through the year.]
Page 18
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
B5 (i) There will be Type I censoring of lives that survive to age 55 years.
There will be random censoring of lives that withdraw before age 55 years.
(ii) The calculations are shown in the table below, where durations are measured
in years since the 50th birthday.
Using the convention that, when deaths and withdrawals are observed at the
same duration, deaths occur first:
tj Nj dj cj dj /Nj ˆ =
Λ t ∑ (d j / N j )
t j ≤t
0 12
0.25 12 1 1 0.0833 0.0833
1.00 10 1 2 0.1000 0.1833
2.75 7 1 2 0.1429 0.3262
4.25 4 1 3 0.25 0.5762
Since Sˆ (t ) = exp(−Λ
ˆ )
t
t Sˆ (t )
Page 19
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
(iii)
11
0.9
0.9
0.8
0.8
0.7
0.7
0.6
0.6
S(t)
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.2
0.1
0.1
00
0 0.5 1 1.5 2 2.5 33 3.5 44 4.5 55
Duration since 50th birthday
Duration
μ x = (polynomial(1)) + exp(polynomial(2)) ,
α 0 + α1x + α 2 x 2 + ...
β0 + β1x + β2 x 2 + ....
αi = 0 for all i
and
βi = 0 for i = 2, 3, ….
Page 20
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
μ x = exp(β0 + β1x)
2
⎛ 1⎞
The term b2 ⎜ x + ⎟ in the formula can allow for this possibility.
⎝ 2⎠
⎡ ⎤ 1 ⎛ 1⎞
2
of log ⎢μˆ 1 ⎥ on x + and ⎜x+ ⎟ .
⎢⎣ x+ 2 ⎥⎦ 2 ⎝ 2⎠
(iii) (a) The null hypothesis is that there is no difference between the graduated
rates and the underlying rates in the population from which the crude
rates are derived.
∑ z x 2 ∼ χ2m ,
x
Page 21
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
Age x zx zx2
last
birthday
50 -0.12031 0.01447
51 -0.20055 0.04022
52 -0.24749 0.06125
53 0.11341 0.01286
54 -0.79336 0.62942
55 -0.66436 0.44137
56 -0.44369 0.19686
57 -0.35225 0.12408
Sum 1.52053
(b) To test for bias we use EITHER the Signs Test or the Cumulative
Deviations test.
Signs Test
Since this probability > 0.025 (two-tailed test) we do not reject the null
hypothesis.
We conclude that the graduated rates are not biased above or below the
crude rates.
Page 22
Subject CT4 — Models Core Technical — September 2006 — Examiners’ Report
o
∑ (μˆ x+ 1 Ex − μ x+ 12 Ex )
x 2 ~ Normal(0,1) .
o
∑ μ x+ 12 Ex
x
o o
Age x μˆ 1 Ex − μ x+ 1 Ex μ x+ 1 Ex
x+ 2 2
2
last
birthday
50 -0.63 27.63
51 -1.06 28.06
52 -1.32 28.32
53 0.61 29.39
54 -4.13 27.13
55 -3.20 23.20
56 -2.03 21.03
57 -1.72 23.72
(-13.48/√208.48) = -0.9335.
using a two-tailed test, the absolute value of the test statistics is less
than 1.96, so we do not reject the null hypothesis.
We conclude that the graduated rates are not biased above or below the
crude rates.
(c) To test for the existence of individual ages at which the graduated rates
depart greatly from the observed rates we can use the Individual
Standardised Deviations Test.
Therefore we do not reject the null hypothesis and conclude that there
are no outliers.
END OF EXAMINERS’ REPORT
Page 23
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 11 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
© Faculty of Actuaries
CT4 A2007 © Institute of Actuaries
1 (a) Define, in the context of stochastic processes, a:
1. mixed process
2. counting process
(i) (a) Suggest, with reasons, a suitable method of graduation in this case.
(ii) Comment on any further considerations that the company should take into
account before using the graduated rates for premium calculations. [2]
[Total 5]
3 The government of a small country has asked you to construct a model for forecasting
future mortality.
Outline the stages you would go through in identifying an appropriate model. [6]
4 The actuary to a large pension scheme carried out an investigation of the mortality of
the scheme’s pensioners over the two years from 1 January 2005 to 1 January 2007.
(i) List the data required by the actuary for an exact calculation of the central
exposed to risk for lives aged x. [2]
CT4 A2007—2
(ii) (a) Derive an expression that could be used to estimate the central exposed
to risk using the available data. State any assumptions you make.
(b) Use the data to estimate μ65 . State any further assumptions that you
make. [4]
[Total 6]
5 (i) Define the hazard rate, h(t), of a random variable T denoting lifetime. [1]
(ii) An investigation is undertaken into the mortality of men aged between exact
ages 50 and 55 years. A sample of n men is followed from their 50th
birthdays until either they die or they reach their 55th birthdays.
The hazard of death (or force of mortality) between these ages, h(t), is
assumed to have the following form:
h(t ) = α + β t
(a) Derive an expression for the survival function between ages 50 and 55
years.
(c) Comment on the appropriateness of the assumed form of the hazard for
modelling mortality over this age range.
[6]
[Total 7]
A B
An instrument was used to monitor this process, but it was set up incorrectly and only
recorded the state occupied after every two time periods. From these observations the
following two-step transition probabilities have been estimated:
2
PAA = 0.5625
2
PAB = 0.125
2
PBA = 0.475
2
PCC = 0.4
Calculate the one-step transition matrix consistent with these estimates. [8]
CT4 A2007—4
7 Every person has two chromosomes, each being a copy of one of the chromosomes
from one of their parents. There are two types of chromosomes labelled X and Y. A
child born with an X and a Y chromosome is male and a child with two X
chromosomes is female.
A medical researcher wishes to study the progress of the disease through the first born
child in each generation, starting with a female carrier.
• every parent has a equal chance of passing either of their chromosomes to their
children
• the partner of each person in the study does not carry a defective X chromosome;
and
(i) Show that the expected progress of the disease through the generations may be
modelled as a Markov chain and specify carefully:
For those patients who died during the study exact dates of death were recorded as
follows:
F 1 February 2001
G 1 March 2001
H 1 April 2001
I 1 June 2001
J 1 September 2001
K 1 September 2001
L 1 November 2001
For those patients with whom the hospital lost contact before the end of the
investigation:
(i) Explain whether and where each of the following types of censoring is present
in this investigation:
(ii) Calculate the Kaplan-Meier estimate of the survival function for these
patients. State any assumptions that you make. [7]
(iii) Hence estimate the probability that a patient will die within 4 years of surgery.
[1]
[Total 11]
CT4 A2007—6
9 An insurance company is concerned that the ratio between the mortality of its female
and male pensioners is unlike the corresponding ratio among insured pensioners in
general. It conducts an investigation and estimates the mortality of male and female
pensioners, μˆ m x+1/ 2 and μ
f
ˆ x+1/ 2 . It then uses the μ
ˆmx+1/ 2 to calculate what the expected
mortality of its female pensioners would be if the ratio between male and female
mortality rates reflected the corresponding ratio in the PMA92 and PFA92 tables,
S x +1/ 2 , using the formula
μ
xf +1/ 2 = μˆ m
x +1/ 2 S x +1/ 2 .
The table below shows, for a range of ages, the numbers of female deaths actually
observed in the investigation and the number which would be expected from the
μ f
x+1/ 2 .
65 30 28.4
66 20 30.1
67 25 31.2
68 40 33.5
69 45 34.1
70 50 41.8
71 50 46.5
72 45 44.5
(i) Describe and carry out an overall test of the hypothesis that the ratios between
male and female death rates among the company’s pensioners are the same as
those of insured pensioners in general. Clearly state your conclusion. [5]
(ii) Investigate further the possible existence of unusual ratios between male and
female death rates among the company’s pensioners, using two other
appropriate statistical tests. [6]
[Total 11]
The partnership decides to investigate whether this claim is true using a multiple-state
model with three states: (1) fully qualified but not yet a partner, (2) fully qualified and
a partner, (3) working for another partnership. The period of the investigation is to
be 1 January 1997 to 31 December 2006.
(i) (a) Draw and label a state-space diagram depicting the chosen model,
showing possible transitions between the three states.
(b) State any assumptions implied by the diagram you have drawn and
comment on their appropriateness.
[3]
(ii) (a) State what data would be required in order to estimate the transition
intensity of moving from state (1) to state (2) for employees aged 30
years last birthday.
The investigation assumes that all transition intensities are constant within
each year of age. [7]
(iii) Derive, using these data, an estimate for the profession as a whole of the
corresponding transition intensity of becoming a partner among persons aged
30 years last birthday during the period of the investigation. State any
assumptions you make. [5]
[Total 15]
CT4 A2007—8
11 (i) Consider two Poisson processes, one with rate λ and the other with rate μ .
Prove that the sum of events arising from either of these processes is also a
Poisson process with rate ( λ + μ ). [2]
(b) Describe the circumstances in which the outcome of the Markov jump
chain differs from the standard Markov chain with the same transition
matrix. [4]
An airline has N adjacent check-in desks at a particular airport, each of which can
handle any customer from that airline. Arrivals of passengers at the check-in area are
assumed to follow a Poisson process with rate q. The time taken to check-in a
passenger is assumed to follow an exponential distribution with mean 1/a.
(iii) Show that the number of desks occupied, together with the number of
passengers waiting for a desk to become available, can be formulated as a
Markov jump process and specify:
(iv) State the Kolmogorov forward equations for the process, in component form.
[2]
(vi) (a) Set out the transition matrix of the jump chain associated with the
airline check-in process.
(b) Determine the probability that all desks are in use before any passenger
has completed the check-in process, given that no passengers have
arrived at check-in at the outset. [4]
[Total 17]
END OF PAPER
CT4 A2007—9
Faculty of Actuaries Institute of Actuaries
EXAMINATION
April 2007
EXAMINERS’ REPORT
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
M A Stocker
Chairman of the Board of Examiners
June 2007
© Faculty of Actuaries
© Institute of Actuaries
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
Comments
Comments on solutions presented to individual questions for this April 2007 paper are given
below and further comments, where appropriate, are given in the solutions that follow.
Question 3 This was reasonably well answered, although many candidates took no
account of the particular circumstances referred to in the question.
Question 7 Overall this was reasonably well answered, with the stronger candidates
scoring highly.
Question 10 Parts (i) and (ii) were fairly well answered overall, but few candidates scored
well in part (iii).
Page 2
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
1 Mixed process
(a) Is a stochastic process that operates in continuous time, which can also change
value at predetermined discrete instants.
Counting process
(a) Is a process, X, in discrete or continuous time, whose state space is the natural
numbers {0, 1, 2, …}.
There are likely to be existing standard tables which are suitable and
this method is suitable for relatively small data sets.
• Plot the crude rates against qxs from the standard table to identify a
simple relationship.
• Test the graduation for goodness of fit. If the fit is not adequate,
the process should be repeated.
Page 3
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
• The rates will be based on current mortality; the company should also take
into account expected future changes, especially any reductions in
mortality rates.
• Premiums charged by other insurer: if rates are too high the company will
fail to attract business; if too low, it may attract too much, unprofitable
business.
3 Clarify the purpose of the exercise. Why does the government want forecasts of
mortality? What is the period for which the forecast is wanted? Is it short (e.g. 5–10
years) or long (e.g. 50–70 years).
Consult the existing literature on models for forecasting mortality, and speak to
experts in this field of application. Consider using or adapting existing models which
are employed in other countries.
Establish what data are available (e.g. on past mortality trends in the country,
preferably with deaths classified by age and cause of death).
On the basis of what data are available, define the model you propose to use. If the
data are simple and not detailed, then a complex model is not justified. Will a
deterministic or a stochastic model be appropriate in this case?
Identify suitable computer software to implement the model, or, if none exists, write a
bespoke program.
Debug the program or, if existing software is used, check that it performs the
operations you intend it to do.
Run the model and test the reasonableness of the output. Consider, for example, the
forecast values of quantities such as the expectation of life at birth.
Write a report documenting the results and the model and communicate the results
and the output to the government of the small country.
Page 4
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
4 (i) For each pensioner in the investigation, the actuary would need:
(ii) (a) The central exposed to risk of pensioners aged x nearest birthday is
given by
2
E xc = ∫ Px,t
0
1
≈ ∑ 12 ( Px,t + Px,t +1 ) = 1P
2 x ,0
+ Px,1 + 12 Px,2
0
(b) This is a life year rate interval, from age x-½ to x+½. The age in the
middle of the rate interval is x, so μ̂ estimates μ x , assuming a constant
force of mortality over the life year.
d65,2005 + d65,2006
μˆ 65 = c
E65
16 + 23 39
= =
( 1 × 1678 + 1720 + 1 × 1622
2 2 ) 3370
= 0.01157
Page 5
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
1
h(t ) = lim ( Pr[T ≤ t + dt | T > t ]) .
dt →0+ dt
⎛ t ⎞
S (t ) = exp ⎜ − ∫ h( s )ds ⎟ ,
⎜ ⎟
⎝ 0 ⎠
then
⎛ t ⎞ ⎡ βs 2 ⎤
t
⎡ βt 2 ⎤
S (t ) = exp ⎜ − ∫ ( α + β s ) ds ⎟ = exp ⎢ −αs − ⎥ = exp ⎢ −αt − ⎥
⎜ ⎟ ⎢
⎣ 2 ⎦⎥ ⎢
⎣ 2 ⎦⎥
⎝ 0 ⎠ 0
where 0 ≤ t ≤ 5 .
1
0.9
0.8
0.7
0.6
S(t)
0.5
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5
Duration since age 50 years
(c) If both α and β are positive, then the formula implies a force of
mortality which increases with age, which is sensible for this age
range.
Page 6
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
6 Based on the given transition diagram, the one-step transition matrix must be of the
form:
⎛a 0 c⎞
⎜ ⎟
⎜d e f⎟
⎜0 i ⎟⎠
⎝ h
c ⎞ ⎜⎛ a c(a + i ) ⎞
2
⎛a 0 c ⎞ ⎛a 0 ch
⎟
⎜ ⎟ ⎜ ⎟
⎜d e f ⎟*⎜ d e f ⎟ = ⎜ d (a + e) e 2 + fh cd + ef + fi ⎟
⎜0 ⎜ ⎟
⎝ h i ⎟⎠ ⎜⎝ 0 h i ⎟⎠ ⎜ dh h (e + i ) fh + i 2 ⎟
⎝ ⎠
2
PAA = 0.5625 ⇒ a 2 = 0.5625 ⇒ a = 0.75
So, a+c=1
and c = 0.25
2
PAB = 0.125 ⇒ ch = 0.125 ⇒ h = 0.5
h+i=1
so i=0.5
2
PCC = 0.4 ⇒ f × 0.5 + 0.52 = 0.4 ⇒ f = 0.3
2
PBA = 0.475 ⇒ d (0.75 + e) = 0.475
Substitute for e:
Page 7
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
⎛ 0.75 0 0.25 ⎞
⎜ 0.5 0.2 0.3 ⎟
⎜ ⎟
⎜ 0 0.5 0.5 ⎟⎠
⎝
7 (i) Consider the sequence of the status of the first born child in each generation.
Page 8
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
0.5
0.25
0.25
FN FC
0.5 0.5
0.25
0.25
0.5
MN MH
0.5
0.5
Each of the transition probabilities depends only on state currently occupied,
so the process possesses the Markov property.
(ii) (a) The chain is reducible because once it enters states FN or MN it cannot
access FC or MH.
FN FC MN MH
FN (0) 0.5 0 0.5 0
A = FC (1) 0.25 0.25 0.25 0.25
MN (2) 0.5 0 0.5 0
MH (3) 0 0.5 0.5 0
Page 9
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
So,
An alternative solution combines the states FN and MN to give a 3-state model. This
was given credit.
8 (i) (a) Type I censoring is present for those lives still under observation at 31
December 2005 as the censoring times are known in advance.
Page 10
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
(ii) The durations at which lives died or were censored are shown below. Duration
is measured in years and months from the date of surgery.
tj nj dj cj λˆ j = d j / n j
0 15 0 0 0
0.5 15 1 1 1/15
0.833 13 2 0 2/13
1.917 11 1 3 1/11
4.333 7 1 6 1/7
(
The estimated survival function is given by, Sˆ ( t ) = ∏ 1 − λ j . So, )
t j ≤t
t Ŝ ( t )
0.000 ≤ t < 0.500 1.0000
0.500 ≤ t < 0.833 0.9333
0.833 ≤ t < 1.917 0.7897
1.917 ≤ t < 4.333 0.7179
4.333 ≤ t < 5.0 0.6154
Page 11
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
(iii) The probability that a patient will die within 4 years of surgery is estimated
by:
1 − Sˆ ( 4 ) = 1 – 0.7179
= 0.2821
E xcμˆ xf +1/ 2 − E xc μ
xf +1/ 2
zx = .
E xc μ
xf +1/ 2
∑ zx2 = 11.3343.
x
Since 11.3343 < 15.51, we have no reason to reject the null hypothesis that the
sex ratios of death rates among the company’s pensioners are the same as
those prevailing in the PMA92 and PFA92 tables.
Page 12
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
Using the individual standardised deviations test, we note that none of the zxs
exceeds 1.96 in absolute value, so there is no evidence that the sex ratios
among the company’s pensioners are unusual at any specific ages
Signs test
Since this is greater than 0.025 (two-tailed test), the sex ratios of death rates
among the company’s pensioners are not systematically higher or lower than
those derived from the PMA92 and PFA92 tables.
and since |0.875| < 1.96 using a two-tailed test, the sex ratios of death rates
among the company’s pensioners are not systematically higher or lower than
those derived from the PMA92 and PFA92 tables.
Page 13
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
Credit was only given for one of the Signs test and the Cumulative Deviations
test as they both test for bias.
1 7 1 8
z (1) = ∑x
7 1
z = 0.3029 , and z (2)
= ∑ zx = 0.2707
7 2
Sum 4.8231
2 2
Age ⎡ z x − z (1) ⎤ ⎡ z x +1 − z (2) ⎤
⎣ ⎦ ⎣ ⎦
65 0.0000 4.4592
66 4.5962 1.9064
67 1.9963 0.7264
68 0.6726 2.5467
69 2.4450 0.9951
70 0.9320 0.0588
71 0.0442 0.0383
4.8231
r1 = = 0.4503
(10.6863)(10.7310)
Page 14
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
Since 1.27 < 1.645, we conclude that there is no evidence that the sex ratios of
death rates among the company’s pensioners vary with age in a way different
from the ratios derived from PMA92 and PFA92.
Note that the Grouping of Signs test is not appropriate with 8 ages, 6 positive
and 2 negative signs.
μ 21
x +t
1 Fully qualified but 2 Fully qualified
not yet a partner and a partner
μ12
x +t
μ13
x +t μ31
x +t μ32
x +t μ 23
x +t
3 Working for
another company
(b) The chosen model ignores death among persons in the relevant age
groups. Since mortality in this age group among professional people is
likely to be low, this seems reasonable.
This diagram assumes that demotion is possible, i.e. some-one who has
become a partner can return to non-partnership status without leaving
the company.
The assumption is also made that a new employee joining from another
company can do so as a partner.
Credit was given for models based on alternative assumptions, provided these
were reasonable.
• the total waiting time during the calendar years 1997–2006 in state
(1) when aged 30 last birthday
Page 15
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
N
L = ∏ K exp[−(μ13 + μ12 )vi ](μ12 ) di
i =1
where
N
log e L = ∑ log e K − (μ12 + μ13 )vi + di log e μ12
i =1
∂ log e L N ∑ di
12
= −∑ vi + i =1
12
,
∂μ i =1 μ
and setting this equal to zero and solving for μ12 gives
N
∑ di
12 i =1
μˆ = N
.
∑ vi
i =1
∂ 2 log e L
∑ di
i =1
12 2
=− 12 2
which must be negative.
(∂μ ) (μ )
Page 16
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
(iii) The data on becoming a partner are classified by age last birthday, which is
the same classification as used in the company’s own investigation, therefore
the relevant intensities will relate to the same age range.
For the correct exposed to risk we only consider those who are members of the
institute but not yet partners.
Let the number of such members in the census in year t who were born in year
s be Pt ,s .
All persons born in year s would be aged x last birthday on 1 January in year
s+x+1.
Therefore, assuming that the Pt ,s change linearly during each calendar year
the correct exposed to risk for the year 1997 is
1
( P1997,1956 + P1998,1957 )
2
and the exposed to risk for the entire 10-year period of the investigation is
t = 2006
1
∑ ( Pt ,t −31 + Pt +1,t −30 ) .
t =1997 2
t = 2006
∑ θt
t =1997
t = 2006
.
1
∑ 2 ( Pt ,t −31 + Pt +1,t −30 )
t =1997
Page 17
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
(ii) (a) A jump chain is formed by recording the state of a Markov jump
process only at the instant when a transition has just been made.
(b) The outcome of the jump chain can only differ from that of the
standard Markov chain if the jump process enters an absorbing state.
(iii) The possible states are 0 to N desks in use with no passengers queuing, and N
desks in use with 0, 1, 2, ….. passengers in the queue.
When all desks are occupied and there are M passengers in the queue denote
the state as N:M.
Transition diagram:
a 2a Na Na Na
q q q q q
d
P0 (t ) = aP1 (t ) − qP0 (t )
dt
d
Pr (t ) = a(r + 1) Pr +1 (t ) + qPr −1 (t ) − (ar + q ) Pr (t ) r+1 ≤ N
dt
d
PN :0 (t ) = aNPN :1 (t ) + qPN −1 (t ) − (aN + q) PN :0 (t )
dt
d
PN :m (t ) = aNPN :m+1 (t ) + qPN :m−1 (t ) − (aN + q ) PN :m (t ) m≥1
dt
Page 18
Subject CT4 — Models Core Technical — April 2007 — Examiners’ Report
However if the airline operates many flights this may not be an issue.
⎛ 0 1 ⎞
⎜ ⎟
⎜ a 0
q ⎟
⎜a+q a+q ⎟
⎜ ⎟
⎜ 2a q ⎟
0
⎜ 2a + q 2a + q ⎟
⎜ ⎟
⎜ % % % ⎟
⎜ Na q ⎟
⎜ 0 ⎟
⎜ Na + q Na + q ⎟
⎜ Na q ⎟
⎜ 0 ⎟
⎜ Na + q Na + q ⎟
⎜ % % ⎟⎠
⎝
(b) This is the probability that all the first N transitions are to the right in
the transition diagram.
N −1
1
Required probability is therefore q N −1
.∏
i =1
ia + q
Page 19
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 11 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the
Formulae and Tables and your own electronic calculator.
© Faculty of Actuaries
CT4 S2007 © Institute of Actuaries
1 List the factors you would consider when assessing the suitability of an actuarial
model for its purpose. [4]
2 A particular baker’s shop in a small town sells only one product: currant buns. These
currant buns are delicious and customers travel many miles to buy them.
Unfortunately, the buns do not keep fresh and cannot be stored overnight.
The baker’s practice is to bake a certain number of buns, K, before the shop opens
each morning, and then during the day to continue baking c buns per hour. He is
concerned that:
(i) Describe a model which would allow you to estimate the probability that the
baker will run out of buns. State any assumptions you make. [3]
(ii) Determine the relevant expression for the probability that the baker will run
out of buns, in terms of K, c, and Bj, the number of buns bought by the day’s
jth customer. [1]
[Total 4]
3 A no-claims discount system has 3 levels of discount: 0%, 25% and 50%. The rules
for moving between discount levels are:
• After a claim-free year, move up to the next higher level or remain at the 50%
discount level.
• After a year with one or more claims, move down to the next lower level or
remain at the 0% discount level.
The long-run probability that a policyholder is in the maximum discount level is 0.75.
Calculate the probability that a given policyholder has a claim-free year, assuming
that this probability is constant.
[5]
CT4 S2007—2
4 A national mortality investigation was carried out. It was suggested that the mortality
of the male population could be represented by the following graduated rates:
o
μ x + 1 =μ sx + 2 1
2 2
The table below shows the graduated rates for part of the age range, together with the
exposed to risk, expected and actual deaths at each age. The squared standardised
deviations that were calculated are also shown.
⎛ c o ⎞
⎜ θ x − E x ⋅μ x + 12 ⎟
The standardised deviations were calculated as z x = ⎝ ⎠
o
E xc ⋅μ x + 1
2
(ii) List two methods of graduating a set of crude mortality rates and state, for
each method:
Age x lx
58 88,792
62 84,173
(ii) Find the actual value of l60 in the tables and hence comment on the relative
validity of the two assumptions you used in part (i). [3]
[Total 8]
CT4 S2007—4
7 In order to boost sales, a national newspaper in a European country wishes to compile
a “fair play league table” for the country’s leading football clubs. On 1 December it
undertakes a survey of all the players who play for these clubs, in which it collects the
following data:
• number of games played by each player since the beginning of the season (the
football season in this country begins in September); and
• for each player who had been dismissed from the field of play between the
beginning of the season and 1 December (inclusive), the number of games he had
played before the game in which he was first dismissed
The statistic the newspaper proposes to use in order to construct its “fair play league
table” is the probability that a player will not have been dismissed in any of his first
10 games. It plans to calculate this statistic for each of the 20 leading clubs.
The following table shows the data collected for the players of the club which was top
of the league on 1 December.
1 12 0
2 12 0
3 12 1 5
4 12 0
5 12 1 7
6 12 0
7 10 0
8 9 1 0
9 9 1 5
10 8 0
11 6 2 2
12 5 0
13 5 0
14 4 1 0
15 4 0
(i) (a) Explain how the Kaplan-Meier estimator can be used to estimate the
newspaper’s statistic from these data.
(b) Comment on the way in which censoring arises and on the type of
censoring produced. [4]
(ii) Calculate the newspaper’s statistic using the data above. [4]
[Total 8]
22 150 160 20
23 160 155 25
(ii) (a) Estimate the initial exposed to risk at ages 22 and 23.
(b) Hence estimate q22 and q23.
[4]
(iii) Explain how you would calculate the initial exposed-to-risk for lives aged 22
years last birthday. [4]
[Total 10]
9 In a game of tennis, when the score is at “Deuce” the player winning the next point
holds “Advantage”. If a player holding “Advantage” wins the following point that
player wins the game, but if that point is won by the other player the score returns to
“Deuce”.
When Andrew plays tennis against Ben, the probability of Andrew winning any point
is 0.6. Consider a particular game when the score is at “Deuce”.
(i) Show that the subsequent score in the game can be modelled as a Markov
Chain, specifying both:
(iii) Calculate the number of points which must be played before there is more than
a 90% chance of the game having been completed. [3]
(iv) (a) Calculate the probability that Andrew wins the game.
(b) Comment on your answer. [4]
[Total 12]
CT4 S2007—6
10 (i) Compare the advantages and disadvantages of fully parametric models and the
Cox regression model for assessing the impact of covariates on survival. [3]
You have been asked to investigate the impact of a set of covariates, including age,
sex, smoking, region of residence, educational attainment and amount of exercise
undertaken, on the risk of heart attack. Data are available from a prospective study
which followed a set of several thousand persons from an initial interview until their
first heart attack, or until their death from a cause other than a heart attack, or until 10
years had elapsed since the initial interview (whichever of these occurred first).
(ii) State the types of censoring present in this study, and explain how each arises.
[2]
(iii) Describe a criterion which would allow you to select those covariates which
have a statistically significant effect on the risk of heart attack, when
controlling the other covariates of the model. [4]
Suppose your final model is a Cox model which has three covariates: age (measured
in age last birthday minus 50 at the initial interview), sex (male = 0, female = 1) and
smoking (non-smoker = 0, smoker = 1), and that the estimated parameters are:
Age 0.01
Sex -0.4
Smoking 0.5
Sex x smoking -0.25
(iv) Describe the final model’s estimate of the effect of sex and of smoking
behaviour on the risk of heart attack. [3]
(v) Use the results of the model to determine how old a female smoker must be at
the initial interview to have the same risk of heart attack as a male non-smoker
aged 50 years at the initial interview. [3]
[Total 15]
The following additional data in respect of secondary transitions were collected from
observation of the same process.
ABC 42 BCA 38
ABA 68 BCB 7
ACA 85 CAB 64
ACB 4 CAC 56
BAB 50 CBA 8
BAC 30 CBC 7
(iii) State the distribution of the number of transitions from state i to state j, given
the number of transitions out of state i. [1]
(nijk − E ) 2
[Hint: Use the test statistic χ = ∑∑∑
2
where E is the expected
i j k E
number of triplets under the distribution in (iii)]
(v) Identify two other aspects of the appropriateness of the fitted model that could
be tested, stating suitable tests in each case. [2]
(vi) Outline two methods for simulating the Markov jump process, without
performing any calculations. [4]
[Total 20]
END OF PAPER
CT4 S2007—8
Faculty of Actuaries Institute of Actuaries
EXAMINATION
September 2007
EXAMINERS’ REPORT
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
M A Stocker
Chairman of the Board of Examiners
December 2007
© Faculty of Actuaries
© Institute of Actuaries
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
Comments
Comments on solutions presented to individual questions for this September 2007 paper are
given below and further comments may be written within in the solutions that follow.
Q2 This was the most poorly answered question on the examination paper. Very few
candidates recognised that the baker’s problem could be modelled using the
compound Poisson process described in Unit 2 section 3.4 of the Core Reading.
Q3 This was well answered, with many candidates scoring full marks.
Q4 Although most candidates performed the chi-squared test correctly, few realised that
when using this to test a graduation some degrees of freedom are lost, a fact which is
clearly stated in the Core Reading in Unit 12, section 7.3. In part (ii) comments
tended not to be related to the data in the question; rather they focused rather
mechanically on the shortcomings of the chi-squared test.
Q6 Most candidates obtained the correct numerical answers in part (i) of this question,
but answers to part (ii) were rather sketchy and vague.
Q7 This was more demanding than some previous questions on the Kaplan-Meier or
Nelson-Aalen estimators, and the standard of the answers was lower than expected.
Q8 This exposed-to-risk question was easier than many questions on the same topic in
previous papers. Most candidates scored well on parts (i) and (ii), although few
explained that the method relied on the assumption of a uniform distribution of
deaths. Answers to part (iii) were less impressive and tended to lack detail. Some
candidates couched their answers to this part in aggregate terms, despite the question
clearly referring to individual-level data.
Q9 Most candidates scored well on parts (i) and (ii). Common errors included the use of
a three-state model (Deuce, Advantage and Game) which is inappropriate as the
transition out of the state “Advantage” is ill-defined. Few candidates made attempts
at parts (iii) and (iv) and several of these wrongly thought that part (iv) could be
solved by finding the stationary distribution of the chain.
Q10 Parts (i) and (iv) of this question tested knowledge of Unit 7, sections 2, 3 and 5 of the
Core Reading, which has not been tested in previous CT4 examination papers.
Perhaps because of this, many candidates gave very sketchy and vague answers. In
part (ii), while most candidates spotted that Type I censoring was present, only a
small minority also registered the existence of random censoring. In part (iii) few
candidates correctly interpreted the sex x smoking interaction. Part (v) was well
answered by most candidates.
Q11 Many candidates only attempted parts (i) and (ii) of this question. The remainder was
very poorly answered, with few candidates making serious attempts at part (vi),
despite this being bookwork based on Core Reading, Unit 4, section 5.4.
Page 2
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
Not all these factors needed to be mentioned for full marks to be awarded.
2 (a) Assume that, during each day, customers arrive at the shop according to a
Poisson process.
Assume that the numbers of buns bought by each customer, the Bj, are
independent and identically distributed random variables.
Then if Xt is the total number of buns sold between the beginning of the day
and time t, (where t is measured in hours since the shop opens), Xt is a
compound Poisson process defined by
Nt
Xt = ∑ B j ,
j =1
where the number of customers arriving between the shop opening and time t
is Nt .
(b) The probability that the baker will run out of buns is
Nt
Pr[ K + ct − ∑ B j < 0]
j =1
for some t.
Page 3
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
⎛1 − α α ⎞
⎜ ⎟
⎜1 − α α⎟ .
⎜ 1 − α α ⎟⎠
⎝
(I) π1 = (1 − α ) π1 + (1 − α ) π2
(II) π2 = απ1 + (1 − α ) π3
(III) π3 = απ2 + απ3 .
(IV) π1 + π2 + π3 = 1 .
⎛ 1− α ⎞
(III) gives π2 = ⎜ ⎟ π3 .
⎝ α ⎠
2
⎛ 1− α ⎞
Substituting in (I) gives π1 = ⎜ ⎟ π3 ,
⎝ α ⎠
⎛ ⎛ 1 − α ⎞2 ⎛ 1 − α ⎞ ⎞
and so (IV) leads to ⎜ ⎜ + + 1⎟ π = 1 .
⎜ ⎝ α ⎟⎠ ⎜⎝ α ⎟⎠ ⎟ 3
⎝ ⎠
⎛ (1 − α )2 + α (1 − α ) + α 2 ⎞
⎜ ⎟ × 0.75 = 1 ,
⎜ α 2 ⎟
⎝ ⎠
(( ) ( )
⇒ 0.75 1 − 2α + α 2 + α + α 2 + α 2 = α 2 , )
⇒ 0.25α 2 + 0.75α − 0.75 = 0 .
Page 4
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
4 (i) The null hypothesis is that graduated rates are the same as the true underlying
rates in the population.
we reject the null hypothesis and conclude that the graduated rates do not
adhere to the data.
(ii) From the data we can see that the actual deaths are lower than those
expected at all ages.
The graduated rates are too high; the graduation should be revisited.
There appear to be no outliers (all the zxs have absolute values below
1.96).
Page 5
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
By parametric formula:
Depends on a suitable formula being found which fits the data well.
Should be used if a standard table for a class of lives similar to the experience
is available, and the experience we are interested in does not provide much
data.
and provided the function linking the graduated rates to the rates in the
standard table is simple, this smoothness will be “transferred to the graduated
rates”.
Graphical
The graduation should be tested for smoothness using the third differences of
the graduated rates, which should be small in magnitude and progress
regularly with age.
Page 6
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
Therefore
= 86,482.5.
(b) ALTERNATIVE 1
⎛ 4 ⎞
Then we have 4 p58 = exp ⎜ − ∫ μdx ⎟ = e −4μ .
⎜ ⎟
⎝ 0 ⎠
l62 84,173
But 4 p58 = = = 0.94798 .
l58 88, 792
whence μ = 0.01336.
so l60 = 86,452.
ALTERNATIVE 2
⎛ 4 ⎞
Then we have 4 p58 = exp ⎜ − ∫ μdx ⎟ = e −4μ .
⎜ ⎟
⎝ 0 ⎠
l62
But 4 p58 = .
l58
Page 7
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
l62
and, since 2 p58 = e −2μ = e−4μ = ,
l58
l62
l60 = l58 = l58l62 = (88, 792)(84,173)
l58
so l60 = 86,452
This shows that neither assumption is very accurate, but that the uniform
distribution of deaths (UDD) is closer than the constant force of mortality.
The UDD assumption is better than the constant force of mortality assumption
because UDD implies an increasing force of mortality over this age range,
which is biologically more plausible than the assumption of a constant force.
The fact that the actual value of l60 is considerably greater than that implied by
the UDD assumption suggests that the true rate of increase of the force of
mortality over this age range in English Life Table 15 (males) is even greater
than that implied by UDD.
7 (i) (a) If, for player i, Ti is the number of games played before he is
dismissed, and Ci is the total number of games played before
1 December, and di = 1 if the player had been dismissed before
1 December and 0 otherwise.
then
EITHER
min(Ti,Ci)
and di,
e.g. for player 1, min(Ti,Ci) = 12 and di = 0
OR
Page 8
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
Player min(Ti,Ci) di
1 12 0
2 12 0
3 5 1
4 12 0
5 7 1
6 12 0
7 10 0
8 0 1
9 5 1
10 8 0
11 2 1
12 5 0
13 5 0
14 0 1
15 4 0
(b) Censoring in these data arises because not all players have been
dismissed before 1 December. Those players who have yet to be
dismissed on that data are right-censored.
0≤ t < 2 0.8667
2≤ t < 5 0.8000
5≤ t < 7 0.6222
7≤ t < 12 0.5333
^
Therefore the value of the chosen statistic, S (10) is 0.5333.
Page 9
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
0≤ t < 2 0.8667
2≤ t < 5 0.8000
5≤ t < 7 0.6545
7≤ t < 12 0.5610
^
Therefore the value of the chosen statistic, S (10) is 0.5610.
8 (i) The central exposed to risk at age x, Exc , is the observed waiting time in a
multiple-state or a Poisson model. It is the sum of the times spent under
observation by each life at age x.
The initial exposed to risk requires adjustments for those lives who die, whom
we continue observing until the end of the rate interval.
(ii) The age definition used for both deaths and exposed to risk is the same, so no
adjustment is necessary.
Using the census formula, and assuming that the population aged 22 and 23
years changes linearly over the year, we have, for the central exposed to risk:
1
E xc = ∫ Px,t dt ,
0
so that
1
E xc = ( Px,0 + Px,1 ) .
2
Page 10
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
This assumes that deaths are uniformly distributed across each year of age.
1 20
E22 = (150 + 160) + = 165 ,
2 2
and
1 25
E23 = (160 + 155) + = 170 .
2 2
20 25
Hence q22 = = 0.1212 and q23 = = 0.1471 .
165 170
(iii) ALTERNATIVE 1
The central exposed to risk is calculated as ∑ (bi − ai ) , for all lives i for
i
whom bi − ai > 0 ,
where ai and bi are measured in years since the person’s 22nd birthday, and
Page 11
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
ALTERNATIVE 2
Transition matrix:
The chain is Markov because the probability of moving to the next state does
not depend on history prior to entering that state (because the probability of
each player winning a point is constant)
(ii) The chain is reducible because it has two absorbing states Game A and
Game B.
States Game A and Game B are absorbing so have no period. The other three
states each have a period of 2 so the chain is not aperiodic.
Page 12
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
[This can also be obtained by calculating the square of the transition matrix.]
so that
ln 0.1
N> > 3.14 .
ln(0.48)
But the game can only finish every two points so we require 4 cycles, that is 8
points.
(iv) (a) Define AX to be the probability that A ultimately wins the game when
the current state is X.
We require ADeuce.
Similarly:
and
Page 13
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
So,
0.6
ADeuce = AAdv A ,
0.76
0.6
AAdv A = 0.6 + 0.4 × AAdv A ,
0.76
and
AAdv A = 0.8769 ,
and
ADeuce = 0.6923 .
ALTERNATIVELY
= 0.6923
(b) This is higher than 0.6 because Ben has to win at least two points in a
row to win the game.
10 (i) Fully parametric models are good for comparing homogenous groups, as
confidence intervals for the fitted parameters give a test of difference between
the groups which should be better than non-parametric procedures, or semi-
parametric procedures such as the Cox model.
Page 14
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
(ii) Type I censoring, since the investigation ends after a period which is fixed in
advance.
Random censoring, since death from a cause other than a heart attack is a
random variable and may occur at any time.
Suppose we fit a model with p covariates and another model with p+q
covariates which include all the p covariates of the first model.
Then if the maximised log-likelihoods of the two models are Lp and Lp+q, then
the statistic
−2( L p − L p + q )
This statistic can be used either will full likelihoods or with partial likelihoods
in the Cox model
This statistic can be used to test the statistical significance of any set of q
covariates in the presence of any other disjoint set of p covariates.
The relative risks, compared with the baseline category of male non-smokers
are as follows.
Page 15
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
(v) Let the required age for the woman smoker be 50+x.
The hazard for a male non-smoker aged 50 at the initial interview is simply
h0(x), since this is the baseline category.
Thus we have
so that
or
exp(0.01x - 0.15) = 1
so that
0.01x = 0.15
Page 16
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
ALTERNATIVELY
(ii) (a) The estimator [it is the MLE but this need not
be stated] of λi, λ̂ , is the inverse of the average duration of each visit
to state i.
The estimator [it is the MLE but this need not be stated] of rij, rˆij , is
the proportion of observed jumps out of state i to state j.
rˆAB = 11/20
rˆAC = 9/20
rˆBA = 80/125 =16/25
rˆBC = 9/25
rˆCA = 24/27 =8/9
rˆCB = 1/9
⎛ −4 11 9 ⎞
⎜ 5 5 ⎟
⎜16 −5 9 ⎟
⎜ 5 5 ⎟
⎜⎜ 4 ⎟
1 −3 ⎟
⎝ 3 6 2⎠
Page 17
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
(iv) Null hypothesis is that the Markov property applies to successive transitions,
or that the observed triplets are from a Binomial distribution with the
estimated parameters (given the number of transitions to the middle state).
Using test statistic given in the hint, we can draw up the table below.
(nijk − E ) 2
Triplet nijk E=nij rˆjk
E
Under the null hypothesis, the test statistic follows a χ 2 distribution with the
following number of degrees of freedom:
Number of triplets 12
Minus Number of pairs 6
Plus Number of states 3
Minus One 1
8 degrees of freedom
Page 18
Subject CT4 — Models Core Technical — September 2007 — Examiners’ Report
(v) [Refer back to part (i) — the test in (iv) has only tested that there is no
evidence that the destination that the next jump depends on the previous state
occupied. Need to test the other assumptions].
There is no obvious test statistic for doing this. A suitable test would be to
classify jumps as being from short, medium and long holding times and
investigating these graphically.
Divide time into very short intervals, h, such that σij h is much less than 1.
EXACT METHOD
{ }
Once the path Xˆ n : n = 0,1,... has been generated, the holding times
{Tn : n = 0,1,...} are a sequence of independent exponential random variables,
having parameter λ Xˆ .
n
Page 19
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 11 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the Formulae
and Tables and your own electronic calculator from the approved list.
© Faculty of Actuaries
CT4 A2008 © Institute of Actuaries
1 List four factors in respect of which life insurance mortality statistics are often
subdivided. [2]
2 Describe how smoothness is ensured when mortality rates are graduated using each of
the following methods:
Data are also available about the population of never-married women on 1 January
each year, classified by age last birthday.
You have been asked to estimate the intensity, λ x , of first marriage for women
aged x.
(i) State the rate interval implied by the first marriages data. [1]
(iii) Explain to what age x your estimate of λ x applies. State any assumptions
that you make. [2]
[Total 7]
CT4 A2008—2
6 An investigation was carried out into mortality rates among a certain class of female
pensioners. Crude mortality rates were estimated by single years of age from ages
65–89 years last birthday inclusive. The investigators decided to ask an actuary to
compare the crude rates with a standard table. They calculated the relevant
standardised deviations, printed them out and sent them to the actuary.
The null hypothesis to be tested is that the crude mortality rates come from a
population with underlying mortality consistent with that in the standard table.
(i) List two statistical tests of the null hypothesis which the actuary could carry
out on the basis of the information received. [1]
(ii) Carry out both tests. For each test, state what feature of the experience it is
specifically testing, and give your conclusion. [10]
[Total 11]
7 In a certain small country all listed companies are required to have their accounts
audited on an annual basis by one of the three authorised audit firms (A, B and C).
The terms of engagement of each of the audit firms require that a minimum of two
annual audits must be conducted by the newly appointed firm. Whenever a company
is able to choose to change auditors, the likelihood that it will retain its auditors for a
further year is (80%, 70%, 90%) where the current auditor is (A,B,C) respectively. If
changing auditors a company is equally likely to choose either of the alternative firms.
(i) A company has just changed auditors to firm A. Calculate the expected
number of audits which will be undertaken before the company changes
auditors again. [2]
(ii) Formulate a Markov chain which can be used to model the audit firm used by
a company, specifying:
(iii) Calculate the expected proportion of companies using each audit firm in the
long term. [5]
[Total 11]
Instrument
Piano 0
Violin [-0.05,0.19]
Trumpet [0.07,0.21]
Tuition method
Traditional 0
New [-0.15,0.05]
Sex
Male [-0.08,0.12]
Female 0
(i) Write down a general expression for the Cox proportional hazards model,
defining all terms that you use. [3]
(ii) State the regression parameters for the fitted model. [2]
(iii) Describe the class of children to which the baseline hazard applies. [1]
(iv) Discuss the suggestion that the new tuition method has improved the chances
of children continuing to play their instrument. [3]
(v) Calculate, using the results from the model, the probability that a boy will still
be playing the piano after 4 years if provided with the new tuition method,
given that the probability that a girl will still be playing the trumpet after 4
years following the traditional method is 0.7. [3]
[Total 12]
CT4 A2008—4
9 An investigation into the mortality of patients following a specific type of major
operation was undertaken. A sample of 10 patients was followed from the date of the
operation until either they died, or they left the hospital where the operation was
carried out, or a period of 30 days had elapsed (whichever of these events occurred
first). The data on the 10 patients are given in the table below.
1 2 Died
2 6 Died
3 12 Died
4 20 Left hospital
5 24 Left hospital
6 27 Died
7 30 Study ended
8 30 Study ended
9 30 Study ended
10 30 Study ended
(i) State whether the following types of censoring are present in this
investigation. In each case give a reason for your answer.
(a) Type I
(b) Type II
(c) Random [3]
(ii) State, with a reason, whether the censoring in this investigation is likely to be
informative. [1]
(iii) Calculate the value of the Kaplan-Meier estimate of the survival function at
duration 28 days. [5]
(iv) Write down the Kaplan-Meier estimate of the hazard of death at duration 8
days. [1]
The probabilities that the customer is online and offline at time t are PON(t) and
POFF(t) respectively.
(ii) ′ (t ) .
Write down Kolmogorov’s forward equation for POFF [2]
(iii) Solve the equation in part (ii) to obtain a formula for the probability that a
customer is offline at time t, given that they were offline at time 0. [3]
(iv) Calculate the expected proportion of time spent online over the period [0,t].
[HINT: Consider the expected value of an indicator function which takes the
value 1 if offline and 0 otherwise.]
[4]
CT4 A2008—6
11 An investigation was carried out into the relationship between sickness and mortality
in an historical population of working class men. The investigation used a three-state
model with the states:
1 Healthy
2 Sick
3 Dead
Let the probability that a person in state i at time x will be in state j at time x+t be
t p x . Let the transition intensity at time x+t between any two states i and j be μ x +t .
ij ij
(i) Draw a diagram showing the three states and the possible transitions between
them. [2]
∂
t p x = t p x μ x +t + t p x μ x +t .
23 21 13 22 23
[5]
∂t
(iii) Write down the likelihood of the data in the investigation in terms of the
transition rates and the waiting times in the Healthy and Sick states, under the
assumption that the transition rates are constant. [3]
(iv) Derive the maximum likelihood estimator of the transition rate from Sick to
Dead. [3]
(a) the value of the constant transition rate from Sick to Dead
(b) 95 per cent confidence intervals around this transition rate
[4]
[Total 17]
END OF PAPER
CT4 A2008—7
Faculty of Actuaries Institute of Actuaries
EXAMINERS’ REPORT
April 2008
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
M A Stocker
Chairman of the Board of Examiners
June 2008
© Faculty of Actuaries
© Institute of Actuaries
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
Comments
Comments on solutions presented to individual questions for this April 2008 paper are given
below.
Question 2 Answers to this question were disappointing. In part (a) many candidates did
not realise that smoothness is automatically ensured when graduating with a
parametric formula with a small number of parameters. In part (b) many
candidates presented descriptions of the method of graphical graduation,
rather than answering the question which was set.
Question 3 Most candidates scored reasonably well on part (i), but few candidates could
state the conditions required for a compound Poisson process to be a Poisson
process in part (ii).
Question 4 A reasonable attempt was made at this bookwork question by most candidates,
although few made sufficient distinct points to score close to full marks.
Question 5 This exposed-to-risk question was quite well answered by many candidates,
who correctly identified the rate interval and the appropriate census-type
formula. An encouraging number of candidates also recognised the need to
adjust the age definition in order to ensure correspondence between the first
marriages data and the exposed-to-risk data.
Question 6 Many candidates scored well on this question. Common errors were failure to
use (or incorrect use of) the continuity correction in the normal approximation
to the signs test; calculating only the probability of 18 positive signs (rather
than the probability of 18 or more signs) when using the exact binomial
computation of the signs test; and calculating only the probability of 2 positive
runs (rather than the probability of 2 or fewer positive runs) when using the
exact computation of the grouping of signs test.
Question 7 Only a small proportion of candidates correctly answered part (i). In part (ii)
a very large number of candidates adopted a three-state solution to this
problem, with state space {A, B, C}. Partial credit was given for this, and also
for correctly following this three-state solution through in part (iii) to obtain
the steady-state proportions of 3/11, 2/11 and 6/11 using auditors A, B and C
respectively.
Question 8 This question was not as well answered as some others. Some candidates
failed to write the numerical values of the estimated parameters down in part
(ii). There were few correct attempts at part (v). Many candidates simply
calculated the ratio between the two hazards, which is incorrect. Others made
unnecessary assumptions about the form of the baseline hazard (e.g. that it
was constant).
Page 2
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
Question 9 This straightforward calculation of the survival function was very well
answered, apart from part (iv), in which only a handful of candidates realised
that the Kaplan-Meier estimate of the hazard at any duration at which no
event is observed to take place is 0. Given that the Kaplan-Meier estimate of
the hazard is a step function, it is clear than this must be so. It was very
encouraging to see the high proportion of sensible answers to part (ii). Credit
was given in part (ii) to candidates who stated that the censoring was non-
informative provided that the reason given was consistent with this statement.
Question 10 Few candidates scored highly on this question. Many candidates got no
further than part (ii). Although there were a fair number of attempts to solve
the differential equation in part (iii), only a minority of candidates spotted that
PON (t ) + POFF (t ) = 1 .
Question 11 This question was very well answered. Many candidates provided
substantially correct answers to all parts, losing marks only for failure to
include certain details in part (ii) (for example that we need to condition on
the state occupied at time x+t); or for failing to point out that we need to
substitute the estimated values from the data into the formula for the variance
of μ 23 in part (v).
Page 3
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
1 Sex
Age
Type of policy
Smoker/non-smoker
Level of underwriting
Duration in force
Sales channel
Policy size
Known impairments
Occupation
Page 4
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
N0 = 0
Nt has independent increments
Nt has stationary increments
[λ (t − s )]n e−λ (t − s )
P[ Nt − N s = n] = s < t, n = 0, 1, 2…..
n!
OR
N0 = 0
P[ Nt + h − Nt = 1] = λh + o(h)
P[ Nt + h − Nt = 0] = 1 − λh + o(h)
P[ Nt + h − Nt ≠ 0,1] = o(h)
Nt
X t = ∑ Yi
i =1
Page 5
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
4 Benefits
In a model of a complex system we can usually get much better control over the
experimental conditions so that we can reduce the variance of the results output
from the model without upsetting their mean values
Avoids costs and risks of making changes in the real world, so we can study
impact of changing inputs before making decisions.
Limitations
Models can look impressive when run on a computer so that there is a danger that one
gets lulled into a false sense of confidence.
If a model has not passed the tests of validity and verification its impressive
output is a poor substitute for its ability to imitate its corresponding real world
system.
Models rely heavily on the data input. If the data quality is poor or lacks credibility
then the output from the model is likely to be flawed.
It is important that the users of the model understand the model and the uses to which
it can be safely put. There is a danger of using a model as a black box from which it is
assumed that all results are valid without considering the appropriateness of using that
model for the particular data input and the output expected.
It is not possible to include all future events in a model. For example a change in
legislation could invalidate the results of a model, but may be impossible to
predict when the model is constructed.
It may be difficult to interpret some of the outputs of the model. They may only be
valid in relative, rather than absolute, terms. For example comparing the level of risk
of the outputs associated with different inputs.
Page 6
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
A definition of the population data which is compatible with these data on first
marriages is
Px,t = number of lives under observation at time t since the start of the
investigation who were aged x next birthday on the 1 January
immediately preceding t
Since we follow each cohort of lives through each calendar year, this exposed
to risk is
N
E xc = ∫ Px,t dt
0
N −1
∑ 2 ( Px,t + Px+1,t +1)
1
E xc =
0
Px,t = Px-1,t*
K +N
∑
1
E xc = ( Px −1,t * + Px,t +1* ) .
t =K 2
Page 7
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
(iii) The age range at the start of the rate interval is (x–1, x)
exact.
EITHER
17.5 − 12.5 −5
= = −2.00 .
6.25 2.5
Page 8
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
OR
⎛ 25 ⎞ 25
Pr[j positive deviations] = ⎜ ⎟ 0.5 .
⎝j ⎠
This is equal to
= 0.02164.
THEN EITHER
Page 9
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
OR
⎛17 ⎞ ⎛ 8 ⎞
⎜ ⎟⎜ ⎟
Pr[1 positive run] = ⎝ ⎠ ⎝ ⎠ =
0 1 8
= 0.0000166
⎛ 25 ⎞ 480, 700
⎜ ⎟
⎝18 ⎠
⎛17 ⎞ ⎛ 8 ⎞
⎜ ⎟⎜ ⎟
Pr[2 positive runs] = ⎝ ⎠ ⎝ ⎠ =
1 2 (17)(28)
= 0.000990
⎛ 25 ⎞ 480, 700
⎜ ⎟
⎝ 18 ⎠
OR
⎛ (18)(8) [(18)(7)]2 ⎞
N⎜
⎜ 25
, ⎟ = N ( 5.76,1.02 )
⎝ (25)3 ⎟⎠
2 − 5.76
= −3.722 .
1.02
Page 10
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
∞
= ∑ probability ith audit takes place prior to changing auditors
i =1
= 1 + 1 + 0.8 + 0.82+0.83+……..
= 1 + 1/(1-0.8) = 6
Transition matrix A is
AL A BL B CL C
AL 0 1 0 0 0 0
A 0 0.8 0.1 0 0.1 0
BL 0 0 0 1 0 0
B 0.15 0 0 0.7 0.15 0
CL 0 0 0 0 0 1
C 0.05 0 0.05 0 0 0.9
π = πA
Page 11
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
Combining (1) and (2), (3) and (4), and (5) and (6)
π A = 1.5π B
πC = 3π B
∑ πi = 1 implies
i
(1.5 + 0.3 + 1 + 0.3 + 3 + 0.3)π B = 1
⎛ π AL ⎞ ⎛ 0.046875 ⎞
⎜ ⎟
⎜ π A ⎟ ⎜⎜ 0.234375 ⎟⎟
⎜ ⎟ ⎜
π 0.046875 ⎟
So ⎜ BL ⎟ = ⎜
⎜ π ⎟ ⎜ 0.15625 ⎟⎟
⎜ B ⎟ ⎜
⎜ πCL ⎟ ⎜ 0.046875 ⎟⎟
⎜⎜ ⎟⎟ ⎜ 0.46875 ⎟⎠
⎝ πC ⎠ ⎝
Page 12
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
8 (i) h( z, t ) = h0 (t ).exp(β.ziT )
a female,
following traditional tuition method,
playing the piano
However the 95% confidence interval for the parameter spans zero.
So at the 5% significance level it is not possible to conclude that the new
tuition method has improved the chances of children continuing to play their
instrument.
(v) The hazard for a girl being taught the trumpet by the traditional
method giving up is h0 (t ) exp(0.14) .
⎛ 4 ⎞ ⎛ 4
⎞
S female (4) = exp ⎜ − ∫ h0 (t ) exp(0.14)dt ⎟ = exp ⎜ −1.150274 ∫ h0 (t )dt ⎟
⎝ 0 ⎠ ⎝ 0 ⎠
Page 13
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
⎛ 4
⎞
exp ⎜ −1.150274 ∫ h0 (t )dt ⎟ = 0.7, so that
⎝ 0 ⎠
4
log e 0.7 = −1.150274 ∫ h0 (t )dt ,
0
4
log e 0.7
and hence ∫ h0 (t )dt =
0
−1.150274
= 0.310078.
The hazard of giving up for a boy taught the piano by the new
method is h0 (t ) exp(−0.05 + 0.02) = h0 (t ) exp(−0.03).
⎛ 4 ⎞
S male (4) = exp ⎜ − ∫ h0 (t ) exp(−0.03)dt ⎟ = exp [ −0.310078(0.970446) ]
⎝ 0 ⎠
ALTERNATIVELY
The hazard of giving up for a girl being taught the trumpet by the
traditional method is h0 (t ) exp( β 2 ) .
⎛ 4 ⎞ ⎛ 4
⎞
S female (4) = exp ⎜ − ∫ h0 (t ) exp( β 2 )dt ⎟ = exp ⎜ − exp( β 2 ) ∫ h0 (t )dt ⎟
⎝ 0 ⎠ ⎝ 0 ⎠
and hence
4
log e [ S female (4)]
∫ h0 (t )dt =
0
− exp β 2
= − exp(− β 2 ) log e [ S female (4)] .
Page 14
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
⎛ 4
⎞
S male (4) = exp ⎜ − exp( β3 + β 4 ) ∫ h0 (t )dt ⎟ .
⎝ 0 ⎠
4
Substituting for ∫ h (t )dt produces
0
0
(ii) Yes
dj dj dj ^
tj nj dj cj
nj
1-
nj
∏1 − nj
= S (t )
t j ≤t
0 10
2 10 1 0 1/10 9/10 9/10 = 0.9
6 9 1 0 1/9 8/9 8/10 = 0.8
12 8 1 2 1/8 7/8 7/10 = 0.7
27 5 1 4 1/5 4/5 14/25 = 0.56
Page 15
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
1
0.9
0.8
0.7
0.6
S(t)
0.5
0.4
0.3
0.2
0.1
0
0 10 20 30 40
Duration t (days)
PON (t ) + POFF (t ) = 1
′ (t ) + POFF (t ) = 0.8
POFF
so that
Page 16
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
t t
E[Ot POFF (0) = 1] = ∫ E[ I s POFF (0) = 1]ds = ∫ POFF ( s )ds
0 0
t t
−s −t t
∫ POFF (s)ds = ∫ (0.8 + 0.2e )ds = 0.8t − 0.2e = 0.8t + 0.2(1 − e−t )
0
0 0
0.2
Percentage online
0
t=0 Time
Page 17
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
11 (i)
1 Healthy 2 Sick
3 Dead
But dt p33
x +t = 1, so
dt px23+t = μ 23
x +t dt + o(dt ) and dt p13
x +t
=
= μ13
x +t dt + o(dt )
o(dt )
where o(dt ) is defined such that lim =0.
dt →0 dt
t + dt px23 = t px22 [μ 23
x +t dt + o( dt )] + t p x [μ x +t dt + o( dt )] + t p x ,
21 13 23
d t p x − t + dt p x
t p x = lim = t px21μ13
x +t + t p x μ x +t
23 22 23
dt dt →0 dt
Page 18
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
12 21 13 23
exp[(−μ12 − μ13 )v1 ]exp[(−μ 23 − μ 21 )v 2 ](μ12 ) d (μ 21 ) d (μ13 ) d (μ 23 ) d
d log L d 23
= −v 2 + .
d μ 23 μ 23
^
d 23
μ 23 = .
v2
d 2 (log L) d 23
This is a maximum because =−
(d μ 23 ) 2 (μ 23 ) 2
Page 19
Subject CT4 — Models Core Technical — April 2008 — Examiners’ Report
^
Approximating μ 23 by μ 23 and E[V ] by v 2 we
0.2857
estimate for the variance as = 0.00204 .
140
or (0.1972, 0.3742).
Page 20
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 12 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the Formulae
and Tables and your own electronic calculator from the approved list.
© Faculty of Actuaries
CT4 S2008 © Institute of Actuaries
1 You work for a consultancy which has created an actuarial model and is now
preparing documentation for the client.
List the key items you would include in the documentation on the model. [4]
Give four examples, one of each type, of stochastic models which may be used to
model observed processes, and suggest a practical problem to which each model may
be applied. [4]
3 Compare the advantages and disadvantages of the Binomial and the multiple-state
models in the following situations:
4 In the village of Selborne in southern England in the year 1637 the number of babies
born each month was as follows
January 2 July 5
February 1 August 1
March 1 September 0
April 2 October 2
May 1 November 0
June 2 December 3
Data show that over the 20 years before 1637 there was an average of 1.5 births per
month. You may assume that births in the village historically follow a Poisson
process.
An historian has suggested that the large number of births in July 1637 is unusual.
(i) Carry out a test of the historian’s suggestion, stating your conclusion. [4]
(ii) Comment on the assumption that births follow a Poisson process. [1]
[Total 5]
CT4 S2008—2
5 An investigation into the mortality experience of a sample of the male student
population of a large university has been carried out. The university authorities wish
to know whether the mortality of male students at the university is the same as that of
males in the country as a whole. They have drawn up the following table.
18 13 10
19 15 12
20 14 14
21 20 12
22 12 8
23 8 5
Carry out an overall test of the university authorities’ hypothesis, stating your
conclusion. [5]
age last birthday at the most recent policy anniversary prior to the portfolio
transfer + number of premiums received by insurer B.
(i) (a) State the rate interval implied by the above data.
(b) Write down the range of ages at the start of the rate interval. [2]
(ii) Give an expression which can be used to estimate the initial exposed to risk at
age x, Ex, stating any assumptions made. [2]
x dx ∑ Px,t ∑ Px,t +1
39 28 10,536 11,005
40 36 10,965 10,745
41 33 10,421 10,577
Suppose N identical and independent lives are observed from age x exact for one year
or until death if earlier.
Define:
CT4 S2008—4
Ci to be a random variable denoting the number of claims arising from lives with i
policies.
(ii) Derive an expression for the ratio of the variance of the number of claims
arising compared with that if each policy covered an independent life. [4]
(iii) Explain how the expression derived in (ii) could be used in practice. [2]
[Total 8]
8 A No-Claims Discount system operated by a motor insurer has the following four
levels:
Level 1: 0% discount
Level 2: 25% discount
Level 3: 40% discount
Level 4: 60% discount
• Following a year with no claims, move to the next higher level, or remain at
level 4.
• Following a year with one claim, move to the next lower level, or remain at
level 1.
• Following a year with two or more claims, move down two levels, or move to
level 1 (from level 2) or remain at level 1.
For a given policyholder in a given year the probability of no claims is 0.85 and the
probability of making one claim is 0.12.
(i) Write down the transition matrix of this No-Claims Discount process. [1]
(ii) Calculate the probability that a policyholder who is currently at level 2 will be
at level 2 after:
(i) Describe the nature of the state space and time space for this process. [2]
(ii) Draw and label a transition diagram indicating appropriate transitions between
the states. [2]
t p1ix the probability that a life is in state i at age x+t, given they are in state 1 at age
x
(iii) Write down equations which could be used to determine the evolution of t p1ix
(for each i) appropriate for:
CT4 S2008—6
10 In an investigation of reconviction rates among those who have served prison
sentences, let X be a random variable which measures the duration from the date of
release from prison until the ex-prisoner is convicted of a subsequent offence. The
investigation monitored a sample of 100 ex-prisoners (who were all released on the
same date) at one-monthly intervals from their date of release for a period of 6
months. Those who could not be traced in any month were removed from the sample
at that point and not traced in subsequent months. Reconviction was assumed to take
place at the duration that a prisoner was first known to have been reconvicted.
(i) Express the hazard rate at duration x months in terms of probabilities. [1]
The investigation produced the following data for a sample of 100 ex-prisoners.
1 100 0
2 97 0
3 95 4
4 90 3
5 85 5
6 80 0
(iii) Estimate confidence intervals around the integrated hazard using the results
from part (ii) to test the hypothesis that the rate of reconviction has declined
since the previous investigation. [6]
[Total 12]
(i) Write down the state space and transition graph of the sequence Xn. [2]
(a) is aperiodic.
(b) is reducible.
(c) admits a stationary distribution. [3]
(iii) Derive an expression for the number of upward movements in the sequence Xn
between t and (t + m) if Xt = i and Xt+m = j. [2]
(iv) Derive expressions for the m-step transition probabilities pij( m) . [3]
(v) Show how the one-step transition probabilities would alter if Xn was restricted
to non-negative numbers by introducing:
CT4 S2008—8
12 (i) Explain the meaning of the rates of mortality usually denoted qx and mx , and
the relationship between them. [3]
(iii) Calculate mx under each of the assumptions (a) to (c) above. [8]
END OF PAPER
CT4 S2008—9
Faculty of Actuaries Institute of Actuaries
EXAMINERS’ REPORT
September 2008
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
R D Muckart
Chairman of the Board of Examiners
November 2008
© Faculty of Actuaries
© Institute of Actuaries
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
Comments
Comments on solutions presented to individual questions for the September 2008 paper are
given below.
Q1 This standard bookwork question was fairly well answered. Some candidates simply
wrote down a list of steps in the development of the model, rather than answering the
question that was set.
Q2 This straightforward question was well answered. Some candidates were vague
about emphasising that continuous time models are applied to problems which
require continuous monitoring.
Q3 Answers to this question were very poor. Many candidates did not go beyond making
the point that the Binomial model is hard to extend to multiple decrements, whereas
the multiple state model extends quite naturally.
Q4 Only a minority of candidates answered this question using the approach intended.
Many tried to do a chi-squared test comparing the observed and expected numbers of
births. This received some credit, especially when candidates combined the months
into half-years, or thirds of a year, before performing the chi-squared test, so that the
expected values in each cell were greater than 5.
Q5 This straightforward question was very well answered. The most common error was
in reducing the number of degrees of freedom below 6. This is incorrect in this case,
because the comparison is between an observed experience and a pre-existing
experience, not between crude rates and graduated rates.
Q7 Part (ii) of this question was standard bookwork, but was nevertheless answered in a
brief or cursory fashion by many candidates. On the other hand, a good number of
candidates were able to make the points required in part (i).
Q8 This question was very well answered, with many candidates scoring full marks.
Some candidates were penalised in part (iii) for simply calculating the stationary
distribution and not stating explicitly which of the numbers represented the long-run
probability of being in discount level 2.
Q9 Answers to this question were disappointing, especially to part (iii). In part (ii),
candidates who included additional transitions between states 2 and 3, and between
states 2 and 1, were not penalised. However, such candidates were expected to
produce answers to part (iii) which were consistent with the transition diagram they
had drawn in part (ii).
Q10 Part (i) of this question was very disappointingly answered, as the required definition
is in the Core Reading. Most candidates were able to compute the estimated survival
function in part (ii). Some candidates interpreted the question as meaning that the
Page 2
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
numbers contacted at any duration include those reconvicted prior to that duration,
so that those reconvicted must be subtracted from those contacted to obtain the
relevant nj. These candidates were given credit for part (ii). In part (iii) many
candidates correctly calculated the variance of the integrated hazard but then incorrectly
used this variance to compute a confidence interval around the survival function, rather than
first computing the confidence interval around the integrated hazard and then using the
formula S(x) = exp(−Λx) to convert this into a confidence interval around S(6).
Q11 Answers to this question were disappointing. Many candidates were able to answer
parts (i) and (ii) reasonably well, but made little or no attempt at the remaining
sections.
Q12 Answers to this question varied widely, but overall were disappointing. There was a
large variation by centre, with average scores for some centres being several marks
higher than for other centres. Perhaps this is the result of different training and
education materials being used in different locations? While most candidates could
q
write down the formula mx = 1 x and the formulae required to answer part (ii), it
∫ t px dt
0
was clear from the answers to parts (iii) and (iv) that understanding of what these
formulae mean was very shaky.
Page 3
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
General random walk, ARIMA time series model, moving average model
Share price at end of each day
Page 4
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
The Binomial model also does not use all the information
available if exact ages at entry into and exit from
observation are available.
The estimators have the same form and the same statistical
properties as in the classic life table.
4 (i) Suppose that the number of births each month, B, is the outcome of a Poisson
process with a rate λ = 1.5.
Therefore we have
b Pr[ B = b]
0 0.223
1 0.335
2 0.251
3 0.126
4 0.047
5 0.014
6+ 0.004
Page 5
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
But only July has a number of births more than 5, and at the 5% level of
statistical significance we expect 1 month in 20 to have such a large
number, then unless we have a prior expectation that July is unusual, we
should be cautious before accepting the historian’s suggestion.
Age x zx zx2
18 0.9487 0.9
19 0.8660 0.75
20 0 0
21 2.3094 5.3333
22 1.4142 2
23 1.3416 1.8
Page 6
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
4
1 3
(ii) Central exposed to risk E xc = ∫ Px,t dt ≈ ∑ ( Px,t + Px,t +1 )
t =0
2 t =0
1 3 1
Initial exposed to risk Ex ≈ ∑
2 t =0
( Px,t + Px,t +1 ) + d x ,
2
assuming deaths are uniform over the rate interval OR deaths occur on
average half way through the rate interval.
(but NOT deaths are uniform over the “year”, or occur on average half
way through the “year”)
dx
(iii) qˆ x = estimates qx for the average age
Ex
at the start of the rate interval.
28
qˆ39 = = 0.002596
1 1
(10536 + 11005) + * 28
2 2
36
qˆ40 = = 0.003311
1 1
(10965 + 10745) + *36
2 2
Page 7
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
7 (i) Individual life offices are likely to have their systems set up to provide
information on a “by policy” basis.
(ii) If the mortality rate is qx then since the lives are independent the number of
deaths Di will be distributed Binomial ( q x , πi N )
So ∑ Ci = ∑ iDi .
i i
⎡ ⎤ ⎡ ⎤
Hence Var[C] = Var ⎢ ∑ Ci ⎥ = Var ⎢ ∑ iDi ⎥ = ∑ i 2 Var [ Di ]
⎣⎢ i ⎦⎥ ⎣⎢ i ⎦⎥ i
by independence of deaths
= ∑ i 2 πi Nqx (1 − q x )
i
∑ i 2 πi
So the variance is increased by the ratio i
∑ iπi
i
(iii) If the proportions of lives holding i policies were known, the variance ratio
could be allowed for in statistical tests
Page 8
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
⎛ 0.15 0.85 0 0 ⎞
⎜ ⎟
0.15 0 0.85 0 ⎟
P= ⎜
⎜ 0.03 0.12 0 0.85 ⎟
⎜ ⎟
⎝0 0.03 0.12 0.85 ⎠
(ii) (a) For the one year transition, p22 = 0, as can be seen
from above (or is obvious from the statement).
This reads
∑i=1 πi = 1 .
4
Discard the first of these equations and use also the fact that
Page 9
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
17
0.85π2 + 0.12 × π3 = π3 ⇒ π3 = 2.65625π 2
3
1
0.85π1 = π2 − 0.12π3 − 0.03π4 = π3 − 0.12π3 − 0.17π3 = 0.0865π3,
2.65625
or that π3 = 0.13996.
9 (i) The state space is discrete with states as given in the question.
Page 10
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
(ii)
2
μ12
x +t μ 24
x +t
No longer
employed
μ14
x +t
4
1
Currently μ13
x +t Pensioner
employed μ 25
x +t
3
Ill health
μ 45
x +t
μ 35
x +t
μ15
x +t
Dead
∂ p11 12 13 15 11
∂t t x = −(μ x +t + μ x +t + μ x +t ) t px
∂ p12 12 11 25 12
∂t t x = μ x +t . t p x − μ x + t . t p x
∂ p13 13 11 35 13
∂t t x = μ x +t . t p x − μ x +t . t p x
∂ p15 15 11 25 12 35 13
∂t t x = μ x +t . t p x + μ x +t . t p x + μ x +t . t p x
and t p14
x is zero.
(b) For x + t = 65
t p11 12
x and t p x become 0 at x + t = 65+ δ
t +δ p14
x = t −δ p11 12
x + t −δ p x
Page 11
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
t p11 12
x = t px = 0
∂ p13 35 13
∂t t x = −μ x +t . t p x
∂ p14 45 14
∂t t x = −μ x +t . t p x
∂ p15 35 13 45 14
∂t t x = μ x +t . t p x + μ x +t . t p x
10 (i) EITHER
Pr[ X ≤ x + dt | X > x]
lim .
dt →0 dt
OR
OR
1 d
The hazard rate at duration x is given by h( x) = − [ S ( x)] ,
S ( x) dx
dj dj
xj nj dj cj
nj
Λx = ∑ nj
x j ≤x
0 100 0 0 0 0
1 100 0 3 0 0
2 97 0 2 0 0
3 95 4 1 4/95 = 0.0421 0.0421
4 90 3 2 3/90 = 0.0333 0.0754
5 85 5 0 5/85 = 0.0588 0.1343
6 80 0 80 0 0.1343
Page 12
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
x S(x)
0≤ x<3 1.0000
3≤ x < 4 0.9588
4≤ x<5 0.9274
5≤ x 0.8744
~ d j (n j − d j )
Var[ Λ x ] = ∑ n3j
x j ≤x
d j (n j − d j ) d j (n j − d j )
xj nj dj
n3j
∑ n3j
x j ≤x
0 100 0 0 0
1 100 0 0 0
2 97 0 0 0
3 95 4 0.000425 0.000425
4 90 3 0.000358 0.000783
5 85 5 0.000651 0.001434
6 80 0 0 0.001434
^ ^
Λ 6 ± 1.96 var Λ 6
= (0.0601, 0.2085).
Page 13
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
THEN EITHER
^
The estimated survival function, S ( x) is given
^
by exp(− Λ x ) ,
^
so that the 95 per cent confidence interval for S ( x) is
[exp(−0.0601), exp(−0.2085)]
^
Since the 95 per cent confidence interval around S ( x) in the current
^
investigation does not include the value 0.8, and our estimate of S ( x) > 0.8
we conclude that the rate of reconviction has declined since the previous
investigation.
OR
Λ6 = −loge(0.8) = 0.2231.
Since this is higher than the upper limit in the range (0.0601, 0.2085) we
conclude that the rate of reconviction has declined since the previous
investigation.
Page 14
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
Transition graph:
p p p p
-2 -1 0 1 2
1-p 1-p 1-p 1-p
(c) No stationary distribution will exist because the state space is infinite.
and u – (m – u) = j – i
m+ j −i
So u = .
2
⎛m⎞
These can be ordered in ⎜ ⎟ ways.
⎝u⎠
Page 15
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
⎧⎛ m ⎞ u
⎪ p (1 − p) m−u if j − i ≤ m and m + j − i even
pij( m) = ⎨⎜⎝ u ⎟⎠
⎪ 0 otherwise
⎩
(v) EITHER
OR
Reflecting:
⎛ 0 1 0 0 0 ... ⎞
⎜ ⎟
⎜1 − p 0 p 0 0 ... ⎟
⎜ 0 1− p 0 p 0 ... ⎟
⎜ ⎟
⎜ 0 0 1− p 0 p ... ⎟
⎜ 0 0 0 1− p 0 ... ⎟
⎜⎜ ⎟⎟
⎝ : : : : : ⎠
Absorbing:
⎛ 1 0 0 0 0 ... ⎞
⎜ ⎟
⎜1 − p 0 p 0 0 ... ⎟
⎜ 0 1− p 0 p 0 ... ⎟
⎜ ⎟
⎜ 0 0 1− p 0 p ... ⎟
⎜ 0 0 0 1− p 0 ... ⎟
⎜⎜ ⎟⎟
⎝ : : : : : ⎠
OR
Page 16
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
Reflecting
1 p
0 1 2
1-p 1-p
Absorbing:
p
1
0 1 2
1-p 1-p
(vi) In both cases the zero transition probabilities remain
zero as the period is still 2 where relevant.
Otherwise
Page 17
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
12 (i) qx is the probability that a life aged exactly x will die before reaching
exact age x+1, and is called the initial rate of mortality.
mx is called the central rate of mortality and represents the probability that a
life alive between the ages of x and x+1 dies
qx
mx = 1
∫ t px dt
0
t qx = t * qx
t qx = 1 − e −μ*t
1−t q x + t = (1 − t ) * q x
1
1 1 ⎡t2 ⎤
∫t x ∫
p dt = (1 − 0.1t ) dt = 1 − 0.1⎢ ⎥ = 0.95
⎢⎣ 2 ⎥⎦
0 0 0
mx = 0.1/0.95 = 0.105263
(b) CFM
μ given by:
1 − e −μ = 0.1
μ = − ln 0.9 = 0.1053605
Page 18
Subject CT4 — Models Core Technical — September 2008 — Examiners’ Report
EITHER
so mx = 0.1053605
OR
1 1 1
))dt = − ⎡ e ⎤
−μt 1 −μt 1
∫ t px dt = ∫ (1 − (1 − e = (1 − e−μ ) = 0.949122
μ⎣ ⎦0 μ
0 0
mx = 0.1/0.949122=0.1053605
(c) Balducci
So
1 px 0.9 0.9
t p x= = =
1−t p x +t 1 − 1−t q x +t 0.9 + 0.1t
1 1
0.9 0.9
∫ t px dt = ∫ 0.9 + 0.1t dt = 0.1 [ln(0.9 + 0.1t )]0 = −9 ln 0.9 = 0.9482446
1
0 0
So mx = 0.1/0.9482446=0.1054580
Page 19
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 12 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the Formulae
and Tables and your own electronic calculator from the approved list.
© Faculty of Actuaries
CT4 A2009 © Institute of Actuaries
1 A life insurance company has a small group of policies written on impaired lives and
has conducted an investigation into the mortality of these policyholders. It is
proposed that the crude mortality rates be graduated for use in future premium
calculations.
Discuss the suitability of two methods of graduation that the insurance company could
use. [3]
⎛1 − a a ⎞
⎜ ⎟
⎝ b 1− b ⎠
(ii) Explain the range of values that a and b can take which result in this being a
valid Markov chain which is:
(a) irreducible
(b) periodic [3]
[Total 5]
30 98,617
40 97,952
(i) Calculate 5 q30 under each of the two following alternative assumptions:
(ii) Calculate the number of survivors to exact age 35 years out of 100,000 births
under each of the assumptions in (i) above. [1]
English Life Table 15 (females) was originally calculated using data classified by
single years of age. The number of survivors to exact age 35 years was 98,359.
CT4 A2009—2
5 Explain the basis underlying the grouping of signs test, and derive the formula for the
probability of exactly t positive groups by considering the possible arrangements of a
set of positive and negative signs. [5]
• number of persons who recovered from the operation during the month (defined as
being discharged from the hospital) classified by the month of their operation.
On the first day of each month from January 2008 to January 2009, the hospital listed
all in-patients who were yet to recover from this operation, classified according to the
length of time elapsing since their operation, to the nearest month.
(i) (a) Write down an expression which will enable the hospital to calculate
rates of recovery, rx, during 2008 at various durations x since the
operation using the available data.
(b) Derive a formula for the exposed to risk based on the information in
the hospital’s monthly lists of in-patients which corresponds to the data
on recovery from the operation.
[5]
(ii) Determine the value of f such that the expression in (i)(a) applies to an actual
duration x + f since the operation. [2]
[Total 7]
7 (i) Explain how the classification of stochastic processes according to the nature
of their state space and time space leads to a four way classification. [2]
• a food retailer
• a general insurance company
[6]
[Total 8]
(i) If the number of cats currently infected is x, explain why the number of
possible pairings of cats which could result in a new flea infection is x(10 – x).
[1]
(ii) Show how the number of infected cats at any time, X(t), can be formulated as
a Markov jump process, specifying:
(iii) State the distribution of the holding times of the Markov jump process. [2]
(iv) Calculate the expected time until all the cats have fleas, starting from a single
flea-infected cat. [2]
[Total 9]
9 (i) Prove that, under Gompertz’s Law, the probability of survival from age x to
age x + t, t p x , is given by:
c x ( ct −1)
⎡ ⎛ − B ⎞⎤
t p x = ⎢ exp ⎜ ⎟⎥ . [3]
⎣ ⎝ ln c ⎠ ⎦
1 p50 = 0.995
2 p50 = 0.989 .
(iii) Comment on the calculation performed in (ii) compared with the usual process
for estimating the parameters from a set of crude mortality rates. [3]
[Total 9]
CT4 A2009—4
10 Let Tx be a random variable denoting future lifetime after age x, and let T be
another random variable denoting the lifetime of a new-born person.
( )
S x (t ) = exp −(λt )β ,
(iii) Derive an expression for the Weibull force of mortality in terms of λ and β.
[3]
(iv) Sketch, on the same graph, the Weibull force of mortality for 0 ≤ t ≤ 5 for the
following pairs of values of λ and β:
λ = 1, β = 0.5
λ = 1, β = 1.0
λ = 1, β = 1.5
[4]
[Total 10]
1 Alive
μ12
x +t
μ14
x +t
μ13
x +t
∂ 12 12 11
t p x = μ x +t t p x . [5]
∂t
The investigation was carried out separately for each year of age, and the transition
intensities were assumed to be constant within each single year of age.
(ii) (a) Write down, defining all the terms you use, the likelihood for the
transition intensities.
The investigation produced the following data for persons aged 64 last birthday:
(iii) (a) Calculate the maximum likelihood estimate (MLE) of the force of
mortality from heart disease at age 64 last birthday.
(b) Estimate an approximate 95% confidence interval for the MLE of the
force of mortality from heart disease at age 64 last birthday. [3]
(iv) Discuss how you might use this model to analyse the impact of risk factors on
the death rate from heart disease and suggest, giving reasons, a suitable
alternative model. [3]
[Total 16]
CT4 A2009—6
12 A motor insurer operates a no claims discount system with the following levels of
discount {0%, 25%, 50%, 60%}.
The rules governing a policyholder’s discount level, based upon the number of claims
made in the previous year, are as follows:
• Following a year with no claims, the policyholder moves up one discount level, or
remains at the 60% level.
• Following a year with one claim, the policyholder moves down one discount level,
or remains at 0% level.
• Following a year with two or more claims, the policyholder moves down two
discount levels (subject to a limit of the 0% discount level).
(i) Determine the transition matrix for the no claims discount system. [3]
(iii) Calculate the expected average long term level of discount. [1]
The following data shows the number of the insurer’s 130,200 policyholders in the
portfolio classified by the number of claims each policyholder made in the last year.
This information was used to estimate the mean of 0.30.
No claims 96,632
One claim 28,648
Two claims 4,400
Three claims 476
Four claims 36
Five claims 8
(iv) Test the goodness of fit of these data to a Poisson distribution with mean 0.30.
[5]
(v) Comment on the implications of your conclusion in (iv) for the average level
of discount applied. [2]
[Total 16]
END OF PAPER
CT4 A2009—7
Faculty of Actuaries Institute of Actuaries
EXAMINERS’ REPORT
April 2009
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
R D Muckart
Chairman of the Board of Examiners
June 2009
© Faculty of Actuaries
© Institute of Actuaries
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
Comments
Comments on solutions presented to individual questions for the April 2009 paper are given
below.
Q1 Answers to this question were satisfactory. Most candidates realised that graduation by
reference to a standard table was potentially appropriate, and that graphical graduation
might have to be used as a last resort. Credit was given for sensible points other than those
mentioned in the specimen solution below.
Q2 In part (ii) some explanation of the correct possible values of a and b was required for
full credit. A common error in part (ii) (a) was to write 0 < a < 1 and 0 < b < 1, ignoring
the possibility that a and b could equal 1.
Q3 This bookwork question was well answered by many candidates. Credit was given for
sensible points other than those mentioned in the specimen solution below.
Q4 Answers to parts (i) and (ii) were generally good, with a substantial proportion of
candidates scoring full marks. Part (iii) was much less convincingly answered. Although not
all the points mentioned in the specimen solutions below were required for full credit, many
candidates only included the briefest of comments, and consequently scored few marks.
Q5 Most candidates simply wrote down the formula for Pr[G=t] (which is given in the book
of Formulae and Tables) and then explained what each bracketed expression in the formula
meant. Few candidates gave more than the briefest explanation of why the test is useful, and
what it is designed to achieve, and still fewer gave any indication of how the test was to be
performed.
Q6 Answers to this question were very disappointing. Although this was slightly more
demanding than some exposed-to-risk questions in the past, many candidates seemed to have
little notion of how to approximate the central exposed to risk.
Q7 This question was generally well answered, although part (ii) was less well answered
than similar questions on previous papers in which examples relevant to actuarial work were
asked for. Marks were deducted in part (ii) for problems which seemed trivial, or where
essentially the same examples were given for more than one class of models.
Q8 Few candidates made a serious attempt at this question. Many answers consisted of an
attempt at part (i) followed by a description of the state space in part (ii)(a), the general
expression for the Kolmogorov equations, and a statement in part (iii) that the distribution of
holding times was exponential. Few candidates attempted to write down the matrix in part
(ii). Note that credit was given in part (iv) for errors carried forward from incorrect
matrices in part (ii).
Q9 Part (i) of this question was well answered by a good proportion of candidates. Fewer
managed to calculate the values of B and c in part (ii), partly due to algebraic errors. Credit
was given for the calculation of B to candidates who calculated an incorrect value for c but
then correctly computed the value of B which corresponded to their value of c. Part (iii) was
poorly answered, with many candidates offering no comments at all.
Page 2
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
Q10 Answers to this question were very disappointing. Parts (i) and (ii) were bookwork
based on the Core Reading, yet many candidates seemed not to understand what was
required. Part (iii) was rather better answered. Candidates who derived an incorrect hazard
function in part (iii) could score full credit in part (iv) for correct sketches of these incorrect
hazards. Indeed, of the relatively small number of candidates who scored highly for the
sketches in part (iv), some did indeed produce correct plots of incorrect (and sometimes
much more complicated) hazard functions.
Q11 This question was well answered by many candidates. The only general weaknesses
were steps missing in part (i) and the lack of explanation of where the approximate variance
came from in part (iii)(b). In part (iv), an encouraging number of candidates realised that
the Cox model was an obvious alternative model, though few made any further comments on
how it might be applied to the problem mentioned in the question.
Q12 This question was also well answered by the majority of candidates. Many scored full
marks on parts (i), (ii) and (iii), and made a good attempt at part (iv). The comments asked
for in part (v) were, however, much less convincingly made. In part (iv), several candidates
combined the two categories “4 claims” and “5 claims” because the expected value was
small. Full credit was given for this if the chi-squared statistic was computed correctly, and
the number of degrees of freedom was correct for this alternative. However, candidates who
performed the test on the reduced number of categories “0 claims”, “1 claim” and “2 or
more claims” were penalised.
Page 3
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
However the fact that the company insures non-standard lives makes it unlikely that a
suitable standard table would exist.
However it is a last resort as it is difficult to obtain results which are smooth and
which adhere to the data.
EITHER
OR
If the l-step transition probabilities are dependent only on the time lag, the
chain is said to be time-homogeneous.
(ii) (a) In this case the chain is irreducible if the transition probability
out of each state is non-zero (or, equivalently, if it is possible to
reach the other state from both states)
(b) The chain is only periodic if the chain must alternate between
the states.
So a = 1 and b = 1.
Page 4
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
3 Benefits
In models of complex systems we can control the experimental conditions and thus
reduce the variance of the results without upsetting their mean values.
Often models are the only practicable means of answering actuarial questions.
Systems with a long time-frame can be studied and results obtainedrelatively quickly.
Limitations
In a stochastic model, many independent runs of the model are needed to obtain
results for a given set of inputs.
Models can look impressive and there is a danger this results in false sense of
confidence.
Users need to understand the model and the uses to which it can safely be put — the
model is not a “black box”.
It is not possible to include all future events in a model (e.g. change in legislation).
Models are better for comparing the impact of input variations than for optimising
outputs.
Page 5
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
4 (i) (a) Under UDD the number of deaths between exact ages 30 and 35 years
is half the number of deaths between exact ages 30 and 40 years.
332.5
and 5 q30 = = 0.0033716 .
98, 617
⎛ t ⎞
Then, since t p x = exp ⎜ − ∫ μ x + s ds ⎟ ,
⎜ ⎟
⎝ 0 ⎠
10 p30 = exp ( −10μ )
so
= 1 − exp[(−5)(0.0006766)] = 0.0033773 .
(ii) EITHER
OR
Page 6
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
(iii) The actual number of survivors to exact age 35 years is higher (or,
equivalently, mortality is lighter) than that under either the UDD or the
constant force assumptions.
The actual number of survivors implies that there were 258 deaths between
ages 30 and 35 years and 407 deaths between ages 35 and 40 years.
The actual data reveal that the force of mortality is higher between ages 35 and
40 years than it is between ages 30 and 35 years for females in English Life
Table 15, which suggests that the force of mortality is increasing over this age
range.
The actual force of mortality seems to be increasing even faster than is implied
by UDD.
5 Suppose we have a set of n crude mortality rates for a given age range x to x + n − 1,
and we wish to compare them to a standard set of n mortality rates for the same age
range.
If the mortality underlying the crude rates is the same as that of the standard set of
rates (the null hypothesis), then we should expect the difference between the two sets
of rates to be due only to sampling variability.
The grouping of signs test tests the null hypothesis by examining the number of
groups of consecutive positive deviations among the n ages, where a positive
deviation occurs when the crude rate exceeds the corresponding rate in the standard
set.
Then the number of possible ways to arrange t positive groups among n – m negative
⎛ n − m + 1⎞
deviations is ⎜ ⎟.
⎝t ⎠
Page 7
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
⎛ m − 1⎞
There are ⎜ ⎟ ways to arrange m positive signs into t positive groups.
⎝ t −1 ⎠
⎛n ⎞
There are ⎜ ⎟ ways to arrange m positive and n – m negative signs.
⎝m⎠
⎛ n − m + 1⎞ ⎛ m − 1⎞
⎜ ⎟⎜ ⎟
⎝ t ⎠ ⎝ t −1 ⎠
Pr[G = t ] =
⎛n ⎞
⎜ ⎟
⎝m⎠
The grouping of signs test then evaluates Pr[t ≤ G ] under the null hypothesis.
If this is less than 0.05 we reject the null hypothesis at the 5% level.
dx
rx = , x = 0, 1, 2, ... months
Exc
(b) We need to ensure that the Exc correspond to the data on persons
recovering
The hospital’s data imply a calendar month rate interval for the
recoveries, running from the first day of each month until the last day
of each month.
Page 8
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
12 11
2(
E xc = ∫ P *x,t dt = ∑ 1 P *x,t + P *x +1,t +1 ) ,
0 0
We assume all months are the same length, and that the numbers in the
hospital vary linearly across each month.
(ii) At the start of the rate interval, durations since the operation range from x − 1
to x months, so the average duration is x − ½, assuming operations take place
evenly across the month.
7 (i) Processes can be classified, first, according to whether their state space (i.e.
the range of states they can possibly occupy) is discrete or continuous
For processes operating in both discrete and continuous state space the time
domain can either be discrete or continuous
EITHER
OR
Discrete Discrete
Discrete Continuous
Continuous Discrete
Continuous Continuous
Page 9
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
(ii)
Flea transmission requires one of the x infected cats to meet one of the (10 − x)
uninfected cats.
⎛10 ⎞
(ii) The total number of pairings of cats is ⎜ ⎟ = 45.
⎝2⎠
As this depends only on the number of cats currently infected, and meetings
occur according to a Poisson process, the number of infected cats over time
follows a Markov jump process.
Page 10
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
⎛0 0 ⎞
⎜ ⎟
⎜ −9 9 ⎟
⎜ −16 16 ⎟
⎜ ⎟
⎜ −21 21 ⎟
⎜ −24 24 ⎟
μ⎜ ⎟
A= ⎜ −25 25 ⎟
90 ⎜ ⎟
−24 24
⎜ ⎟
⎜ −21 21 ⎟
⎜ −16 16 ⎟
⎜ ⎟
⎜ −9 9 ⎟
⎜ ⎟
⎝ 0⎠
Kolmogorov’s equations:
EITHER
d
forward form P(t ) = P(t ) A
dt
OR
d
backward form P(t ) = AP(t )
dt
90 μx(10 − x)
With mean OR parameter .
μx(10 − x) 90
(iv) Total expected time is the sum of the mean holding times.
90 9 90 ⎛ 1 1 1 1 1 1⎞
∑
1 1 1 1
= = ⎜ + + + + + + + + ⎟
μ x =1 x(10 − x) μ ⎝ 9 16 21 24 25 24 21 16 9 ⎠
= 50.92/µ
Page 11
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
μ x = Bc x .
Since
⎛ t ⎞
⎜ ∫ x+ w ⎟
p
t x = exp ⎜ − μ dw ⎟,
⎝ 0 ⎠
⎛ t ⎞ ⎛ t⎞
Bc x c w
we have t p x = exp ⎜ − ∫ Bc x + w dw ⎟ = exp ⎜ − ⎟,
⎜ ⎟ ⎜ ln c ⎟
⎝ 0 ⎠ ⎝ 0⎠
⎛ ⎡ Bc x ct − Bc x ⎤ ⎞ c x ( ct −1)
⎡ ⎛ − B ⎞⎤
which is exp ⎜ − ⎣ ⎦⎟=
⎢ exp ⎜ ln c ⎟ ⎥ .
⎜ ln c ⎟ ⎣ ⎝ ⎠⎦
⎝ ⎠
c50
⎡ ⎛ − B ⎞⎤
(ii) Define Q = ⎢ exp ⎜ ⎟⎥
⎣ ⎝ ln c ⎠ ⎦
ln 0.995 = (c − 1) ln Q
ln 0.989 = (c2 − 1) ln Q
(c 2 − 1) (c − 1)(c + 1)
= = 2.20665
(c − 1) (c − 1)
c = 1.20665
Therefore Q = 0.976036128
1.2066550
⎡ ⎛ −B ⎞⎤
⎢exp ⎜ ln1.20665 ⎟ ⎥ = 0.976036128
⎣ ⎝ ⎠⎦
B = 3.797*10−7.
(iii) In this example, only two observations are provided so there is an analytical
solution to the Gompertz model.
Page 12
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
The more general graduation process allows the fitting of more complex
models from the Gompertz-Makeham family which have the form
μ x = polynomial(1) + exp(polynomial(2))
(b) EITHER
Pr[T > x + t ]
Since Pr[Tx > t ] = Pr[T > x + t | T > x] =
Pr[T > x]
S (x + t)
then S x (t ) = .
S ( x)
OR
x +tp0 S ( x + t )
Therefore t px = S x (t ) = = .
x p0 S ( x)
(ii) EITHER
1 d
μ x +t = − [Pr(Tx > t )]
Pr[Tx > t ] dt
OR
1
μ x +t = lim
+
( Pr[Tx ≤ t + h | Tx > t )
h →0 h
(iii) EITHER
d
f x (t ) = S x (t )μ x +t = − S x (t )
dt
Page 13
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
1 d
Therefore μ x +t = − S x (t )
S x (t ) dt
( )
If S x (t ) = exp −(λt )β , therefore, we have
μ x +t = −
1 d
(
exp −(λt )β )
(
exp −(λt )β ) dt
μ x +t = −
1
exp ( −(λt ) )
β ( exp ( −(λt ) )) ( −λ βt ) = λ βt
β β β−1 β β−1
OR
⎡ t ⎤
S x (t ) = exp ⎢ − ∫ μ x + s ds ⎥ = exp ⎡ −(λt )β ⎤ .
⎢⎣ 0 ⎥⎦ ⎣ ⎦
So
d ⎡ ⎤
t
d
⎢ ∫ μ x + s ds ⎥ = μ x +t = ⎡ (λt )β ⎤ ,
dt ⎢ ⎥⎦ dt ⎣ ⎦
⎣0
and hence
μ x +t = βλβt β−1 .
(iv)
Page 14
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
We have
x = t p x dt p x +t + t p x dt p x +t .
p12 11 12 12 22
t + dt
Also, dt px22+t = 1,
x +t = μ x +t dt + o( dt )
p12 12
dt
where o(dt) is the probability that a life makes two or more transitions in the
time interval dt, and
o(dt )
lim =0.
dt →0 dt
x = μ x +t t p x dt + t p x + o( dt )
p12 12 11 12
t + dt
Thus
x − t p x = μ x +t t p x dt + o(dt )
p12 12 12 11
t + dt
and
∂ 12 x − t px
p12 12
t + dt
t p x = lim + = μ12 11
x +t t p x
∂t dt →0 dt
(ii) (a) Suppose we observe d12 deaths from heart disease, d13 deaths from
cancer and d14 deaths from other causes.
Suppose also that we observe the waiting time for each life, and that
the total observed waiting time is V, being the sum of the waiting times
for each life.
Page 15
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
( )
L ∝ exp ⎡ − μ12 + μ13 + μ14 V ⎤ (μ12 ) d (μ13 ) d (μ14 ) d .
12 13 14
⎣ ⎦
log L = −(μ12 + μ13 + μ14 )V + d 12 log μ12 + d 13 log μ13 + d 14 log μ14 + K
(where K is a constant )
∂ log L d 12
= −V + ,
∂μ12 μ12
and setting the partial derivative equal to zero leads to the solution
d 12
μˆ 12 = .
V
∂ 2 log L d 12
Since =− , the second derivative is always negative
(∂μ12 ) 2 (μ12 )2
and so we have a maximum.
(iii) (a) The maximum likelihood estimate of the force of mortality from heart
disease is 34/1,065 = 0.0319249
0.0319249
= 0.000029976 .
1, 065
Page 16
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
0.000029976 = 0.00547507 .
(iv) Using the four state model, the lives in the investigation would have to be
stratified according to the risk factors and the transition intensities estimated
separately for each stratum.
Using a Cox regression model with death from heart disease as the event of
interest and the risk factors as covariates would avoid this problem.
Lives who died from other causes could be treated as censored at the durations
when they died.
⎛ 0.259182 0.740818 0 0 ⎞
⎜ ⎟
⎜ 0.259182 0 0.740818 0 ⎟
⎜ 0.036936 0.222245 0 0.740818 ⎟
⎜ ⎟
⎝ 0 0.036936 0.222245 0.740818 ⎠
(ii) π = πP
Page 17
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
Using (4)
In (3)
Then in (2)
So
π4 = 1/ (1+0.349859+0.17226+0.07771)=0.625067
π3 = 0.218685
π2 = 0.107674
π1 = 0.048574
60%*0.625067+50%*0.218685+25%*0.107674 = 51.13%
Null hypothesis: the data come from a source where the underlying
distribution of number of claims follows a Poisson distribution with mean
0.30.
This is a one-tailed test, and the upper 5% point of the chi-squared distribution
with 5 degrees of freedom is 11.07.
Page 18
Subject CT4 — Models Core Technical — April 2009 — Examiners’ Report
(v) As the goodness of test fails, the discount level calculated assuming the
Poisson distribution may be incorrect.
The average discount level calculated from the data could usefully be
compared with that estimated using the Poisson distribution.
Page 19
Faculty of Actuaries Institute of Actuaries
EXAMINATION
1. Enter all the candidate and examination details as requested on the front of your answer
booklet.
2. You must not start writing your answers in the booklet until instructed to do so by the
supervisor.
4. Attempt all 11 questions, beginning your answer to each question on a separate sheet.
Hand in BOTH your answer booklet, with any additional sheets firmly attached, and this
question paper.
In addition to this paper you should have available the 2002 edition of the Formulae
and Tables and your own electronic calculator from the approved list.
© Faculty of Actuaries
CT4 S2009 © Institute of Actuaries
1 Describe the difference between the following assumptions about mortality between
any two ages, x and y (y > x):
In your answer, explain the shape of the survival function between ages x and y under
each of the two assumptions. [2]
2 (i) List the key steps in constructing a new actuarial model. [4]
You work for an actuarial consultancy which is taking over responsibility for a
modelling process which has previously been conducted in house by a client.
(ii) Discuss the extent to which the steps required for this task differ from those
listed in your answer to (i). [2]
[Total 6]
3 (i) List the data needed for the exact calculation of a central exposed to risk
depending on age. [2]
An investigation studied the mortality of persons aged between exact ages 40 and 41
years. The investigation began on 1 January 2008 and ended on 31 December 2008.
The following table gives details of 10 lives involved in the investigation.
Persons with no date of death given were still alive when the investigation ended.
(ii) Calculate a central exposed to risk using the data for the 10 lives in the
sample. [3]
(iii) (a) Calculate the maximum likelihood estimate of the hazard of death at
age 40 last birthday.
CT4 S2009—2
4 (i) In the context of mortality investigations describe the principle of
correspondence and give an example of a situation in which it may be hard to
adhere to this principle. [2]
The government of the country wishes to obtain estimates of the force of mortality,
μ x , by single years of age x for the period between 1 January 2005 and 1 January
2008. Annual population censuses have been taken on 30 June each year since 2004,
which classify the population by age last birthday. However the only copy of the data
from the population census of 30 June 2006 was lost when the computer disc on
which it was stored was being transferred between government departments.
Let the population aged x last birthday on 30 June in year t be denoted by the symbol
Px,t , and the number of deaths during the period of investigation of persons aged x be
denoted by the symbol dx.
(ii) Derive an expression in terms of Px,t and dx which may be used to estimate
μx . [6]
[Total 8]
(ii) Prove that if the process has independent increments it satisfies the Markov
property. [3]
(iii) (a) Describe the difference between a Markov chain and a Markov jump
process.
An actuarial student can see the office lift (elevator) from his desk. The lift has an
indicator which displays on which of the office’s five floors it is at any point in time.
For light relief the student decides to construct a model to predict the movements of
the lift.
(a) irreducible
(b) has the Markov property
[3]
[Total 9]
The company models the arrival of complaints using a Poisson process with rate 1.25
per working day.
(i) List the assumptions underlying the Poisson process model. [2]
(iii) Define a state space under which the number of outstanding complaints can be
modelled as a Markov jump process. [2]
(iv) Discuss the appropriateness of using the model for this purpose, with reference
to the assumptions being made. [3]
[Total 9]
7 A firm rents cars and operates from three locations — the Airport, the Beach and the
City. Customers may return vehicles to any of the three locations.
The company estimates that the probability of a car being returned to each location is
as follows:
Car returned to
Car hired from Airport Beach City
It is suggested that the cars should be based at each location in proportion to the
stationary distribution.
CT4 S2009—4
(iv) Sketch, using your answers to parts (i) and (ii), a graph showing the
probability that a car currently located at the Airport is subsequently at the
Airport, Beach or City against the number of times the car has been rented. [3]
[Total 10]
8 A researcher is studying a certain incurable disease. The disease can be fatal, but
often sufferers survive with the condition for a number of years. The researcher
wishes to project the number of deaths caused by the disease by using a multiple state
model with state space:
{H – Healthy, I – Infected, D(from disease) – Dead (caused by the disease), D(not from disease)
– Dead (not caused by the disease)}.
(i) Draw a transition diagram for the multiple state model. [2]
(iii) Determine integral expressions, in terms of the transition rates and any
expressions previously determined, for:
(a) PHH(x, x + t)
(b) PHI(x, x + t)
When the test was complete, the sub-contractor reported that he had terminated the
test after 150 days. He further reported that:
However, he reported that he was only able to return 11 batteries, as one had exploded
after 110 days, and he had treated this battery as censored at that duration when
working out the Kaplan-Meier estimate of the survival function.
(i) State, with reasons, the forms of censoring present in this study. [2]
(ii) Calculate the Kaplan-Meier estimate of the survival function based on the
information supplied by the sub-contractor. [5]
In his report, the sub-contractor claimed that the Kaplan-Meier estimate of the
survival function at the duration when the investigation was terminated was 0.2727.
CT4 S2009—6
10 An investigation into the mortality of men engaged in a hazardous occupation was
carried out. The following is an extract from the results.
30 950 12 0.0126
31 1,200 14 0.0117
32 1,200 16 0.0133
33 900 9 0.0100
34 1,000 11 0.0110
35 1,100 15 0.0136
36 800 10 0.0125
37 1,250 16 0.0128
38 1,400 17 0.0121
It was decided to graduate the results with reference to English Life Table 15 (males).
o
The formula used for the graduation was q x = 10qxs .
(i) Using a test of the overall fit of the graduated rates to the data, test the
hypothesis that the underlying mortality of men in the hazardous occupation is
in accordance with the graduation formula given above. [6]
(ii) Test the graduation using two other tests which detect different features of the
graduation. For each test you apply:
The study investigated the impact of age, sex and educational qualifications on the
hazard of returning to work using the following covariates:
The investigators decided to use a Cox proportional hazards regression model for the
study.
(ii) Explain why the Cox model is a popular model for the analysis of survival
data. [3]
(iii) (a) Write down the equation of the model that was estimated, defining
the terms you use (other than those defined above).
(b) List the characteristics of the young person to whom the baseline
hazard applies. [3]
• The hazard of resuming work for males who started claiming benefit aged 17
years exact and who had passed the mathematics examination was 1.5 times the
hazard for males who started claiming benefit aged 16 years exact but who had
not passed the mathematics examination.
• Females who had passed the mathematics examination were twice as likely to take
up a new job as were males of the same age who had failed the mathematics
examination.
• Females who started claiming benefit aged 20 years exact and who had passed the
mathematics examination were twice as likely to resume work as were males who
started claiming benefit aged 16 years exact and who had also passed the
mathematics examination.
END OF PAPER
CT4 S2009—8
Faculty of Actuaries Institute of Actuaries
EXAMINERS REPORT
Introduction
The attached subject report has been written by the Principal Examiner with the aim of
helping candidates. The questions and comments are based around Core Reading as the
interpretation of the syllabus to which the examiners are working. They have however given
credit for any alternative approach or interpretation which they consider to be reasonable.
R D Muckart
Chairman of the Board of Examiners
December 2009
Comments for individual questions are given with the solutions that follow.
Faculty of Actuaries
Institute of Actuaries
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
Examiners’ Comments
Comments on solutions presented to individual questions for the September 2009 paper are
given below. In general, those using this report should be aware that in the case of non-
numerical answers full credit could often be obtained for rather less than is given in the
solutions which follow. The solutions are meant as a guide to the various points which could
have been made and considered relevant.
1
A uniform distribution of deaths means
EITHER
that deaths are evenly spaced between the ages x and y.
OR
that t qx tqx (t y x)
OR
that t px x t is constant for t y x .
It also means that the survival function decreases linearly between ages x and y. The
assumption of a constant force of mortality between any two ages means
EITHER
that the hazard does not change with age over this age range.
OR
that t px ( px )t .
This implies that the survival function decreases exponentially between ages x and y.
2
(i) Define objectives of modelling process.
Plan the modelling process and how it will be validated.
Collect and validate the data required.
Page 2
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
(ii) Whilst in theory all steps are still required, some may take the form of
reviewing the appropriateness of existing decisions made, such as how the
form of the model was determined.
Extent of work will depend on whether the existing model is to be used,
adapted or superseded.
An understanding of how results compare with those previously used by the
company will be required.
Process maps for the existing approach, or discussions with the people running
the process about what they do, may be helpful.
The scope needs to be tightly defined up front to ensure it is clear what is
expected of the consultancy.
Data sources may already be established.
3
(i) For each life we need
EITHER date of birth OR exact age at entry into observation OR exact age at
exit from observation
Date of entry into observation
Date of exit from observation
Page 3
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
[Alternatives were given full credit, provided the information given allowed the
calculation of the date of entry into and exit from observation and the life’s age]
(ii) The contribution of each life to the central exposed to risk is the number of
months between STARTDATE and ENDDATE, where STARTDATE is the
latest of date of 40th birthday 1 January 2008 and ENDDATE is the earliest of
date of 41st birthday date of death 31 December 2008
Summing the number of months over the 10 lives gives a total of 53 months,
which is 4.42 years, which is the central exposed to risk.
(iii)
a. The total number of deaths during the period of observation is 2. So the
maximum likelihood estimate of the hazard of death is 2/4.42 =
0.4528.
b. ALTERNATIVE 1
If the hazard of death at age 40 years is 40 , then
q40 1 p40 1 exp(40 )
= 1 exp(0.4528) 1 0.6358 0.3642.
ALTERNATIVE 2
If the central exposed to risk is E40c , then if we work in years
d 40
q40
E 0.5d 40
c
40
Page 4
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
2 2
= 0.3690.
4.42 1 5.42
This was well answered. A common error was to count 3 deaths rather than 2.
Although 3 deaths are mentioned in the data given in the question, one of these
occurred after the life’s 41st birthday and so should not be included in the estimation
of μ40. Another common error was to forget that exposure ends at exact age 41 years.
Each of these errors was only penalised once, so that calculations which followed
through correctly in (iii) were awarded full marks for part (iii). Note also that
candidates who made BOTH the above errors were only penalised for one, as if
exposure is assumed to continue past exact age 41 years, it is consistent to count 3
deaths!
4
(i) The principle of correspondence states that a life alive at time t should be
included in the exposure at age x at time t if and only if, were that life to die
immediately, he or she would be counted in the deaths data at age x. Problems
in adhering to this can arise when the deaths data and the exposed-to-risk data
come from two different sources. These may classify lives differently.
(ii) Since deaths are classified by age last birthday at date of death, a central
exposed to risk which corresponds to the deaths data is given by
t 3
Exc Px,t
t 0
where Px,t is the population aged x last birthday at time t, and t is measured in
years since 1 January 2005. We have censuses on 30 June 2004, 30 June 2005,
30 June 2007 and 30 June 2008.
Assuming that the population varies linearly across the period between each
successive census for which we have data the population aged x last birthday
on 1 January 2005 is equal to
1 (P Px,30 / 6 / 2005 )
2 x,30 / 6 / 2004
and the population aged x last birthday on 1 January 2008 is equal to
1 (P Px,30 / 6 / 2008 ) .
2 x,30 / 6 / 2007
Dividing the period of the investigation into three sub-periods
from 1 January 2005 to 30 June 2005
from 30 June 2005 to 30 June 2007
from 30 June 2007 to 1 January 2008
and applying the trapezium rule to each sub-period produces the following
exposed to risk for persons aged x last birthday
For the sub-period between 1 January 2005 and 30 June 2005
Page 5
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
1 1 (P Px,30/6/2005 )
2 2 x,1/1/2005
1 1 ( 1 ( Px,30/6/2004 Px,30/6/2005 ) Px,30/6/2005 )
2 2 2
For the sub-period between 30 June 2005 and 30 June 2007
2 1 ( Px,30/6/2005 Px,30/6/2007 )
2
For the sub-period between 30 June 2007 and 1 January 2008
1 1 (P Px,1/1/2008 )
2 2 x,30/6/2007
1 1 ( Px,30/6/2007 1 ( Px,30/6/2007 Px,30/6/2008 ))
2 2 2
Summing these gives
Exc 1 Px,30/6/2004 1 Px,30/6/2005 1 Px,30/6/2005 Px,30/6/2005
8 8 4
Px,30/6/2007 1 Px,30/6/2007 1 Px,30/6/2007 1 Px,30/6/2008
4 8 8
which simplifies to
Exc 1 Px,30/6/2004 11 Px,30/6/2005 11 Px,30/6/2007 1 Px,30/6/2008 .
8 8 8 8
The force of mortality may be estimated using the formula
dx
x ,
E xc
where d x denotes deaths to persons aged x last birthday when they died.
This was very poorly answered. It was perhaps rather more difficult than some
exposed-to-risk questions in previous examination papers, but nevertheless the
standard of most attempts was disappointing. In part (ii) credit was given for various
alternative approximations provided that they were explained clearly.
5
(i) The Markov property states that the future development of a process can be
predicted from its present state alone without reference to its past history.
(ii) Formally, for times s1 s2 ... sn s t and for states x1, x2 ,..., xn , x in the
state space S and all subsets A of S, the Markov property can be written
Pr[ X (t ) A | X (s1) x1, X (s2 ) x2 ,...., X ( sn ) xn , X ( s) x] Pr[ X t A | X ( s) x]
For independent increments we can write
Pr[ X (t ) A | X ( s1 ) x1, X ( s2 ) x2 ,...., X ( sn ) xn , X ( s) x]
Pr[ X (t ) X ( s ) x A | X ( s1 ) x1, X ( s2 ) x2 ,...., X ( sn ) xn , X ( s ) x]
Pr[ X (t ) X ( s ) x A | X ( s ) x]
Pr[ X (t ) A | X ( s ) x]
Page 6
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
(iii)
a. A Markov chain is a stochastic process with the Markov property
which has a discrete time set with a discrete state space. A Markov
jump process is a stochastic process with the Markov property which
has a continuous time set with a discrete state space.
b.A Markov chain is irreducible if any state can be reached from any
other state.
(iv)
a. A lift could not serve its purpose unless it could return to each of the
floors which it serves. This means an irreducible model would be
appropriate.
b.Suppose, for example, the lift is currently at the third floor, with its last
two states being the fourth floor and the fifth floor. In such a case the
lift is more likely to be heading downwards than upwards. So the past
history is likely to provide information on the likely future movement
of the lift, unless the state space is very complicated (involving a
number of past floors as well as the current floor). Therefore a Markov
model is unlikely to be appropriate.
This question was generally well answered, apart from section (iv)(b) in which few
candidates spotted the point that the direction of travel of the lift as well as its current
floor will influence its next location.
6
(i) A Poisson process is a continuous-time integer valued process
Nt, t 0 with
N0 = 0
independent increments
EITHER
increments follow a Poisson distribution
OR
[(t s)]n exp[(t s)]
P[ Nt N s n] , for s < t, n = 0, 1, 2, ....
n!
(ii) Average work created by a complaint is
Page 7
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
but do need to know how many of each type. If cases are allocated randomly
rather than in order, then the state space consists of (in terms of complaints not
resolved):
r – straightforward,
s – medium,
t – complicated.
where r = 0,1,2,3,4,5,….
s = 0,1,2,3,4,5,……
t = 0,1,2,3,4,5,…..
(iv) EITHER The model will only give an approximation.
OR The model is not suitable for this purpose.
The model could not be used to do this without extending the state space to
consider the time the complaint has been in the queue. There are only two
employees, so holidays and sickness are important factors not taken into
account.
The model assumes complaints are time-homogeneous. We do not know the
nature of the business, but for some industries complaints would be seasonal
e.g. holiday companies.
The model assumes that complaint arrivals are independent, but more
complaints might be expected if the company has had a quality control
problem at a particular time. If struggling to meet the service standard, action
would be. Taken, such as overtime, or prioritising easy cases. Staff may be
able to deal with complaints which are similar to other recent complaints very
quickly, using standard „template‟ responses.
The memoryless property is unlikely to be realistic as the work required to
complete the case could be assessed and then worked through to a schedule.
The Markov jump process could be used to estimate the probability that a
complaint is responded to within a given number of days of receipt.
So the model could be used to estimate the probability of a complaint not
being responded to in the stated time, that is the failure to meet the service
standard.
7
(i) Two step transition matrix
Page 8
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
3 0.251 0.53
and 1 2 3 1
1 23
2 33
1 1
3
2 1
2
3 1
6
(iii)The stationary distribution gives the long run probability that a particular car
will be at each location. However this does not take into account the demand
for hiring vehicles at each location, or the amount of space available at each
location. These factors are likely to be more important in determining how
many cars to base at each site.
1.2
0.8
Probabilty at location
Airport
0.6 Beach
City
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9
Page 9 Number of rentals
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
The calculations in parts (i) and (ii) were, as is usually the case in CT4 examinations,
successfully completed by the vast majority of candidates. However only a minority
made the point that, whereas the stationary distribution gives the long run probability
that cars will be returned to each location, the company would be better advised to
position cars at the three locations to reflect the demand for rentals. In part (iv), some
candidates drew a set of histograms. Credit was given for this, provided that
histograms were presented for 1 rental, 2 rentals, and the long run distribution,
together with a statement that at 0 rentals the car must be at the Airport.
8
(i)
( x) Dead
( x)
(from
Healthy Infected disease)
( x) ( x) ( x)
Dead (not
from
disease)
d
(ii) P( x) P( x) A( x) where with order of state space
dt
{Healthy, Infected, Dead (not disease), Dead(from disease)}
( x) ( x) ( x) ( x) 0
0 ( x) ( x) ( x) ( x)
A(x)=
0 0 0 0
0 0 0 0
Page 10
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
(iii)
t
a. PHH(x, x+t)= exp[ (( x w) ( x w))dw]
w0
t t
b. PHI(x, x+t)= PHH ( x, x w).( x w).exp[ ( x u) du].dw
w0 u w
c. EITHER
t
PHD(from disease)(x, x+t)= PHI ( x, x w).(x w).dw
w0
OR (backwards alternative)
PHD(from disease)(x, x+t)
t
=
w0
P HH ( x w). ( x w).PID ( fromdisease ) ( x w, x t ).dw .
s
and PII ( x w, x s ) exp ( x u )du .
u w
So PHD(from disease)(x, x+t)
t
s t
PHH ( x w). ( x w). exp ( x u )du . ( x s).ds.dw
w 0 sw u w
This question was considerably better answered than were similar questions in
previous examinations. In particular, the proportion of candidates making serious
attempts at part (iii) was greater than has been the case for similar questions in the
past.
9
(i) Type II censoring as the study was terminated after a pre-determined number
of failures. Random censoring of the device which exploded.
(ii) According to the information supplied by the sub-contractor, the Kaplan-
Meier estimate of the survival function should be calculated as follows:
j tj Nj dj cj dj/Nj 1 – dj/Nj
0 0 12
1 97 12 2 1 2/12 10/12
Page 11
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
0 0 11
1 97 11 2 0 2/11 9/11
2 120 9 3 0 3/9 6/9
3 141 6 2 0 2/6 4/6
4 150 4 1 0 1/4 3/4
+½ +½
The Kaplan-Meier estimate is then
dj
Sˆ (t ) 1
Nj
t j t
so we have
t Sˆ (t )
Page 12
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
0 t 97 1
97 t 120 9/11
120 t 141 6/11
141 t 150 4/11
150 t 3/11 = 0.2727
Many candidates scored highly on this question. Credit was given in part (i) for other
types of censoring provided that a sensible reason was given. In part (iii), for full
credit some kind of calculation of an alternative survival function was needed,
together with an explanation of why this provided evidence to support the suggestion
that the sub-contractor has stolen the battery.
10
(i) The chi-squared test is for the overall fit of the graduated rates to the data
The test statistic is z x2 , where
o
x Ex qx
zx .
o o
E x q x (1 q x )
o
The calculations are shown in the table below (since qx is
o
x Ex qx
small we use the approximation z x .
o
Ex qx
o o
Age x x qx Ex qx zx z x2
∑ 4.808
The test statistic has a chi-squared distribution with degrees of freedom (d.f.)
given by number of ages
o
– 1 (for parameter of function linking qx and q xs )
– some d.f. for constraints imposed by choice of standard table
The critical value of the chi-squared distribution is
11.07 with 5 d.f.
12.59 with 6 d.f.
14.07 with 7 d.f.
15.51 with 8 d.f.
16.92 with 9 d.f. at the 5% level (from tables)
Since 4.808 < 11.07 (or 12.59 etc.) there is no evidence to reject the null
hypothesis that the graduated rates are the true rates underlying the crude
rates.
(ii) EITHER
Signs test
a. The Signs test looks for overall bias.
b. The number of positive signs among the z x s
is distributed Binomial (9, 0.5).
We observe 6 positive signs.
The probability of obtaining 6 or more positive signs is
(from tables)
1 – 0.7461 = 0.2539.
a. When applied over the whole age range, the Cumulative Deviations
test looks for overall bias
b. The test statistic is
o
x Ex qx
x
Normal(0,1)
o
Ex qx
x
o o
Age x x Ex qx x Ex qx
30 12 8.645 3.355
31 14 11.28 2.72
32 16 11.64 4.36
33 9 8.91 0.09
34 11 10.60 0.40
35 15 12.76 2.24
36 10 10.16 -0.16
37 16 17.25 -1.25
38 17 20.86 -3.86
∑ 112.105 7.895
7.895
So the value of the test statistic is 0.7457
112.105
Using a 5% level of significance, we see that
1.96 < 0.7457 < 1.96
c. We accept the null hypothesis at the 5% level of significance and
conclude there is no overall bias in the graduation.
Grouping of Signs test
a. The Grouping of Signs test looks for runs or clumps of deviations of
the same sign OR the grouping of signs test tests for overgraduation.
b. We have:
9 ages in total
6 positive deviations
3 negative deviations
We have 1 positive run
Pr[1 positive run] is therefore equal to
Page 15
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
5 4
0 1 4
4
0.0476
9 9.8.7 84
6 3.2
Since this is less than 0.05 (using a one-tailed
test)
c. We reject the null hypothesis that the graduated rates are the true rates
underlying the crude rates (OR we conclude that the graduation is
unsatisfactory OR there is evidence of over-graduation).
Individual Standardised Deviations test
a. The Individual Standardised Deviations tests looks for individual large
deviations at particular ages.
b. If the graduated rates were the true rates underlying the observed rates
we would expect the individual deviations to be distributed Normal
(0,1) and therefore only 1 in 20 z x s should have absolute magnitudes
greater than 1.96. Looking at the z x s we see that the largest individual
deviation is 1.278. Since this is less in absolute magnitude than 1.96
c. we cannot reject the null hypothesis that the graduated rates are the
true rates underlying the crude rates.
11
(i) A proportional hazards (PH) model is a model which allows investigators to
assess the impact of risk factors, or covariates, on the hazard of experiencing
an event.
In a PH model the hazard is assumed to be the product of two terms, one
which depends only on duration, and the other which depends only on the
values of the covariates.
Under a PH model, the hazards of different lives with covariate vectors z1 and
z2 are in the same proportion at all times:
for example in the Cox model
(t ; z1 ) exp( z1T )
.
(t ; z2 ) exp( z2T )
Page 16
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
(ii) Cox‟s model ensures that the hazard is always positive. Standard software
packages often include Cox‟s model.
Cox‟s model allows the general “shape” of the hazard function for all
individuals to be determined by the data, giving a high degree of flexibility
while an exponential term accounts for differences between individuals.
This means that if we are not primarily concerned with the precise form of the
hazard, we can ignore the shape of the baseline hazard and estimate the effects
of the covariates from the data directly.
(iii)
a. (t ) 0 (t ) exp( A A E E S S ) , where (t ) is the estimated
hazard and 0 (t ) is the baseline hazard.
b. A female aged exactly 16 years when she first claimed benefit who had
not passed the school mathematics examination.
(iv) “The hazard of resuming work for males aged 17 years who had passed the
mathematics examination was 1.5 times the hazard for males aged 16 years
who had not passed the mathematics examination” implies that
exp[( A *1) S E ]
exp( A E )
exp(S )
exp( A ) exp( E ) 1.5
“Females who had passed the examination were twice as likely to take up a
new job as were males of the same age who had failed” implies that
exp( E )
2
exp(S )
since the age terms cancel out.
“Females aged 20 years who had passed the examination were twice as likely
to resume work as were males aged 16 years who had also passed the
examination” implies that
exp( A * 4)
2.
exp(S )
Page 17
Subject CT4 — Models Core Technical — September 2009 — Examiners’ Report
exp(5 A ) 1.5
log e 1.5
A 0.0811
5
From (1) then, we obtain
exp( E ) exp(0.0811) 1.5
E 0.0811 0.4055
E 0.3244 .
Finally, from (2) we obtain
exp(0.3244)
2
exp(S )
0.3244 S log e 2 0.6931
S 0.3688
This was satisfactorily answered by many candidates. Although it is still the case
than only a minority of candidates seem to understand the essential feature of a
proportional hazards model that the hazard can be factorised into one part depending
on duration and another part depending on the values of covariates, many candidates
could list some advantages of the Cox model in part (ii). In part (iii)(b) very few
candidates spotted that the baseline person was aged 16 years when first claiming
benefit. In part (iv) candidates who failed to write down the correct equations
implied by the three statements in the question were given some credit for correctly
solving the equations they did produce.
Page 18