Tutorials › Bayesian Clinical Trials › A Regulatory Evidence Package for a Bayesian Trial

Case Studies · Tutorial 696 of 696

A Regulatory Evidence Package for a Bayesian Trial

Build a transparent Bayesian trial evidence package and distinguish what the observed-data posterior says from how the decision rule performs across repeated trials.

Advanced 12 min read

What You'll Learn

  • Specify a clinical estimand in terms of population, treatment conditions, outcome, time point, and summary measure.
  • Explain a prior distribution on the clinical effect scale and examine alternative priors.
  • Derive a normal posterior and calculate a posterior probability against a clinically meaningful threshold.
  • Convert a posterior decision rule into an observed-estimate cutoff without changing an inclusive inequality.
  • Calculate operating characteristics under fixed true treatment effects.
  • Separate posterior evidence, repeated-trial decision performance, and sensitivity to modeling choices.

Why Assemble an Evidence Package?

A regulatory conclusion should not depend on a single probability quoted without its clinical target, assumptions, or decision context. A useful Bayesian evidence package connects the question the trial was designed to answer with the prior information, the posterior estimate, the prespecified decision rule, and the behavior of that rule across possible repeated trials. These components answer different questions: the posterior summarizes uncertainty about the effect given the observed data and model, while operating characteristics describe how often a decision rule would be triggered under specified true effects.

The worked case below is hypothetical. It uses a simple normal model so that every calculation can be checked directly. Real submissions may require more elaborate outcome, missing-data, safety, and borrowing models, but the same reporting logic applies. The estimand must remain clinically interpretable, and the statistical model and decision rule must be stated precisely enough to reproduce the results.

Start With the Estimand and Model

Definition: An estimand is the precise target of estimation, described by the population, treatment conditions being compared, outcome variable, handling of events that occur after treatment assignment (the intercurrent-event strategy), and population-level summary measure. In this example, the estimand is the treatment-minus-control difference in mean change from baseline in a symptom score at week 12 among all randomized patients, analyzed according to randomized assignment; outcomes after treatment discontinuation are included.

Suppose higher symptom-score change means greater improvement. The treatment effect \(\delta\) is therefore measured in score points, with positive values favoring treatment. A value of \(\delta=2\) points is stipulated in this hypothetical case as the minimum clinically important benefit. The decision target is not merely whether \(\delta\) is above zero, but whether the evidence supports a benefit exceeding 2 points.

Assume a randomized, parallel-group trial with 20 patients per group. For a deliberately transparent illustration, suppose the outcome is normally distributed within each group with known standard deviation \(\sqrt{20}\) points, and observations are independent across patients. Let \(\hat{\delta}\) be the observed treatment-minus-control difference in mean change. Under this model,

$$ \hat{\delta}\mid\delta\sim N(\delta,2), \qquad SE(\hat{\delta})=\sqrt{20\left(\frac{1}{20}+\frac{1}{20}\right)}=\sqrt{2}. $$

The variance 2 is the sampling variance of the estimated difference, not the variance of an individual outcome. The model assumes the stated outcome scale and week-12 measurement are used consistently, the randomized groups are analyzed as assigned, the normal likelihood and known variance are adequate, and the data are complete for this simplified example. A real analysis must address whether these assumptions are credible, including how missing week-12 outcomes are handled. The known-variance assumption is a simplifying device here, not a general recommendation.

Worked Example: State the Target and Check the Sampling Error

Suppose the hypothetical trial reports a mean change of 5.1 points in the treatment group and 0.9 points in the control group. Confirm the treatment-effect estimate and its standard error under the model.

The estimand is the mean difference in week-12 change for all randomized patients under assignment to treatment versus assignment to control, including observed outcomes after treatment discontinuation. The estimate, in score points, is

$$ \hat{\delta}=5.1-0.9=4.2\text{ points}. $$

With 20 patients in each group and individual outcome variance 20 points squared, the estimated difference has sampling variance \(20/20+20/20=2\) points squared. Thus \(SE=\sqrt{2}=1.414\) points. This calculation depends on the independent-group, equal-known-variance model; it would not be justified if the groups were paired or if the assumed variance model were inappropriate.

Prior Rationale and Posterior Evidence

Use a normal prior \(\delta\sim N(0,25)\), where the second parameter is the variance and the prior standard deviation is 5 points. Centering at zero expresses no expected direction of benefit before seeing the trial data. A standard deviation of 5 allows substantial effects in either direction: approximately 95% of the prior distribution lies between \(-9.8\) and \(9.8\) points. This is not a “noninformative” prior; it is a stated distribution whose scale should be judged against clinical knowledge and the outcome scale. In a regulatory package, the rationale should explain its source, clinical plausibility, and whether it uses external data. Here it is an illustrative, non-borrowing prior.

For the normal likelihood with known variance and a normal prior, the posterior is normal. Its precision is the sum of prior precision and data precision. With prior mean zero, the posterior variance and mean are as follows.

$$ \operatorname{Var}(\delta\mid\hat{\delta}) =\left(\frac{1}{25}+\frac{1}{2}\right)^{-1} =\frac{50}{27}, \qquad E(\delta\mid\hat{\delta}) =\frac{50}{27}\left(\frac{\hat{\delta}}{2}\right) =\frac{25}{27}\hat{\delta}. $$

The posterior standard deviation is \(\sqrt{50/27}=1.36083\) points. Substituting the observed estimate of 4.2 points gives posterior mean \( (25/27)(4.2)=3.88889\) points. The posterior distribution is therefore \(N(3.88889,50/27)\), with the variance—not the standard deviation—as its second parameter.

Worked Example: Calculate the Probability of Clinically Important Benefit

Calculate the posterior probability that the treatment effect exceeds the 2-point clinical threshold, and apply the prespecified rule that this probability must be at least 0.90.

For a normal posterior, standardize the threshold using the posterior mean and standard deviation:

$$ P(\delta>2\mid\hat{\delta}=4.2) = 1-\Phi\left(\frac{2-3.88889}{1.36083}\right) = \Phi\left(\frac{3.88889-2}{1.36083}\right) = \Phi(1.388) \approx 0.9175. $$

The posterior probability is about 0.9175, or 91.75%. Since \(0.9175\geq0.90\), the stated Bayesian decision rule is met. In context, conditional on this model, prior, and observed data, the posterior probability that the mean improvement benefit exceeds 2 points is about 91.75%. It is not the probability that the decision rule is correct, nor a probability statement about repeated trials. A complete report also gives the posterior location and spread in clinical units, as above, rather than only the threshold probability.

Justify the Decision Cutoff

The posterior probability rule can be expressed as a cutoff on the observed estimate. This is useful for checking the decision and for evaluating its repeated-trial behavior. Because the posterior is normal, its probability of exceeding 2 is at least 0.90 exactly when its standardized distance above 2 is at least the 90th percentile of the standard normal distribution.

$$ P(\delta>2\mid\hat{\delta})\geq0.90 \iff \frac{(25/27)\hat{\delta}-2}{\sqrt{50/27}} \geq z_{0.90} \iff \hat{\delta}\geq \frac{2+z_{0.90}\sqrt{50/27}}{25/27}. $$

Here \(z_{0.90}=\Phi^{-1}(0.90)\approx1.2815516\). Substituting without rounding intermediate quantities gives

$$ \hat{\delta}\geq \frac{2+\operatorname{norm.ppf}(0.9)\sqrt{50/27}}{25/27} \approx 4.04349\text{ points}. $$

The inequality is inclusive: a value exactly at the cutoff gives posterior probability 0.90 and meets a rule requiring “at least 0.90.” The observed estimate of 4.2 points exceeds the cutoff, consistent with the direct posterior-probability calculation. The derivation relies on this specific prior, likelihood, clinical threshold, and probability criterion; a different choice changes the cutoff.

Operating Characteristics Are a Separate Result

Operating characteristics describe the probability that a decision rule leads to a specified action over hypothetical repetitions of the trial, conditional on a fixed true treatment effect. They are not posterior probabilities. Under the assumed sampling model, \(\hat{\delta}\mid\delta\sim N(\delta,2)\), and the decision is made when \(\hat{\delta}\geq c\), where \(c=4.04349\). Therefore the probability of meeting the rule at a fixed true effect is

$$ P(\text{rule met}\mid\delta) = P(\hat{\delta}\geq c\mid\delta) = 1-\Phi\left(\frac{c-\delta}{\sqrt{2}}\right) = \Phi\left(\frac{\delta-c}{\sqrt{2}}\right). $$

This calculation assumes the sampling model and decision rule are fixed and that the true effect is the stated value. It is a conditional repeated-sampling probability, not an estimate of the probability that the true effect takes that value. For example, at the clinical boundary \(\delta=2\), it describes the chance of meeting the rule when the true benefit is exactly 2 points. If “clinically important benefit” means \(\delta>2\), then the rule’s probability of a positive decision at \(\delta=2\) is also a boundary false-positive rate for the composite null \(\delta\leq2\); monotonicity makes it the largest such probability within that null.

Worked Example: Evaluate Repeated-Trial Performance

Calculate the probability of meeting the rule if the true treatment effects are 2 points and 5 points. Interpret both results.

At the boundary \(\delta=2\), using \(c=4.04349\) and \(SE=\sqrt{2}\),

$$ P(\text{rule met}\mid\delta=2) = \Phi\left(\frac{2-4.04349}{1.41421}\right) = \Phi(-1.445) \approx0.0743. $$

Thus, if the true effect were exactly the 2-point threshold, about 7.43% of repeated trials under this model would meet the rule. At a true effect of 5 points,

$$ P(\text{rule met}\mid\delta=5) = \Phi\left(\frac{5-4.04349}{1.41421}\right) = \Phi(0.6763) \approx0.7505. $$

The rule would be met in about 75.05% of repeated trials at a true effect of 5 points. Neither result is the posterior probability for the observed trial. In particular, the observed data yielded a posterior probability of about 91.75%, while the repeated-trial chance of meeting the rule at a fixed true effect of 4.2 points is only \(\Phi((4.2-4.04349)/1.41421)\approx0.5441\). These answer different questions.

Sensitivity, Diagnostics, and Alternatives

A regulatory evidence package should show whether the conclusion depends materially on modeling choices, especially the prior. Repeating the analysis with reasonable alternatives is a direct sensitivity analysis. For an alternative normal prior \(N(m_0,v_0)\), the posterior variance is \((1/v_0+1/2)^{-1}\), and the posterior mean is that variance multiplied by \(m_0/v_0+\hat{\delta}/2\). The likelihood and observed estimate remain fixed in this comparison.

Worked Example: Compare Two Alternative Priors

Using the same observed estimate and likelihood, compare a more concentrated prior \(N(0,4)\) with a wider prior \(N(0,100)\). Determine whether each meets the 0.90 rule for \(\delta>2\).

For \(N(0,4)\), posterior variance is \((1/4+1/2)^{-1}=4/3\), posterior mean is \((4/3)(4.2/2)=2.8\), and posterior standard deviation is 1.15470. Therefore,

$$ P(\delta>2\mid\hat{\delta}=4.2,N(0,4)) = \Phi\left(\frac{2.8-2}{1.15470}\right) = \Phi(0.6928) \approx0.7558. $$

This more concentrated prior does not meet the rule. For \(N(0,100)\), posterior variance is \((1/100+1/2)^{-1}=100/51\), posterior mean is \((100/51)(4.2/2)=4.11765\), and posterior standard deviation is 1.40028. Thus,

$$ P(\delta>2\mid\hat{\delta}=4.2,N(0,100)) = \Phi\left(\frac{4.11765-2}{1.40028}\right) = \Phi(1.512) \approx0.9348. $$

The wider prior meets the rule; the more concentrated prior does not. This difference is substantive evidence that the decision depends on prior scale. It should be reported and clinically interpreted, not hidden by presenting only the primary analysis.

Diagnostics should examine whether the model represents the data and estimand adequately. For a simple normal model, inspect group outcome distributions, influential observations, variance assumptions, and missingness patterns; verify that the analysis includes outcomes in a way consistent with the treatment-policy estimand. For more complex models, posterior predictive checks can identify systematic discrepancies between replicated data and observed data. A check that reveals poor fit is a reason to investigate the model, not a mechanical instruction to accept or reject the trial result.

The main alternative summaries answer related but distinct questions. A frequentist confidence interval describes the long-run coverage behavior of an interval procedure under its assumptions; it is not a posterior probability statement. A Bayesian credible interval gives posterior probability for an interval under the specified model and prior, but alone does not encode the particular clinical threshold or regulatory action. A posterior threshold rule makes that action explicit, while operating-characteristic calculations show its repeated-trial behavior. These approaches can be reported together, provided their interpretations are not conflated.

Important limitations remain. The example treats the sampling variance as known, uses a single continuous outcome, and assumes complete data and a correct normal model. It does not address multiplicity, interim adaptation, missing outcomes, treatment-effect heterogeneity, or joint benefit-risk evidence. Sensitivity to prior scale does not prove that any one prior is correct. The decision threshold and posterior-probability requirement should be clinically justified and specified in advance; operating characteristics should be evaluated for a range of plausible true effects, not selected after seeing the result. A clear package states the estimand, model, prior, posterior summary, decision rule, operating characteristics, and sensitivity results together so readers can see both the evidence and the conditions on which it rests.

Key takeaway: A posterior probability is conditional evidence about the estimand under a specified model and prior; an operating characteristic is a repeated-trial probability under a fixed true effect. A regulatory summary should report both when relevant, explain the prior and decision threshold, and show whether reasonable sensitivity analyses change the conclusion.

Check Your Understanding

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

  1. In this case, what is the estimand, and why does including outcomes after treatment discontinuation matter to its interpretation?
  2. Why is the posterior probability rule written with \(\geq\) rather than \(>\), and what is the cutoff on \(\hat{\delta}\) to three decimal places?
  3. At a true effect of 2 points, what does the operating characteristic of approximately 0.0743 mean?
  4. Does the \(N(0,4)\) sensitivity analysis meet the decision rule? Give the posterior probability and explain what this result says about robustness.
  5. Explain why the observed-data posterior probability of about 0.9175 is not the same as the repeated-trial probability of meeting the rule at a true effect of 4.2 points.

Solutions. (1) It is the treatment-minus-control difference in mean week-12 change for all randomized patients under assignment to treatment or control, including outcomes after discontinuation. Including those outcomes specifies the treatment-policy strategy and makes the target an effect of assignment that incorporates what happens after discontinuation, rather than an effect restricted to continued treatment. (2) “At least 0.90” includes equality, so the condition is \(\geq\). The cutoff is \((2+z_{0.90}\sqrt{50/27})/(25/27)\approx4.04349\), or 4.043 points to three decimals; at the unrounded cutoff the posterior probability equals 0.90. (3) Under the assumed sampling model, if the true effect is exactly 2 points, about 7.43% of repeated trials would produce an estimate at or above the decision cutoff. It is not the probability that the true effect is 2 points. (4) No. Its posterior probability is approximately 0.7558, below 0.90. The difference from the primary analysis shows that this decision is sensitive to a more concentrated prior centered at zero. (5) The posterior probability conditions on the observed data and model to quantify uncertainty about \(\delta\). The repeated-trial probability conditions on a fixed true \(\delta=4.2\) and quantifies how often new trial data would cross the rule. They are probabilities over different sources of uncertainty and answer different questions.

Clinical Biostats