4.2 Poisson Regression & Exposure Offsets

The Poisson Distribution for Insurance Claim Counts

In property and casualty insurance, claim frequency represents the rate at which insured losses occur over a defined period of coverage. The classical starting point for modeling claim counts is the Poisson distribution.

A discrete random variable $Y$ representing the number of claims reported on a policy follows a Poisson distribution with parameter $\lambda > 0$ if its probability mass function is given by:

$$P(Y = y) = \frac{e^{-\lambda} \lambda^y}{y!}, \quad y \in \{0, 1, 2, \dots\}$$

The defining characteristic of the Poisson distribution is equidispersion, meaning its expected value and variance are mathematically identical:

$$\mathbb{E}[Y] = \lambda$$ $$\text{Var}(Y) = \lambda$$

In an insurance portfolio, policyholders do not all share the same exposure to risk. A commercial fleet policy covering 100 power units over an entire policy year has a vastly higher expected claim count than a personal auto policy covering a single vehicle for one month. Therefore, Poisson modeling must explicitly incorporate the volume of earned exposure.

---

Incorporating Exposure: The Offset Derivation

Let:

• $Y_i$ denote the observed claim count for policy $i$ ($y_i \in \{0, 1, 2, \dots\}$).
• $e_i$ denote the earned exposure for policy $i$ (measured in earned car-years, policyholder years, house-years, or 100 units of payroll).
• $\mu_i$ denote the underlying expected claim frequency per unit of exposure (the annualized hazard rate).

The expected number of claims $\lambda_i = \mathbb{E}[Y_i]$ is directly proportional to earned exposure:

$$\mathbb{E}[Y_i] = \lambda_i = e_i \cdot \mu_i$$

To relate frequency $\mu_i$ to rating variables $\mathbf{x}_i = [1, x_{i1}, \dots, x_{ik}]^T$, we specify a Generalized Linear Model with the standard log link function:

$$\ln(\mu_i) = \mathbf{x}_i^T \boldsymbol{\beta} = \beta_0 + \beta_1 x_{i1} + \dots + \beta_k x_{ik}$$

Substituting this into the expected claim count expression and taking the natural logarithm:

$$\ln(\mathbb{E}[Y_i]) = \ln(e_i \cdot \mu_i) = \ln(e_i) + \ln(\mu_i)$$ $$\ln(\mathbb{E}[Y_i]) = \ln(e_i) + \mathbf{x}_i^T \boldsymbol{\beta}$$ $$\mathbb{E}[Y_i] = \exp\left( \ln(e_i) + \mathbf{x}_i^T \boldsymbol{\beta} \right) = e_i \cdot \exp(\mathbf{x}_i^T \boldsymbol{\beta})$$

Definition of an Offset

In statistical regression, an offset is defined as an explanatory variable whose regression coefficient is not estimated from data, but is constrained a priori to equal exactly 1.0.

Here, $\ln(e_i)$ acts as an offset. It enters the linear predictor with a fixed parameter of 1.0:

$$\eta_i = 1.0 \cdot \ln(e_i) + \beta_0 + \beta_1 x_{i1} + \dots + \beta_k x_{ik}$$

Why the Coefficient Must Equal 1.0: Exposure Additivity

Constraining the coefficient of $\ln(e_i)$ to 1.0 enforces strict exposure proportionality (homogeneity of degree 1):

1. If exposure doubles from $1.0$ car-year to $2.0$ car-years, the expected claim count must double exactly: $\mathbb{E}[Y] = 2.0 \cdot \mu$.
2. If a policyholder purchases insurance for two consecutive 6-month terms ($e_1 = 0.5, e_2 = 0.5$), the total expected claim count across both terms equals the expected claim count of a single 12-month policy ($e = 1.0$):

$$\mathbb{E}[Y_1] + \mathbb{E}[Y_2] = 0.5\mu + 0.5\mu = 1.0\mu = \mathbb{E}[Y_{\text{full year}}]$$

If an actuary mistakenly includes $\ln(e_i)$ as an unconstrained covariate and the regression estimates a coefficient $\hat{\beta}_e = 0.85$, the expected claim count scales as $e^{0.85}$. A 6-month policy would have an expected count factor of $0.5^{0.85} \approx 0.5548$. Two consecutive 6-month policies would produce an expected factor of $0.5548 + 0.5548 = 1.1096$—charging the policyholder roughly $11\%$ more than an identical 12-month policy. This violates fundamental actuarial principles of cost-causation and exposure additivity.

---

Frequency Modeling Alternatives: Count Offset vs. Weighted Frequency

In actuarial practice, there are two distinct ways to structure a Poisson GLM in statistical software. Both approaches are mathematically equivalent and yield identical parameter estimates, standard errors, and deviance statistics.

Approach 1: Raw Count with Offset

• Target variable: Raw claim counts $Y_i \in \{0, 1, 2, \dots\}$
• Distribution family: Poisson with log link
• Exposure specification: $\text{offset} = \ln(e_i)$

The log-likelihood function is:

$$\ell(\boldsymbol{\beta}) = \sum_{i=1}^n \left[ y_i \ln(e_i \mu_i) - e_i \mu_i - \ln(y_i!) \right]$$ $$\ell(\boldsymbol{\beta}) = \sum_{i=1}^n \left[ y_i \ln(e_i) + y_i \mathbf{x}_i^T \boldsymbol{\beta} - e_i \exp(\mathbf{x}_i^T \boldsymbol{\beta}) - \ln(y_i!) \right]$$

Differentiating with respect to the regression parameter $\beta_j$ gives the score equation:

$$\frac{\partial \ell}{\partial \beta_j} = \sum_{i=1}^n \left[ y_i x_{ij} - e_i x_{ij} \exp(\mathbf{x}_i^T \boldsymbol{\beta}) \right] = \sum_{i=1}^n \left( y_i - e_i \mu_i \right) x_{ij} = 0$$

Approach 2: Empirical Frequency with Prior Weights

• Target variable: Empirical claim frequency $F_i = \frac{Y_i}{e_i}$
• Distribution family: Poisson with log link
• Prior weight: $\omega_i = e_i$ (no offset)

In the Exponential Dispersion Family, the dispersion function with prior weights is $a(\phi) = \frac{\phi}{\omega_i} = \frac{1}{e_i}$ (since $\phi = 1$ for Poisson). The weighted quasi-log-likelihood for frequency $f_i$ is:

$$\ell_i(\boldsymbol{\beta}) = \frac{f_i \theta_i - b(\theta_i)}{a(\phi)} = e_i \left[ f_i \ln(\mu_i) - \mu_i \right] = e_i \left[ \left(\frac{y_i}{e_i}\right) \mathbf{x}_i^T \boldsymbol{\beta} - \exp(\mathbf{x}_i^T \boldsymbol{\beta}) \right]$$ $$\ell(\boldsymbol{\beta}) = \sum_{i=1}^n \left[ y_i \mathbf{x}_i^T \boldsymbol{\beta} - e_i \exp(\mathbf{x}_i^T \boldsymbol{\beta}) \right]$$

Differentiating with respect to $\beta_j$ gives the exact same score equation:

$$\frac{\partial \ell}{\partial \beta_j} = \sum_{i=1}^n \left( y_i - e_i \mu_i \right) x_{ij} = 0$$

Because the score equations and the negative Hessian matrices (Fisher Information) are identical up to constants independent of $\boldsymbol{\beta}$, both methods produce identical coefficients $\hat{\boldsymbol{\beta}}$ and covariance matrices $\widehat{\text{Var}}(\hat{\boldsymbol{\beta}})$.

Software Implementation Syntax Comparison

                               POISSON MODELING PATHS
                                         │
                  ┌──────────────────────┴──────────────────────┐
                  ▼                                             ▼
       Approach 1: Count + Offset                   Approach 2: Weighted Frequency
       Response: Y (Claims)                         Response: F = Y / Exposure
       Offset: ln(Exposure)                         Weights: Exposure
       Link: Log                                    Link: Log
                  │                                             │
                  └──────────────────────┬──────────────────────┘
                                         ▼
                            Score Equation: Σ (y_i - e_i * mu_i) * x_ij = 0
                            Identical Coefficients: beta_hat
                            Identical Relativities: exp(beta_hat)

---

Underlying Assumptions of the Poisson Regression Model

When applying a Poisson GLM to insurance frequency data, the actuary implicitly relies on four key statistical assumptions:

1. Independence of Events: The occurrence of a claim on policy $i$ is independent of claims on all other policies, and independent of prior claims on the same policy. This assumption is violated during catastrophic weather events (such as hail storms or freezes that damage thousands of homes simultaneously) or multi-claimant vehicle accidents.
2. Homogeneous Intensity Over Time: Within any policy's exposure period, the instantaneous rate of claim occurrence $\mu_i$ remains constant. In reality, claim intensity fluctuates with seasonal driving patterns, weather changes, and economic cycles.
3. Equidispersion: The conditional variance of claim counts equals the conditional mean ($\text{Var}(Y_i \mid \mathbf{x}_i) = \mathbb{E}[Y_i \mid \mathbf{x}_i]$). In empirical insurance data, this assumption almost never holds; variance typically exceeds the mean, a condition known as overdispersion.
4. Infinitesimal Probability: Over an infinitesimally small time interval $\Delta t$, the probability of exactly one claim is approximately $\mu_i \Delta t$, and the probability of two or more claims is negligible ($o(\Delta t)$).

---

Actuarial Exam Traps & Implementation Pitfalls