Why the Covariate Scale Matters
In a basket trial, patient-level covariates can improve adjustment for baseline risk, but their units determine what model coefficients mean. If age is entered in years, a coefficient describes a one-year change in log odds. If age is centered at a clinically meaningful value and scaled in decades, the coefficient instead describes a ten-year change, and the intercept refers to a patient at the chosen age. These changes can make parameters easier to interpret and often put them on more manageable numerical scales.
The transformation does not, by itself, change the information in the data or the fitted probabilities. It does change the parameterization. In a Bayesian model, that distinction matters: using the same-looking priors after changing scales will generally specify a different prior model.
Model and Definitions
Let \(Y_{ij}\) be a binary response for patient \(i\) in basket \(j\), \(T_{ij}\) indicate treatment assignment, and \(x_{ij}\) be a continuous patient-level covariate. For illustration, consider the logistic model
Here \(\alpha_j\) is basket \(j\)'s control log odds at \(x_{ij}=0\), \(\theta_j\) is its treatment log odds ratio, and \(\gamma\) is the common covariate coefficient. A hierarchical specification might model \(\theta_j\) as arising from a distribution shared across baskets; the covariate transformation discussed here does not change that hierarchy in this no-interaction model.
This formulation assumes conditionally independent Bernoulli outcomes given the model parameters and covariates, a suitable logistic mean structure, and a common covariate slope across baskets. It also assumes no treatment-by-covariate interaction. The last condition is important: without an interaction, the treatment log odds ratio \(\theta_j\) is constant across values of \(x\). The model is an example, not a universal prescription; a continuous endpoint or a different link would require the corresponding likelihood and interpretation.
How Coefficients Change
Substitute \(x_{ij}=c+s z_{ij}\) into the raw-scale model. Collecting the constant and covariate terms gives
Thus the transformed model has intercept \(a_j=\alpha_j+\gamma c\) and covariate coefficient \(b=\gamma s\). The inverse transformation is \(\alpha_j=a_j-bc/s\) and \(\gamma=b/s\). These relations establish that the linear predictor, and hence each fitted probability, can remain exactly the same under the change of coordinates.
The transformed intercept \(a_j\) is the control log odds for a patient with \(x=c\), not for a patient with \(x=0\). The coefficient \(b\) is the change in control log odds for an increase of \(s\) units in the original covariate. If \(c\) is the pooled mean age and \(s\) is 10 years, for example, \(a_j\) refers to a patient at the pooled mean age and \(b\) describes a 10-year increase.
Worked Example: Re-expressing an Age Effect
Hypothetical model values. Suppose a basket's control linear predictor is \(-2.4+0.04x\), where age \(x\) is measured in years. Re-express it using \(c=60\) years and \(s=10\) years.
Define \(z=(x-60)/10\), so \(x=60+10z\). Substitution gives
The centered intercept is \(0\), the control log odds at age 60, and the scaled slope is \(0.40\), the log-odds change per 10 years. At age 70, \(z=1\), so the modeled response probability is \(1/(1+\exp(-0.40))\approx0.5987\). At age 60 it is \(1/(1+\exp(0))=0.5000\). The raw intercept \(-2.4\) referred to age zero, which is not a useful reference for this adult-patient model.
Centering and Treatment Effects
With no treatment-by-covariate interaction, centering a covariate changes the interpretation of the intercept but not the treatment contrast. In the model above, the treatment log odds ratio in basket \(j\) is \(\theta_j\) at every covariate value. Adjusting for a prognostic covariate may change its posterior estimate through the likelihood, but re-expressing that covariate alone does not make \(\theta_j\) conditional on a particular age.
If effect modification is scientifically relevant, include an interaction explicitly. For example,
Here \(\theta_j\) is the treatment log odds ratio at \(x=c\), and \(\delta\) is the change in that log odds ratio for each \(s\)-unit increase in \(x\). Centering therefore gives the treatment effect a reference value that can be stated and defended. The interaction is not created by centering; it must be part of the model if treatment effects are allowed to vary with the covariate.
Worked Example: Treatment Effects at Different Ages
Hypothetical model values. Let \(c=60\) years, \(s=10\) years, and suppose one basket has control linear predictor \(-1+0.30z\), treatment log odds ratio \(0.50\) at age 60, and interaction coefficient \(\delta=0.20\). The treatment log odds ratio is \(0.50+0.20z\).
At age 50, \(z=(50-60)/10=-1\), so the treatment log odds ratio is \(0.50+0.20(-1)=0.30\). The corresponding odds ratio is \(\exp(0.30)\approx1.350\). The control linear predictor is \(-1+0.30(-1)=-1.30\), giving probability \(1/(1+\exp(1.30))\approx0.2142\). Under treatment, the linear predictor is \(-1.30+0.30=-1.00\), giving probability \(1/(1+\exp(1.00))\approx0.2689\).
At age 70, \(z=1\), so the treatment log odds ratio is \(0.50+0.20=0.70\), with odds ratio \(\exp(0.70)\approx2.014\). The control linear predictor is \(-1+0.30=-0.70\), giving probability \(0.3318\); the treated linear predictor is \(-0.70+0.70=0\), giving probability \(0.5000\). These are model-based contrasts at two covariate values, not evidence that age causes a different response. The example shows why an interaction model's treatment coefficient needs a clearly stated reference age.
Choosing the Center and Scale
A common center across baskets is usually useful when intercepts or treatment effects are to be compared. A pooled mean, a clinically meaningful threshold, or a prespecified reference value can serve, provided the choice is documented. If each basket is centered at its own mean, each intercept refers to a different covariate profile. Those intercepts are not directly comparable as baseline risks without translating them to a common reference.
Centering at the pooled mean is not automatically best. If the mean lies outside a clinically relevant range for a basket, its intercept may describe a patient unlike those in that basket. A meaningful reference value can be preferable. Scaling by a standard deviation makes coefficient magnitudes comparable across covariates, while scaling by a clinically meaningful increment may yield more direct interpretation. Scaling does not require the covariate to be normally distributed.
For a binary covariate, centering at zero or one can make the intercept refer to a clearly defined group. Standardizing a binary variable is possible, but may make the coefficient less clinically transparent. The choice should follow the intended interpretation and numerical needs rather than a rule that every covariate must be standardized.
Worked Example: Transforming the Prior Correctly
Hypothetical prior specification. On the raw age scale, suppose \(\alpha_j\sim N(-1,0.5^2)\) and \(\gamma\sim N(0,0.01^2)\), independently. Let \(c=60\) years and \(s=10\) years. The transformed parameters are \(a_j=\alpha_j+60\gamma\) and \(b=10\gamma\).
Their prior means are \(E(a_j)=-1\) and \(E(b)=0\). Using independence of \(\alpha_j\) and \(\gamma\),
and
The transformed priors are therefore correlated, even though the raw-scale priors were independent. Their correlation is \(0.06/(\sqrt{0.61}\times0.10)\approx0.768\). Assigning independent priors \(a_j\sim N(-1,0.5^2)\) and \(b\sim N(0,0.1^2)\) would not reproduce the original prior: in particular, it omits the covariance and gives \(a_j\) a different variance. Such a choice may be reasonable, but it is a different prior specification and should be justified on the centered scale.
Justification, Checks, and Limitations
The algebraic substitution shows that centering and scaling preserve the likelihood when the parameter transformation is one-to-one and the coefficients are transformed consistently. In Bayesian analysis, preserving the full model also requires transforming the prior distribution consistently. Equivalently, one may define priors directly on the new parameters, but then those priors express the intended beliefs on the new scale rather than automatically inheriting the old ones. If transforming a density explicitly, its change-of-variables Jacobian must also be included.
Before fitting, check that the chosen \(c\) is a meaningful reference and that \(s>0\) is stated in the original units. Confirm that the resulting \(z\) values are of a numerically reasonable magnitude and that the model's implied probabilities are plausible over the observed covariate range. Prior predictive simulation is useful for detecting priors that imply near-certain response probabilities or implausibly steep covariate effects. After fitting, inspect whether the covariate has enough variation within baskets to support the intended adjustment; a common slope cannot be learned well from little within-basket variation.
Centering does not correct a misspecified functional form, omitted interactions, measurement error, or lack of covariate overlap between treatment groups. Scaling does not make a coefficient clinically important, and it cannot create information where covariate values are sparse. If the model includes nonlinear terms or interactions, transform and interpret those terms explicitly; changing the main-effect reference can alter their coefficient meanings.
A common mistake is to say that centering changes only notation, then retain all raw-scale priors unchanged as independent priors on the new coefficients. The likelihood may be algebraically equivalent, but the prior distribution is not, as the prior example demonstrates. Another mistake is to interpret a centered intercept as the risk for an average patient regardless of the model. In a logistic model, a prediction at mean covariates is not generally the same as the average predicted risk across patients. State the reference covariate value and treatment condition precisely.
Using the raw covariate scale keeps coefficients in familiar units but can leave an intercept at an irrelevant value and create awkward coefficient magnitudes. Centering alone improves the intercept reference while retaining the original slope units. Scaling alone changes the slope units but leaves the intercept tied to \(x=0\). Centering and scaling together can improve both interpretability and computation, but only when the chosen reference, scale, and prior implications are made explicit.
Check Your Understanding
Answer each question, then compare with the solutions.
- A model uses age in years with \(c=55\) and \(s=5\). What value of \(z\) corresponds to age 65, and what age change does a one-unit increase in \(z\) represent?
- In a logistic model without a treatment-by-age interaction, does centering age change the treatment log odds ratio? Explain.
- In a model with treatment log odds ratio \(\theta+\delta z\), what does \(\theta\) represent when \(z=(x-c)/s\)?
- Why can basket-specific centering make basket intercepts difficult to compare?
- If independent priors are assigned to a raw intercept and raw slope, must their centered and scaled versions also be independent? Explain.
Solutions. (1) \(z=(65-55)/5=2\); one unit represents a five-year increase. (2) No. With no interaction, the treatment contrast is constant over age; centering changes the intercept reference, not the treatment contrast's algebraic meaning. (3) \(\theta\) is the treatment log odds ratio at \(x=c\), because \(z=0\) there. (4) Each intercept then refers to its own basket's mean covariate value, so differences combine basket effects with differences in reference profiles. (5) No. Since the centered intercept contains a multiple of the slope, the transformed parameters are generally dependent; the transformed joint prior must account for that dependence to preserve the original prior model.