Tutorials › Bayesian Clinical Trials › Worked Elicitation for a Binary Trial Endpoint

Choosing and Eliciting Priors · Tutorial 67 of 696

Worked Elicitation for a Binary Trial Endpoint

Learn to translate an expert’s expected event rate and uncertainty into a beta prior, check whether the judgments are coherent, and interpret the resulting parameters.

Intermediate 11 min read

What You'll Learn

  • State the binary endpoint model and the assumptions behind a beta prior for its event probability.
  • Convert an elicited prior mean and standard deviation into beta shape parameters.
  • Check whether an expert’s uncertainty is mathematically compatible with the elicited mean.
  • Interpret the beta concentration without treating it as literal observed data.
  • Recognize when a beta prior may poorly represent the expert’s beliefs.

From an Expert Judgment to a Prior

For a binary clinical endpoint, each participant either experiences the event or does not. Before trial data are observed, an expert may be able to estimate the event probability for a clearly specified population, treatment condition, and follow-up period, and to describe how uncertain that estimate is. The task in this tutorial is to turn those judgments into a prior distribution for the event probability.

The endpoint must be defined before asking for a number. “The event rate” is not sufficiently precise if it could mean events by 30 days or by one year, or if the population or treatment condition is unclear. The elicited probability must refer to the same binary outcome and time horizon that the trial will analyze. The setup below also assumes that participants have the same event probability, conditional on the parameter, and that their binary outcomes are independent.

Definition: Let \(Y_i=1\) if participant \(i\) experiences the specified event by the specified follow-up time, and \(Y_i=0\) otherwise. In a single-arm model, conditional on the event probability \(p\), assume \(Y_i\mid p\sim\operatorname{Bernoulli}(p)\), independently for \(i=1,\ldots,n\). Equivalently, the event count \(X=\sum_{i=1}^n Y_i\) satisfies \(X\mid p\sim\operatorname{Binomial}(n,p)\). A beta prior, written \(p\sim\operatorname{Beta}(\alpha,\beta)\), has density proportional to \(p^{\alpha-1}(1-p)^{\beta-1}\) for \(0<p<1\), where \(\alpha>0\) and \(\beta>0\).

The parameter \(p\) is unknown but fixed for the modeled population and conditions; the prior describes uncertainty about its value. This is distinct from the Bernoulli variation in individual outcomes conditional on \(p\). It is also distinct from variation in event probabilities between participants. If meaningful patient-level heterogeneity is part of the trial model, a single common \(p\) may not be adequate.

The beta distribution is useful because it is defined on probabilities and can represent a range of shapes, from concentrated and approximately symmetric to strongly skewed. Its two parameters are not themselves the expert’s expected probability and uncertainty. Those are more naturally elicited on the probability scale, then translated into \(\alpha\) and \(\beta\).

Moment Matching for a Beta Prior

Suppose the expert’s judgments are summarized by a prior mean \(m\) and prior standard deviation \(s_p\), both on the probability scale. The notation \(s_p\) here means the standard deviation of the prior distribution for \(p\), not a sample standard deviation. Let \(v=s_p^2\) be the elicited prior variance. For \(p\sim\operatorname{Beta}(\alpha,\beta)\), the mean and variance are

$$ E(p)=\frac{\alpha}{\alpha+\beta}, \qquad \operatorname{Var}(p)=\frac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}. $$

These moment formulas can be used in reverse. Define the total concentration \(\kappa=\alpha+\beta\). Matching the beta mean to \(m\) gives \(\alpha=m\kappa\) and \(\beta=(1-m)\kappa\). Substituting those expressions into the variance formula gives \(v=m(1-m)/(\kappa+1)\), and hence the elicitation formulas below.

$$ \kappa=\frac{m(1-m)}{v}-1, \qquad \alpha=m\kappa, \qquad \beta=(1-m)\kappa. $$
Key result: To match an expert’s mean \(m\) and variance \(v\) with a beta prior, calculate \(\kappa=m(1-m)/v-1\), then set \(\alpha=m\kappa\) and \(\beta=(1-m)\kappa\). The required conditions are \(0<m<1\) and \(0<v<m(1-m)\), which ensure \(\alpha>0\) and \(\beta>0\).

The variance condition is a coherence check, not an optional technical detail. A probability-valued parameter with mean \(m\) cannot have variance greater than \(m(1-m)\). If the elicited variance equals that upper bound, the moment-matching formula gives \(\kappa=0\), which is not a proper beta distribution. A variance greater than the bound gives negative concentration. Either result means the stated mean and uncertainty cannot define a beta prior; it does not mean that a negative beta parameter should be used.

The concentration \(\kappa\) describes how tightly the prior is distributed around its mean: holding \(m\) fixed, larger \(\kappa\) means smaller prior variance. Although \(\alpha\) and \(\beta\) are sometimes described as prior event and non-event counts, they are not literal observed participants. That shorthand must not replace the elicitation or the model assumptions.

Worked Example: Eliciting a Prior for a 20% Event Probability

Suppose clinicians estimate that the event probability by 90 days in a specified hypothetical trial population is \(m=0.20\). They describe their uncertainty as a prior standard deviation of 0.08, or 8 percentage points. Assume the common-probability, independent Bernoulli model stated above.

First square the elicited standard deviation to obtain the variance: \(v=0.08^2=0.0064\). The feasibility bound is \(m(1-m)=0.20(0.80)=0.16\), so \(0<v<0.16\) and the judgments are compatible with a proper beta prior.

The concentration is

$$ \kappa=\frac{0.20(1-0.20)}{0.0064}-1 =\frac{0.16}{0.0064}-1 =25-1=24. $$

Thus \(\alpha=0.20(24)=4.8\) and \(\beta=0.80(24)=19.2\), giving \(p\sim\operatorname{Beta}(4.8,19.2)\). Its mean is \(4.8/(4.8+19.2)=0.20\), and its variance is \(4.8(19.2)/(24^2\cdot25)=92.16/14400=0.0064\), so its standard deviation is 0.08. In context, this prior centers uncertainty about the 90-day event probability at 20%, with prior standard deviation 8 percentage points.

Check the Elicitation Before Accepting the Parameters

A calculation that returns positive parameters is necessary but not sufficient. The elicitation should be reviewed for clarity, and the resulting distribution should be checked against the judgments that produced it. In particular, ask the expert to distinguish a central estimate from a likely range, and to state uncertainty about the underlying event probability rather than uncertainty about the outcome of one participant.

The beta family can be asymmetric, especially for probabilities near zero or one. Thus, matching a mean and standard deviation does not guarantee that the prior’s median, mode, or tails match an expert’s intended beliefs. If the expert gives quantiles or a range instead of a standard deviation, those statements should be fitted or checked directly; they should not be silently treated as a mean plus or minus a fixed number of standard deviations. A numerical fit can match selected quantiles, but it may not match the mean and variance at the same time.

Worked Example: A Low Event Probability With Asymmetric Uncertainty

In another hypothetical setting, an expert gives a mean event probability of \(m=0.05\) and a prior standard deviation of 0.03. The uncertainty is 3 percentage points, not 3 percent of the mean. Assume the endpoint, population, and follow-up time are fixed, and use the same common-\(p\) Bernoulli model.

The variance is \(v=0.03^2=0.0009\). The maximum feasible variance for this mean is \(m(1-m)=0.05(0.95)=0.0475\), so \(0.0009<0.0475\). Then

$$ \kappa=\frac{0.0475}{0.0009}-1 =52.7778-1 =51.7778, $$

and therefore

$$ \alpha=0.05(51.7778)=2.5889, \qquad \beta=0.95(51.7778)=49.1889. $$

The resulting prior is \(\operatorname{Beta}(2.5889,49.1889)\). Its mean is 0.05, and its variance is \(0.0475/(51.7778+1)=0.0009\), giving standard deviation 0.03. The beta density is not required to be symmetric around its mean; this is useful when uncertainty concerns a probability near a boundary. The calculation matches the two elicited moments, but the expert should still review whether the prior’s shape, especially its lower and upper tails, reflects the intended beliefs.

Worked Example: Expressing Uncertainty as a Concentration

An expert’s mean event probability for a hypothetical endpoint is \(m=0.30\). After discussing uncertainty, the elicitation team considers a beta prior with concentration \(\kappa=12\) to be an appropriate summary. This is a way of specifying uncertainty through concentration rather than asking directly for a standard deviation. Assume the same fixed endpoint and follow-up conditions.

The parameters are \(\alpha=m\kappa=0.30(12)=3.6\) and \(\beta=(1-m)\kappa=0.70(12)=8.4\). Thus \(p\sim\operatorname{Beta}(3.6,8.4)\). Its variance and standard deviation are

$$ \operatorname{Var}(p)=\frac{0.30(0.70)}{12+1} =\frac{0.21}{13} =0.01615, \qquad s_p=\sqrt{0.01615}=0.1271. $$

The implied prior standard deviation is about 12.7 percentage points. This calculation makes the uncertainty implied by the chosen concentration explicit, so it can be discussed with the expert rather than treated as an arbitrary default. If that spread is inconsistent with the expert’s expressed uncertainty, revise the concentration or elicit the uncertainty directly.

Failure Modes and Alternatives

Several problems can be detected before the prior is used. First, check that the event definition, population, and follow-up time are unambiguous. A probability for a different time horizon is not a harmless change of units. Second, calculate the feasibility bound \(m(1-m)\) before interpreting the parameters. Third, verify the fitted moments by substituting the calculated \(\alpha\) and \(\beta\) back into the beta mean and variance formulas. Finally, examine the shape for compatibility with what the expert said, not just whether the parameters are positive.

A common mistake is to set \(\alpha=m\) and \(\beta=1-m\), reasoning that these parameters must encode the desired mean. That does give the correct mean, because their ratio is \(m\), but it fixes \(\alpha+\beta=1\) and therefore imposes a particular, generally very diffuse prior. For example, if \(m=0.20\), this choice gives \(\operatorname{Beta}(0.20,0.80)\), not a prior whose uncertainty has been elicited. The correction is to obtain an uncertainty statement as well and use both moments to determine the concentration and the two shape parameters.

The main alternatives serve different purposes. A uniform prior is a specific choice, \(\operatorname{Beta}(1,1)\), not a neutral substitute for eliciting uncertainty. A fixed probability omits uncertainty about \(p\). A beta prior is convenient for a probability parameter, but a single beta distribution may be unsuitable if expert beliefs are distinctly multimodal or cannot be represented adequately by two moments. Direct elicitation of quantiles, followed by fitting and checking an appropriate distribution, can be preferable in those cases. The choice should be documented with the endpoint definition, elicited judgments, parameterization, and checks, rather than justified only by mathematical convenience.

Key takeaway: Elicit a well-defined event probability and its uncertainty, then match both moments and check feasibility. The resulting beta prior is a model for uncertainty about the event probability, not a record of participants already observed.

Check Your Understanding

Try each calculation, then compare with the solutions below.

  1. An expert gives a mean event probability of 0.25 and a prior standard deviation of 0.10. Calculate the beta parameters.
  2. For a mean event probability of 0.10, an expert states a prior standard deviation of 0.35. Can these judgments define a proper beta prior by moment matching? Explain.
  3. For a beta prior with \(\alpha=6\) and \(\beta=14\), calculate its mean and standard deviation.
  4. Why is \(\operatorname{Beta}(m,1-m)\) not, by itself, a complete elicitation from a stated mean \(m\)?

Solutions. (1) \(v=0.10^2=0.01\), and \(m(1-m)=0.25(0.75)=0.1875\). Thus \(\kappa=0.1875/0.01-1=17.75\), \(\alpha=0.25(17.75)=4.4375\), and \(\beta=0.75(17.75)=13.3125\). (2) The variance is \(0.35^2=0.1225\), but \(m(1-m)=0.10(0.90)=0.09\). Since \(0.1225>0.09\), the stated variance is infeasible for a probability with that mean, so the judgments must be clarified or revised. (3) The mean is \(6/(6+14)=0.30\). The variance is \(6(14)/(20^2\cdot21)=84/8400=0.01\), so the standard deviation is 0.10. (4) It fixes the concentration at \(\alpha+\beta=1\), and therefore fixes the variance as well as the mean. It cannot represent an independently elicited level of uncertainty unless that imposed variance happens to be intended.

Clinical Biostats