What a Continuous-Trial Audit Examines
A prior predictive audit asks what trial data a proposed model and prior could generate before any trial outcomes are observed. For a continuous endpoint, that means looking beyond a prior distribution for a treatment effect: the model must also generate patient values and arm summaries. This makes it possible to detect priors that imply implausible outcome levels, excessive variability, or treatment differences that are too large or too small for the clinical setting.
The audit below uses a parallel-group trial with a continuous endpoint measured at a fixed follow-up time. Suppose lower scores indicate better health. Define the treatment effect as treatment-arm mean minus control-arm mean, so a negative effect favours treatment. The numerical ranges and clinical thresholds are hypothetical and must be adapted to the endpoint being studied.
To make the audit concrete, use \(\mu_C\sim N(50,10^2)\), \(\delta\sim N(0,5^2)\), independently, and initially take \(\sigma=8\) as fixed. These priors are illustrative, not recommendations. The control-arm mean prior represents uncertainty about the typical control outcome; the treatment-effect prior represents uncertainty about the difference in arm means. A complete audit should also assess the prior for \(\sigma\) if it is estimated rather than fixed.
From Parameter Draws to Trial Data
Each prior predictive replicate represents one possible trial under the specified model. Draw the parameters once for the replicate, then use those same values to generate every patient’s outcome in the corresponding arms. Drawing a new treatment effect for each patient would simulate patient-level treatment-effect variation, a different model from the one stated above.
Here \(r\) indexes a replicate, and \(\pi\) denotes a prior distribution. For each replicate, calculate the summaries that would matter in the planned analysis: for example, arm means, arm standard deviations, the mean difference \(\bar{Y}_T-\bar{Y}_C\), and the minimum and maximum patient values. If the design has unequal arm sizes, use those planned sizes in every replicate.
Worked Example: Summarizing One Generated Trial
Consider one hypothetical prior predictive replicate with four patients per arm. Suppose the generated control values are 42, 49, 51, and 58 points, and the generated treatment values are 39, 47, 48, and 52 points. These values are a small illustration of the calculations, not results from a real trial.
The control mean is \(\bar{Y}_C=(42+49+51+58)/4=200/4=50\) points. The treatment mean is \(\bar{Y}_T=(39+47+48+52)/4=186/4=46.5\) points. Thus the generated treatment-minus-control contrast is \(46.5-50=-3.5\) points. With lower scores indicating better health, this replicate has a treatment-arm mean 3.5 points lower than the control-arm mean.
The control sample standard deviation is \(s_C=\sqrt{[(42-50)^2+(49-50)^2+(51-50)^2+(58-50)^2]/(4-1)}=\sqrt{130/3}=6.583\) points. The treatment sample standard deviation is \(s_T=\sqrt{[(39-46.5)^2+(47-46.5)^2+(48-46.5)^2+(52-46.5)^2]/3=\sqrt{89/3}=5.447\) points. The generated range is 42 to 58 points in control and 39 to 52 points in treatment. A single replicate does not establish that the prior is plausible; it shows what summaries to inspect across many replicates.
Audit Treatment Differences Against Clinical Meaning
A useful check is whether prior predictive treatment differences have a reasonable frequency and magnitude relative to a prespecified clinically meaningful difference. The distinction between a parameter effect and an observed arm difference matters: the generated sample contrast varies both because \(\delta\) varies across prior draws and because patient outcomes vary within each trial.
For equal arm sizes \(n\), under the normal model with fixed \(\sigma\), the sample contrast \(D=\bar{Y}_T-\bar{Y}_C\), conditional on \(\delta\), has mean \(\delta\) and variance \(2\sigma^2/n\). With the independent normal prior \(\delta\sim N(0,5^2)\), its prior predictive distribution is normal with variance equal to the prior variance of \(\delta\) plus the sampling variance. This follows by writing \(D=\delta+\varepsilon\), where \(\varepsilon\) is independent of \(\delta\), has mean zero, and has variance \(2\sigma^2/n\).
Worked Example: Probability of a Clinically Large Mean Difference
Suppose the planned trial has 40 patients per arm, and a treatment-minus-control difference of \(-5\) points or less is considered a clinically large benefit. Under the illustrative priors and fixed standard deviation above, the prior predictive standard deviation of \(D\) is \(\sqrt{25+2(64)/40}=\sqrt{28.2}=5.310\) points. Therefore,
Thus the model assigns about a 17.3% prior predictive probability to an observed treatment-minus-control mean difference of \(-5\) points or less in a trial of this size. This is a probability about generated data under the specified model, not the probability that a treatment effect is true or that the null hypothesis is true. For comparison, the parameter-level probability is \(P(\delta\leq-5)=\Phi(-5/5)=\Phi(-1)=0.1587\). The difference between these probabilities reflects sampling variation in \(D\).
In a simulation with 10,000 replicates, one might count the replicates with \(D\leq-5\) and divide by 10,000. For example, 1,734 such replicates give an estimated probability \(1,734/10,000=0.1734\), close to the analytic calculation. The simulation estimate will vary across runs because of Monte Carlo error.
The audit should consider both tails and a range of effect sizes, rather than focusing only on one threshold. A prior that makes clinically important benefit almost impossible may be too concentrated near no effect; a prior that frequently generates implausibly large benefit or harm may be too diffuse or poorly centred. The appropriate judgement depends on clinical knowledge and the decision the trial is intended to support.
Check Outcome Scale and Extreme Values
For continuous endpoints, a prior predictive check must inspect patient-level values as well as arm averages. Arm means can look reasonable even when the model routinely generates individual outcomes outside the scale’s possible or clinically credible range. Conversely, an unusual individual value may be unsurprising under a reasonable heavy-tailed model, so an extreme-value flag should prompt investigation rather than automatic rejection.
Worked Example: A Scale Prior That Generates Extreme Scores
Suppose the endpoint is a score restricted to the interval 0 to 100, and the control mean is 50. Consider a proposed prior \(\sigma\sim\text{Half-Normal}(0,20^2)\). Its 95th percentile is approximately \(1.96(20)=39.2\) points. At the illustrative parameter value \(\sigma=39.2\), a normally distributed patient outcome has probability
So, conditional on this mean and standard deviation, about 20.2% of patient values fall outside the scale. With 40 independent patients in an arm, the probability that at least one value falls outside the scale is \(1-(1-0.202)^{40}=1-(0.798)^{40}\approx0.9999\). This calculation is conditional on \(\sigma=39.2\); it is not the overall prior predictive probability under the half-normal prior. A simulation that draws \(\sigma\) from that prior would estimate the overall frequency directly.
If the endpoint truly cannot be below 0 or above 100, frequent generated values outside that range signal a mismatch between the outcome model and its measurement scale, or a scale prior that places too much mass on large variability. Possible responses include revisiting the prior for \(\sigma\), using a model suited to bounded outcomes, or modelling a suitable transformation. The audit identifies the issue; it does not decide which change is scientifically appropriate.
A Practical Audit Sequence
Record the endpoint units, possible range, clinically credible values, and meaningful treatment differences. State which direction favours treatment.
Specify the outcome distribution, arm structure, sample sizes, prior distributions, and any shared or arm-specific parameters. Check whether independence, common variance, and support assumptions are defensible.
Draw parameters from their priors, then generate all planned patient outcomes conditional on those draws. Preserve the trial’s actual allocation and sample sizes.
Inspect patient-level values and arm summaries, including means, standard deviations, contrasts, and extremes. Compare their distributions with clinical benchmarks and investigate problematic patterns.
Useful diagnostics include histograms or quantiles of generated patient values and arm means, the prior predictive distribution of \(D\), and the frequency of values outside known bounds. For each diagnostic, state what is being counted and its denominator. A probability based on a patient value is not interchangeable with the proportion of trials having at least one such value. Also check whether rare extreme outcomes occur in a small number of replicates or repeatedly across the simulation.
Monte Carlo summaries are estimates, not exact probabilities. If an estimated frequency is important to the conclusion, increase the number of replicates or report its Monte Carlo uncertainty. The earlier material on Monte Carlo error in prior predictive summaries gives the general calculation. The earlier material on sensitivity of prior predictive checks to prior parameters is also relevant when the diagnosis changes substantially under plausible alternative priors.
Alternatives, Limitations, and a Common Pitfall
Simulation is useful because it can carry the entire model through to patient values and complicated trial summaries. For a simple normal model, analytic calculations can independently verify some results, as in the distribution of \(D\) above. Analytic checks are often faster and exact under their assumptions, but become less convenient with more complex models or summaries. Graphical comparisons with clinical benchmarks are valuable for communicating scale and tail behaviour; a numerical discrepancy statistic can condense a prespecified feature into one number, but cannot replace clinical interpretation. These approaches complement rather than substitute for a coherent outcome model.
The normal model is not automatically appropriate for every continuous endpoint. It may assign positive probability to impossible values, fail to represent skewness, or understate tail behaviour. A common standard deviation may also be unsuitable when variability differs by arm. The audit cannot establish that the model is true; it reveals consequences of the specified priors and model assumptions. If generated outcomes are systematically implausible, revise the assumptions or priors for substantive reasons and repeat the check.
Check Your Understanding
Answer each question, then compare your reasoning with the brief solutions.
- In the stated model, why are parameter draws made once per replicate rather than once per patient?
- For equal arm sizes of 25 and fixed \(\sigma=6\), what is the sampling variance of the arm-mean contrast conditional on \(\delta\)?
- Explain why \(P(D\leq-5)\) need not equal \(P(\delta\leq-5)\).
- If 80 out of 2,000 simulated trials have at least one score outside the endpoint’s allowed range, what prior predictive frequency does this estimate?
- Name one limitation of a normal model for a continuous endpoint restricted to 0–100.
Solutions: (1) A parameter draw defines the trial-level means and variability shared by patients in that replicate; drawing a new effect per patient changes the model. (2) The conditional variance is \(2\sigma^2/n=2(36)/25=2.88\) score-squared units. (3) \(D\) includes both the prior variation in \(\delta\) and within-trial sampling variation, so its distribution is wider than the distribution of \(\delta\) alone. (4) The estimated frequency is \(80/2,000=0.04\), or 4% of simulated trials. (5) A normal distribution can generate values below 0 or above 100, which are impossible on that scale.