Why Examine Prior Sensitivity?
A Bayesian design ranks candidate sampling schedules by averaging their performance over a prior distribution for the unknown model parameters. The resulting ranking can be useful, but it is conditional on the prior used. If plausible prior assumptions place more probability on different parameter values, the preferred schedule may change. Prior sensitivity analysis makes this dependence visible before a schedule is selected.
The question is not whether a prior is “correct” in an absolute sense. Rather, it is whether the design recommendation remains similar across a defensible range of prior locations, scales, and dependence structures. A design that wins under all such alternatives has a different evidential basis from one that wins only under a narrow prior specification.
Let \(d\) denote a candidate sampling design, \(\theta\) the unknown model parameters, \(U(d,\theta)\) the utility of design \(d\) at parameter value \(\theta\), and \(\pi(\theta)\) a prior density. The Bayesian expected utility is
For a discrete prior with parameter scenarios \(\theta_j\) and probabilities \(w_j\), the same quantity is \(\operatorname{EU}(d;\pi)=\sum_j w_j U(d,\theta_j)\), where \(w_j\geq 0\) and \(\sum_j w_j=1\). The preferred design maximizes expected utility when larger utility is better. In a sensitivity analysis, the candidate designs, utility definition, model, and constraints should be held fixed while the prior is varied, unless the purpose is explicitly to examine several changes together.
A Controlled Comparison
A useful analysis changes one prior feature at a time where possible. To examine location, shift the prior center while retaining its scale and dependence structure. To examine scale, widen or narrow the prior while keeping its center fixed. To examine correlation, vary the joint dependence while preserving the marginal distributions. These controlled comparisons help identify which assumption is responsible for a change in ranking.
Prior alternatives should be scientifically defensible, not selected merely to make a preferred design win. Plausible parameter ranges, previous knowledge, and elicited uncertainty can inform the alternatives. The earlier material on selecting plausible parameter ranges provides a basis for choosing such alternatives; the sensitivity analysis then asks how the design behaves across them.
For two designs \(d_A\) and \(d_B\), define their expected-utility difference under prior \(\pi_k\) as \(\Delta_k=\operatorname{EU}(d_A;\pi_k)-\operatorname{EU}(d_B;\pi_k)\). A positive value favors \(d_A\), a negative value favors \(d_B\), and a value near zero indicates little separation under that prior. The sign and magnitude of this difference are more informative than a rank label alone.
Worked Example: Changing Prior Location
Suppose a hypothetical model calculation gives two candidate schedules utility scores, in arbitrary information-score units, as functions of a standardized parameter \(X\): \(U_A(X)=12-(X-1)^2\) and \(U_B(X)=11-(X+1)^2\). Consider a family of discrete priors placing probability \(1/2\) at each of \(m-0.5\) and \(m+0.5\). Thus, each prior has location \(m\) and the same variance, \(0.25\).
For this prior, \(E[(X-a)^2]=(m-a)^2+0.25\): the mean squared deviation is the squared distance from the prior mean plus the prior variance. Therefore,
At \(m=0\), the calculations are \(12-(1+0.25)=10.75\) for \(A\) and \(11-(1+0.25)=9.75\) for \(B\); \(A\) ranks higher. At \(m=1\), they are \(12-(0+0.25)=11.75\) and \(11-(4+0.25)=6.75\), again favoring \(A\). At \(m=-1\), they are \(12-(4+0.25)=7.75\) and \(11-(0+0.25)=10.75\), so \(B\) ranks higher.
The candidate schedules and their utility functions have not changed; only the prior location has. The ranking reversal shows that the design recommendation depends on where prior probability is placed. The scores are hypothetical and illustrate the calculation, not results from an actual study.
Scale and Correlation Sensitivity
A prior’s scale represents the degree of uncertainty about parameter values, not simply the location of its center. If a design’s utility rises or falls sharply as a parameter moves away from the center, widening the prior can change its average utility substantially. Correlation matters when joint parameter combinations affect utility: priors with identical marginal distributions can assign different probabilities to combinations such as high clearance with high volume.
Worked Example: Changing Prior Scale
Let \(X\sim N(0,s^2)\) represent a hypothetical standardized parameter deviation. Suppose two candidate schedules have utility functions \(U_A(X)=12-2X^2\) and \(U_B(X)=10\), in information-score units. The prior mean remains zero while its standard deviation \(s\) changes. The normal prior has \(E[X^2]=s^2\), so
At \(s=0.5\), \(\operatorname{EU}(A)=12-2(0.5^2)=12-0.5=11.5\), compared with \(10\) for \(B\); \(A\) is preferred. At \(s=1.0\), the scores are \(12-2(1)=10\) and \(10\), a tie. At \(s=1.2\), they are \(12-2(1.44)=9.12\) and \(10\), so \(B\) is preferred. The ranking is therefore scale-sensitive around \(s=1\). These are hypothetical utility functions and values.
This result follows from the prior expectation, not from a claim that a wider prior always favors one type of schedule. The direction of change depends on each design’s utility across the parameter space.
Worked Example: Changing Prior Correlation
Consider two standardized parameters \(X\) and \(Y\), each taking values \(-1\) and \(+1\) with marginal probability \(1/2\). Construct a joint prior with correlation \(\rho\): the total probability of equal signs is \((1+\rho)/2\), split equally between \((+1,+1)\) and \((-1,-1)\); the total probability of opposite signs is \((1-\rho)/2\), split equally between \((+1,-1)\) and \((-1,+1)\). This is a valid joint distribution for \(-1\leq\rho\leq1\); its marginals remain unchanged and \(E[XY]=\rho\).
Suppose hypothetical utilities are \(U_A(X,Y)=12+2XY\) and \(U_B(X,Y)=12\). Then
For \(\rho=0.6\), \(A\) has expected utility \(12+2(0.6)=13.2\), while \(B\) has \(12\), so \(A\) ranks higher. For \(\rho=-0.6\), \(A\) has \(12+2(-0.6)=10.8\), and \(B\) remains at \(12\), so \(B\) ranks higher. At \(\rho=0\), they tie. The ranking changes even though both parameter marginals are held fixed: the joint prior dependence matters for this utility. This discrete illustration is a sensitivity calculation, not a general model for pharmacokinetic parameters.
Calculation, Diagnostics, and Failure Modes
For complex models, expected utilities are often computed by numerical integration or Monte Carlo sampling. If \(\theta^{(r)}\) are draws from prior \(\pi_k\), an estimator is \(\widehat{\operatorname{EU}}(d;\pi_k)=R^{-1}\sum_{r=1}^R U(d,\theta^{(r)})\). To compare designs reliably, use the same parameter draws for every candidate schedule under a given prior. This paired calculation makes the estimated utility difference less noisy because variation shared across designs can cancel.
When comparing several priors, record their location, scale, and correlation assumptions explicitly. Recompute the expected utility for every candidate design under every prior, then report the utility differences and rankings. A table is often useful: it shows both ranking changes and near-ties that a rank-only summary would conceal.
Fix the candidate schedules, model, utility definition, and design constraints. State which prior feature is being varied.
Give their locations, scales, and dependence assumptions, and verify that each is a valid probability distribution.
Evaluate each candidate’s expected utility under each prior, then calculate pairwise utility differences and rankings.
Increase integration accuracy or Monte Carlo sample size, use shared draws for competing designs, and verify that apparent ranking changes exceed numerical error.
Identify which prior changes alter the recommended schedule, and distinguish a strong reversal from a small difference near a tie.
Several diagnostics help identify a misleading result. First, check that prior weights sum to one, that covariance matrices are valid, and that parameter scales use the intended units or parameterization. Second, inspect whether a small utility difference is being treated as a decisive ranking when numerical uncertainty is of comparable size. Third, confirm that the same utility definition and candidate set were used throughout. If a comparison changes both prior scale and utility criterion, it cannot isolate the effect of prior scale.
A particularly common mistake is to vary prior correlation while accidentally changing the marginal distributions as well. The resulting ranking change cannot then be attributed to correlation alone. The correction is to construct alternatives with the same marginal distributions and different valid dependence structures, or to state clearly that the sensitivity comparison changes more than one feature.
How to Interpret the Results
Prior sensitivity analysis complements rather than replaces other design comparisons. A single-prior Bayesian expected-utility analysis answers which candidate has the highest average utility under that specified prior. Scenario-based robust design compares performance across selected parameter scenarios, while prior sensitivity asks whether the Bayesian ranking persists when the prior itself is changed. These approaches answer related but distinct questions.
A stable ranking does not prove that the prior assumptions are correct, nor does it guarantee good performance under every possible parameter value. It indicates only that the ranking was stable over the alternatives examined and for the chosen utility. Conversely, a ranking reversal does not automatically invalidate Bayesian design: it signals that the recommendation depends on assumptions that should be made transparent and, where possible, investigated further.
Report the actual prior alternatives and utility differences, not just “robust” or “sensitive.” Mention the parameterization, units, and what was held fixed. For example, a conclusion might state that schedule \(A\) ranked first under prior centers from one range, while schedule \(B\) ranked first under centers from another, with the differences measured in a specified information-score unit. Avoid implying that the probability a design is best is given by its rank under a prior; that is a different quantity and requires a separate probabilistic definition.
The main limitation is that sensitivity analysis explores only the alternatives selected. It cannot protect against implausible or omitted prior structures, misspecified utility functions, or errors in the model used to calculate utility. It can also be computationally expensive when each utility evaluation involves a complex model. These limitations make transparent prior specification and numerical checks essential.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- Two designs have expected utilities of 14 and 12 under one prior, but 10 and 11 under another. What changes, and what should be reported?
- If \(X\sim N(0,s^2)\) and a design has utility \(U(X)=9-X^2\), what is its expected utility when \(s=0.8\)?
- Why should a correlation sensitivity analysis hold the marginal distributions fixed when the aim is to isolate correlation?
- In a Monte Carlo comparison, why can using the same prior draws for two candidate designs improve the estimate of their utility difference?
Solutions. (1) The first design ranks higher under the first prior, \(14>12\), while the second ranks higher under the second, \(11>10\). Report both prior specifications, all four utilities, and the reversal; do not report a single unconditional winner. (2) Since \(E[X^2]=s^2\), the expected utility is \(9-s^2=9-0.8^2=9-0.64=8.36\) utility units. (3) Otherwise a changed ranking could reflect changes in marginal location or scale as well as dependence, so the effect of correlation would not be isolated. (4) Shared draws pair the design evaluations at the same parameter values. Variation common to both estimates can cancel in their difference, often reducing Monte Carlo noise; numerical stability should still be checked.