Tutorials › Population Pharmacokinetics › Worked Example: From Forest Plot to Clinical Conclusion

Clinical Relevance of Covariate Effects · Tutorial 322 of 599

Worked Example: From Forest Plot to Clinical Conclusion

Work from forest-plot estimates and confidence intervals to covariate-specific clinical conclusions, while keeping statistical evidence, effect magnitude, and uncertainty distinct.

Advanced 12 min read

What You'll Learn

  • Identify the reference value, contrast, effect scale, and null value for each forest-plot row.
  • Convert a continuous-covariate coefficient and its confidence limits into a clinically selected contrast.
  • Interpret categorical covariate ratios relative to their stated reference group.
  • Compare effect estimates and confidence intervals with a prespecified clinical relevance range.
  • Write conclusions that communicate magnitude and uncertainty without treating statistical significance as clinical importance.

From a Forest Plot to a Clinical Claim

A forest plot can make several covariate effects look directly comparable: each row has an estimate and an interval, often on a common ratio scale. But a plotted number is not yet a clinical conclusion. The meaning depends on the model, the contrast represented by that row, the reference values or categories, and the range of effects considered clinically important. A continuous covariate’s effect also depends on the size of the change being considered; a per-unit estimate may be numerically small while a realistic contrast is not.

This tutorial uses a hypothetical population pharmacokinetic analysis in which covariates affect typical clearance. The forest plot reports ratios of typical clearance for stated contrasts. All numerical values below are invented for teaching. The clinical relevance range used in the examples is also hypothetical: a clearance ratio from 0.80 to 1.25 is treated as a range of effects not considered large enough, by itself, to be clinically important for the decision at hand. Such a range must be justified for the particular drug and clinical context; it is not a universal rule.

Definition: A forest-plot estimate is a model-based estimate of an effect for a specified contrast, commonly displayed with a confidence interval. For a clearance ratio, the null value is 1: a ratio of 1 indicates equal typical clearance between the two stated covariate values or categories. The estimate and interval have meaning only with their contrast, reference, scale, and model assumptions stated.

Read the Contrast Before Reading the Number

Suppose typical clearance is described by a log-linear covariate model, \(CL_{\mathrm{typ}}(x)=CL_{\mathrm{ref}}\exp\{\beta(x-x_{\mathrm{ref}})\}\), for a continuous covariate \(x\). Here \(CL_{\mathrm{ref}}\) is typical clearance at the reference value \(x_{\mathrm{ref}}\), and \(\beta\) is the change in log typical clearance per unit of \(x\). Comparing values \(x_1\) and \(x_0\) gives a ratio, not a difference:

$$ \frac{CL_{\mathrm{typ}}(x_1)}{CL_{\mathrm{typ}}(x_0)} =\exp\{\beta(x_1-x_0)\}. $$

Thus the coefficient alone does not specify a clinically useful effect until a contrast \(x_1-x_0\) is chosen. For a categorical covariate with a reference category, a coefficient \(\gamma\) on the log-clearance scale gives the ratio \(\exp(\gamma)\) for the named category compared with that reference. The forest-plot row should make the comparison direction explicit. “Female versus male,” for example, is not interchangeable with “male versus female”; reversing the contrast takes the reciprocal of the ratio.

A forest plot may place these estimates and intervals on a logarithmic horizontal axis. On a ratio scale, the null line is at 1. A confidence interval entirely on one side of 1 is evidence, under the model and interval procedure, that the ratio differs from 1 at the corresponding confidence level. That fact alone does not say whether the effect is clinically important. Clinical relevance requires comparing the effect’s magnitude and uncertainty with an appropriate clinical benchmark.

Worked Example: Orienting the Rows of a Forest Plot

Consider this hypothetical plot of adjusted typical-clearance ratios. Each interval is a 95% confidence interval. Continuous contrasts and categorical reference groups are shown in the row labels.

Covariate contrastTypical-clearance ratio95% confidence interval
Weight, 90 kg versus 70 kg1.12751.0243 to 1.2411
eGFR, 60 versus 90 mL/min/1.73 m²0.8200.680 to 0.990
Female versus male0.9100.780 to 1.060

The first row compares two specified weights; it does not mean that clearance rises by 12.75% for every kilogram. The second row compares the lower eGFR value with the higher one, so a ratio below 1 means lower typical clearance at 60 than at 90 mL/min/1.73 m². The third row compares females with males, with males as the reference. Its estimate of 0.910 corresponds to an estimated 9.0% lower typical clearance in females, not a 9.0% increase in males relative to females. Reversing that comparison would give \(1/0.910\), approximately 1.099.

All three intervals are compared with the null ratio 1. The weight and eGFR intervals exclude 1; the female-versus-male interval includes 1. These are statements about uncertainty relative to no ratio difference. They do not, by themselves, establish whether any difference is large enough to matter clinically.

Calculate a Clinically Relevant Continuous Contrast

For the weight row, suppose the fitted log-linear model estimates \(\hat{\beta}_W=0.00600\) per kg, with a 95% confidence interval from 0.00120 to 0.0108 per kg. The chosen comparison is 90 kg versus 70 kg, so the contrast is 20 kg. Exponentiating the coefficient multiplied by the contrast gives the estimated ratio:

$$ \widehat{R}_{90:70} =\exp\{(90-70)\hat{\beta}_W\} =\exp(20\times 0.00600) =\exp(0.120) =1.12749685\ldots \approx 1.1275. $$

The corresponding percentage difference is calculated from the ratio, not by treating the ratio itself as a percentage:

$$ (1.12749685\ldots-1)\times 100 =12.749685\ldots\% \approx 12.75\%. $$

Under this model, the estimated typical clearance at 90 kg is about 12.75% higher than at 70 kg. The result is an adjusted model contrast, conditional on the other covariates and model structure, not a causal claim that increasing a person’s weight would cause clearance to rise.

If the coefficient interval is \([0.00120,0.0108]\), the contrast interval follows by multiplying both endpoints by 20 and exponentiating. Because the exponential function is increasing, the order of the endpoints is preserved:

$$ \left[ \exp(20\times0.00120),\, \exp(20\times0.0108) \right] = [\exp(0.024),\,\exp(0.216)] \approx [1.0243,\,1.2411]. $$

The ratio interval corresponds to an estimated range from about 2.43% higher to 24.11% higher typical clearance for 90 kg versus 70 kg. This transformation is valid when the coefficient interval is for the same model coefficient and the contrast is treated as fixed. It does not add uncertainty about the model form, covariate measurement, or whether the selected contrast represents the intended clinical population.

Justification: If a confidence interval for \(\beta\) is \([L,U]\), then for a fixed contrast \(d\), an interval for the ratio \(\exp(d\beta)\) is obtained by transforming the endpoints. For \(d>0\), this is \([\exp(dL),\exp(dU)]\), because multiplication by \(d\) and exponentiation preserve order. For \(d<0\), the order reverses after multiplication, so the transformed endpoints must be arranged from smaller to larger. The resulting interval has the same nominal coverage as the coefficient interval under this monotone transformation.

Compare Magnitude and Uncertainty With Clinical Relevance

For illustration, suppose ratios from 0.80 to 1.25 have been specified in advance as the range not considered clinically important for this hypothetical decision. This benchmark is centered approximately symmetrically on a log scale: 0.80 is the reciprocal of 1.25. A useful assessment considers both the point estimate and the confidence interval against this range. The benchmark should be scientifically justified and chosen for the decision, rather than selected after seeing which boundary makes an estimate appear important.

For weight, the point estimate 1.1275 is within the 0.80-to-1.25 range, and the full 95% interval, 1.0243 to 1.2411, is also inside it. The estimate suggests a difference from the null ratio, but the entire interval remains within the illustrative range of effects not considered clinically important. Under this benchmark, the plot does not suggest a large weight-associated clearance contrast for 90 versus 70 kg, even though the interval excludes 1.

For eGFR, the estimate is 0.820, or an estimated 18.0% lower typical clearance at 60 than at 90 mL/min/1.73 m². The interval from 0.680 to 0.990 excludes 1, but it extends below 0.80. The data are compatible both with a relatively small decrease near 1 and with a decrease exceeding the illustrative clinical boundary. The point estimate alone is just inside that boundary; uncertainty prevents a confident claim that the contrast is clinically important under the stated criterion.

For female versus male, the point estimate is 0.910, an estimated 9.0% lower typical clearance in females. The interval, 0.780 to 1.060, includes 1 and extends below 0.80. It therefore includes no difference, modest differences in either direction, and a potentially larger decrease. It would be inaccurate to conclude that the categories have identical clearance, just as it would be inaccurate to claim a precisely established clinically important decrease. The interval indicates that the effect is uncertain.

Worked Example: Writing a Conclusion for a Categorical Effect

Use the hypothetical female-versus-male row: ratio 0.910, 95% confidence interval 0.780 to 1.060. The ratio is female typical clearance divided by male typical clearance. Convert the estimate to a percentage difference:

$$ (0.910-1)\times100=-9.0\%. $$

For the interval, the lower ratio corresponds to \((0.780-1)\times100=-22.0\%\), and the upper ratio corresponds to \((1.060-1)\times100=6.0\%\). Thus the interval ranges from a 22.0% lower clearance in females to a 6.0% higher clearance in females, relative to males. It contains the null value 1 and crosses both sides of it.

A defensible conclusion is: “The model estimated 9.0% lower typical clearance in females than males, but the 95% confidence interval ranged from 22.0% lower to 6.0% higher and included no difference. The estimate is uncertain; this plot alone does not establish a clinically important sex-related difference.” The conclusion reports direction, size, and uncertainty, without turning an interval that includes 1 into proof of no effect.

Diagnostics, Alternatives, and Limits

Before interpreting a plot, check that each row is on the same effect scale, that the direction of every contrast is explicit, and that each interval is paired with the right estimate. For continuous covariates, verify that the displayed contrast is clinically interpretable and that the model’s assumed functional form is appropriate over that range. For categorical covariates, check the stated reference category and whether the model includes relevant interactions. A single ratio may not summarize the effect if the covariate relationship changes across another covariate.

Look for signs that the plot is being asked to support more than the analysis can: a wide interval, an estimate near a clinical boundary, sparse representation of one category, or a contrast extending beyond the covariate values well represented in the data. Such issues can make an apparently precise plotted summary misleading. The plot also does not show all model uncertainty, such as uncertainty in the chosen structural model, residual error specification, or covariate functional form.

There are several reasonable ways to express an effect. A continuous coefficient per unit is useful for describing the fitted relationship, but a selected contrast is usually easier to interpret clinically. A ratio is natural for multiplicative effects; a difference in clearance units may be preferable when absolute clearance changes matter to the decision. Either representation should be derived from the same fitted model and should state its contrast and units. On a ratio scale, a log-scale plot can display reciprocal changes symmetrically; the clinical interpretation, however, still needs to be stated in words.

Finally, a confidence interval is not a probability statement about the fixed parameter for this particular interval. In repeated applications of the same procedure under its assumptions, 95% of the intervals constructed this way would contain the parameter. The interval’s width describes estimation uncertainty under the fitted model and interval method; it does not guarantee that the model is correct or that every clinically relevant source of uncertainty has been represented.

Common Mistakes

  • Calling a confidence interval that excludes 1 clinically important. Exclusion of 1 concerns evidence for a non-null ratio; clinical importance depends on a clinical benchmark and the magnitude of the effect.
  • Reporting a continuous coefficient without its contrast. A coefficient per kg does not state the effect from 70 to 90 kg. Multiply by the contrast on the log scale, then exponentiate.
  • Reversing a categorical comparison without changing the ratio. Female versus male and male versus female have reciprocal ratios and different verbal interpretations.
  • Treating an interval that includes 1 as proof of no difference. It indicates that no difference is among the values compatible with the interval; the interval may also include effects that matter clinically.
  • Rounding intermediate values too early. Carry sufficient precision through exponentiation and percentage calculations, then round the final result. For example, \(\exp(0.120)=1.12749685\ldots\), which rounds to 1.1275 at four decimal places or 1.127 to three decimal places.
Key takeaway: A sound forest-plot conclusion names the covariate contrast and reference, reports the effect on an interpretable scale, and describes its uncertainty relative to a justified clinical benchmark. Statistical evidence against a null effect and evidence of clinical importance are related but distinct judgments.

Check Your Understanding

Answer each question from the definitions and hypothetical results above. Brief solutions follow the questions.

  1. A continuous-covariate coefficient is \(0.0040\) per unit. What ratio contrast does the model imply for a 25-unit increase?
  2. A categorical contrast has ratio 1.20 for group A versus group B. What is the ratio for group B versus group A?
  3. A ratio estimate is 0.84 with a 95% confidence interval of 0.72 to 0.98. Does the interval exclude the null ratio, and does it establish a clinically important effect if the illustrative boundary is 0.80?
  4. Why is “the covariate has no effect” not justified solely because a confidence interval includes 1?

Solutions. (1) The ratio is \(\exp(25\times0.0040)=\exp(0.100)=1.1052\), approximately 1.105, or about 10.52% higher for the stated increase. (2) Reverse the comparison by taking the reciprocal: \(1/1.20=0.8333\). (3) The interval excludes 1, but it crosses the 0.80 clinical boundary, so it does not establish that the effect exceeds that boundary; the interval includes both ratios below 0.80 and ratios from 0.80 to 0.98. (4) An interval containing 1 is compatible with no difference, but it may also contain non-null effects, including clinically important ones. Its width and endpoints determine how much uncertainty remains.

Clinical Biostats