Tutorials › CAS Exam 8 - Advanced Ratemaking › Uncertainty in Validation Performance Comparisons

CAS Exam 8 · Out-of-Sample Validation · Tutorial 264 of 989

Uncertainty in Validation Performance Comparisons

Use bootstrap resampling to measure uncertainty in validation performance, compare paired predictions fairly, and account for dependent observations within policies or risk clusters.

Intermediate 11 min read

What You'll Learn

  • Define a validation score statistic and its bootstrap standard error
  • Construct an approximate percentile interval from bootstrap replicates
  • Use paired resampling to estimate uncertainty in score differences
  • Resample whole policies or risk clusters when records are dependent
  • Distinguish validation-sample uncertainty from model-fitting and temporal uncertainty

Why Validation Scores Need an Uncertainty Measure

A validation score summarizes how a model performed on a particular set of observations. It is an estimate, not a fixed property of the model: a different validation sample from the same target population would generally produce a different score. That variation matters when a score is used to assess predictive performance or to describe a difference between two models. A small difference in observed scores may be less informative than a larger difference if the validation estimate is highly variable.

Bootstrap resampling provides a practical way to estimate this variability from the observed validation data. The key design choice is the resampling unit. If validation observations are independent policies, resample policies. If several records belong to the same policy or risk and are dependent, resample the whole policy or risk cluster, retaining all its records together. Resampling dependent records one at a time can understate uncertainty.

Definition: Let \(T\) be a validation performance statistic computed from the validation sample, such as an average loss or an exposure-weighted score. A bootstrap replicate is the value of that statistic after drawing a new sample with replacement from the observed independent validation units and recomputing the statistic. The standard deviation of many bootstrap replicates estimates the sampling uncertainty of \(T\), conditional on the fitted model or models.

Bootstrap Procedure for a Validation Score

Suppose there are \(n\) independent validation units and unit \(i\) has loss score \(L_i\) and, when appropriate, weight \(w_i\). An exposure-weighted aggregate score is

$$ T=\frac{\sum_{i=1}^{n} w_i L_i}{\sum_{i=1}^{n} w_i}. $$

The weights and the score definition must match the intended performance measure. For example, an average policy-level score gives each policy equal weight, whereas an exposure-weighted score gives more influence to units with more exposure. The bootstrap must recompute the same statistic, including its weights, in every replicate.

For bootstrap replicate \(b\), draw \(n\) units with replacement from the observed units, and calculate \(T^{*(b)}\). If \(B\) replicates are generated, the bootstrap standard error is the sample standard deviation of those replicate statistics:

$$ \widehat{\operatorname{SE}}_{\mathrm{boot}}(T) = \sqrt{\frac{1}{B-1}\sum_{b=1}^{B} \left(T^{*(b)}-\bar{T}^{*}\right)^2}, \qquad \bar{T}^{*}=\frac{1}{B}\sum_{b=1}^{B}T^{*(b)}. $$

This standard error describes the spread of bootstrap scores. It is not the standard deviation of the bootstrap scores divided by \(\sqrt{B}\); that latter quantity concerns Monte Carlo error in estimating the bootstrap mean from a finite number of replicates. In practice, use enough replicates that Monte Carlo noise is small relative to the uncertainty being reported.

1
Specify the target statistic.
Define the validation score, its direction of preference, and any weights before resampling.
2
Identify independent units.
Use individual observations only when they are independent; otherwise identify the policy or risk clusters that must remain intact.
3
Generate bootstrap replicates.
Sample independent units with replacement and recompute the complete statistic for each resample.
4
Summarize uncertainty.
Use the replicate standard deviation as the bootstrap standard error. An approximate percentile interval uses the lower and upper quantiles of the bootstrap replicate distribution, commonly the 2.5th and 97.5th percentiles for a 95% interval.

A percentile interval is an approximate interval, not a guarantee of exact coverage. It can be affected by bias, skewness, a small number of independent units, or a score with unusual boundary behavior. A large number of bootstrap replicates reduces simulation noise but does not create more independent information than the original validation sample contains.

Why the Bootstrap Standard Deviation Estimates Uncertainty

The bootstrap treats the observed validation units as an empirical approximation to the population from which they were sampled. Resampling from that empirical distribution mimics the act of drawing another validation sample. The resulting variation in \(T^{*}\) approximates the variation that \(T\) would have across repeated samples.

The variance formula follows directly from the usual sample variance applied to the bootstrap statistics. If the empirical bootstrap distribution assigns probability \(1/B\) to each replicate, its center is \(\bar{T}^{*}\); the squared deviations from that center measure how much the statistic changes across resamples. Dividing their sum by \(B-1\) gives the sample variance of the \(B\) simulated values, and taking its square root gives their standard deviation. As \(B\) increases, this calculation more stably estimates the spread of the bootstrap distribution. It does not remove uncertainty due to having observed only a finite validation sample.

Paired Uncertainty in Score Differences

When two models are scored on the same validation units, uncertainty in their difference should preserve that pairing. Let \(L_{Ai}\) and \(L_{Bi}\) be the losses for models A and B on unit \(i\). Define \(\Delta=T_A-T_B\). If lower scores are better, a positive \(\Delta\) favors model B. For each bootstrap replicate, draw one resample of validation units and calculate both model scores on those same units; then subtract the scores. The standard deviation of those replicate differences estimates uncertainty in \(\Delta\).

The reason pairing matters is that the two scores on the same unit are often related. For simple averages, the variance of their difference is

$$ \operatorname{Var}(T_A-T_B) = \operatorname{Var}(T_A)+\operatorname{Var}(T_B) -2\operatorname{Cov}(T_A,T_B). $$

Using common resamples retains this covariance. Resampling separately for each model discards the pairing and estimates a different, generally less appropriate quantity. The paired bootstrap estimates uncertainty in the difference on common validation units; it does not by itself establish that the models or their validation periods are otherwise comparable.

Worked Example: Bootstrap Uncertainty for an Average Score

A hypothetical insurer evaluates a model using one loss score per independent validation policy. The four observed scores are 4, 6, 8, and 10, with lower scores preferred. The observed average is \((4+6+8+10)/4=7\). For a hand-checkable illustration, four bootstrap samples select policy indices \((1,1,2,3)\), \((1,2,3,4)\), \((2,3,4,4)\), and \((1,3,3,4)\). Estimate the bootstrap standard error from these replicates.

Method. Resample policies with replacement and recompute the average score. The four replicate averages are:

  • Indices \((1,1,2,3)\): \((4+4+6+8)/4=22/4=5.5\).
  • Indices \((1,2,3,4)\): \((4+6+8+10)/4=28/4=7.0\).
  • Indices \((2,3,4,4)\): \((6+8+10+10)/4=34/4=8.5\).
  • Indices \((1,3,3,4)\): \((4+8+8+10)/4=30/4=7.5\).

The mean of the four replicates is \((5.5+7.0+8.5+7.5)/4=28.5/4=7.125\). The squared deviations from that mean are \(2.640625\), \(0.015625\), \(1.890625\), and \(0.140625\); their sum is \(4.6875\). Thus the bootstrap variance estimate is \(4.6875/(4-1)=1.5625\), and the bootstrap standard error is \(\sqrt{1.5625}=1.25\) score units.

The result is a demonstration using only four replicates, not a reliable final estimate for a real analysis. A practical calculation would use many more replicates. The common mistake is to divide the replicate standard deviation by \(\sqrt{B}\) and report that as the bootstrap standard error of the original validation score. That would instead estimate Monte Carlo error in the replicate mean.

Worked Example: Paired Bootstrap for Two Models

Two hypothetical rate models are evaluated on the same four policies. Their loss scores are A: 9, 12, 8, 11 and B: 7, 13, 5, 11. Use the paired differences to estimate the uncertainty in the average score difference, given resamples with indices \((1,1,2,4)\), \((1,2,3,4)\), \((2,2,3,4)\), and \((1,3,3,4)\). Positive differences \(A-B\) favor B.

Method. The observed differences \(A-B\) by policy are \(2,-1,3,0\), so the observed average difference is \((2-1+3+0)/4=1\). Resample the same policy indices for both models, which is equivalent to resampling these paired differences:

  • \((1,1,2,4)\): \((2+2-1+0)/4=3/4=0.75\).
  • \((1,2,3,4)\): \((2-1+3+0)/4=4/4=1.00\).
  • \((2,2,3,4)\): \((-1-1+3+0)/4=1/4=0.25\).
  • \((1,3,3,4)\): \((2+3+3+0)/4=8/4=2.00\).

The replicate mean is \((0.75+1.00+0.25+2.00)/4=1.00\). The squared deviations are \(0.0625\), \(0\), \(0.5625\), and \(1.0000\), summing to \(1.625\). The sample variance is \(1.625/3=0.5416667\), so the paired bootstrap standard error is \(\sqrt{0.5416667}=0.7360\), rounded to four decimals. The observed difference of 1.00 favors B, but this small illustrative set is not enough to establish a precise population difference.

The mistake to avoid is independently resampling policies for A and B. That would compare scores from different sets of policies and discard the within-policy relationship that the paired difference is intended to retain.

Resampling Policies or Risk Clusters

A policy can contribute multiple records, such as multiple coverages, locations, or periods. If those records are dependent, the independent unit for the bootstrap is the policy, not the record. Draw policies with replacement and include all records belonging to each selected policy. When a policy is selected more than once, its full set of records appears that many times in the replicate. Recompute the score from the resampled records, including any exposure weighting.

The same principle applies when the independent unit is a broader risk cluster, such as an insured account whose policies share common experience. The resampling unit should reflect the dependence structure relevant to the validation sample. If the data have important time dependence or a sampling design that is not captured by simple independent clusters, an ordinary cluster bootstrap may not be adequate; the resampling method must be adapted to that design.

Worked Example: Exposure-Weighted Score With Cluster Resampling

A hypothetical validation set contains three independent policies. Their total score numerators and exposures are: policy 1, 10 score units and 2 exposure units; policy 2, 24 score units and 3 exposure units; policy 3, 25 score units and 5 exposure units. The target statistic is total score numerator divided by total exposure. Calculate the observed statistic and the bootstrap standard error for four specified cluster resamples: \((1,1,2)\), \((1,2,3)\), \((2,3,3)\), and \((1,3,3)\).

Method. Keep each selected policy's numerator and exposure together, then recompute the ratio. In the original data, the numerator is \(10+24+25=59\) and exposure is \(2+3+5=10\), so the observed score is \(59/10=5.9\).

  • Resample \((1,1,2)\): numerator \(10+10+24=44\), exposure \(2+2+3=7\), score \(44/7=6.2857\).
  • Resample \((1,2,3)\): numerator \(59\), exposure \(10\), score \(59/10=5.9000\).
  • Resample \((2,3,3)\): numerator \(24+25+25=74\), exposure \(3+5+5=13\), score \(74/13=5.6923\).
  • Resample \((1,3,3)\): numerator \(10+25+25=60\), exposure \(2+5+5=12\), score \(60/12=5.0000\).

The mean of the four replicate scores is approximately \((6.2857+5.9000+5.6923+5.0000)/4=5.7195\). Using the unrounded ratios, the sum of squared deviations from this mean is approximately \(0.8716\). The sample variance is \(0.8716/3=0.2905\), giving a bootstrap standard error of \(\sqrt{0.2905}=0.5390\) score units, rounded to four decimals.

The common mistake is to resample score records independently or to attach a selected policy's numerator to a different exposure. Either changes the dependence structure or the defined statistic. Here the numerator and exposure are resampled together as one policy-level unit.

Interpretation and Practical Limits

A bootstrap standard error describes estimated sampling variation under the chosen resampling scheme. It does not automatically include uncertainty from fitting the model. If the goal is to include training-sample variability, each replicate must reproduce the relevant model-fitting process as well as the validation calculation, with care to preserve the separation between training and validation information. If the fitted models are held fixed and only validation units are resampled, the estimate is conditional on those fitted models.

The bootstrap also relies on the validation units being representative of the population to which the score is intended to apply. It cannot correct an unrepresentative validation sample, immature outcomes, or a mismatched score definition. For rates or other statistics with strong skewness or few independent clusters, inspect the replicate distribution rather than relying on a standard error alone. A percentile interval can communicate asymmetry, but its interpretation remains approximate.

  • State whether lower or higher scores are preferred.
  • Use the same score definition and weighting in the original calculation and every replicate.
  • For model differences on common units, calculate paired replicate differences.
  • Resample whole independent clusters whenever records within a cluster are dependent.
  • Report the resampling unit, number of replicates, and whether model fitting was repeated.
Key takeaway: Bootstrap uncertainty is determined as much by the resampling unit and statistic as by the number of replicates. Preserve common units when comparing scores, preserve clusters when observations are dependent, and be explicit about whether uncertainty is conditional on fitted models.

Check Your Understanding

For each problem, identify the appropriate calculation or resampling unit before checking the final answer.

  1. A validation score has four bootstrap replicate values 3, 5, 7, and 5. Calculate their mean and their sample standard deviation. (Answer: Mean \(=5\); sample standard deviation \(=\sqrt{8/3}=1.6330\).)
  2. For four policies, paired model score differences \(A-B\) are \(1,-2,3,0\). Calculate the observed average difference and say which model it favors when lower scores are better. (Answer: Average difference \(=0.5\); it favors model B.)
  3. Three independent policy clusters have score numerators and exposures \((6,2)\), \((20,4)\), and \((18,3)\). In a bootstrap replicate, clusters 1, 2, and 2 are selected. Calculate the replicate's exposure-weighted score. (Answer: \((6+20+20)/(2+4+4)=46/10=4.6\).)
  4. Two models are scored on the same policies, and multiple records from each policy are dependent. State the resampling unit and how the two model scores should be paired in each replicate. (Answer: Resample whole policies with all their records intact, and calculate both model scores on the same sampled policies.)
  5. A report defines the variance of four bootstrap replicate scores using the denominator \(B=4\), and the sum of squared deviations from their mean is 15.5. Calculate the Monte Carlo standard error of the replicate mean using that variance estimate. (Answer: \(\sqrt{(15.5/4)/4}=\sqrt{0.96875}=0.9843\); this is Monte Carlo error of the replicate mean, not the bootstrap standard error of the original validation statistic.)
Clinical Biostats