Why Baseline Risk and Placebo Change Need Separate Terms
In a binary PK/PD analysis, an event probability under active treatment reflects more than the drug’s concentration effect. It also depends on the event risk before treatment and on changes that may occur under placebo, such as the natural course of symptoms, study contact, or expectations. If these components are merged, a treatment contrast can be calculated against the wrong reference and an exposure effect can be misinterpreted as a change from baseline.
The model below separates these components on the logit scale. Its purpose is not to claim that every placebo-related change is psychological; the placebo term represents the modeled change shared by the placebo condition, whatever its causes. The drug exposure term is then an incremental effect relative to that condition.
\(\operatorname{logit}\{P(Y=1\mid C,A)\}=\operatorname{logit}(p_B)+\delta_P A_P+E(C)A_D.\)
Here \(A_P=1\) when the modeled placebo-related change applies and \(0\) otherwise; \(A_D=1\) when the exposure effect applies and \(0\) otherwise; \(\delta_P\) is the placebo-related log-odds change; and \(E(C)\) is the drug effect at concentration \(C\), constrained so that \(E(0)=0\). In a design where active treatment and placebo share the same study-period conditions, both indicators are 1 in the active arm; in the placebo arm, \(A_P=1\) and \(A_D=0\).
The baseline reference must be specified in context: for example, the event probability at study entry for a participant with specified covariates, or a reference-group baseline probability. It is not automatically the observed proportion in either the placebo or active arm. Similarly, \(\delta_P\) is a model parameter for the change from that baseline reference under placebo conditions, not a probability-point change.
From Model Parameters to Treatment Contrasts
Write \(\operatorname{expit}(\eta)\) as the probability corresponding to log odds \(\eta\). Under the model, the baseline, placebo, and active probabilities at concentration \(C\) are
In the first expression, \(p_B\) is the defined baseline probability itself. In the second and third, the same baseline log odds are shifted by the placebo term and, for active treatment, by the drug exposure term. The active-versus-placebo contrast at concentration \(C\) is therefore based on \(p_A(C)\) and \(p_P\), not on \(p_A(C)\) and \(p_B\).
Thus \(E(C)\) is the active-versus-placebo log odds ratio under this model, while \(p_A(C)-p_P\) is the corresponding probability difference. These are distinct scales. The probability difference depends on the reference probability, even when the log odds ratio is held fixed.
Why the Reference Probability Changes the Risk Difference
The probability-scale contrast follows directly from the inverse-logit transformation. Let the placebo log odds be \(\eta_P=\operatorname{logit}(p_B)+\delta_P\), and let \(r=\exp\{E(C)\}\) be the active-versus-placebo odds ratio. Placebo odds are \(o_P=\exp(\eta_P)\), and active odds are \(r o_P\). Converting each set of odds to a probability gives
The last equality follows by substituting \(o_P=p_P/(1-p_P)\) and multiplying numerator and denominator by \(1-p_P\). Consequently, a fixed odds ratio does not imply a fixed risk difference: the transformation also depends on \(p_P\). This is why a treatment effect should be reported against a stated reference probability and, when useful, on both the odds and probability scales.
Worked Example: Baseline Risk, Placebo Change, and a Protective Drug Effect
Suppose a hypothetical binary endpoint has baseline probability \(p_B=0.20\). Under placebo, the odds are assumed to increase by a factor of 1.50, so \(\delta_P=\log(1.50)\). At a specified exposure, suppose \(E(C)=-0.80\). Calculate the placebo and active probabilities and the active-versus-placebo risk difference.
Baseline odds are \(0.20/(1-0.20)=0.25\). The placebo odds are \(0.25(1.50)=0.375\), giving
The active-to-placebo odds ratio is \(\exp(-0.80)=0.4493\). Active odds are therefore \(0.375(0.4493)=0.1685\), and
The active-versus-placebo risk difference is \(0.1442-0.2727=-0.1285\), or about 12.85 percentage points fewer events under active treatment at this exposure. The active-versus-baseline difference is \(0.1442-0.20=-0.0558\), which answers a different question because it includes the placebo-related increase as well as the drug effect.
Placebo Change Is Not the Drug Effect
The placebo parameter can move risk up or down. Its interpretation depends on the coding and reference period: \(\delta_P>0\) means higher event odds under the modeled placebo condition than at the defined baseline; \(\delta_P<0\) means lower odds. The drug effect \(E(C)\) describes the additional log-odds shift from placebo at concentration \(C\). A negative \(E(C)\) is protective only when \(Y=1\) is an undesirable event; the sign alone does not establish clinical benefit.
For a nonlinear exposure effect, one possible specification is an Emax function on the logit scale. For example, with a signed maximum effect \(E_{\max}\), half-maximal concentration \(EC_{50}>0\), and \(C\geq0\),
This function has \(E(0)=0\), approaches \(E_{\max}\) as concentration increases, and gives \(E(EC_{50})=E_{\max}/2\). The baseline and placebo terms remain separate from this concentration-response component. The specific exposure function should be chosen and justified for the application; the contrast calculations do not require this particular form.
Worked Example: A Placebo Increase and an Incremental Drug Increase
Consider a second hypothetical endpoint with \(p_B=0.40\), placebo log-odds change \(\delta_P=0.30\), and drug effect \(E(C)=0.60\) at the concentration of interest. Calculate both probabilities and the active-versus-placebo contrast.
Baseline log odds are \(\operatorname{logit}(0.40)=\log(0.40/0.60)=-0.4055\). Placebo log odds are \(-0.4055+0.30=-0.1055\), so
Active log odds are \(-0.4055+0.30+0.60=0.4945\), giving \(p_A(C)=\operatorname{expit}(0.4945)=0.6212\). The drug odds ratio versus placebo is \(\exp(0.60)=1.822\), and the risk difference is \(0.6212-0.4737=0.1475\), or 14.75 percentage points more events under active treatment. The increase from baseline to active treatment is larger, \(0.6212-0.40=0.2212\), because it includes both the placebo-related shift and the drug effect.
The arithmetic does not imply that the drug caused every active-versus-baseline change. The active-versus-placebo contrast isolates the modeled incremental exposure effect, provided the design and model assumptions support that comparison.
Worked Example: The Same Odds Ratio at Two Baseline Risks
Suppose, hypothetically, that placebo change is zero and the drug effect is \(E(C)=-0.60\), so the active-versus-placebo odds ratio is \(\exp(-0.60)=0.5488\). Compare the predicted risk differences for two baseline risks, \(p_B=0.10\) and \(p_B=0.50\).
For \(p_B=0.10\), baseline odds are \(0.10/0.90=0.1111\). Active odds are \(0.1111(0.5488)=0.06098\), so \(p_A=0.06098/(1+0.06098)=0.0575\). The risk difference is \(0.0575-0.10=-0.0425\), about 4.25 percentage points fewer events.
For \(p_B=0.50\), baseline odds are 1. Active odds are \(1(0.5488)=0.5488\), so \(p_A=0.5488/(1+0.5488)=0.3543\). The risk difference is \(0.3543-0.50=-0.1457\), about 14.57 percentage points fewer events.
The odds ratio is the same in both calculations, but the probability differences are not. The second risk difference is larger in magnitude because the reference event risk differs. This matters when translating a model coefficient into an expected number of events prevented or increased.
Assumptions, Diagnostics, and Alternatives
The interpretation above relies on several assumptions. First, the binary outcome is correctly represented as \(Y=1\) or \(0\), with the event definition and observation window held fixed across comparisons. Second, the baseline probability refers to a clearly specified population and covariate profile. Third, the logit-scale terms have the stated meanings: the placebo shift is shared by the relevant treatment conditions, and the drug term is incremental to placebo and equals zero at zero exposure. Finally, the chosen exposure-response function and any covariate structure must adequately describe the conditional event probabilities.
The Bernoulli logistic formulation treats observations as conditionally independent given the predictors in the model. If outcomes are repeated within individuals or otherwise clustered, conditional independence may fail; a model that accounts for the dependence may be needed. Exposure must also be appropriately measured and aligned with the outcome interval. If the drug concentration is not known without error, treating an uncertain value as exact can distort exposure-effect estimates.
Useful warning signs include predicted probabilities outside the plausible range for the clinical setting, a placebo shift inconsistent with the study’s baseline-to-placebo pattern, or a concentration-response shape that does not track the observed response pattern. Sparse observations at some concentrations or in some treatment conditions can make estimates unstable. Such patterns prompt investigation of data, design, and model specification; they do not by themselves identify the cause of a problem. A full assessment should also examine uncertainty in the predicted probabilities and contrasts.
A categorical-arm logistic model can estimate separate log odds for baseline, placebo, and active groups without imposing a particular decomposition. That can be useful when exposure is not central or when a compact descriptive comparison is the aim, but it does not by itself provide a continuous concentration-response relationship. A risk-difference model targets probability changes more directly, but must still ensure fitted probabilities are valid and specify how effects vary with baseline risk. A probit model offers a different link-scale parameterization; its coefficients are not log odds ratios. These approaches answer related but not identical modeling questions.
Common Mistakes and Practical Interpretation
- Comparing active treatment with baseline when asked for the drug effect. The active-versus-baseline difference includes placebo change. For an incremental treatment effect, compare active treatment with the corresponding placebo probability.
- Calling a log-odds coefficient a probability change. \(E(C)\) is a log odds ratio; exponentiating it gives an odds ratio. A probability difference requires computing both probabilities at the stated reference risk.
- Assuming the placebo term is necessarily a placebo mechanism. It represents change under the modeled placebo condition. Natural history, time, and study context may contribute, and the model alone does not distinguish their causes.
- Reporting a probability contrast without its reference. State the baseline or placebo risk and the concentration at which the prediction is made. Without them, the probability difference is not fully interpretable.
In practice, report the definition of the event, baseline reference, placebo condition, concentration, and contrast scale. An odds ratio can summarize the relative odds effect, while predicted event probabilities and their difference communicate the corresponding absolute effect for a specified reference risk.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- If \(p_B=0.25\) and \(\delta_P=\log(2)\), what is the placebo event probability?
- If \(E(C)=\log(0.70)\), what is the active-versus-placebo odds ratio, and what does its value mean?
- Why can the same active-versus-placebo odds ratio yield different risk differences in two populations?
- Suppose active-versus-baseline risk difference is reported as the drug effect, although placebo risk differs from baseline. What is the problem?
- For the Emax exposure function, what is \(E(C)\) at \(C=EC_{50}\), and what is its limiting value as \(C\) becomes large?
Solutions. (1) Baseline odds are \(0.25/0.75=1/3\). Placebo odds are \((1/3)(2)=2/3\), giving \(p_P=(2/3)/(1+2/3)=0.400\). (2) The odds ratio is \(\exp\{\log(0.70)\}=0.70\): active-treatment event odds are 30% lower than placebo event odds at that exposure. This is not a 30% reduction in event probability. (3) Converting odds to probability depends on the placebo reference probability, so a common odds ratio can map to different absolute probability changes. (4) That contrast combines the placebo-related change with the incremental drug effect; for the drug effect under this parameterization, compare active treatment with placebo. (5) Substitution gives \(E(EC_{50})=E_{\max}/2\), and as \(C\) increases without bound, \(E(C)\) approaches \(E_{\max}\).