Tutorials › Statistical Methods for Single-Cell RNA Sequencing › Defining Target Comparisons and Estimands

Multi-Sample and Multi-Condition Study Designs · Tutorial 567 of 700

Defining Target Comparisons and Estimands

Specify exactly which population quantity a multi-sample single-cell study is designed to estimate, and distinguish that estimand from its sample estimate and causal interpretation.

Advanced 10 min read

What You'll Learn

  • Identify the population, conditions, outcome summary, and contrast that jointly define an estimand.
  • Distinguish a donor-average outcome from a cell-weighted outcome.
  • Interpret differences on proportion and log-expression scales with appropriate units.
  • Use target-population weights to define marginal contrasts across subgroups.
  • Separate descriptive condition contrasts from causal effects.
  • Diagnose estimands that are ambiguous, unsupported by the sampling frame, or mismatched to the observed assay.

Why Define the Target Before Choosing an Analysis?

A multi-sample single-cell study can produce many summaries from the same assay: cell-type proportions, expression values within a cell type, or contrasts on a log-ratio scale, for example. These summaries answer different questions. Even after choosing an outcome, a comparison can change meaning depending on which donors the study targets and how subgroup results are combined. Defining the target quantity first makes the scientific question explicit; the analysis can then be judged by whether it estimates that quantity.

The central object is the estimand: the precisely defined population quantity a study seeks to learn. It is not the numerical result calculated from the observed donors, and it is not a test statistic or p-value. The earlier material on experimental units and sampling frames motivates specifying who contributes independent biological information; here the focus is the target comparison itself.

Definition: An estimand for a multi-sample comparison specifies (1) the target population, (2) the conditions being compared, (3) the outcome and its summary scale, and (4) the contrast between conditions. For a donor-level outcome \(Y_i(c)\) under condition \(c\), a population mean is \(\theta_c=E_P[Y_i(c)]\), where \(P\) denotes the target population. A difference estimand is \(\Delta=\theta_b-\theta_a\).

The notation \(Y_i(c)\) denotes the outcome for donor \(i\) under condition \(c\). It is a potential-outcome notation: for each donor and condition there is a conceptually defined outcome, though typically only the outcome under the donor’s observed condition is available. The population \(P\) must be described in scientific terms, such as eligible adults with a specified disease, rather than merely as “the donors in the data.” The outcome definition should specify the cell type, assay-derived quantity, and scale whenever these matter.

An estimand may be descriptive or causal. A descriptive contrast compares outcome distributions or averages among populations observed under different conditions. A causal contrast compares potential outcomes under alternative conditions for the same target population. Calling a contrast causal requires a defensible interpretation of the conditions as interventions and assumptions that connect observed data to those potential outcomes, including appropriate assignment or confounding control. A condition label alone does not establish causation.

Choosing the Outcome Summary and Contrast

In single-cell studies, one essential decision is the level at which the outcome is summarized. A donor-level cell-type proportion gives each donor a single outcome, so the population mean describes the average donor’s proportion. Pooling all cells instead gives more influence to donors contributing more captured cells; that is a different target quantity. Similarly, a mean expression summary within a cell type can be defined donor by donor, then averaged over donors, or calculated by pooling cells. The intended weighting must be explicit.

The contrast scale also matters. If \(\theta_a\) and \(\theta_b\) are mean proportions, \(\theta_b-\theta_a\) is an absolute difference in proportion, often reported in percentage points. A ratio \(\theta_b/\theta_a\) is relative. For log-expression outcomes, a difference of mean log2 values is measured in log2 units; exponentiating that difference gives a ratio of geometric means under the stated log-scale summary. These are not interchangeable descriptions of one number.

A useful specification states the outcome and contrast in one sentence. For example: “Among eligible adults with condition X, estimate the difference between conditions B and A in the population mean donor-level proportion of captured cells assigned to cell type T.” This wording identifies a population quantity and an assay-defined outcome. It does not claim to estimate absolute cell numbers unless the measurement procedure supports that claim.

Worked Example: Mean Donor-Level Proportion Difference

Suppose a hypothetical study defines its target as the mean donor-level proportion of captured cells assigned to a specified cell type among eligible adults. Four donors are observed in each condition. Their proportions are 0.08, 0.10, 0.12, and 0.14 under A, and 0.10, 0.13, 0.15, and 0.18 under B. Treat these values as illustrative data, not results from a real study.

The sample means are \(\bar{x}_A=(0.08+0.10+0.12+0.14)/4=0.11\) and \(\bar{x}_B=(0.10+0.13+0.15+0.18)/4=0.14\). The estimated difference, B minus A, is

$$ \widehat{\Delta}=0.14-0.11=0.03. $$

Thus the estimated difference is 0.03 in proportion, or 3 percentage points. The corresponding population estimand is the difference between the target-population mean donor proportions, not necessarily this observed value. Because the outcome is a proportion of captured cells, the result concerns observed composition under the assay; it is not by itself a difference in absolute cell number. Each donor contributes one proportion to the average, rather than contributing in proportion to the number of cells captured.

Population Weighting and Marginal Contrasts

When a target population contains meaningful subgroups, an overall estimand depends on their weights. For example, the mean under a condition in a population comprising age groups \(s\) can be written as a weighted mean of subgroup means. The weights should describe the target population, not be chosen simply because they match the observed sample.

$$ \theta_c=\sum_s w_s E_P[Y_i(c)\mid S_i=s], \qquad \sum_s w_s=1, \qquad \Delta=\sum_s w_s\left\{E_P[Y_i(b)\mid S_i=s]-E_P[Y_i(a)\mid S_i=s]\right\}. $$

Here \(S_i\) is a subgroup label and \(w_s=P(S_i=s)\) is the target-population proportion in subgroup \(s\). The formula is a direct application of the law of total expectation. It clarifies why the same subgroup-specific results can yield different overall contrasts for different target populations.

Justification: Partition the target population into mutually exclusive subgroups \(S=s\). The population mean under condition \(c\) is the subgroup mean averaged over the probability of belonging to each subgroup: \(E_P[Y_i(c)]=\sum_s P(S_i=s)E_P[Y_i(c)\mid S_i=s]\). Subtracting this expression for conditions \(a\) and \(b\) gives the weighted contrast above. The weights must refer to the intended target population.

Worked Example: A Contrast Standardized to a Target Population

Consider a hypothetical outcome measured as each donor’s cell-type proportion. Suppose the target population is 75% subgroup 1 and 25% subgroup 2. The estimated B-minus-A differences are 0.02 in subgroup 1 and 0.08 in subgroup 2. The target-population standardized contrast is

$$ \widehat{\Delta}_{\mathrm{target}} =0.75(0.02)+0.25(0.08) =0.015+0.020 =0.035. $$

The standardized estimate is 0.035, or 3.5 percentage points. If the observed sample instead contained equal numbers from the two subgroups, its equally weighted subgroup contrast would be \(0.50(0.02)+0.50(0.08)=0.05\), or 5 percentage points. These differ because they target different subgroup mixtures. The calculation assumes the subgroup-specific estimates are relevant to the target subgroups and that both conditions are represented sufficiently to estimate each contrast; it does not establish that the condition contrast is causal.

Scale, Interpretation, and Assumptions

An estimand is defined before a particular estimator or statistical model is chosen. Nevertheless, it should be possible to identify the assumptions needed for the planned estimate to inform that target. At minimum, the target population and outcome must be operationally clear, the observed donors must provide information relevant to that population, and the comparison must be supported by the data. A causal interpretation additionally requires conditions and assignment assumptions that justify interpreting observed contrasts as contrasts between potential outcomes.

For expression outcomes, the summary scale should be stated precisely. Suppose each donor has a normalized log2 expression summary for a chosen cell type. A mean difference of 0.5 log2 units is not an increase of 0.5 expression units on the original scale. If the estimand is a difference in mean log2 expression, exponentiating gives \(2^{0.5}\approx1.414\), a ratio of geometric means under that definition. It is not generally the ratio of arithmetic mean expression values.

Observed single-cell outcomes also depend on measurement. Capture and classification can affect observed cell-type proportions; library normalization and filtering affect expression summaries. The estimand should name what the assay measures, and its interpretation should not silently replace an observed assay quantity with a biological quantity the assay does not identify.

Worked Example: Difference in Mean Log2 Expression

Suppose a hypothetical study defines the outcome as donor-level normalized log2 expression within a specified cell type. Three donors in condition A have values 4.2, 4.5, and 4.3; three in condition B have values 4.7, 4.9, and 4.8.

The sample mean under A is \((4.2+4.5+4.3)/3=13.0/3=4.333\). Under B it is \((4.7+4.9+4.8)/3=14.4/3=4.800\). Therefore the estimated B-minus-A contrast is \(4.800-4.333=0.467\) log2 units. On the geometric-mean ratio scale, this corresponds to \(2^{0.467}\approx1.382\). The estimand here is a difference in population mean donor-level log2 expression, so the ratio interpretation is a ratio of geometric means, not a ratio of arithmetic means on the original expression scale.

Diagnostics, Alternatives, and Common Mistakes

A well-defined estimand makes design and analysis problems easier to recognize. Before fitting a model, check whether the condition labels and outcome summary match the scientific question; whether the target-population weights are stated when subgroup aggregation matters; and whether the observed data include relevant donors under the conditions being contrasted. Large differences between donor-weighted and cell-weighted summaries are a warning that the two analyses answer different questions, not simply a technical disagreement. Very sparse subgroups, absent condition observations in a subgroup, or an outcome definition that changes across conditions can also prevent the intended contrast from being supported.

Common alternatives answer distinct questions. A donor-average contrast gives equal weight to donors; a cell-weighted contrast gives greater weight to donors with more measured cells. A marginal contrast averages over a specified target population; a conditional contrast reports comparisons within specified subgroups. A difference gives an absolute change, while a ratio or log-ratio gives a relative change. Choosing among them is a scientific decision about the target, not a way to select whichever result is most favorable.

A frequent mistake is to write, “The study estimates the effect of condition B,” before specifying the population, outcome, and scale. The correction is to state the complete contrast—for example, the difference in target-population mean donor-level proportion between B and A—and then separately justify any causal language. Another mistake is to treat a large number of cells as if it automatically defines a population of cells as the target. If the intended question concerns the average donor, the donor is the unit being averaged; cell count does not change that estimand.

No estimand can repair a mismatch between the target and the data. A carefully computed estimate may still be imprecise, biased for the target, or based on an assay quantity that only partially represents the biological quantity of interest. Define the target transparently, report the operational outcome and contrast scale, and state the limits of interpretation.

Check Your Understanding

For each question, identify the relevant feature of the estimand. Answers follow the questions.

  1. A study reports a 0.04 difference in the mean donor-level proportion of captured cells in a cell type. What does 0.04 mean on the proportion scale, and what does it not establish?
  2. Two subgroups have condition contrasts of 0.01 and 0.09. Why can an overall contrast not be specified without further information?
  3. Explain the difference between averaging one outcome per donor and pooling all captured cells before calculating a proportion.
  4. A difference in mean log2 expression is 1. What ratio does this imply under a geometric-mean interpretation?

Solutions. (1) It is a 4-percentage-point difference in the specified donor-level captured-cell proportion; by itself it does not establish a difference in absolute cell number or a causal effect. (2) The subgroup weights in the target population are needed to combine the subgroup contrasts into a marginal contrast. (3) The first gives each donor equal weight, whereas pooling cells gives more influence to donors with more captured cells; these are different summaries. (4) The geometric-mean ratio is \(2^1=2\), meaning a twofold ratio on that scale.

Clinical Biostats