Why Anchor Samples Matter
Suppose one cohort has measurements on modality \(X\) and a second cohort has measurements on modality \(Y\), but few or no subjects have both. Separate summaries of the two modalities do not, by themselves, reveal which values of \(X\) tend to occur with which values of \(Y\). A group of subjects measured on both modalities can provide that missing link. The usefulness of this link depends not only on how many such subjects are measured, but also on which subjects and cohorts they represent.
An anchor sample is not a mathematical device for creating subject-level pairings between unrelated people. Its paired records supply direct information about cross-block dependence. That information can help connect modality-specific cohorts, provided the relationship measured among anchors is relevant to the target population. The earlier material on Estimating Cross-Omics Associations From Partial Overlap addresses association estimation with partial overlap; here the emphasis is on how the number and placement of jointly measured subjects govern the bridge.
What the Anchor Sample Estimates
For scalar features \(X\) and \(Y\), the target cross-block covariance is \(\operatorname{Cov}(X,Y)=E[(X-\mu_X)(Y-\mu_Y)]\), where \(\mu_X=E(X)\) and \(\mu_Y=E(Y)\) refer to the target population. For \(n\) anchor subjects, the sample covariance and correlation are
These quantities use the within-subject pairing. Reordering the \(Y_i\) values relative to the \(X_i\) values generally changes both estimates, even though the separate modality-specific distributions stay exactly the same. For feature vectors, the same idea yields a cross-covariance matrix: each entry pairs one feature from one modality with one feature from another modality on the same anchor subjects.
Using an anchor estimate for a broader cohort requires assumptions beyond having paired observations. The anchors must be measured comparably to the target data, their sampling and missingness must not select subjects in a way that distorts the relationship of interest (or that selection must be accounted for), and the cross-modality relationship must transport to the target population. If the target spans different cohorts or subgroups, anchor coverage should reflect those differences. A large anchor sample from one narrow subgroup cannot establish that the same relationship holds everywhere else.
Worked Example: Calculating a Cross-Block Association From Paired Anchors
Consider five hypothetical anchor subjects with paired measurements \(X=(1,2,3,4,5)\) and \(Y=(2,1,4,3,5)\), in arbitrary assay units. Assume subjects are independent draws from the population of interest, and both measurements are available and comparable for every listed subject.
The sample means are \(\bar{X}=3\) and \(\bar{Y}=3\). The paired deviations are \((-2,-1,0,1,2)\) for \(X\) and \((-1,-2,1,0,2)\) for \(Y\). Their product sum is \(2+2+0+0+4=8\); the sum of squared deviations for each variable is \(10\). Thus
Therefore \(r_{XY}=2/(\sqrt{2.5}\sqrt{2.5})=2/2.5=0.800\). In these five hypothetical anchors, the measured features have a positive sample association. This does not establish causation, nor does a sample of five precisely determine the population association. If the \(Y\) values had been paired differently, the marginal means and standard deviations would remain unchanged, but the cross-block estimate could differ.
Anchor Count and Precision
More independent anchors generally provide more information about a cross-block relationship, but the gain is not a guarantee of low uncertainty if the sample is unrepresentative or the measurements are noisy. Under independent, identically distributed bivariate normal observations, an approximate interval for a population correlation \(\rho\) can be constructed on the Fisher-transformed scale. Define \(z_r=\frac{1}{2}\log\{(1+r)/(1-r)\}\). For moderate or large \(n\), its approximate standard error is \(1/\sqrt{n-3}\). A 95% interval on that scale is \(z_r\pm1.96/\sqrt{n-3}\); applying the inverse transformation \(\tanh(\cdot)\) gives an approximate interval for \(\rho\). This is an approximation, not an exact small-sample result, and relies on the stated distributional and sampling conditions.
The formula makes the role of anchor count explicit: increasing \(n\) decreases the approximate standard error. But count alone is not enough. Dependence among subjects, influential observations, restricted ranges, and subgroup imbalance can make the nominal precision misleading or make an estimate irrelevant to the target.
Worked Example: Comparing Two Anchor-Sample Sizes
Suppose two hypothetical studies observe the same sample correlation, \(r=0.600\), under independent bivariate normal sampling. Study A has \(n=12\) anchors and Study B has \(n=48\). The Fisher transform of the observed correlation is
For Study A, \(SE_z=1/\sqrt{12-3}=1/3=0.3333\). Its approximate 95% interval on the transformed scale is \(0.6931\pm1.96(0.3333)\), or \([0.040,\ 1.346]\). Applying \(\tanh\) to the endpoints gives an approximate interval for \(\rho\) of \([0.040,\ 0.873]\).
For Study B, \(SE_z=1/\sqrt{48-3}=1/\sqrt{45}\approx0.1491\). The transformed interval is \(0.6931\pm1.96(0.1491)=[0.401,\ 0.985]\). Transforming back gives approximately \([0.381,\ 0.755]\). Both estimates are \(r=0.600\), but the larger anchor sample produces a narrower approximate interval. These intervals describe the uncertainty of a procedure under its assumptions; they do not mean there is a 95% probability that the fixed population correlation lies in either calculated interval.
Placement: Which Subjects Should Be Anchors?
Placement has at least two meanings. First, anchors should cover the cohorts, strata, or covariate ranges to which the cross-block estimate will be applied. Second, when estimating a relationship such as a regression slope, anchors need useful variation in the predictor. A sample concentrated near one predictor value may estimate a relationship poorly even if its total size appears adequate.
For a simple linear model \(Y_i=\beta_0+\beta_1X_i+\varepsilon_i\), assume the conditional mean is linear, errors are independent with mean zero and constant variance \(\sigma^2\), and the observed \(X_i\) values are treated as fixed (or are independent of the errors). The least-squares slope is \(\hat{\beta}_1=\sum_i(X_i-\bar{X})(Y_i-\bar{Y})/\sum_i(X_i-\bar{X})^2\). Given the \(X_i\) values, its variance is
This follows by writing \(\hat{\beta}_1-\beta_1=\sum_i(X_i-\bar{X})\varepsilon_i/S_{XX}\), where \(S_{XX}=\sum_i(X_i-\bar{X})^2\). Independence and equal error variance give \(\operatorname{Var}(\hat{\beta}_1\mid X)=\sum_i(X_i-\bar{X})^2\sigma^2/S_{XX}^2=\sigma^2/S_{XX}\). Thus, greater spread in anchor \(X\) values reduces slope variance under this model. It does not justify extrapolating a linear relationship beyond the range or populations represented by the anchors.
Worked Example: Same Anchor Count, Different Predictor Spread
Compare two hypothetical sets of eight anchors for estimating the association of an assay feature \(Y\) with a feature \(X\). Assume the linear model above is appropriate and the error standard deviation is known to be \(\sigma=1\) assay unit. In Design A, four anchors have \(X=0\) and four have \(X=1\); in Design B, four have \(X=0\) and four have \(X=4\). For both designs, let the corresponding \(Y\) values be the group mean plus residuals \((0,1,-1,0)\): the means are 10 at \(X=0\), 10.5 at \(X=1\) in A, and 12 at \(X=4\) in B. This gives the same observed slope, \((10.5-10)/(1-0)=0.500\) in A and \((12-10)/(4-0)=0.500\) in B.
For Design A, \(\bar{X}=0.5\), so \(S_{XX}=4(0-0.5)^2+4(1-0.5)^2=2\). Its slope standard error is \(1/\sqrt{2}\approx0.707\) assay units of \(Y\) per unit of \(X\). For Design B, \(\bar{X}=2\), so \(S_{XX}=4(0-2)^2+4(4-2)^2=32\). Its slope standard error is \(1/\sqrt{32}\approx0.177\) assay units of \(Y\) per unit of \(X\). Under the assumed model, the wider predictor spread gives Design B greater precision despite the same number of anchors and the same fitted slope. In practice, the wider design is only preferable if the linear relationship remains appropriate over that range and the extreme values are credible and measured comparably.
Bridging Cohorts and Choosing an Anchor Allocation
When modality-specific cohorts differ, the target relationship must be specified before deciding where to place anchors. Let \(G\) denote cohort membership. The overall covariance decomposes as
The first term summarizes within-cohort cross-modality covariance, averaged over the target cohort distribution. The second captures association between cohort-specific means. Consequently, a pooled anchor covariance can reflect both within-cohort relationships and between-cohort differences. If anchors come only from one cohort, they can estimate that cohort's relationship but cannot, without a transport assumption, reveal whether the within-cohort relationship is the same elsewhere or determine the target's between-cohort contribution.
If the goal is a target-population estimate, allocate anchors to the cohorts and important subgroups that define that target. If the goal is a within-cohort relationship, ensure enough anchors within the cohort for that estimate. When a cohort is rare in the target but scientifically important, disproportionate anchor sampling may be useful, but the resulting estimate must be combined using target-population weights rather than treated as a simple unweighted average. If the relationship is believed invariant across cohorts, anchors can be pooled to improve precision, but that invariance is an assumption to scrutinize, not a consequence of having many paired subjects.
Alternatives, Diagnostics, and Limits
Complete overlap, in which every subject has every modality, gives direct paired information throughout the study but may be costly or infeasible. A deliberately chosen anchor subset is more economical, though it trades measurement cost for dependence on how representative and well placed the subset is. With no anchors, approaches such as model-based imputation or other partial-overlap methods may still estimate cross-block quantities under additional modeling and identification assumptions; separate marginal data alone do not identify an arbitrary cross-block covariance. Anchors contribute empirical paired evidence, but do not remove assumptions about selection, measurement comparability, or transport.
Before using an anchor estimate, check the pairing and data alignment, compare anchor and target cohort composition, inspect the observed range of each feature, and look for influential measurements. Examine whether the estimated relationship differs meaningfully across cohorts or key strata when enough anchors exist to make that comparison informative. Sparse coverage in a subgroup is not evidence that the relationship is identical there. A linear model should also be checked for systematic departures from linearity and changing residual spread; wider predictor coverage cannot repair a misspecified relationship.
A common mistake is to claim that a large anchor count guarantees a valid bridge. The reasoning is wrong because sample size mainly controls sampling variability under the assumed model; it does not fix selection bias, poor cohort coverage, assay incompatibility, or an incorrect transport assumption. The correction is to assess both precision (how much information the anchors provide) and relevance (whether those anchors represent the target relationship). Anchor estimates can be precise for the wrong population.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- Two cohorts have measurements on \(X\) and \(Y\), respectively, but share no measured subjects. Can their separate means and variances identify \(\operatorname{Cov}(X,Y)\)? Why?
- Under the Fisher approximation, what is the transformed-scale standard error for \(n=28\) independent bivariate normal anchors?
- For a linear regression with known error standard deviation \(\sigma=2\), calculate the slope standard error when \(S_{XX}=16\).
- Why might an anchor sample of 100 subjects from one cohort be less useful for a multi-cohort target than 50 anchors drawn across the target cohorts?
Solutions. (1) No. Separate marginal distributions do not supply the subject-level pairings needed for covariance; additional assumptions or paired data are required. (2) \(SE_z=1/\sqrt{28-3}=1/5=0.200\). (3) \(SE(\hat{\beta}_1)=\sigma/\sqrt{S_{XX}}=2/\sqrt{16}=0.500\) units of \(Y\) per unit of \(X\). (4) The 100 anchors may estimate one cohort's relationship precisely but provide little evidence about relationships or mean differences in the other target cohorts. The 50 distributed anchors may better cover the target, although their precision depends on the within-cohort counts and the intended target estimand.