Why Recompute an Estimate Without Each Donor?
A population-level result can look precise while depending heavily on one donor. This can happen when donor responses differ, when a donor has an unusual estimate, or when an analysis gives some donors more influence than others. A direct way to examine that dependence is to repeat the analysis after excluding each donor in turn.
This procedure is a sensitivity analysis: it asks how much the estimated result changes under a specified set of omissions. It does not identify a donor as erroneous, prove that a conclusion generalizes, or replace the primary analysis. Its value is diagnostic. It makes dependence on individual observed donors visible and helps readers judge whether the substantive conclusion is stable to their inclusion.
The phrase “repeating the analysis” matters. For a summary defined as the equally weighted mean of donor-level estimates, recalculate that mean over the remaining donors. For a fitted model, remove the donor’s observations and refit the same model with the same outcome, covariates, contrast, and estimation procedure. Simply deleting a donor’s displayed estimate from a table is not enough if the original estimate came from a model whose parameters need refitting.
Define the Target and the Recalculation
Before calculating anything, specify the target estimate. If \(\hat{\theta}_i\) is the estimate for donor \(i\), one common target is the equally weighted mean across donors:
This target gives each donor equal weight, regardless of how many cells were sampled from that donor. It is appropriate when the scientific summary is the average donor-level response. Other analyses may use a different estimator, such as a regression coefficient or a model-based contrast. In those cases, the leave-one-donor-out quantity must be defined by refitting that same estimator after omission. The resulting estimates need not equal an average of the individual donor estimates.
A useful comparison is the change from the full-data estimate:
The sign of \(\Delta_i\) indicates the direction of movement after donor \(i\) is removed; its magnitude gives the change in the estimate’s units. For an estimate of a log fold change, for example, \(\Delta_i\) is also in log-fold-change units. A percentage change can sometimes help describe scale, but is unhelpful when the full estimate is near zero or changes sign. Report the estimates and their units rather than relying on a percentage alone.
A Useful Identity for an Equally Weighted Mean
For the equally weighted mean, the relationship between each donor estimate and its leave-one-out change can be derived directly. Since the sum of all donor estimates is \(D\hat{\theta}\), omitting donor \(i\) gives:
Subtract the full-data mean and simplify:
Thus, for this estimator, omitting a donor whose estimate exceeds the overall mean moves the mean downward, and omitting a donor whose estimate is below the mean moves it upward. The size of the shift grows with the donor’s distance from the overall mean and shrinks as the number of donors increases. This identity is a calculation check, not a general formula for a mixed-model coefficient or another refitted estimator.
Worked Example: Leave-One-Out Changes in a Mean Donor Contrast
Suppose five hypothetical donors have donor-level condition contrasts of \(0.18\), \(0.22\), \(0.25\), \(0.27\), and \(0.58\) log2 units. Assume the intended summary is the equally weighted mean donor contrast. The donors are the independent analysis units, each contrast is defined in the same way, and all five estimates are available.
The full-data mean is:
The leave-one-donor-out estimates are calculated by removing each value in turn:
The largest contrast belongs to donor 5. Excluding that donor reduces the mean from \(0.3000\) to \(0.2300\) log2 units, a shift of \(-0.0700\) log2 units. The other four omissions produce estimates between \(0.3075\) and \(0.3300\). This pattern indicates that donor 5 has comparatively strong influence on the magnitude of the mean. It does not establish that the donor is an outlier or should be removed: the value could reflect genuine biological variation, measurement features, or other factors requiring investigation.
Reading the Pattern, Not Just the Largest Change
A useful review considers the full set of recalculated estimates rather than focusing only on the most extreme donor. Ask whether the scientific interpretation is preserved across omissions: does the estimate keep its direction, remain similar in magnitude, or cross a scientifically important threshold? For a fitted model, also check whether standard errors, intervals, or other reported uncertainty summaries change materially. A stable point estimate does not guarantee stable uncertainty, and a stable sign does not by itself establish a robust biological conclusion.
Interpretation depends on the pattern. If one omission produces a large shift and the other omissions have little effect, the result may be sensitive to that particular donor. If many omissions produce substantial changes, the conclusion may be broadly dependent on the small donor set. If the estimates remain similar, that is evidence of stability to these particular single-donor omissions, not proof of stability to every plausible data or modeling change.
This analysis complements, rather than duplicates, influential-donor diagnostics. An influence measure may summarize how unusual a donor’s contribution is under a particular analysis; leave-one-donor-out estimates show the result of an explicit refit or recalculation. In a donor-level analysis, both perspectives can be informative. Neither licenses removal merely because the result becomes more convenient without a donor.
Worked Example: A Paired Contrast with a Directional Pattern
Consider six hypothetical paired donor contrasts, in a common expression-scale unit: \(-0.10\), \(0.05\), \(0.10\), \(0.15\), \(0.20\), and \(0.80\). Each donor’s contrast is computed comparably, and the target is their equally weighted mean. The sum is \(1.20\), so:
Omitting the donor with contrast \(0.80\) gives:
For comparison, omitting the donor with contrast \(-0.10\) gives \((1.20-(-0.10))/5=0.2600\). Omitting the donor with contrast \(0.05\) gives \(1.15/5=0.2300\); omitting \(0.10\) gives \(1.10/5=0.2200\); omitting \(0.15\) gives \(1.05/5=0.2100\); and omitting \(0.20\) gives \(1.00/5=0.2000\). Thus, the full leave-one-out range is \(0.0800\) to \(0.2600\), and every estimate remains positive.
The direction of the average is stable under these omissions, but its magnitude is sensitive to the \(0.80\) contrast: removing it reduces the mean by \(0.1200\) units. This is a nuanced result. It would be misleading to say either that the estimate is wholly robust or that the positive direction depends entirely on that donor. These recalculations do not provide a new p-value; inference requires an appropriate donor-level procedure and its assumptions.
Assumptions and Practical Diagnostics
Leave-one-donor-out sensitivity analysis is a recalculation, not a test with a universal null distribution. Its interpretation nonetheless relies on a clearly specified and coherent analysis:
- Donor is the omission unit. All observations from a donor are excluded together. Treating cells from the same donor as separate omissions would not assess sensitivity to biological replication.
- The estimate has the same meaning after each omission. Use the same outcome definition, condition contrast, covariates, and estimator. A change in the estimand makes the comparison difficult to interpret.
- The remaining data support the recalculation. Check that each refit converges, relevant groups or conditions remain represented, and the estimate is not undefined. If a model fails after an omission, report that failure rather than quietly replacing the procedure.
- Donor estimates are comparable when averaged. An equally weighted mean requires a meaningful donor-level estimate on a common scale for every included donor. If precision weighting or a fitted model defines the primary estimator, repeat that estimator rather than switching to an unweighted mean for the sensitivity analysis.
A compact diagnostic table should identify the omitted donor, the recalculated estimate, the change from the full estimate, and—when relevant—the corresponding uncertainty summary or model-fit status. Examine the range and direction of changes. Also inspect whether influence aligns with donor-level characteristics or data-quality concerns, while avoiding post hoc exclusions based only on the desired result.
Worked Example: Distinguishing Donor Weighting from Cell Weighting
Suppose three hypothetical donors have estimated responses of \(1\), \(2\), and \(5\) units, based on \(10\), \(10\), and \(80\) sampled cells, respectively. The question is the average response across donors, so each donor—not each cell—is the target unit. The equally weighted donor estimate is:
The leave-one-donor-out estimates are:
Excluding the donor with 80 cells produces the largest shift, from \(2.667\) to \(1.500\) units. That donor is influential for the donor-average estimate because its response is high, not because it has more cells. In contrast, pooling the donor responses with weights proportional to cell counts would yield \((10(1)+10(2)+80(5))/100=4.300\) units. That is a different weighting target and gives the high-cell-count donor much greater weight. More cells can improve measurement of a donor’s response, but they do not automatically make that donor represent more donors.
Limitations and Common Mistakes
Each leave-one-out analysis still uses almost the same data as the full analysis, so the recalculated estimates are not independent replications. They also cannot reveal the effects of donors who were never observed. With few donors, a single omission may leave a very small analysis set, making estimates unstable or a model impossible to fit. These limitations are reasons to report the donor count and all recalculated results, not reasons to suppress them.
A common mistake is to interpret a large leave-one-out shift as proof that the omitted donor is bad data. The reasoning is wrong because influence describes how the estimate responds to an observation; it does not diagnose the observation’s validity. The correction is to check data quality and the scientific context, report the sensitivity, and retain the donor unless there is a defensible, documented reason to exclude it.
Another mistake is to call a conclusion “robust” merely because every recalculated estimate has the same sign. Sign stability is useful, but the magnitude or uncertainty may still change substantially. Conversely, a small change in the point estimate does not ensure stable uncertainty. State separately what happened to the estimate, its uncertainty, and the substantive interpretation.
Check Your Understanding
Answer each question, then compare with the solutions.
- Four hypothetical donor estimates are \(1\), \(2\), \(3\), and \(6\) units. Find the equally weighted full mean and each leave-one-donor-out mean.
- For an equally weighted mean across \(D\) donors, explain why omitting a donor whose estimate exceeds the full mean lowers the leave-one-out mean.
- A sensitivity analysis gives estimates \(0.12\), \(0.11\), \(0.10\), and \(-0.02\) after four different omissions. What can be concluded about sign stability, and what cannot be concluded?
- A donor has many more sampled cells than the other donors. Does that fact alone mean the donor should receive more weight in an equally weighted donor-average estimand?
Solutions. (1) The full mean is \((1+2+3+6)/4=3.000\) units. Omitting \(1\) gives \((2+3+6)/3=3.667\); omitting \(2\) gives \((1+3+6)/3=3.333\); omitting \(3\) gives \((1+2+6)/3=3.000\); omitting \(6\) gives \((1+2+3)/3=2.000\) units. (2) The identity \(\Delta_i=(\hat{\theta}-\hat{\theta}_i)/(D-1)\) is negative when \(\hat{\theta}_i>\hat{\theta}\), so omission lowers the mean. (3) The estimates do not all share a sign, so sign stability fails under these omissions. The pattern alone does not establish which donor is invalid, whether a population effect exists, or whether a formal inferential conclusion holds. (4) No. If the target is an equally weighted mean across donors, each donor receives equal weight; cell counts affect the amount of within-donor information, not the target’s donor weighting.