Why Specify a Baseline Model?
A continuous pharmacodynamic (PD) endpoint may differ among individuals before treatment, even when everyone is observed under the same treatment conditions. A model with one common baseline assigns the same expected starting response to everyone; a covariate-dependent baseline allows that expectation to vary with measured characteristics. The choice affects predictions and can affect estimates of placebo, disease-progression, and drug-response components when baseline variation is not represented adequately.
Here, “baseline” means the model component describing expected response at the reference conditions specified for the other components. It is not necessarily the same as every individual’s measured baseline value. The focus is the form of that component for a continuous endpoint, not the separate question of whether baseline should be treated as an initial condition or as an estimated parameter.
For an additive model on the endpoint scale, a compact representation is
where \(Y_{it}\) is the continuous PD response for individual \(i\) at time \(t\), and \(g_{it}\) collects the other modeled components under those conditions. The earlier material on decomposing baseline, placebo, and natural-history response motivates separating these contributions. The baseline model specifies how the baseline component varies; it does not, by itself, specify a placebo or drug effect.
Constant and Covariate-Dependent Forms
The simplest baseline specification is
This assumes that, at the chosen reference conditions, all individuals share the same expected baseline response. Their observed responses can still differ because of residual variation or other modeled components. A constant baseline is plausible when meaningful baseline predictors are absent, when the study population is sufficiently homogeneous for the intended purpose, or when a simpler model is supported by the data.
A common covariate-dependent form is linear in a continuous predictor \(z_i\):
Here \(z_{\mathrm{ref}}\) is a chosen reference value. The intercept \(\beta_{B0}\) is the expected baseline response for an individual whose covariate equals that reference, and \(\beta_{B1}\) is the expected change in baseline response for a one-unit increase in \(z_i\). The slope’s units are endpoint units per covariate unit. Centering does not change fitted values when the intercept is adjusted accordingly; it changes the interpretation and often improves numerical stability.
For several predictors, an additive linear baseline can be written as
Each coefficient describes the baseline association with its covariate while the other included covariates are held fixed. This is a conditional model statement, not evidence that changing a covariate would cause a change in response. Categorical predictors can also be represented with indicator variables, with the reference category and resulting coefficient contrasts stated explicitly.
A covariate main effect on baseline is not automatically a covariate effect on the drug response. For example, a model may allow older individuals to have a higher expected baseline while imposing the same drug-effect function at every age. A covariate-by-concentration or covariate-by-treatment interaction is a separate modeling choice and should be included only when the scientific question and data support it.
Worked Examples
Worked Example: Predicting Response With a Constant Baseline
Suppose a hypothetical continuous endpoint is measured in units. At the reference time, a model uses a constant baseline of \(B_0=18\) units. Under a particular treatment condition, the modeled placebo contribution is \(+1.5\) units and the drug contribution is \(-4.0\) units. Assume these components add on the endpoint scale.
The expected response is
Substituting the hypothetical values gives
Thus the model predicts 15.5 units under those conditions, for every individual to whom the common baseline applies. If one hypothetical observed response is 15.1 units, its residual relative to this prediction is \(15.1-15.5=-0.4\) units. That residual is not a new baseline estimate; it is the difference between the observation and the model prediction.
This calculation relies on the stated additive scale and on the modeled components being defined at compatible reference conditions. It does not establish that the baseline is truly identical across individuals; that is an assumption to assess against the study design and data.
Worked Example: A Centered Covariate Changes Expected Baseline
In a hypothetical analysis, baseline response is modeled as a function of age. Let the reference age be 50 years, with estimated baseline parameters \(\beta_{B0}=20\) units and \(\beta_{B1}=0.30\) units per year. The model is
For a 40-year-old individual, the predicted baseline is
For a 60-year-old individual, it is
The model therefore predicts a difference of \(23.0-17.0=6.0\) units between these ages, or \(0.30\times20=6.0\) units. This is a modeled baseline association across a 20-year age difference; it is not a causal claim about aging. At age 50, the covariate term is zero, so the intercept of 20 units has a direct reference-population interpretation.
For illustration, suppose the same treatment condition contributes \(-4.0\) units for both individuals and no other component is present. Their respective expected responses are \(17.0-4.0=13.0\) units and \(23.0-4.0=19.0\) units. This model allows different baseline levels while keeping the treatment contribution constant.
Worked Example: Combining Two Baseline Predictors
Consider a hypothetical model for a continuous endpoint in which both age and a centered severity score predict baseline. Let the reference age be 50 years and the reference severity score be 10 points. The fitted baseline function is
For an individual aged 60 years with a severity score of 12 points, substitute both covariates:
For a 50-year-old individual with a severity score of 8 points, the prediction is
The difference is \(29.0-21.0=8.0\) units. It can also be checked from the coefficient contributions: the first individual is 10 years older, contributing \(0.20\times10=2.00\) units, and has a severity score 4 points higher, contributing \(1.50\times4=6.00\) units. The total is \(2.00+6.00=8.00\) units. This interpretation is conditional on the other modeled predictor and on the linear additive form being adequate over the covariate range.
Why Centering and Model Choice Matter
Centering makes the intercept correspond to a specified, interpretable covariate profile. Without centering, the intercept in a linear model is the predicted baseline at covariate value zero. If zero is impossible or far outside the study population, that intercept may have little scientific meaning. Centering does not force the covariate effects to be small or make a poor functional form correct.
The following justification shows why centering preserves predictions. Start with an uncentered one-covariate model \(b(z)=a_0+a_1z\). Add and subtract the reference value inside the predictor:
Therefore, the centered form has slope \(\beta_{B1}=a_1\) and intercept \(\beta_{B0}=a_0+a_1z_{\mathrm{ref}}\). The predicted value at every \(z\) is unchanged, while the intercept now equals the model prediction at \(z_{\mathrm{ref}}\). This is an algebraic reparameterization, not a change in the model’s substantive predictions.
Diagnostics, Alternatives, and Limitations
A constant baseline is simpler and uses fewer parameters, but it can leave systematic baseline patterns unexplained. A covariate-dependent model can account for observed trends, but its extra parameters need support from adequate covariate variation and sample information. A practical comparison is not merely whether one model has a lower residual sum of squares: a more flexible model will generally fit the observed data at least as well. Consider whether the added terms are scientifically motivated, whether fit improves meaningfully, and whether predictions remain credible.
| Approach | What it assumes | Useful consideration |
|---|---|---|
| Constant baseline | One expected baseline level across individuals | Parsimonious, but may miss systematic covariate trends |
| Covariate-dependent baseline | Expected baseline varies according to a specified function of covariates | Can explain structured variation, but depends on functional form and data support |
| Separate baseline by strata | Each prespecified category has its own baseline level | Flexible across categories, but can use many parameters and gives no smooth trend between categories |
Useful checks include plotting observed baseline values against each continuous predictor, plotting residuals against predictors and fitted values, and checking whether residuals retain a trend or curvature. Compare observed and predicted values across clinically meaningful covariate ranges. For a linear baseline term, curvature in such plots may indicate that the linear relationship is inadequate. Sparse observations at the extremes make apparent trends and extrapolated predictions especially uncertain.
A covariate-dependent baseline is not a cure for every form of heterogeneity. Important predictors may be unmeasured, relationships may be nonlinear, and the residual-error assumptions may still be inappropriate. A linear model also predicts without bound as a covariate increases or decreases; that can be implausible outside the range supporting the fit. Transformations or nonlinear functions may be considered when scientifically justified, but should not be selected solely to chase noise.
The additive equations above are defined on the endpoint scale. If a model is instead specified on a transformed scale, such as the logarithm of the response, its baseline terms and coefficients describe that transformed scale. Predictions on the original scale require the appropriate back-transformation, and the back-transformed model prediction need not equal the arithmetic mean response. Scale and units must therefore be stated when interpreting the baseline function.
- Pitfall: treating a baseline covariate association as a drug-effect interaction. A positive age coefficient in the baseline function means that predicted baseline response rises with age, conditional on the other modeled covariates. It does not mean that the drug works better or worse with age. A drug-effect interaction requires an explicit term modifying the drug-response component.
- Pitfall: interpreting the intercept at an irrelevant covariate value. If an uncentered model has intercept \(a_0\), that is the prediction at \(z=0\), not automatically the average participant’s baseline. State the reference value and interpret the intercept there.
- Pitfall: assuming covariate inclusion proves causation. Regression coefficients describe the fitted conditional relationship under the model. They do not show that intervening on a predictor would change baseline response.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A constant baseline model has \(B_0=12\) units. What is its predicted baseline for an individual whose age differs from the sample mean?
- A centered linear baseline model has intercept 30 units and slope \(-0.40\) units per year, with age centered at 60 years. Calculate the predicted baseline at age 65.
- In a model, age predicts baseline response, while the drug-effect function has no age term. Does this model assert that drug effect changes with age? Explain.
- Why can the intercept from an uncentered age model be hard to interpret if age zero is outside the study population?
- Name one diagnostic pattern that could indicate that a linear covariate-dependent baseline is inadequate.
Solutions. (1) It remains 12 units: the constant model assigns everyone the same expected baseline, regardless of age. (2) Substitute into the model: \(30+(-0.40)(65-60)=30-2.00=28.0\) units. (3) No. Age changes the predicted baseline, but the stated model keeps the drug-effect function the same across ages; age modification of drug effect would require an explicit interaction or equivalent change to that component. (4) The uncentered intercept is the prediction at age zero, so it may be far outside the observed data and lack a useful interpretation for the study population. (5) A curved trend in residuals against age, or systematic curvature in observed-versus-predicted baseline values across age, suggests the linear function may not capture the relationship.