Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Decision Robustness to Likelihood Misspecification

Bayesian Multi-Arm Approaches · Tutorial 421 of 499

Decision Robustness to Likelihood Misspecification

You will learn to compare posterior arm decisions under plausible alternative likelihoods and distinguish a stable recommendation from one driven by modeling assumptions.

Advanced 12 min read

What You'll Learn

  • Define likelihood misspecification and decision robustness for Bayesian arm selection.
  • Compare posterior decisions while holding the prior, data, and utilities fixed.
  • Calculate how heavy tails can change posterior odds and the selected arm.
  • Assess a binary decision’s sensitivity to overdispersion using an approximate likelihood.
  • Summarize stability using winner changes, utility ranges, and regret.
  • Recognize diagnostics and limitations of likelihood sensitivity analyses.

Why Examine the Likelihood?

A Bayesian arm decision depends on the observed data, the prior distribution, the likelihood, and the utility assigned to each possible action. A recommendation can therefore be sensitive to the assumed likelihood even when the data, prior, and decision rule are unchanged. This matters in a multi-arm, multi-stage trial: a model that treats observations as independent when they are clustered, or that gives extreme observations too much influence, can change which arm appears preferable at an interim analysis.

Here, likelihood misspecification means that the probability model used for the observed outcomes does not adequately describe their sampling distribution conditional on the model parameters. Decision robustness is not a claim that one likelihood is certainly correct. It is an assessment of whether the arm decision remains useful and stable across a prespecified set of scientifically plausible likelihoods.

Definition: A Bayesian arm decision is robust to likelihood misspecification over a model set if its recommendation, or its decision-relevant utility, changes little when the likelihood is replaced by plausible alternatives while the observed data, prior, utility specification, and decision rule are held fixed. Robustness is relative to the models examined; it does not establish that any model is correct.

The comparison should be controlled. If the prior or utility changes at the same time as the likelihood, a changed recommendation cannot be attributed to likelihood sensitivity alone. Begin with a primary likelihood and specify alternatives that address a concrete concern, such as heavy tails, excess binary-response variation, or clustering. Fit each model to the same data and apply the same decision rule, such as selecting the action with greatest posterior expected utility.

Posterior Decisions Across Likelihoods

Let \(a\) denote an available action, such as selecting an experimental arm or retaining control, and let \(U(a,\theta)\) be its utility at the unknown state \(\theta\). Under likelihood model \(M_k\), the posterior expected utility is

$$ \operatorname{EU}_k(a)=E_k\{U(a,\theta)\mid y\} =\int U(a,\theta)\,\pi_k(\theta\mid y)\,d\theta, $$

where \(y\) is the observed data and \(\pi_k(\theta\mid y)\) is the posterior under model \(M_k\). The decision under that model is an action maximizing \(\operatorname{EU}_k(a)\). A convenient first summary is whether the maximizing action changes. Also report utility differences: a winner change between two nearly tied actions is less consequential than a large loss in expected utility.

For two hypotheses or states, posterior odds show directly why a likelihood change can matter. If the prior odds are fixed, then

$$ \frac{P(H_1\mid y,M_k)}{P(H_0\mid y,M_k)} = \frac{P(H_1)}{P(H_0)} \frac{P(y\mid H_1,M_k)}{P(y\mid H_0,M_k)}. $$

The likelihood ratio can change sharply when a model assigns different relative probabilities to unusual observations. In a multi-arm decision, the same principle applies through the full posterior distribution and the utility integral: changes to the likelihood can alter posterior location, uncertainty, or both.

The comparison assumes that the candidate models describe the same outcome and estimand, that their likelihoods are well-defined for the observed data, and that prior and utility choices are held fixed. If an alternative model changes the estimand—for example, from an individual-level response to a cluster-level response—then the decision comparison needs an explicit mapping to a common utility scale.

Worked Example: A Heavy-Tailed Likelihood Changes the Decision

Suppose an arm’s treatment effect \(\theta\), measured as a benefit difference in outcome units, has a discrete prior: \(P(\theta=1)=P(\theta=-1)=0.5\). The action is to select the arm if \(P(\theta>0\mid y)>0.5\); otherwise select control. Three hypothetical independent treatment-control contrasts are \(1,1,-5\). Compare a normal likelihood \(d_i\mid\theta\sim N(\theta,1)\) with a Student \(t\) likelihood with 3 degrees of freedom and scale 1, independently for each contrast. The prior and decision rule are identical under both models.

Under the normal model, the likelihood ratio for \(\theta=1\) against \(\theta=-1\) is

$$ \frac{L(1)}{L(-1)} =\exp\left\{\sum_{i=1}^3 \frac{(d_i+1)^2-(d_i-1)^2}{2}\right\} =\exp\{2(1+1-5)\} =\exp(-6)\approx 0.002479. $$

The prior odds are 1, so the posterior odds are also approximately 0.002479. Thus \(P(\theta=1\mid y)=0.002479/(1+0.002479)\approx0.002473\). The rule selects control.

For the \(t_3\) likelihood, its density is proportional to \(\{1+(d_i-\theta)^2/3\}^{-2}\). For each observation equal to 1, the likelihood ratio is \((7/3)^2=49/9\). For the observation equal to \(-5\), it is

$$ \frac{1+( -5+1)^2/3}{1+( -5-1)^2/3} \text{ squared} = \left(\frac{19/3}{13}\right)^2 =\left(\frac{19}{39}\right)^2. $$

Consequently, the likelihood ratio for all three observations is

$$ \frac{L(1)}{L(-1)} =\left(\frac{49}{9}\right)^2\left(\frac{19}{39}\right)^2 \approx 7.035. $$

With equal prior odds, \(P(\theta=1\mid y)=7.035/(1+7.035)\approx0.8755\). The rule selects the experimental arm. This is a genuine decision reversal: the normal likelihood lets the single extreme negative contrast dominate, whereas the heavier-tailed model gives that observation less relative influence. Neither calculation establishes which model is appropriate; the reversal is a signal that the decision is model-sensitive.

Binary Overdispersion and Decision Sensitivity

For binary outcomes, the ordinary binomial likelihood assumes conditionally independent responses with a common event probability within each arm. Clustering, repeated measurements, or unmodeled heterogeneity can make outcomes more variable than that likelihood allows. A beta-binomial model is one possible likelihood alternative when cluster-level variation is scientifically plausible. A design-effect adjustment can also serve as a transparent approximate sensitivity calculation, but it is not a substitute for fitting a coherent overdispersed likelihood when one is required.

For clusters of size \(m\) with within-cluster correlation \(\rho\), the usual design effect is \(D=1+(m-1)\rho\). It approximates the variance inflation relative to independent observations. The effective sample size approximation is \(n_{\mathrm{eff}}=n/D\). It is useful for illustrating sensitivity, provided its approximation and assumptions are stated rather than treating fractional effective counts as actual observed data.

Worked Example: Overdispersion Reverses a Utility Choice

Suppose a hypothetical binary-outcome trial observes response proportions of \(0.80\) in an experimental arm and \(0.50\) in control, with 10 participants per arm. The estimated risk difference is \(\hat{\delta}=0.80-0.50=0.30\). For a sensitivity calculation, use the normal approximation \(\hat{\delta}\mid\delta\sim N(\delta,V)\), a prior \(\delta\sim N(0,0.5^2)\), and utility \(U(\text{experimental},\delta)=\delta-0.245\) versus \(U(\text{control},\delta)=0\). Select the experimental arm when its posterior expected utility is positive.

Under independent binomial sampling, the plug-in variance is

$$ V_{\mathrm{ind}}=\frac{0.80(0.20)}{10}+\frac{0.50(0.50)}{10} =0.016+0.025=0.041. $$

For this normal likelihood and normal prior, the posterior mean is \(\{0.25/(0.25+V)\}\hat{\delta}\). Thus the posterior mean risk difference is \(0.25(0.30)/(0.25+0.041)\approx0.2577\). Expected utility for the experimental arm is \(0.2577-0.245=0.0127\), so this model selects the experimental arm.

Now suppose responses are clustered in groups of size \(m=5\), with a plausible within-cluster correlation \(\rho=0.20\). The design effect is \(D=1+(5-1)(0.20)=1.8\). Inflating the variance gives \(V_{\mathrm{over}}=1.8(0.041)=0.0738\). The posterior mean becomes \(0.25(0.30)/(0.25+0.0738)\approx0.2316\), and expected utility is \(0.2316-0.245=-0.0134\). This approximate overdispersion analysis selects control instead. The finding is not that the adjusted model is necessarily correct; it is that this decision lies close enough to the utility threshold to depend on the assumed sampling variation.

How to Assess and Report Robustness

A useful analysis specifies the model set before interpreting the results. Alternatives should correspond to plausible failure modes, not be selected merely because they produce a preferred arm. Keep the outcome definition, prior, utility, and action set fixed where possible. For each model, record the posterior expected utility of every action, the winning action, and the difference between the best and next-best action. For a small model set, the range of each action’s expected utility and its worst-case regret are informative summaries. Regret under model \(M_k\) is the expected-utility loss from choosing \(a\) rather than the action that would maximize utility under that model.

1
Name the concern.
For example, identify extreme residuals, clustering, or excess variation in binary responses as the reason to consider an alternative likelihood.
2
Specify plausible models.
Choose likelihoods that address the concern and describe their assumptions, such as a Student \(t\) error model or a beta-binomial model.
3
Recompute the same decision.
Use the same data, prior, utility, and decision rule under each model. Check numerical integration or simulation error when posterior quantities are computed approximately.
4
Report stability and consequences.
State whether the selected arm changes and quantify differences in expected utility, not just whether model-specific posterior summaries differ.

Diagnostics help decide which alternatives deserve consideration. For continuous outcomes, inspect residual patterns and whether a few extreme values dominate arm contrasts; compare predictive distributions for tail behavior. For binary data, examine variability across clusters or sites and compare observed variation with that predicted under conditional independence. Posterior predictive checks can reveal that a likelihood fails to reproduce features of the observed data, but agreement with such checks does not prove the model is correct.

In a staged trial, the data available at an interim analysis may be limited, and selected arms can have unusually favorable observed results by chance. Repeating the same decision calculation under plausible likelihoods makes this dependence visible, but it does not remove selection effects or guarantee that an interim recommendation will remain appropriate as more data accrue. The analysis also cannot protect against every form of misspecification: the model set may omit the relevant alternative, and a robust decision across a narrow set can still be fragile outside it.

Common Pitfalls and Alternatives

A common mistake is to call a decision robust because the estimated treatment effect is similar across models. The posterior uncertainty may have changed enough to alter the probability of benefit or the expected utility, even if posterior means are close. Conversely, a winner can change when two actions have nearly equal utility, while the practical cost of switching is negligible. Always pair the identity of the selected action with the size of the utility difference.

Another mistake is to change the prior, likelihood, and utility together, then attribute the resulting decision change to the likelihood. Separate sensitivity analyses clarify which assumption matters. Model averaging is an alternative when probabilities over candidate models can be justified: it combines model-specific posterior distributions using model probabilities. It answers a different question from checking whether recommendations agree model by model. Robust likelihood procedures, such as heavy-tailed errors, can reduce sensitivity to specific data features, but they still rely on assumptions and do not replace checking whether the resulting model fits the scientific setting.

A likelihood sensitivity analysis is not a formal guarantee of good operating characteristics. In a MAMS design, trial-level performance across possible data-generating scenarios requires design evaluation, which is distinct from examining the posterior decision for one observed data set. Report the candidate models, the reason each is plausible, and the decision consequences; avoid presenting one preferred model as validated solely because it yields a stable answer.

Key takeaway: A posterior arm recommendation is conditional on its likelihood. Robustness analysis holds the prior, utility, and decision rule fixed, varies plausible likelihood assumptions, and reports both whether the chosen action changes and how much expected utility is at stake.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. Why must the prior and utility be held fixed when isolating sensitivity to the likelihood?
  2. In the heavy-tailed example, which observation most strongly drives the normal-model decision, and why?
  3. For clusters of size 4 with within-cluster correlation \(0.10\), calculate the design effect and effective sample size for \(n=40\).
  4. If the selected arm changes across plausible likelihoods but the expected-utility difference is very small, what should be reported?

Solutions. (1) Holding them fixed makes the likelihood the only changed component, so any decision difference can be attributed to the sampling model rather than a simultaneous change in beliefs or preferences. (2) The contrast \(-5\) is extreme relative to the two positive contrasts. Under the normal likelihood it has a large influence on the likelihood ratio, while the \(t_3\) model assigns relatively more probability to such a tail observation. (3) \(D=1+(4-1)(0.10)=1.30\), so \(n_{\mathrm{eff}}=40/1.30\approx30.77\) participants. This is an approximation, not a literal replacement data set. (4) Report that the model-specific winner changes, give expected utilities or their differences, and explain that the practical consequence may be small because the actions are nearly tied. Do not describe the recommendation as wholly stable or imply that the models are equivalent.

Clinical Biostats