Tutorials › Quantitative Systems Pharmacology (QSP) › From Observed Summaries To Population Targets

Virtual Patient Generation · Tutorial 689 of 799

From Observed Summaries To Population Targets

Translate observed population summaries into clearly defined targets and uncertainty ranges that can be documented for virtual population generation.

Advanced 12 min read

What You'll Learn

  • Specify the population quantity, estimate, units, and sample size behind a reported summary.
  • Distinguish an estimated target from its confidence interval and from an operational tolerance band.
  • Calculate and interpret uncertainty ranges for means, variances, proportions, and medians.
  • Use sample size and stated assumptions to judge how precisely a summary informs a target.
  • Identify when a reported summary alone is insufficient to define a defensible target range.

From a Reported Number to a Population Target

A published summary is not yet a complete specification for a virtual population. “Mean clearance: 42” does not say whether 42 is a sample mean or a known population value, what units apply, how many observations contributed, or how uncertain the estimate is. Those distinctions matter when the summary is used to guide a population of mechanistic parameter vectors: an unsupported narrow band can exclude plausible patients, while a very wide band may provide little constraint.

The task here is to translate reported summaries into explicit, traceable targets. It is useful to keep three items separate: the population quantity of interest, its estimate from observed data, and the tolerance or uncertainty range used in a decision. A confidence interval describes statistical uncertainty under stated assumptions. An operational tolerance band is a practical criterion chosen for a particular purpose; it may be informed by a confidence interval, but it is not automatically the same thing.

Definition: A population target is a specified population quantity to which a virtual population is intended to correspond, such as a mean \(\mu\), variance \(\sigma^2\), quantile, or prevalence \(p\). A target specification should state the quantity, population or subgroup, units, estimate, source sample size, and any uncertainty or operational tolerance range. The estimate is not itself proof that the corresponding population quantity has that exact value.

A reported \(n\) is the number of observations contributing to a summary, not necessarily the number of virtual patients to generate. It helps determine the uncertainty of the estimate, but does not by itself define the target. The target must also match the intended population: age range, disease status, treatment, measurement time, and other inclusion conditions can change the relevant summary.

Means and Their Uncertainty

Suppose \(x_1,\ldots,x_n\) are independent observations from a population with mean \(\mu\). Their sample mean is \(\bar{x}\), and their sample standard deviation is \(s\). If the population is normal, or if the sample is sufficiently large for the sampling distribution of \(\bar{x}\) to be adequately approximated by a normal distribution, a two-sided \(100(1-\alpha)\%\) confidence interval for \(\mu\) is

$$ \bar{x}\ \pm\ t_{1-\alpha/2,\ n-1}\frac{s}{\sqrt{n}}. $$

Here \(t_{1-\alpha/2,\ n-1}\) is the corresponding quantile of the \(t\) distribution with \(n-1\) degrees of freedom. The standard error \(s/\sqrt{n}\) quantifies the estimated sampling variability of the sample mean. This interval does not describe the spread of individual patients; that spread is represented by quantities such as \(s\), not by the standard error.

Worked Example: A Mean Clearance Target

Suppose a hypothetical report gives \(n=16\) patients, a mean clearance \(\bar{x}=42.0\) L/h, and a sample standard deviation \(s=8.0\) L/h. The observations are independent and are treated as approximately normally distributed, with no strong outliers or severe skew that would make a small-sample \(t\) interval unreliable.

The point estimate for the population mean is \(42.0\) L/h. For a 95% confidence interval, \(df=16-1=15\), and \(t_{0.975,15}\approx2.131\). The estimated standard error is

$$ SE=\frac{8.0}{\sqrt{16}}=2.0\ \mathrm{L/h}. $$

The margin of error is \(2.131(2.0)=4.262\) L/h, so the interval is

$$ 42.0\pm4.262=[37.74,\ 46.26]\ \mathrm{L/h}. $$

A defensible specification records the estimated population-mean target as \(42.0\) L/h and the 95% confidence interval as approximately \(37.74\) to \(46.26\) L/h. Under the stated sampling assumptions, this procedure has 95% long-run coverage for \(\mu\). The interval is not a range containing 95% of individual patients, nor does it mean there is a 95% probability that this fixed \(\mu\) lies in this particular interval.

For a virtual population, the confidence interval can inform how much precision the observed mean supports. A separate operational rule might require the generated population mean to fall within a prespecified band. That band should be justified by the application, measurement precision, and intended population—not silently substituted for the confidence interval.

Variances: A Distributional Assumption Matters

A variance target describes population dispersion, in squared units. For a random sample from a normal population, \((n-1)s^2/\sigma^2\) has a chi-square distribution with \(n-1\) degrees of freedom. Inverting this result gives the following \(100(1-\alpha)\%\) confidence interval for \(\sigma^2\):

$$ \left[ \frac{(n-1)s^2}{\chi^2_{1-\alpha/2,\ n-1}}, \frac{(n-1)s^2}{\chi^2_{\alpha/2,\ n-1}} \right]. $$

The lower endpoint uses the larger chi-square quantile. This interval depends strongly on the normality assumption, particularly for small samples. If the data are markedly skewed or contain influential outliers, the interval may not have its stated coverage.

Worked Example: A Variance Target in Squared Units

In a hypothetical sample of \(n=25\) measurements, suppose \(s^2=144\ (\mathrm{mg/L})^2\). Assume independent observations from a normal population. For a 95% interval, \(df=24\), with chi-square quantiles \(\chi^2_{0.975,24}=39.364\) and \(\chi^2_{0.025,24}=12.401\).

Substitution gives

$$ \left[ \frac{24(144)}{39.364}, \frac{24(144)}{12.401} \right] = [87.79596\ldots,\ 278.6872\ldots]\ (\mathrm{mg/L})^2 \approx[87.80,\ 278.69]\ (\mathrm{mg/L})^2. $$

The point target is \(s^2=144\ (\mathrm{mg/L})^2\), while the confidence interval expresses uncertainty about \(\sigma^2\). Relative to the point estimate, the confidence-based departures are approximately \(-56.20\) and \(+134.69\ (\mathrm{mg/L})^2\): \(87.79596\ldots-144=-56.20404\ldots\), and \(278.6872\ldots-144=134.6872\ldots\). The asymmetry is a feature of the variance interval, not an arithmetic error. A range for individual measurements would instead be expressed in mg/L and requires a different calculation.

The calculation illustrates why a target should include its scale and units. A variance in \((\mathrm{mg/L})^2\) cannot be compared directly with a standard deviation in mg/L. If a report gives only a standard deviation, square it to obtain a variance target; if uncertainty for that target is needed, the sampling assumptions and sample size are also required.

Proportions and Sample Size

For a binary characteristic, let \(p\) be the population proportion and \(\hat{p}=x/n\) the observed proportion, where \(x\) of \(n\) observations have the characteristic. When the sample is independent and both expected counts \(n\hat{p}\) and \(n(1-\hat{p})\) are sufficiently large, a normal-approximation interval is

$$ \hat{p}\ \pm\ z_{1-\alpha/2} \sqrt{\frac{\hat{p}(1-\hat{p})}{n}}. $$

This approximation is most defensible when neither category is rare. It can produce endpoints outside zero to one in small or highly imbalanced samples; in those settings, use an interval designed for proportions rather than treating impossible endpoints as meaningful. The sample size \(n\) is essential: the same observed percentage can imply quite different precision in a small and a large sample.

Worked Example: Prevalence and the Information in n

In a hypothetical sample of \(n=200\) patients, 30 have a specified phenotype. Thus \(\hat{p}=30/200=0.1500\). The expected counts under the fitted proportion are 30 and 170, so the normal approximation is reasonable for this illustration.

The estimated standard error is

$$ SE=\sqrt{\frac{0.15(1-0.15)}{200}} =\sqrt{0.0006375} \approx0.02525. $$

A 95% approximate interval is \(0.1500\pm1.96(0.02525)=0.1500\pm0.04949\), or approximately \([0.1005,\ 0.1995]\). The observed prevalence target is 15.0%; the interval corresponds to about 10.1% to 20.0%. It describes uncertainty about the population prevalence under the stated assumptions, not the range of phenotype values for individuals.

If this prevalence were used to plan a new sample with a desired approximate 95% margin of error \(E=0.05\), using 0.15 as the planning value gives

$$ n\geq\frac{1.96^2(0.15)(0.85)}{0.05^2} =\frac{0.489804}{0.0025} =195.9216. $$

Rounding up gives \(n=196\). This is a planning calculation conditional on the anticipated proportion and normal approximation; it does not alter the original sample size of 200 or make that sample a virtual-population size requirement.

Quantiles Require Rank Information

A quantile target identifies a location in a distribution, such as the median or the 90th percentile. A reported quantile alone does not usually provide enough information to quantify its sampling uncertainty. That uncertainty depends on sample size and on the data near the quantile. With raw observations, a distribution-free interval for a population median can be obtained from order statistics, provided the observations are independent and drawn from a continuous population with a well-defined median.

For a sample of size \(n\), let \(X_{(1)}\leq\cdots\leq X_{(n)}\) be the sorted observations. If the population median is \(m\), the number \(B\) of observations below \(m\) follows a binomial distribution with \(n\) trials and success probability \(1/2\). Consequently, the interval \([X_{(k)},X_{(n-k+1)}]\) has coverage \(P(k\leq B\leq n-k)\). This gives a direct justification for rank-based uncertainty: the interval covers the median when the count on either side places it between the selected order statistics.

Worked Example: A Median Target With Rank-Based Uncertainty

Suppose 20 hypothetical patient measurements, in hours, are sorted, with \(X_{(6)}=5.4\), \(X_{(10)}=6.1\), \(X_{(11)}=6.3\), and \(X_{(15)}=7.1\). The sample median is \((6.1+6.3)/2=6.2\) hours. Assume independent observations from a continuous population.

Using \(k=6\), the interval for the population median is \([X_{(6)},X_{(15)}]=[5.4,7.1]\) hours. Its exact coverage is \(P(6\leq B\leq14)\) for \(B\sim\mathrm{Binomial}(20,0.5)\). By symmetry,

$$ P(6\leq B\leq14) =1-2P(B\leq5) =1-2\left(\frac{\sum_{j=0}^{5}{20\choose j}}{2^{20}}\right) =1-2\left(\frac{21700}{1048576}\right) \approx0.9586. $$

Thus the rank interval has approximately 95.86% coverage for the population median under the assumptions. The point estimate is 6.2 hours, while the interval communicates the substantial uncertainty associated with only 20 observations. Without the order statistics or raw data, the reported median alone would not justify this interval.

For other quantiles, the same principle uses the probability mass below the target quantile rather than one half. In practice, a report should provide a confidence interval or enough underlying data to construct one. A quantile point estimate should not be given an arbitrary symmetric tolerance and presented as statistically derived.

Building a Defensible Specification

A concise target record can include: the population or subgroup; the quantity and units; the estimate and sample size; the assumptions or method used for uncertainty; the interval, if available; and any separate operational acceptance band. This makes clear whether the range reflects sampling uncertainty, a scientific constraint, or a decision rule.

A confidence interval and an operational tolerance answer different questions. The interval asks how much uncertainty the observed sample leaves about a population quantity under a model. A tolerance band states what discrepancy is acceptable for a particular use. For example, a project may choose an allowed deviation from a mean because of measurement precision or a downstream decision threshold. The choice should be documented and justified; it is not licensed solely by the confidence level.

1
Identify the estimand.
Name the population quantity, subgroup, measurement time, and units; do not substitute an individual-level range for a population summary.
2
Record the evidence.
Capture the reported estimate, sample size, and any spread or raw-data information needed to assess uncertainty.
3
Quantify uncertainty under stated assumptions.
Use a method appropriate to the summary and check independence, distributional assumptions, and approximation conditions.
4
State the decision tolerance separately.
If a practical acceptance band is needed, specify its rationale and do not label it a confidence interval unless it was constructed as one.

A common mistake is to treat a reported sample mean as the exact population target and then choose a narrow band around it without considering \(n\) or variability. The correction is to separate the point estimate from uncertainty: for a mean, for example, the estimated standard error is \(s/\sqrt{n}\), so precision depends on both spread and sample size. Another mistake is to read a 95% confidence interval as containing 95% of patients. It concerns uncertainty in a population summary; patient-level variation is a different feature of the distribution.

Summaries also have limits. Marginal means, variances, proportions, and quantiles do not determine the joint distribution of patient traits, nor do they guarantee physiological feasibility. A reported mean may conceal skewness or subgroups, and a sample may not represent the intended target population. When essential information is absent—such as sample size, units, or the definition of a subgroup—the honest specification records that limitation rather than manufacturing precision.

Key takeaway: Every population target should distinguish the quantity of interest, its estimate, the evidence supporting its precision, and any separately chosen tolerance. A number without population definition, units, sample size, and uncertainty context is an incomplete target.

Check Your Understanding

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

  1. A study reports \(n=36\), \(\bar{x}=10\) mg/L, and \(s=3\) mg/L. What is the estimated standard error of the mean, and does it describe individual variability?
  2. A sample prevalence is 20% with \(n=50\), while another sample also has prevalence 20% with \(n=500\). Which estimate is generally more precise, and why?
  3. Why can a sample median alone not generally determine a defensible confidence interval for the population median?
  4. Is a project-defined acceptance band automatically a confidence interval? Explain the difference.
  5. A variance interval is reported in \((\mathrm{mg/L})^2\). What is the principal units error in comparing its endpoints directly with a standard deviation in mg/L?

Solutions. (1) \(SE=s/\sqrt{n}=3/\sqrt{36}=0.5\) mg/L. It describes estimated sampling variability of the mean, not the spread of individual measurements, which is summarized here by \(s=3\) mg/L. (2) The estimate with \(n=500\) is generally more precise, all else equal, because the standard error of a proportion decreases approximately as \(1/\sqrt{n}\); both samples still require suitable sampling assumptions. (3) Quantile uncertainty depends on sample size and ranks or observations around the quantile, information not contained in the median alone. (4) No. A confidence interval is constructed to quantify sampling uncertainty under specified assumptions; an operational band is a decision criterion whose width needs its own rationale. (5) The units differ: variance is in squared units and standard deviation is in the original units, so they are not directly comparable.

Clinical Biostats