Why Support Gaps Need Their Own Diagnostic
A generated virtual population can have a plausible mean, variance, and overall range while still omitting a meaningful part of the target population. For example, its minimum and maximum values may bracket the target values, yet almost no virtual patients may lie in a clinically relevant middle range. Such an omission can matter when model predictions depend on that region or when a planned analysis concerns patients with those characteristics.
A support-gap diagnostic asks a focused question: does the candidate or selected population contain patients in target regions that should be represented? It complements distribution summaries and distance measures rather than replacing them. The word “support” requires care: a finite sample cannot establish every point where an underlying continuous population has positive probability. In practice, gaps are assessed over specified intervals or regions at a resolution chosen for the scientific purpose.
The diagnostic depends on a clear target definition. A target can come from a reference sample, population summaries translated into intervals, or explicit eligibility and physiological constraints. Reference and virtual-population values must use comparable definitions, units, transformations, and inclusion rules. If a reference sample is used to represent a broader target population, it should be sufficiently representative for the intended comparison. Without that condition, the plots describe the samples but may not describe the populations of interest.
Empirical Cumulative Distributions
For a single characteristic, the empirical cumulative distribution function (ECDF) records the fraction of observations at or below each value. Let \(T_1,\ldots,T_{n_T}\) denote values in the target reference and \(V_1,\ldots,V_{n_V}\) values in the virtual population. Then:
Here \(I(\cdot)\) is 1 when its condition is true and 0 otherwise. For an interval \((a,b]\), the empirical fraction in the interval is the ECDF increment \(F(b)-F(a)\). Thus, if the target ECDF rises across an interval but the virtual-population ECDF is flat there, that interval contains target observations and no virtual observations. The interval boundaries should be stated explicitly because ECDFs count observations at or below a boundary.
Worked Example: Finding a Missing Middle Interval With ECDFs
Suppose a hypothetical target reference has values \(2,4,8,12,14,18,22,24,28,34\) units, and a candidate virtual population has values \(2,4,8,22,24,28,34,36\) units. The intended admissible range is 0 to 40 units. Check coverage of the target interval \((10,20]\).
The target has three of ten observations in the interval, so its empirical fraction is \(3/10=0.300\). The candidate has no observations there, so its fraction is \(0/8=0\). This can also be read from the ECDFs:
The candidate spans the target’s overall range: both have a minimum of 2 units, and their maxima are 34 and 36 units. Both also lie within the admissible range. Those range checks do not reveal the missing middle interval. The ECDF increments do: at the chosen resolution, \((10,20]\) is an empirical support gap in the candidate. This is a descriptive finding, not evidence that every underlying target population must have exactly 30.0% in that interval.
Density Plots and the Resolution of a Gap
Histograms and density plots can make gaps easier to see across a wider range of values. A histogram displays counts or proportions in chosen bins. A density plot, commonly made by smoothing observations, displays an estimated distribution over the value scale. Both require judgment about resolution. Narrow bins can produce many empty bins in a modest sample, while wide bins can merge a meaningful gap with covered values. Kernel smoothing can draw a continuous-looking curve across a region with no observations; that smooth curve is not itself evidence that virtual patients occupy the region.
A practical check is to plot target and virtual-population values on the same scale, state the bin widths or smoothing bandwidth, and mark important target intervals. Inspect the observed values or bin counts alongside the curve. If a density curve appears to bridge a proposed gap, verify whether any actual observations lie there. Changing bandwidth or bin width is a useful sensitivity check, but it should not replace a scientifically justified region definition.
Worked Example: Checking a Density Plot Against Observed Counts
In a hypothetical comparison, 20 target-reference measurements have five observations in each of four bins: \([0,10)\), \([10,20)\), \([20,30)\), and \([30,40]\) units. Ten candidate measurements have five in the first bin, none in either middle bin, and five in the last bin. Assess coverage of the target region \([10,30)\).
The target fraction in the region is \((5+5)/20=0.500\). The candidate fraction is \((0+0)/10=0\). A histogram using these bins displays two empty candidate bins where the target has 10 observations. The regional proportions confirm that the difference is not just a visual impression.
A smoothed candidate density may nevertheless show a low curve between the clusters near the two ends. Smoothing assigns influence beyond each observed value; it does not create candidate observations in the middle bins. If a different bin width is chosen, the displayed pattern may look less stark, so the bin definition and the direct interval count should be reported together. At the specified resolution, the candidate has a support gap in \([10,30)\).
Range Checks and Multidimensional Regions
A range check compares observed minima and maxima with admissible bounds or target ranges. It is simple and often useful for detecting values outside physiological or design constraints. However, checking only endpoints cannot show how observations are distributed between them. A population can pass a range check and still omit a large interior region. Conversely, a single observation in an interval establishes sample presence, but not adequate representation.
For multiple patient characteristics, checking each marginal range separately is not enough to show that their combinations are covered. A target region may require a particular combination of values. Define that region directly and count patients that satisfy all its conditions. This is a localized diagnostic; it does not by itself establish that the full joint distribution is well represented.
Worked Example: Marginal Ranges Pass but a Joint Region Is Missing
Suppose a hypothetical target reference has eight patients, including two in a target region defined by \(8\leq A\leq12\) units and \(4\leq B\leq6\) units. Its overall observed ranges are \(5\) to \(15\) units for \(A\), and \(2\) to \(8\) units for \(B\). A four-patient candidate has values \((A,B)=(9,2),(11,8),(5,5),(15,5)\). Check the joint target region.
The candidate minimum and maximum for \(A\) are 5 and 15, and for \(B\) they are 2 and 8, matching the target’s overall observed ranges. Yet neither of the first two candidate patients meets the \(B\) condition, while neither of the last two meets the \(A\) condition. Therefore, no candidate patient satisfies both conditions. The target fraction is \(2/8=0.250\), and the candidate fraction is \(0/4=0\).
The separate marginal ranges pass, but the joint region is empty in the candidate. The result identifies a region-specific coverage concern; it does not diagnose every aspect of dependence between \(A\) and \(B\). The region definition and units make the finding interpretable and reproducible.
A Practical Diagnostic Sequence
A defensible support-gap assessment follows a sequence that makes the target, resolution, and evidence explicit. The checks below can be applied to candidate populations before selection and to selected populations after filtering or resampling.
Use scientific or design rationale to define intervals or combinations of traits, including boundaries and units. Do not choose regions solely because the plotted candidate happens to have observations there.
Confirm that reference and virtual values use the same variable definitions, transformations, units, eligibility criteria, and missing-data rules.
Compare ECDFs and, where useful, histograms or density plots. Calculate target and virtual counts or fractions in each important region; do not rely on visual smoothness alone.
Compare minima and maxima with admissible limits, then evaluate any specified multidimensional regions using joint membership rather than separate marginal ranges.
State region boundaries, sample sizes, observed counts, and whether the result is a complete empirical gap or a lower-than-target representation. Note sensitivity to binning or smoothing when relevant.
Comparisons, Failure Modes, and Limitations
These tools answer related but different questions. An ECDF shows cumulative distribution differences and can localize intervals where one population accumulates observations and another does not. A histogram or density plot gives a visual summary whose apparent resolution depends on bins or smoothing. A range check is a fast boundary check but is weak evidence about interior coverage. A global distance metric summarizes distributional discrepancy in a chosen way, but a single overall value may not identify which target region is empty. Combining regional counts with plots and bounds is often more informative than treating any one diagnostic as sufficient.
- Choosing regions after seeing the result. This can make an ordinary fluctuation look like a meaningful gap. Define important regions from the target question where possible, and label exploratory regions as exploratory.
- Treating a smoothed curve as observed coverage. A density estimate can be positive where there are no observations. Check the raw values or bin counts before concluding that a region is represented.
- Equating one observed patient with adequate coverage. Presence is not the same as adequate representation. Report counts and fractions, and interpret them against the intended use rather than declaring sufficiency from a nonzero count alone.
- Inferring a population-level zero from a sample gap. A finite candidate may omit an interval by chance. The diagnostic establishes what is absent from that generated sample, not the exact probability in an unlimited population.
- Checking only marginal ranges. Plausible values for each trait do not guarantee that target combinations occur together. Count membership in relevant joint regions.
Small samples make observed gaps especially sensitive to sample size and the chosen resolution. Repeating generation can show whether a gap recurs, but a recurring gap may still reflect the generation or selection procedure rather than a property of the target. These plots and counts are diagnostics, not automatic acceptance rules. Their role is to reveal where representation deserves investigation and to make the evidence transparent.
Check Your Understanding
For each question, identify the diagnostic and interpret the evidence at the stated resolution.
- A target sample has 4 of 20 observations in \((30,40]\), while a virtual population has 0 of 25. What are the two empirical fractions, and what does their comparison show?
- A candidate ECDF is flat from 5 to 10 units, but a kernel density plot is slightly above zero there. What should be checked before claiming the candidate covers that interval?
- A candidate minimum and maximum fall within the target’s overall range. Does that establish coverage of every target interval? Explain.
- Each of two traits has values within its target marginal range, but no candidate patient satisfies a prespecified combination of both traits. What kind of gap is indicated?
Solutions. (1) The target fraction is \(4/20=0.200\), and the virtual fraction is \(0/25=0\). There is an empirical support gap in the stated interval for this virtual sample. (2) Check the observed values or interval count; smoothing can produce a positive curve where no observations occur. (3) No. The endpoints do not reveal interior coverage, so the interval counts or ECDF increments must also be examined. (4) The candidate has an empirical gap in the specified joint region, even though its separate marginal ranges may be acceptable.