01.06.2015 Views

Actuarial Modelling of Claim Counts Risk Classification, Credibility ...

Actuarial Modelling of Claim Counts Risk Classification, Credibility ...

Actuarial Modelling of Claim Counts Risk Classification, Credibility ...

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.

42 <strong>Actuarial</strong> <strong>Modelling</strong> <strong>of</strong> <strong>Claim</strong> <strong>Counts</strong><br />

the time, the test is performed with j = 1, so that we compare the full model to a simpler<br />

one with one parameter less.<br />

Let ˜ be the maximum likelihood estimator under H 0 , and let ̂ be the maximum likelihood<br />

estimator under H 1 . The likelihood ratio test is based on the ratio <strong>of</strong> the likelihoods between<br />

a full and a restricted (or reduced) nested model with fewer parameters. The restricted model<br />

must be nested within (i.e., be a subset <strong>of</strong>) the full model. The likelihood ratio test statistic is<br />

T = 2ln ̂ (<br />

)<br />

˜ = 2 L̂ − L˜ <br />

The evidence against H 0 will be strong when T is large.<br />

The Chi-square distribution plays a prominent role in likelihood ratio tests. Recall that the<br />

Gamma distribution with = /2 and = 1/2 for some positive integer is known as the<br />

Chi-square distribution with degrees <strong>of</strong> freedom (which is denoted as 2 ), with associated<br />

probability density function<br />

fx = x/2−1 exp−x/2<br />

( ) <br />

<br />

2 2<br />

/2<br />

x>0<br />

If X ∼ 2 then its mean is , and its variance 2. It is useful to recall that the 2 distribution<br />

is closely related to the Normal distribution. Specifically, the 2 arises as the distribution <strong>of</strong><br />

the sum <strong>of</strong> independent squared or0 1 random variables.<br />

Under H 0 , the test statistic T is approximatively Chi-square distributed with degrees<br />

<strong>of</strong> freedom equal to the number <strong>of</strong> parameters in the full model minus the number <strong>of</strong><br />

parameters in the restricted model (that is, with j degrees <strong>of</strong> freedom) when the sample size<br />

n is sufficiently large (and additional mild regularity conditions are fulfilled). Note that the<br />

likelihood ratio test requires us to perform two maximum likelihood estimations, one under<br />

H 0 and another one under H 1 . When the largest model H 1 is misspecified (that is, the data<br />

have not been generated by this probability model), the likelihood ratio statistic is no longer<br />

necessarily Chi-square distributed under H 0 .<br />

Unfortunately, there are cases where regularity conditions do not hold for T to be<br />

approximately j<br />

2 distributed under H 0 . In particular this happens when a constrained<br />

parameter is on the boundary <strong>of</strong> the parameter space, e.g., testing Poisson versus Negative<br />

Binomial. Here Poisson is a particular case <strong>of</strong> Negative Binomial when the latter has a<br />

parameter on its boundary space. In this case, the limiting distribution <strong>of</strong> the statistic T<br />

becomes a mixture <strong>of</strong> Chi-square distributions. We refer the reader to Titterington ET AL.<br />

(1985) for more details about these situations.<br />

Wald Tests<br />

The Wald test provides an alternative to the likelihood ratio test that requires the estimation<br />

<strong>of</strong> only the full model, not the restricted model. The logic <strong>of</strong> the Wald test is that if the<br />

restrictions are correct then the unrestricted parameter estimates should be close to the value<br />

hypothesized under the restricted model.<br />

The Wald test is based on the distribution <strong>of</strong> a quadratic form <strong>of</strong> the weighted sum <strong>of</strong><br />

squared Normal deviates, a form that is known to be Chi-square distributed. Specifically,<br />

using (1.49), we can test H 0 = 0 versus H 1 ≠ 0 with the statistic

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

Saved successfully!

Ooh no, something went wrong!