Why the Representation of Ancestry Matters
An ancestry proportion and a population label describe related but different information. A proportion records gradations: for example, an estimated genome-wide ancestry fraction of 0.35 differs from 0.65. A label assigns an individual to a category, such as a cluster or a study-defined population group. Replacing proportions with labels may make a model easier to describe, but it discards differences among people assigned to the same category.
That loss matters in association analysis when ancestry is related both to genotype dosage and to the outcome. A label may control broad differences between groups while leaving a residual association between dosage and ancestry within each group. Conversely, a continuous linear term may be too restrictive if the relationship with ancestry is nonlinear. Neither representation is automatically correct: the choice depends on the scientific meaning of the measure, its relationship with the variables, and whether the model fits adequately.
Two Regression Representations
Let \(Y_i\) be an outcome for individual \(i\), \(X_i\) a genetic dosage, and \(A_i\) a continuous ancestry proportion between 0 and 1. For a quantitative outcome, a simple continuous-ancestry model is
Here, \(\beta_X\) is the expected change in outcome units per one-unit increase in dosage, conditional on the ancestry proportion. The coefficient \(\gamma\) describes the fitted outcome difference for a change of 1.0 in the ancestry proportion, holding dosage fixed. This is a linear conditional-mean model: it assumes the association of the outcome with \(A_i\) can be represented by a straight line, and that the dosage slope is constant across ancestry values. Other prespecified covariates can be added.
For \(K\) population labels, a label-adjusted model can instead be written as
where \(L_i\) is the label and \(I(L_i=k)\) equals 1 when individual \(i\) has label \(k\), and 0 otherwise. The omitted label is the reference. The \(\delta_k\) terms permit different group intercepts. This specification assumes a common dosage slope across labels and treats individuals within each label as having the same fitted ancestry contribution to the mean outcome.
For standard coefficient tests and confidence intervals, a common set of regression conditions is needed: the conditional mean is correctly specified, observations are independent (or dependence is handled in the variance estimate), and there is adequate variation in dosage after adjustment. Normal errors with constant variance support exact small-sample \(t\)-inference; with larger samples, robust standard errors can address some variance misspecification, but do not repair a wrong conditional mean or dependence that was ignored. Ancestry proportions are estimates, so substantial measurement error can also weaken or distort adjustment.
What Label Adjustment Removes—and What It Leaves
The distinction can be made precise by considering a model with one continuous ancestry measure. For clarity, suppose the models include an intercept and no other covariates. Let \(S_{UV}=\sum_i(U_i-\bar U)(V_i-\bar V)\) denote a centered cross-product. In the continuous model, the dosage coefficient is
The numerator and denominator are the cross-product and sum of squares remaining after each variable is linearly adjusted for ancestry. With additional covariates, the same interpretation applies after first residualizing the relevant variables with respect to those covariates.
In the label model, the corresponding calculation first subtracts each label-specific mean. Write \(X_i-\bar X_{L_i}\) for the dosage deviation from its label mean, and define the analogous deviations for \(Y_i\) and \(A_i\). The label-adjusted dosage slope is the within-label cross-product of dosage and outcome divided by the within-label dosage sum of squares. If the conditional mean contains a continuous ancestry term, then within labels the outcome has a component \(\gamma(A_i-\bar A_{L_i})\). Thus the label-adjusted slope contains the term
This term vanishes only if the within-label dosage and ancestry deviations have zero cross-product (or if \(\gamma=0\)). A label can therefore remove between-label differences without removing within-label ancestry confounding.
Worked Example: A Label-Adjusted Slope Retains Within-Label Variation
Consider eight hypothetical individuals, with dosage \(X\), ancestry proportion \(A\), two labels, and outcome \(Y=2X+3A\). The outcome is defined exactly by this equation for illustration; there is no random error in these values.
| Label | \(X\) | \(A\) | \(Y\) |
|---|---|---|---|
| Low | -1 | 0.10 | -1.70 |
| Low | -1 | 0.20 | -1.40 |
| Low | 1 | 0.30 | 2.90 |
| Low | 1 | 0.40 | 3.20 |
| High | -1 | 0.60 | -0.20 |
| High | -1 | 0.70 | 0.10 |
| High | 1 | 0.80 | 4.40 |
| High | 1 | 0.90 | 4.70 |
Within each label, the dosage mean is 0. The ancestry deviations from the label mean are \(-0.15,-0.05,0.05,0.15\), and the corresponding dosage values are \(-1,-1,1,1\). Thus, for each label,
Across both labels, the within-label dosage sum of squares is \(8\), and the within-label dosage-ancestry cross-product is \(0.40+0.40=0.80\). Since \(Y=2X+3A\), the within-label dosage-outcome cross-product is
Consequently, the label-adjusted slope is \(18.4/8=2.3\) outcome units per dosage unit. It is not the dosage coefficient 2 in the generating equation: the remaining 0.3 is \(3(0.80/8)\), the contribution of ancestry variation correlated with dosage within labels.
For continuous adjustment, the centered sums over all eight individuals are \(S_{XX}=8\), \(S_{XA}=0.80\), \(S_{AA}=0.60\), \(S_{XY}=18.4\), and \(S_{AY}=3.4\). For example, \(S_{XY}=2S_{XX}+3S_{XA}=16+2.4=18.4\). Substitution gives
This recovers the specified dosage coefficient because the outcome is exactly linear in \(X\) and \(A\), and there is enough variation to estimate both terms. In an actual sample, random error would generally make the estimate differ from 2.
Worked Comparison: Categorizing a Smooth Gradient
Worked Example: The Same Label Contains Different Ancestry Values
Suppose four hypothetical individuals follow \(Y=10+4A\), with outcome measured in a specified quantitative scale. A study assigns the first two individuals to a “lower” label and the last two to a “higher” label.
| Individual | \(A\) | Label | \(Y=10+4A\) |
|---|---|---|---|
| 1 | 0.10 | Lower | 10.4 |
| 2 | 0.30 | Lower | 11.2 |
| 3 | 0.60 | Higher | 12.4 |
| 4 | 0.90 | Higher | 13.6 |
The lower-label mean outcome is \((10.4+11.2)/2=10.8\), and the higher-label mean is \((12.4+13.6)/2=13.0\). A label-only model predicts 10.8 for either lower-label individual and 13.0 for either higher-label individual. Its residuals, observed minus fitted, are \(-0.4,0.4,-0.6,0.6\). The label model captures the broad difference between groups, but not the gradient within either group.
The continuous model has an intercept of 10 and a slope of 4 outcome units per one-unit increase in \(A\). For a 0.10 increase in ancestry proportion, its fitted outcome difference is \(4(0.10)=0.40\) units. This interpretation is conditional on the assumed linear relationship; it does not establish that changing ancestry causes an outcome change.
Choosing, Checking, and Comparing the Models
A useful comparison begins with what the variable represents. A continuous proportion is appropriate when the scientific question concerns variation along an ancestry gradient and the estimates are sufficiently informative. A label may be appropriate when the categories correspond to the study's sampling strata or when the estimand is explicitly a within-category association. A label is not made more informative simply by being familiar or easy to communicate.
For multiple ancestry proportions \(A_{i1},\ldots,A_{iK}\) that sum to 1, including all \(K\) proportions with an intercept causes exact collinearity. Use \(K-1\) proportions with a reference component, or another identifiable parameterization. This is distinct from replacing proportions by a label: a \(K\)-category label also needs a reference category, but it represents group membership rather than a full quantitative composition.
Diagnostics should target the assumptions that distinguish the representations. For a continuous model, inspect outcome residuals across the range of each ancestry proportion and check whether a straight-line conditional mean is plausible. For a label model, examine ancestry-proportion distributions within each label and check for residual outcome patterns across those within-label values. In either model, assess dosage variation within ancestry strata, influential observations, and whether the fitted design matrix is nearly collinear. Very little dosage variation after adjustment makes the dosage coefficient imprecise even if the model is otherwise reasonable.
State whether the covariate is an estimated continuous proportion, a categorical label, or both; explain the reference category and scale.
Use continuous terms for a defensible gradient model, labels for category-specific intercept adjustment, and include other covariates required by the study design.
Check whether dosage, ancestry, and outcome still co-vary within labels and whether continuous-model residuals show curvature or systematic structure.
Describe the dosage coefficient as conditional on the covariates and representation used, and report uncertainty and units without treating adjustment as proof of causality.
A practical alternative is to include both ancestry proportions and label indicators, if the design and data support estimating their separate contributions. This can accommodate broad category differences while adjusting for within-label gradients, but it can also increase collinearity and uncertainty. A flexible nonlinear function of the proportion may be preferable to arbitrary categorization when a linear trend is implausible and data cover the range adequately. Each extension requires enough observations and dosage variation to estimate its parameters.
Common Pitfall and Limitations
- Wrong reasoning: “The model includes population labels, so ancestry confounding is removed.” Correction: Labels adjust between-label means; they do not automatically adjust ancestry variation within a label. Check within-label covariance and residual patterns.
- Wrong reasoning: “A continuous proportion must be better because it contains more numerical detail.” Correction: A misspecified straight-line relationship, uncertain proportions, or poor overlap can make a simple continuous model inadequate. Evaluate fit and measurement quality.
- Wrong reasoning: “The adjusted dosage coefficient is a causal genetic effect.” Correction: Regression adjustment defines a conditional association under the model. Causal interpretation requires additional design and identification assumptions.
Both representations have limits. Labels depend on how groups were defined and can hide diversity within categories; continuous proportions depend on estimation and can be difficult to interpret as a single gradient when ancestry is multidimensional. Either model may have weak overlap in ancestry or dosage across the sample, limiting what comparisons the data support. Report the ancestry representation, coding, adjustment set, and diagnostics so that the coefficient's meaning is explicit.
Check Your Understanding
Answer each question, then compare your reasoning with the brief solutions.
- In the first worked example, why does label adjustment produce a dosage slope larger than 2?
- What does a coefficient of 4 mean when \(A\) is a proportion on the 0-to-1 scale?
- Why can a label-only model show a residual pattern even when its group means are correctly estimated?
- What coding issue arises if all \(K\) ancestry proportions and an intercept are included in the same linear model?
Solutions. (1) Within labels, \(X\) and \(A\) have cross-product 0.80 over dosage sum of squares 8. The omitted within-label ancestry contribution is \(3(0.80/8)=0.30\), giving \(2+0.30=2.30\). (2) It is the fitted outcome difference for a 1.0 increase in proportion, holding other model terms fixed; a 0.10 increase corresponds to \(4(0.10)=0.40\) outcome units. (3) Group means capture between-label differences, but individuals in a label can still differ in ancestry and outcome systematically. (4) Because the proportions sum to 1, their columns are exactly linearly dependent with the intercept; omit a component or use an identifiable alternative coding.