Actuarium
SOA
ASA
SRM
PA

Fitting a Poisson GLM to Territory / Vehicle-Age Claim Frequency

Fit a log-link Poisson GLM to a small personal-auto rating grid, interpret coefficients, and use deviance to test significance.

GLM structure
Log link & offsets
Deviance & goodness of fit
Relativities

Background

Meridian Mutual's pricing team is refreshing personal-auto claim-frequency relativities for two rating variables: Territory (Urban vs. Rural) and Vehicle Age (New vs. Old, where "Old" means model year 6+). A pricing analyst pulled one accident year of earned exposure and reported claim counts for the four cells formed by crossing these two variables:

TerritoryVehicle AgeEarned exposureClaim countRaw frequency
UrbanNew1,2001080.0900
UrbanOld9001170.1300
RuralNew800400.0500
RuralOld700490.0700

Portfolio-wide frequency is fˉ=364/3600=0.0872\bar f = 364/3600 = 0.0872 claims per exposure. The raw cell frequencies suggest Rural business runs lighter than Urban, and Old vehicles run heavier than New, but the analyst wants a multiplicative structure that can be extended to many more variables without a combinatorial explosion of cells — the textbook motivation for a generalized linear model.

Model specification

The analyst fits a Poisson GLM with a log link and an exposure offset, so that expected claim counts are

E[Ni]=eiβ‹…exp⁑(Ξ²0+Ξ²1β‹…1[Rural]i+Ξ²2β‹…1[Old]i),E[N_i] = e_i \cdot \exp(\beta_0 + \beta_1 \cdot \mathbb 1[\text{Rural}]_i + \beta_2 \cdot \mathbb 1[\text{Old}]_i),

where eie_i is the exposure in cell ii and Urban/New is the baseline level absorbed into Ξ²0\beta_0. This is equivalent to log⁑(E[Ni]/ei)=Ξ²0+Ξ²11[Rural]i+Ξ²21[Old]i\log(E[N_i]/e_i) = \beta_0 + \beta_1\mathbb 1[\text{Rural}]_i + \beta_2\mathbb 1[\text{Old}]_i β€” the offset log⁑ei\log e_i is added to the linear predictor with a coefficient fixed at 1, which is the standard GLM device for modeling a rate while keeping the response a Poisson count.

Fitting this saturated model (4 cells, 3 free parameters plus the fixed offset β€” one fewer than the number of cells, since a full interaction term would make it exactly saturated) by iteratively reweighted least squares gives:

ParameterEstimateeΞ²^e^{\hat\beta} (relativity)
Intercept (Ξ²0\beta_0)βˆ’2.4034-2.4034β€” (Urban/New base frequency =0.0904=0.0904)
Rural (Ξ²1\beta_1)βˆ’0.6048-0.60480.54620.5462
Old (Ξ²2\beta_2)0.35890.35891.43181.4318

The fitted expected counts are:

CellFitted countActual count
Urban/New108.49108
Urban/Old116.51117
Rural/New39.5140
Rural/Old49.4949

The fitted relativities say a Rural risk has 0.5462Γ—0.5462\times the frequency of an otherwise-identical Urban risk (a 46% discount), and an Old vehicle has 1.4318Γ—1.4318\times the frequency of an otherwise-identical New vehicle (a 43% surcharge), holding the other factor fixed β€” this is exactly the "holding all else equal" interpretation that a GLM provides and a one-way (univariate) analysis cannot, because a univariate Rural/Urban split would be confounded by the fact that Rural business happens to skew towards newer vehicles in this book.

Deviance and goodness of fit

The model's residual deviance is D=0.0154D = 0.0154 on 4βˆ’3=14 - 3 = 1 degree of freedom (this model has one fewer parameter than cells, so it is nearly saturated and fits almost perfectly by construction). The null deviance (intercept-only model, i.e. every cell predicted at the overall frequency) is D0=34.08D_0 = 34.08 on 33 degrees of freedom. The drop in deviance from adding Territory and Rural,

Ξ”D=D0βˆ’D=34.08βˆ’0.0154=34.07Β onΒ 3βˆ’1=2Β d.f.,\Delta D = D_0 - D = 34.08 - 0.0154 = 34.07 \text{ on } 3-1=2 \text{ d.f.},

is large relative to a Ο‡22\chi^2_2 reference distribution (the 99th percentile of Ο‡22\chi^2_2 is about 9.219.21), so both rating variables are statistically significant contributors to explaining claim frequency in this data β€” consistent with, but not proof of, causal risk differentiation (omitted-variable confounding with e.g. annual mileage or driver age is always possible with only two rating variables).

Your task

You are reviewing this GLM output for the rate filing. Treat the fitted coefficients, relativities and deviance figures above as given (do not refit). Answer the questions below using GLM theory and the specific numbers in this case.

Part A β€” Multiple choice (6 Γ— 1 point)

Q1
Log link & offsets

In the Poisson GLM E[Ni]=eiexp⁑(β0+β1X1i+β2X2i)E[N_i] = e_i\exp(\beta_0+\beta_1 X_{1i}+\beta_2 X_{2i}), what role does eie_i (earned exposure) play?

Q2
Relativities

The fitted Rural relativity is eΞ²^1=0.5462e^{\hat\beta_1}=0.5462. What is the correct interpretation?

Q3
GLM structure

Why is the Poisson distribution (rather than Normal/OLS) the natural choice for modeling claim counts here?

Q4
Deviance & goodness of fit

The drop in deviance from the null model to the fitted model is 34.08βˆ’0.0154=34.0734.08-0.0154=34.07 on 22 degrees of freedom. What does comparing this to a Ο‡22\chi^2_2 distribution test?

Q5
Relativities

A one-way univariate analysis (ignoring vehicle age) shows Rural frequency 0.05930.0593 vs Urban 0.10710.1071, a raw ratio of 0.5540.554 β€” close to, but not identical to, the GLM's multivariate relativity of 0.54620.5462. What actuarial concept explains the small remaining difference?

Q6
GLM structure

If the analyst instead fit this model with a log link but Gamma error distribution, what would typically change about its suitability here?

Part B β€” Written response

W1
Relativities
6 points

Explain, using the specific numbers in this case, why a multivariate GLM relativity for Territory can differ from a univariate one-way relativity, and why an actuary should generally prefer the multivariate estimate for ratemaking.

0 words
W2
Deviance & goodness of fit
6 points

The fitted model in this case has 3 parameters for 4 cells (one degree of freedom of residual deviance). Discuss the practical limitation this creates for assessing goodness-of-fit here, and describe how a real ratemaking GLM with dozens of rating variables and thousands of cells avoids this problem.

0 words
Answer every question to submit.
Ask the tutor