Tutorials › Bayesian Clinical Trials › Continuous-Outcome Estimands for Means and Quantiles

Bayesian Estimands and Regulatory Interpretation · Tutorial 283 of 696

Continuous-Outcome Estimands for Means and Quantiles

Learn to define and interpret posterior treatment contrasts in means and quantiles, especially when clinical outcomes are skewed.

Advanced 10 min read

What You'll Learn

  • Define marginal mean-difference and quantile-difference estimands for a continuous outcome.
  • Distinguish a population quantile contrast from an effect on particular individuals.
  • Calculate mean and quantile contrasts for right-skewed lognormal and exponential outcomes.
  • Construct posterior contrasts draw by draw and interpret their uncertainty in context.
  • Diagnose when model fit, tail information, or the choice of summary makes an estimand unreliable.

Why the Summary Scale Matters

A continuous clinical outcome may be strongly skewed: most patients might have modest symptom improvement while a smaller group has very large improvement, or most hospital stays may be short while a few are prolonged. In such settings, a treatment contrast in means answers a different question from a contrast in medians or upper quantiles. Both can be useful, but neither is a universal summary of the entire outcome distribution.

The estimand must first specify the outcome, target population, treatment conditions, time point, and strategy for intercurrent events. Those choices determine which outcome distribution is being summarized. Here, the focus is on defining and interpreting mean and quantile contrasts from the resulting treatment-specific marginal distributions. The distributions are marginal over the target population, rather than conditional on a particular covariate profile.

Definition: Let \(Y^a\) denote the continuous outcome under treatment condition \(a\), as defined by the chosen estimand, and let \(F_a(y)=P(Y^a\le y)\) be its distribution in the target population. If the mean exists, \(\mu_a=E(Y^a)\). The mean-difference estimand, using active treatment \(A\) minus control \(C\), is \(\Delta_\mu=\mu_A-\mu_C\). For \(0<p<1\), define the \(p\)-quantile as \(q_{a,p}=\inf\{y:F_a(y)\ge p\}\). The quantile-difference estimand is \(\Delta_{q,p}=q_{A,p}-q_{C,p}\). Both contrasts are expressed in the outcome’s units.

A negative contrast means the active-treatment mean or selected quantile is lower than the corresponding control value. Whether that is desirable depends on the outcome: lower symptom scores might be better, whereas a higher functional score might be better. A quantile contrast compares two locations in two population distributions. It does not, in general, describe the effect for an identifiable group of individuals ranked at the \(p\)-th position.

Conditions for a Meaningful Contrast

The distributions \(F_A\) and \(F_C\) must refer to the same target population and to well-defined treatment conditions. In a randomized trial, randomization supports comparison of the assigned treatment groups; consistency connects observed outcomes to the specified treatment conditions. If the target is a broader population or the analysis adjusts for covariates, the model and standardization must also support that target. The outcome handling after intercurrent events must follow the chosen estimand strategy; changing that strategy changes the outcome distributions and therefore changes the estimand.

For a mean contrast, each treatment-specific mean must be finite. This is more than a technical detail: a distribution with a sufficiently heavy tail may have no finite mean, or may yield a mean so sensitive to a small number of extreme values that it is not a stable clinical summary. Quantiles can still be defined when a mean is not finite, although their interpretation and estimation require care. For a quantile to be uniquely and stably located, it is helpful for the distribution to be continuous with positive density near that quantile. Sparse information in a tail can make an upper quantile especially uncertain.

In Bayesian analysis, a model represents uncertainty about the treatment-specific distributions. For each posterior draw, calculate the two treatment-specific means or quantiles and then subtract control from active treatment. The resulting draws form the posterior distribution of the estimand. The treatment distributions and their dependence within each draw matter: uncertainty in a contrast should be calculated from the paired treatment-specific quantities in that draw, not by treating their marginal posterior summaries as independent.

How Skewness Separates Means and Quantiles

A short calculation shows why the mean and quantile can tell different stories. Suppose \(Y\) is lognormal: \(\log(Y)=Z\), where \(Z\sim N(\mu,\sigma^2)\). Since exponentiation is increasing, the \(p\)-quantile of \(Y\) is the exponential of the \(p\)-quantile of \(Z\). Also, the normal moment-generating function gives the mean of \(Y\).

$$ q_p=\exp\{\mu+\sigma\Phi^{-1}(p)\}, \qquad E(Y)=E\{\exp(Z)\}=\exp(\mu+\sigma^2/2). $$

The \(\sigma^2/2\) term in the mean, absent from the median, shows that a larger spread on the log scale raises the arithmetic mean even when the median is unchanged. Thus, a mean contrast can reflect both location and the distribution’s upper tail. A quantile contrast focuses on a specified point in the distribution, but says nothing by itself about the rest of that distribution.

Worked Example: Mean and Median Differences Under Skewness

Suppose a hypothetical trial models a positive outcome measured in days as lognormal in each treatment group. Let the control distribution have log-scale parameters \(\mu_C=\ln(15)\) and \(\sigma_C=1\), and the active-treatment distribution have \(\mu_A=\ln(10)\) and \(\sigma_A=0.5\). These parameters describe marginal outcome distributions in the same target population. Calculate the mean difference and the median difference.

Means. Applying \(E(Y)=\exp(\mu+\sigma^2/2)\), the active-treatment mean is \(10\exp(0.5^2/2)=10\exp(0.125)\approx11.33\) days. The control mean is \(15\exp(1^2/2)=15\exp(0.5)\approx24.73\) days. Therefore, the mean difference is \(11.33-24.73\approx-13.40\) days. Using the unrounded means gives \(10\exp(0.125)-15\exp(0.5)\approx-13.3993\) days.

Medians. The median is the \(p=0.5\) quantile, and \(\Phi^{-1}(0.5)=0\). The active median is \(\exp(\ln 10)=10\) days and the control median is \(\exp(\ln 15)=15\) days. The median difference is \(10-15=-5\) days.

The mean difference is much more negative than the median difference because the control distribution has substantially greater spread and a longer upper tail. These are two different population summaries, not competing calculations of one uniquely defined treatment effect. The lognormal assumption is an illustrative model assumption here; in an actual analysis, its fit would need to be checked.

Posterior Contrasts: Calculate Within Each Draw

Suppose the posterior distribution is represented by \(M\) draws. In draw \(m\), let \(\mu_A^{(m)}\) and \(\mu_C^{(m)}\) be the model-implied marginal means. The draw for the mean-difference estimand is \(\Delta_\mu^{(m)}=\mu_A^{(m)}-\mu_C^{(m)}\). For a quantile estimand, calculate \(q_{A,p}^{(m)}-q_{C,p}^{(m)}\) in each draw. Posterior probabilities about a contrast are then proportions of draws satisfying the stated condition. For example, the estimated posterior probability that the mean difference is below zero is the number of draws with \(\Delta_\mu^{(m)}<0\), divided by \(M\).

Worked Example: Comparing Posterior Probabilities on Two Scales

In a hypothetical randomized trial, lower symptom-score values are preferable. A fitted Bayesian model gives \(10{,}000\) paired posterior draws of the marginal active and control distributions. In one draw, the means are 102 and 108 points, and the 75th percentiles are 91 and 95 points. Across all draws, 9,400 mean differences are below zero and 7,700 75th-percentile differences are below zero. The posterior mean difference is \(-4.2\) points with a 95% credible interval from \(-8.1\) to \(-0.3\) points; the posterior mean 75th-percentile difference is \(-2.0\) points with a 95% credible interval from \(-6.0\) to \(2.1\) points.

For the displayed draw, the mean contrast is \(102-108=-6\) points and the 75th-percentile contrast is \(91-95=-4\) points. Across all draws, the posterior probability of a negative mean difference is \(9{,}400/10{,}000=0.94\). The posterior probability of a negative 75th-percentile difference is \(7{,}700/10{,}000=0.77\).

The first 95% credible interval lies entirely below zero, while the second includes zero. In context, under this model and the observed data, the posterior probability that the active-treatment population mean symptom score is lower is 0.94. The corresponding probability for the 75th percentile is 0.77. These statements concern posterior uncertainty about specified population contrasts; they are not probabilities that a null hypothesis is true. Nor does the higher probability for the mean establish that every patient benefits.

Quantiles and the Shape of the Distribution

Quantiles make it possible to ask targeted questions about different parts of a skewed distribution. The median describes its center in the sense of the 50th percentile; a 90th-percentile contrast describes the difference at a high outcome value. These contrasts need not be similar. A treatment may have a modest median difference but a larger upper-tail difference, which could matter when unusually high outcomes carry substantial clinical burden.

The quantile level should be selected for a clinical reason, not because it looks favorable after inspecting results. The same applies to the mean: it is often valuable because it incorporates all outcome values and can connect to total burden, but it can be sensitive to a long tail. Neither summary replaces examination of the outcome distributions themselves.

Worked Example: Mean and Upper-Tail Contrasts With Exponential Outcomes

Consider a second hypothetical illustration in which hospital stay, measured in days, follows an exponential distribution in each group. Assume the active-treatment distribution has mean \(\theta_A=14\) days and the control distribution has mean \(\theta_C=20\) days. This model is strongly restrictive: it assumes a particular right-skewed distribution, so the calculations illustrate interpretation rather than establish an appropriate model for a real trial.

For an exponential distribution with mean \(\theta\), the cumulative distribution function is \(F(y)=1-\exp(-y/\theta)\) for \(y\ge0\). Solving \(F(q_p)=p\) gives \(q_p=-\theta\ln(1-p)\). The mean difference is \(14-20=-6\) days.

For the median, \(p=0.5\), so the quantile difference is \(-14\ln(0.5)-[-20\ln(0.5)]=(14-20)\ln 2\approx-4.159\) days. For the 90th percentile, it is \((20-14)\ln(0.1)\approx-13.82\) days. Thus the model implies an active-treatment median stay about 4.16 days shorter and a 90th-percentile stay about 13.82 days shorter. The differences vary by quantile because the model scales the entire distribution as its mean changes; the mean difference alone does not communicate that pattern.

Diagnostics, Alternatives, and Limitations

Before relying on a posterior contrast, check whether the fitted model represents the parts of the outcome distribution that the estimand uses. Posterior predictive checks should compare relevant features of replicated and observed outcomes, including skewness, extreme values, and upper-tail behavior. For a high quantile, inspect whether the data contain enough information near and beyond that part of the distribution. A credible interval can be precise under a restrictive model and still fail to reflect uncertainty about model misspecification.

For outcomes with many extreme values, examine whether a small number of observations dominate the mean. For quantiles, check sensitivity to the assumed distribution and to sparse tail data; a high quantile is often less stable than the median. Also verify that reported contrasts use the intended marginal target population. A conditional model coefficient is not automatically the marginal mean or quantile contrast: if the scientific question is marginal, obtain the treatment-specific marginal distributions by averaging over the target population before calculating the contrasts.

The main alternative is not a single summary but a different estimand. A mean difference is natural when average outcome burden is the clinical target and the mean is well-defined. A quantile difference is useful when a particular distributional location is clinically important or when the mean is unrepresentative of a skewed distribution. Ratios or contrasts on a transformed scale may answer other questions, but they change the scale and interpretation and should be specified as such. Reporting several summaries can describe different features, but selecting whichever one appears most favorable after seeing the results creates a misleading account.

Finally, marginal quantile differences should not be interpreted as effects on the same individuals at a particular rank. The patient at the active-treatment 75th percentile need not be the same patient who would occupy the control 75th percentile. Likewise, a posterior probability that a contrast is negative is conditional on the likelihood, prior, model, and data used; it does not remove uncertainty about whether those choices adequately represent the trial and its target population.

Key takeaway: Define the target population and outcome distribution first, then choose the summary that answers the clinical question. A mean difference summarizes average outcomes; a quantile difference compares a specified location in two distributions. Under skewness, their magnitudes and posterior uncertainties can differ substantially.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. In a hypothetical outcome analysis, the active-treatment mean is 8 units and the control mean is 11 units. State the mean-difference estimand and its interpretation.
  2. Why does a 90th-percentile difference not generally describe the effect on the same individuals in both treatment conditions?
  3. In \(8{,}000\) posterior draws, \(6{,}800\) mean differences are below zero. What posterior probability does this estimate?
  4. For exponential outcomes with active and control means of 12 and 18 days, calculate the mean difference and the median difference.

Solutions. (1) The active-minus-control mean difference is \(8-11=-3\) units: the active-treatment population mean is 3 units lower. (2) The two quantiles are defined separately from two population distributions; the individuals at the same percentile need not be the same people under the two treatment conditions. (3) The estimated posterior probability of a negative mean difference is \(6{,}800/8{,}000=0.85\). (4) The mean difference is \(12-18=-6\) days. For exponential outcomes the median is \(\theta\ln2\), so the median difference is \((12-18)\ln2\approx-4.159\) days.

Clinical Biostats