01.08.2014 Views

Estimation, Evaluation, and Selection of Actuarial Models

Estimation, Evaluation, and Selection of Actuarial Models

Estimation, Evaluation, and Selection of Actuarial Models

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

127<br />

The total <strong>of</strong> the last column is the test statistic <strong>of</strong> 2.85. With three degrees <strong>of</strong> freedom (four<br />

groups less one less zero estimated parameters) the critical value is 7.81 <strong>and</strong> the null hypothesis <strong>of</strong><br />

a Poisson distribution cannot be rejected.<br />

Exercise 86 These are discrete data from a discrete population, so the normal <strong>and</strong> gamma models<br />

are not appropriate. There are two ways to distinguish among the three discrete options. One<br />

is to look at successive values <strong>of</strong> kn k /n k−1 . They are 2.67, 2.33, 2.01, 1.67, 1.32, <strong>and</strong> 1.04. The<br />

sequence is linear <strong>and</strong> decreasing, indicating a binomial distribution is appropriate. An alternative<br />

is to compute the sample mean <strong>and</strong> variance. They are 2 <strong>and</strong> 1.494 respectively. The variance is<br />

considerably less than the mean, indicating a binomial model.<br />

Exercise 87 The various tests for the three data sets produce the following results. For Data<br />

Set B truncated at 50, the estimates are ˆα =0.40982, ˆτ =1.24069, <strong>and</strong>ˆθ =1, 642.31. For Data<br />

Set B censored at 1,000 there is no maximum likelihood estimate. For Data Set C, the maximum<br />

likelihood estimate is ˆα =4.50624, ˆτ =0.28154, <strong>and</strong>ˆθ =71.6242.<br />

B truncated at 50 B censored at 1,000<br />

Criterion Exponential Weibull Trans gam Exponential Weibull Trans gam<br />

K-S 0.1340 0.0887 0.0775 0.0991 0.0991 N/A<br />

A-D 0.4292 0.1631 0.1649 0.1713 0.1712 N/A<br />

χ 2 1.4034 0.3615 0.5169 0.5951 0.5947 N/A<br />

p-value 0.8436 0.9481 0.7723 0.8976 0.7428 N/A<br />

loglikelihood −146.063 −145.683 −145.661 −113.647 −113.647 N/A<br />

SBC −147.535 −148.628 −150.078 −115.145 −116.643 N/A<br />

C<br />

Criterion Exponential Weibull Trans gam<br />

K-S N/A N/A N/A<br />

A-D N/A N/A N/A<br />

χ 2 61.913 0.3698 0.3148<br />

p-value 10 −12 0.9464 0.8544<br />

loglikelihood −214.924 −202.077 −202.046<br />

SBC −217.350 −206.929 −209.324<br />

For Data Set B truncated at 50, there is no reason to use a three-parameter distribution. For<br />

Data Set C, the transformed gamma distribution does not provide sufficient improvement to drop<br />

the Weibull as the model <strong>of</strong> choice.<br />

Exercise 88 The loglikelihood values for the two models are −385.9 for the Poisson <strong>and</strong> −382.4 for<br />

the negative binomial. The test statistic is 2(−382.4 + 385.9) = 7.0. There is one degree <strong>of</strong> freedom<br />

(two parameters minus one parameter) <strong>and</strong> so the critical value is 3.84. The null hypothesis is<br />

rejected <strong>and</strong> so the data favors the negative binomial distribution.<br />

Exercise 89 The penalty function subtracts ln(100)/2 =2.3 for each additional parameter. For<br />

the five models, the penalized loglikelihoods are – generalized Pareto: −219.1 − 6.9 =−226.0,<br />

Burr: −219.2 − 6.9 =−226.1, Pareto: −221.2 − 4.6 =−225.8, lognormal: −221.4 − 4.6 =−226.0,<br />

<strong>and</strong> inverse exponential: −224.3 − 2.3 =−226.6. The largest value is for the Pareto distribution.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!