Tutorials › Bayesian Clinical Trials › Standardized Effects And Plausible Trial Differences

Prior Predictive Checks · Tutorial 81 of 696

Standardized Effects And Plausible Trial Differences

Translate standardized treatment effects into endpoint units, then use prior predictive calculations to judge whether the implied trial differences are credible.

Intermediate 11 min read

What You'll Learn

  • Define a standardized treatment effect for a continuous endpoint and relate it to the raw-scale mean difference.
  • Distinguish uncertainty in the treatment-effect parameter from sampling variability in an observed trial difference.
  • Calculate raw-scale implications when the endpoint standard deviation is fixed or uncertain.
  • Simulate prior predictive differences while preserving the model’s assumptions and units.
  • Diagnose implausible prior implications and compare standardized-effect priors with raw-scale alternatives.

Why Standardized Effects Need a Clinical Scale

A prior on a standardized treatment effect can be easier to elicit or compare across settings than a prior stated in endpoint units. But standardization does not make the clinical scale irrelevant. A standardized effect of \(-0.3\) corresponds to different mean differences when the endpoint standard deviation is 5 units versus 20 units. Prior predictive checking therefore asks a concrete question: what raw-scale treatment differences does this prior imply for this endpoint?

This tutorial considers a continuous endpoint measured at a specified follow-up time and a two-arm trial. The earlier material on predictive means under Normal endpoint models covers the basic construction of predictive arm means. Here the focus is translating standardized effects into raw units and judging the resulting differences against the endpoint’s measurement scale and clinical context.

Definition: Let \(\mu_C\) and \(\mu_T\) be the control- and treatment-arm means for a continuous endpoint, and let \(\sigma\) be their common within-arm standard deviation. The raw-scale treatment difference is \(\Delta=\mu_T-\mu_C\). The standardized treatment effect is \(\delta=\Delta/\sigma\), so \(\Delta=\delta\sigma\). The units of \(\delta\) are standard deviations; the units of \(\Delta\) are the endpoint’s original units.

The sign convention here is treatment minus control. Thus, if lower endpoint values are better, negative values of \(\Delta\) and \(\delta\) favor treatment; if higher values are better, positive values favor treatment. State the direction explicitly in an application. Standardization changes the scale, not the clinical meaning of the endpoint.

Model and Prior Predictive Translation

For the basic model, suppose that each participant contributes one endpoint measurement, participants are independent within and across randomized arms, and the endpoint has a Normal distribution with a common standard deviation:

$$ Y_{C,i}\mid\mu_C,\sigma\sim N(\mu_C,\sigma^2),\qquad Y_{T,i}\mid\mu_T,\sigma\sim N(\mu_T,\sigma^2),\qquad \mu_T=\mu_C+\delta\sigma. $$

The \(N(\mu,\sigma^2)\) notation uses the second argument as the variance. The model assumes the same endpoint definition, measurement units, and follow-up time in both arms; a common within-arm spread; and no unmodeled clustering or dependence among participants. The prior must describe plausible values of the parameters in this model. In particular, a prior on \(\delta\) does not by itself determine \(\Delta\): the scale \(\sigma\), fixed or drawn from a prior, is also needed.

In a prior predictive simulation, draw the parameters from their joint prior and calculate \(\Delta=\delta\sigma\) for each draw. If trial data are also simulated, draw participant outcomes under the model and calculate the observed arm-mean difference. These are related but distinct quantities: \(\Delta\) is the parameter contrast, while the observed difference varies around it because of sampling.

Key result: If \(\sigma\) is fixed and \(\delta\sim N(m_\delta,s_\delta^2)\), then \(\Delta=\sigma\delta\sim N(\sigma m_\delta,\sigma^2s_\delta^2)\). If \(\sigma\) is uncertain, the distribution of \(\Delta\) is generally a mixture induced by the joint prior on \((\delta,\sigma)\); it need not be Normal.

The fixed-scale result follows directly from the rules for a linear transformation of a Normal random variable. With uncertain scale, the product \(\delta\sigma\) must be evaluated under the joint prior. If \(\delta\) and \(\sigma\) are independent, their independent draws can be paired in the simulation. If they are dependent, the simulation must preserve that dependence rather than sample from separate marginal priors.

Worked Example: Translating a Prior for a Blood-Pressure Endpoint

Suppose a hypothetical trial measures systolic blood pressure in mmHg at a fixed follow-up. Lower values are better. For this illustration, take the common standard deviation as fixed at \(\sigma=12\) mmHg and assign \(\delta\sim N(-0.25,0.20^2)\). Translate the prior into raw units.

Using \(\Delta=\delta\sigma\), its prior mean is \(12(-0.25)=-3.00\) mmHg and its prior standard deviation is \(12(0.20)=2.40\) mmHg. Thus \(\Delta\sim N(-3.00,2.40^2)\) mmHg. A central 95% interval under this Normal prior is approximately \[ -3.00\pm1.96(2.40)=[-7.704,\ 1.704]\text{ mmHg}. \]

The interval describes prior uncertainty about the treatment-arm mean minus the control-arm mean, not an interval calculated from trial data. It includes negative differences favoring treatment and positive differences favoring control. Whether its range is reasonable must be judged against the blood-pressure endpoint and the intended clinical setting; the standardization alone cannot settle that question.

Why the Observed Trial Difference Is More Variable

For a trial with \(n_C\) and \(n_T\) participants in the two arms, let \(D=\bar{Y}_T-\bar{Y}_C\) be the observed difference in sample means. Under the Normal model and conditional on the parameters, independent observations give \[ D\mid\delta,\sigma\sim N\!\left(\delta\sigma,\ \sigma^2\left(\frac{1}{n_C}+\frac{1}{n_T}\right)\right). \] The conditional mean is the parameter difference; the conditional variance is sampling variance. This distinction matters when a prior predictive check concerns trial-level summaries rather than only the clinical effect parameter.

With fixed \(\sigma\), a Normal prior on \(\delta\) also gives a Normal prior predictive distribution for \(D\). The derivation uses the law of total expectation and the law of total variance:

$$ E(D)=E\{E(D\mid\delta)\}=\sigma m_\delta, $$
$$ \begin{aligned} \operatorname{Var}(D) &=E\{\operatorname{Var}(D\mid\delta)\} +\operatorname{Var}\{E(D\mid\delta)\}\\ &=\sigma^2\left(\frac{1}{n_C}+\frac{1}{n_T}\right) +\sigma^2s_\delta^2. \end{aligned} $$

The first variance term is trial sampling variability; the second is prior uncertainty about the treatment effect. This result assumes the specified Normal model, fixed common \(\sigma\), and a Normal prior for \(\delta\). With an uncertain \(\sigma\), the predictive distribution generally requires integrating over its prior or simulating; the fixed-scale Normal formula should not be applied as if the scale were known.

Worked Example: Separating Effect Uncertainty From Trial Noise

Consider a hypothetical trial of a continuous endpoint where lower values are better, with \(n_C=n_T=100\), fixed \(\sigma=8\) endpoint units, and \(\delta\sim N(-0.20,0.25^2)\). Calculate the prior predictive distribution of \(D\), the observed treatment-minus-control sample-mean difference.

The prior mean of the parameter difference is \(8(-0.20)=-1.60\) units, and its prior standard deviation is \(8(0.25)=2.00\) units. The sampling variance of \(D\) conditional on the parameters is \[ 8^2\left(\frac{1}{100}+\frac{1}{100}\right)=64(0.02)=1.28\text{ units}^2. \] The prior variance of the parameter difference is \(2.00^2=4.00\text{ units}^2\). Therefore, \[ D\sim N(-1.60,\ 1.28+4.00)=N(-1.60,\ 5.28), \] where 5.28 is the variance in squared endpoint units. Its predictive standard deviation is \(\sqrt{5.28}=2.298\) units.

The observed difference has greater spread than the prior distribution of the parameter difference because it includes sampling noise. For example, a simulated trial difference can depart from its simulated parameter value even when the model and prior are correct. Comparing the spread of \(D\) directly with the spread of \(\Delta\) without accounting for this additional term confuses clinical-effect uncertainty with trial precision.

When the Endpoint Scale Is Uncertain

The scale may itself be uncertain, particularly before a trial has collected endpoint data. In that case, simulation makes the translation explicit: draw \(\sigma\) and \(\delta\) jointly from their prior, calculate \(\Delta=\delta\sigma\), and inspect the resulting raw-unit differences. For trial-level checks, then simulate observations under the endpoint model and calculate \(D\). The output should retain the endpoint’s units so it can be judged against relevant clinical knowledge.

Worked Example: Simulating Differences With an Uncertain Scale

Suppose a hypothetical HbA1c endpoint is measured in percentage points, with lower values better. Use an illustrative independent prior in which \(\delta\sim N(-0.40,0.30^2)\), and \(\sigma\) takes values 0.8, 1.2, or 1.6 percentage points with probabilities 0.25, 0.50, and 0.25, respectively. These values define an example prior, not a universal model for HbA1c.

The mean scale is \(0.25(0.8)+0.50(1.2)+0.25(1.6)=1.2\) percentage points. Its variance is \[ 0.25(0.8-1.2)^2+0.50(1.2-1.2)^2+0.25(1.6-1.2)^2=0.08. \] For independent \(\delta\) and \(\sigma\), \(E(\Delta)=E(\delta)E(\sigma)=(-0.40)(1.2)=-0.48\) percentage points. Also \(E(\sigma^2)=0.25(0.64)+0.50(1.44)+0.25(2.56)=1.52\), so \[ \operatorname{Var}(\Delta) =E(\sigma^2)\operatorname{Var}(\delta) +\operatorname{Var}(\sigma)\{E(\delta)\}^2 =1.52(0.09)+0.08(0.16)=0.1496. \] The prior standard deviation of \(\Delta\) is \(\sqrt{0.1496}=0.3868\) percentage points.

Individual simulation draws show how the scale changes the raw effect: \(\delta=-0.40,\sigma=1.2\) gives \(\Delta=-0.48\) percentage points; \(\delta=0.20,\sigma=1.6\) gives \(+0.32\) percentage points; and \(\delta=-1.40,\sigma=1.6\) gives \(-2.24\) percentage points. A clinical team should judge such draws against a justified endpoint-specific range, not a generic rule about what a standardized effect is allowed to be. A draw outside that range is a prompt to review the prior and assumptions; it is not, by itself, evidence that a simulated trial is impossible.

Diagnostics, Alternatives, and Limitations

A useful check is to compare the simulated raw differences with quantities that have meaning on the endpoint’s scale: the instrument’s possible range, known measurement resolution, relevant historical variability, and clinically credible differences for the specified population and follow-up. The basis for any comparison range should be stated. For example, a difference larger than the instrument can record signals a support or scale problem; a difference that conflicts with clinical knowledge may instead signal an implausible prior or an inappropriate model. Check both directions when effects can favor either arm.

Inspect the scale draws as well as the effect draws. A modest standardized effect multiplied by an unexpectedly large \(\sigma\) can yield a large raw difference. Conversely, an apparently substantial standardized effect may imply a small change in original units when \(\sigma\) is small. If an endpoint has natural bounds, a Normal model may also generate impossible individual outcomes even when the mean difference looks reasonable; this tutorial’s contrast calculation does not remove that model concern.

A prior directly on \(\Delta\) is a reasonable alternative when elicitation is naturally expressed in endpoint units or when clinical decisions are tied to a raw-unit threshold. A prior on \(\delta\) can aid communication across scales, but it requires a relevant \(\sigma\) model and can obscure the implications when that scale is uncertain. Neither parameterization is automatically preferable: compare the resulting prior predictive implications in the endpoint’s units. For unequal arm variances, a single common \(\sigma\) no longer defines the model’s standardized effect in the same way; the standardization convention and variance model must be specified rather than silently applying the common-SD formula.

Common mistake: “The prior mean standardized effect is \(-0.2\), so the treatment is expected to lower the endpoint by 0.2 units.” This treats standard deviations as endpoint units. The correction is to multiply by the relevant scale: if \(\sigma=8\) endpoint units, then the corresponding mean raw difference is \((-0.2)(8)=-1.6\) endpoint units. If \(\sigma\) is uncertain, use the joint prior to calculate or simulate the implied raw differences.

The central diagnostic is not whether a standardized effect looks familiar in isolation. It is whether the specified joint prior and endpoint model generate parameter differences, and if needed trial-level differences, that are credible for this particular measurement scale and study setting. This check cannot establish that the prior is correct; it can expose consequences that deserve reconsideration before the data are analyzed.

Check Your Understanding

Answer each question using the stated sign convention and units.

  1. A lower-is-better endpoint has \(\sigma=5\) units and a prior mean standardized effect of \(-0.30\). What is the prior mean raw difference?
  2. With fixed \(\sigma\), why is the prior predictive variance of an observed arm-mean difference larger than the prior variance of the parameter difference?
  3. For \(\delta=-0.5\) and \(\sigma=4\) mmHg, calculate \(\Delta\) and state which arm has the lower mean.
  4. Why should a simulation preserve dependence between \(\delta\) and \(\sigma\) when their prior is dependent?

Solutions. (1) \(E(\Delta)=5(-0.30)=-1.5\) endpoint units. (2) The observed difference includes sampling variance \(\sigma^2(1/n_C+1/n_T)\) in addition to prior variance in the parameter contrast. (3) \(\Delta=(-0.5)(4)=-2\) mmHg, so the treatment-arm mean is lower by 2 mmHg. (4) Independent sampling would change the joint prior and therefore could change the distribution of the product \(\delta\sigma\), giving incorrect raw-scale implications.

Clinical Biostats