Tutorials › Statistical Methods for Single-Cell RNA Sequencing › Diagnosing Imbalance Across Donors

Biological Replication and Donor-Level Inference · Tutorial 461 of 700

Diagnosing Imbalance Across Donors

You will learn to summarize donor-level count imbalance and assess how unequal cell sampling relates to outcome summaries and their precision.

Advanced 10 min read

What You'll Learn

  • Summarize donor cell-count imbalance with ratios, shares, and concentration measures.
  • Compare equally weighted donor summaries with summaries weighted by cell counts.
  • Distinguish the number of cells measured from the number of biological replicates.
  • Use within-donor correlation to assess how counts affect donor-summary precision.
  • Recognize diagnostic warnings without treating them as proof of bias or donor influence.

Why Diagnose Donor Imbalance?

A single-cell data set may contain many cells but very different numbers of cells from different donors. That imbalance matters because a donor with few cells may have a less precisely measured outcome, while a donor with many cells may contribute a large share of a cell-weighted summary. Neither fact means that the donor is less biologically important or that the donor should automatically be excluded.

The aim here is a disciplined diagnostic: describe how many cells contribute for each donor, summarize the donor-level outcomes, and check whether the count distribution is associated with the outcome summaries or the conclusions one might draw from them. Counts are not independent biological replicates. For comparisons about donors or populations, the donor remains the biological replication unit.

Definition: Let \(J\) be the number of donors and \(n_j\) the number of eligible cells measured for donor \(j\), for \(j=1,\ldots,J\). Let \(\bar{y}_j\) be a donor-level summary of the outcome, such as a mean expression value or a proportion, calculated from that donor's eligible cells. Donor-count imbalance means that the \(n_j\) values differ substantially across donors. Outcome imbalance describes differences in the \(\bar{y}_j\) values; it is not the same thing as count imbalance.

Counts must refer to a comparable set of observations. For example, if the question concerns one cell type, use counts for that cell type rather than each donor's total across all cell types. For a condition comparison, record counts separately by donor and condition. A count can be small because of biology, sampling, filtering, or technical quality; the count alone does not identify the cause.

Summarizing the Count Distribution

Start with a table or plot of \(n_j\) by donor. Useful descriptive summaries include the minimum and maximum count, the maximum-to-minimum ratio, each donor's share \(w_j=n_j/\sum_{k=1}^{J}n_k\), and the coefficient of variation \(s_n/\bar{n}\), where \(\bar{n}\) and \(s_n\) are the sample mean and sample standard deviation of the donor counts. These summaries describe different features: a ratio highlights the range, shares show how concentrated the total is, and the coefficient of variation describes spread relative to the mean.

A concentration summary can make the count distribution easier to compare across data sets. Define \(N_{\mathrm{conc}}=1/\sum_j w_j^2\). This is the number of equally contributing donors that would produce the same sum of squared shares. It is a descriptive equivalent number, not the number of biological replicates and not a replacement for \(J\).

Worked Example: Finding Count Concentration

Suppose five hypothetical donors contribute 40, 50, 60, 70, and 280 eligible cells. The total is \(500\), and the mean count is \(500/5=100\) cells per donor. The maximum-to-minimum ratio is \(280/40=7\). The largest donor contributes \(280/500=0.56\), or 56% of all eligible cells.

The sample standard deviation of the counts is \(s_n=\sqrt{[(40-100)^2+(50-100)^2+(60-100)^2+(70-100)^2+(280-100)^2]/4}\) \(=\sqrt{41000/4}=\sqrt{10250}=101.2\) cells. Thus the coefficient of variation is \(101.2/100=1.012\).

The donor shares are \(0.08, 0.10, 0.12, 0.14,\) and \(0.56\). Therefore \(N_{\mathrm{conc}}=1/(0.08^2+0.10^2+0.12^2+0.14^2+0.56^2)=1/0.364=2.747\). The data contain five donors, but their cell-count concentration is equivalent to about 2.75 equally contributing donors for this particular descriptive measure. That does not mean there are only 2.75 biological replicates; all five donors remain donors.

No single threshold defines unacceptable imbalance in every study. A maximum-to-minimum ratio can look dramatic when the minimum is very small, while a moderate ratio can conceal one donor contributing a large share when there are many donors. Report the distribution and relevant summaries rather than relying on one cutoff.

Compare Donor Summaries With Cell-Weighted Summaries

An equally weighted donor summary gives each donor the same contribution: \(\bar{y}_{D}=J^{-1}\sum_{j=1}^{J}\bar{y}_j\). A cell-count-weighted summary instead gives donor \(j\) weight \(w_j=n_j/\sum_k n_k\): \(\bar{y}_{C}=\sum_j w_j\bar{y}_j\). These answer different descriptive questions. The first is the average of the observed donor summaries, giving each donor equal weight. The second is the average over the observed cells, so donors with more eligible cells contribute more. The appropriate target depends on the scientific question; neither summary is automatically the correct inferential estimator.

$$ \bar{y}_{C}-\bar{y}_{D} =\sum_{j=1}^{J}\left(w_j-\frac{1}{J}\right)(\bar{y}_j-\bar{y}_{D}) $$

To see this identity, write \(\bar{y}_{C}=\sum_j w_j(\bar{y}_j-\bar{y}_{D})+\bar{y}_{D}\sum_jw_j\). Since the weights sum to one, subtracting \(\bar{y}_{D}\) gives the displayed expression. Thus the difference between the summaries depends on whether donors with above-average counts tend to have above- or below-average outcomes. If all counts are equal, every \(w_j=1/J\), and the summaries agree. A large difference is a useful warning that the weighting choice is consequential; by itself, it does not establish bias or causation.

Worked Example: How Counts Change a Summary

Four hypothetical donors have cell counts \(20, 40, 80,\) and \(160\), and donor-level outcome means \(2, 3, 3,\) and \(8\) score units. The equally weighted donor mean is \((2+3+3+8)/4=4.000\) score units.

The total count is \(300\). The cell-weighted mean is \((20\times2+40\times3+80\times3+160\times8)/300\) \(=(40+120+240+1280)/300=5.600\) score units.

The donor with the largest outcome also has the largest count, so cell weighting raises the summary by \(5.600-4.000=1.600\) score units. The calculation does not show that the cell-weighted summary is wrong. It shows that it describes a different weighting of these observed donors and that the result is sensitive to whether donors or cells receive equal weight.

Counts and the Precision of Donor Outcomes

A low count can make a donor's cell-based outcome summary noisy, but the amount of additional information from more cells depends on dependence among cells from the same donor. For a simple diagnostic, suppose the cell outcomes within a donor have common marginal variance \(\sigma^2\) and common pairwise correlation \(\rho\), with \(0\leq\rho\leq1\). Under this exchangeable-correlation model, the variance of a donor's cell mean based on \(n_j\) cells is \(\sigma^2[\rho+(1-\rho)/n_j]\). This calculation assumes the cells are sampled comparably within donor and that the variance and correlation adequately describe their dependence. It is a precision model for a donor summary, not a model for between-donor variation.

The expression also explains why raw cell counts should not be read as independent sample sizes. When \(\rho>0\), the variance approaches \(\sigma^2\rho\), rather than zero, as \(n_j\) grows. More cells can improve measurement of a donor, but they do not create more donors. The earlier material on intraclass correlation and cell sampling uncertainty develops these dependence and precision ideas further.

Worked Example: Comparing Donor-Mean Precision

Suppose, for illustration, that cell outcomes have marginal variance \(\sigma^2=4\) score-units squared and within-donor correlation \(\rho=0.20\). Compare donors measured on 10 and 100 cells. The model gives \(\operatorname{Var}(\bar{Y}_j)=4[0.20+0.80/n_j]\).

For 10 cells, the variance is \(4[0.20+0.80/10]=4(0.28)=1.120\) score-units squared, giving \(SE=\sqrt{1.120}=1.058\) score units. For 100 cells, the variance is \(4[0.20+0.80/100]=4(0.208)=0.8320\), giving \(SE=\sqrt{0.8320}=0.9121\) score units.

The donor with 100 cells has a smaller model-based standard error, but not one-tenth the standard error. This comparison presumes the same \(\sigma^2\) and \(\rho\) for both donors and treats them as known for illustration. In practice, uncertainty in these quantities and differences in donor biology also matter. The calculation is a prompt to assess precision, not a justification for automatically weighting donors by their cell counts.

A Practical Diagnostic Sequence

Use a consistent set of summaries before fitting or interpreting donor-level comparisons. The following sequence keeps count imbalance, outcome variation, and inferential choices distinct.

1
Define the counting unit.
State the cell type, condition, and inclusion rules that define an eligible cell. Record counts by donor for that exact set.
2
Inspect the count distribution.
Report the number of donors, the range, donor shares, and a spread measure. Plot counts by donor when possible; a sorted dot plot makes concentration visible.
3
Pair counts with outcomes.
Show each donor's outcome summary alongside its count. If available, include a measure of within-donor sampling uncertainty, while making clear that it does not represent uncertainty across donors.
4
Check weighting sensitivity.
Compare an equally weighted donor summary with a cell-weighted descriptive summary when both are meaningful. Explain what population or collection of observations each represents.
5
Investigate the source and report limits.
Check whether count differences align with sampling, quality filters, cell availability, or study design. Counts alone cannot identify which explanation is responsible.

A table can be particularly useful when outcomes are condition-specific: include one row per donor-condition combination, with its cell count and outcome summary. This reveals whether an apparently balanced total count hides a sparse donor in one condition. In paired designs, both conditions need usable information for the donor-level contrast; a large total across conditions does not guarantee adequate measurement in each one.

Comparisons, Failure Modes, and Limits

Diagnostic or summaryWhat it describesWhat it does not establish
Count range, shares, or concentrationHow cell contributions are distributed among observed donorsWhether a donor is biologically atypical or should be excluded
Equally weighted donor summaryThe average of the observed donor-level summariesThe average across all observed cells
Cell-weighted summaryThe outcome averaged over the observed cellsAn equally weighted donor-level population effect
Model-based within-donor standard errorSampling precision under its stated within-donor assumptionsBetween-donor uncertainty or the number of independent donors

A common mistake is to see one donor with many more cells and conclude that the analysis has many more independent replicates for that donor. The correction is to distinguish cell-level measurement from donor-level replication: extra cells may refine a donor summary, but they do not add independent donors. Another mistake is to equal-weight donors in every setting without asking what the target is. Equal weighting is informative for a donor-average description, while cell weighting describes the observed-cell average; the scientific estimand and analysis model must guide the choice.

Other warning signs include comparing counts from different cell populations, hiding a few low-count donors behind a large total, or interpreting a count-outcome association as evidence that count imbalance caused the outcome pattern. The diagnostic summaries are descriptive. They do not resolve whether count differences reflect biological composition, technical sampling, or design, and they do not replace donor-aware uncertainty estimation. Nor should a large count or an extreme outcome alone be treated as a formal finding of donor influence; that requires a separate influence assessment.

The practical conclusion should be calibrated: state how uneven the counts are, whether the donor summaries appear to depend on that unevenness, and which assumptions limit interpretation. When imbalance is substantial, report it transparently and show how the chosen donor-aware analysis handles unequal measurement. Do not claim that a larger count automatically deserves greater biological weight.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. Three donors contribute 10, 30, and 60 cells. What share of cells comes from the largest donor, and what is the maximum-to-minimum count ratio?
  2. Two donors have outcome means of 2 and 6 units and counts of 90 and 10. Calculate the equally weighted donor mean and the cell-weighted mean. What does their difference show?
  3. Why does a donor with ten times as many cells not necessarily provide ten times as much independent biological information?
  4. Name two diagnostic summaries that describe count imbalance, and state one limitation of using a single summary.

Solutions. (1) The largest donor's share is \(60/(10+30+60)=0.600\), or 60%; the maximum-to-minimum ratio is \(60/10=6\). (2) The equally weighted mean is \((2+6)/2=4.000\) units. The cell-weighted mean is \((90\times2+10\times6)/100=240/100=2.400\) units. The gap shows that the observed-cell summary gives much greater weight to the low-outcome donor; it does not by itself determine which summary answers the scientific question. (3) Cells from one donor can be correlated, so additional cells refine measurement but do not create additional independent donors. (4) Examples include the maximum-to-minimum ratio, donor shares, the coefficient of variation, and the concentration equivalent number. Each captures a different feature, so one statistic can miss an important aspect of the distribution.

Clinical Biostats