Tutorials › Population Pharmacokinetics › Power-Function Models for Body-Size Effects

Body Size, Age and Maturation · Tutorial 202 of 599

Power-Function Models for Body-Size Effects

Formulate and interpret a body-weight power model, including its reference parameter, exponent, predictions, assumptions, and limits.

Intermediate 10 min read

What You'll Learn

  • Write a positive PK parameter as a power function of weight relative to a reference weight.
  • Interpret the reference parameter and exponent in context and with units.
  • Calculate model predictions and finite proportional changes across body weights.
  • Estimate an exponent from two parameter values and their corresponding weights.
  • Identify assumptions, diagnostic concerns, and risks of extrapolating a power model.

Why Use a Power Function for Body Weight?

A pharmacokinetic parameter may change with body weight in a way that is better represented as a proportional relationship than as a fixed additive change. A power function provides one compact model for that pattern. It separates the parameter's value at a chosen reference weight from the strength and direction of its weight relationship. The model is a useful form to consider; its use does not, by itself, establish that the relationship is correct for a particular drug or population.

The reference weight makes the parameter interpretable in its original units, while the exponent describes how the parameter changes across weights. Both pieces matter: an exponent alone does not give a predicted parameter value, and the reference parameter alone does not describe how values change away from the reference.

Definition: A body-weight power-function model for a positive pharmacokinetic parameter \(P\) is
$$ P(W)=P_{\mathrm{ref}}\left(\frac{W}{W_{\mathrm{ref}}}\right)^\theta. $$
Here \(W\) is body weight, \(W_{\mathrm{ref}}\) is a specified positive reference weight, \(P_{\mathrm{ref}}\) is the value of the parameter at \(W=W_{\mathrm{ref}}\), and \(\theta\) is the power exponent. \(W\) and \(W_{\mathrm{ref}}\) must use the same units. The ratio \(W/W_{\mathrm{ref}}\) is dimensionless, so \(P(W)\) has the same units as \(P_{\mathrm{ref}}\).

What the Reference Value and Exponent Mean

Substituting \(W=W_{\mathrm{ref}}\) into the model gives \(P(W_{\mathrm{ref}})=P_{\mathrm{ref}}(1)^\theta=P_{\mathrm{ref}}\). Thus, the reference parameter is the model-predicted parameter value at the reference weight. It is not necessarily the parameter for every person at that weight, nor does it have to be an observed value from a particular participant.

The exponent describes the relationship on a proportional scale. For example, when weight is multiplied by a factor \(c\), the predicted parameter is multiplied by \(c^\theta\). If \(\theta\) is positive, the prediction increases as weight increases; if it is negative, the prediction decreases; and if it is zero, the prediction is constant with respect to weight. The absolute change depends on both the exponent and the starting parameter value.

A useful way to justify this interpretation is to take logarithms of the model:

$$ \ln P(W)=\ln P_{\mathrm{ref}}+\theta\ln\left(\frac{W}{W_{\mathrm{ref}}}\right). $$

This is a straight-line relationship between \(\ln P(W)\) and \(\ln(W/W_{\mathrm{ref}})\), with slope \(\theta\) and intercept \(\ln P_{\mathrm{ref}}\). Equivalently, differentiating with respect to log weight gives the proportional sensitivity:

$$ \frac{d\ln P}{d\ln W}=\theta. $$

So \(\theta\) is the elasticity of the parameter with respect to weight: for a small relative change in weight, the relative change in the predicted parameter is approximately \(\theta\) times as large. For a finite weight change, use the power function rather than treating that approximation as exact:

$$ \frac{P(W_2)}{P(W_1)}=\left(\frac{W_2}{W_1}\right)^\theta. $$

This ratio also shows that the relative change between two weights does not depend on \(P_{\mathrm{ref}}\), provided the same exponent applies at both weights. The absolute difference does depend on the parameter's scale.

Worked Example: Predicting a Parameter Above the Reference Weight

Suppose a hypothetical model for clearance has \(P_{\mathrm{ref}}=8.00\ \text{L/h}\) at \(W_{\mathrm{ref}}=70\ \text{kg}\), with \(\theta=0.75\). Predict clearance at \(84\ \text{kg}\).

The weight ratio is \(84/70=1.20\). Substitution gives:

$$ P(84)=8.00\left(\frac{84}{70}\right)^{0.75} =8.00(1.20)^{0.75} \approx8.00(1.14653) \approx9.172\ \text{L/h}. $$

The model predicts about \(9.172\ \text{L/h}\), approximately \(1.172\ \text{L/h}\) above the reference value. The multiplier is about \(1.14653\), so the predicted clearance is about \(14.65\%\) higher than at \(70\ \text{kg}\). This is a model-based comparison, not a claim that changing an individual's weight would necessarily cause that change in clearance.

Worked Example: Interpreting a Negative Exponent

Consider a hypothetical positive parameter with \(P_{\mathrm{ref}}=30.0\ \text{L}\) at \(W_{\mathrm{ref}}=60\ \text{kg}\), and \(\theta=-0.30\). Find the predicted value at \(75\ \text{kg}\) and describe the proportional change.

The weight ratio is \(75/60=1.25\). Therefore:

$$ P(75)=30.0\left(\frac{75}{60}\right)^{-0.30} =30.0(1.25)^{-0.30} \approx30.0(0.93525) \approx28.06\ \text{L}. $$

The prediction is about \(28.06\ \text{L}\). Relative to \(30.0\ \text{L}\), the model-based decrease is \(1-0.93525=0.06475\), or about \(6.475\%\) (approximately \(6.48\%\)). The negative exponent means the predicted parameter is lower at the higher weight under this model. It does not mean that the parameter falls by \(0.30\) units for each kilogram.

Obtaining an Exponent from Two Values

If two positive parameter values are assumed to follow the same power relationship, their ratio can be used to solve for the exponent. Starting from the model at weights \(W_1\) and \(W_2\), divide one equation by the other. The reference parameter cancels:

$$ \frac{P(W_2)}{P(W_1)} =\left(\frac{W_2}{W_1}\right)^\theta. $$

Taking natural logarithms and rearranging yields:

$$ \theta= \frac{\ln\left(P(W_2)/P(W_1)\right)} {\ln\left(W_2/W_1\right)}. $$

This calculation requires distinct weights, positive parameter values, and a common power exponent across the two weights. It gives the exponent that exactly connects those two values under the assumed model. With noisy observed or estimated parameter values, that two-point calculation is not by itself a statistical fitting procedure.

Worked Example: Calculating an Exponent from Two Parameter Values

Suppose hypothetical parameter values are \(5.00\ \text{L/h}\) at \(50\ \text{kg}\) and \(7.20\ \text{L/h}\) at \(80\ \text{kg}\). Find the exponent that connects them with a power function.

Both ratios are dimensionless, and the weights differ, so the two-point formula applies:

$$ \theta= \frac{\ln(7.20/5.00)}{\ln(80/50)} =\frac{\ln(1.44)}{\ln(1.60)} \approx\frac{0.36464}{0.47000} \approx0.7758. $$

The connecting exponent is approximately \(0.7758\). Its positive sign indicates an increasing parameter across these weights. For a weight multiplier of \(1.60\), the parameter multiplier is \(1.60^{0.7758}\), approximately \(1.44\), which agrees with \(7.20/5.00=1.44\). These two values determine a line on the log-log scale, but cannot establish whether that relationship adequately describes weights between or beyond \(50\) and \(80\ \text{kg}\).

Assumptions, Diagnostics, and Alternatives

The power-function form has simple mathematical requirements and stronger scientific assumptions. The parameter and weights must be positive for the model and its logarithmic interpretation to apply. The reference and measured weights must be expressed in the same units. More substantively, the model assumes that one exponent is an adequate description over the weight range being considered, with other modeled influences held fixed. It also assumes that the selected parameter is meaningfully represented by this positive, multiplicative relationship.

A fitted power relationship should be evaluated rather than accepted solely because its formula is convenient. Useful checks include whether predictions are sensible across the observed weight range, whether residual patterns suggest curvature or distinct subgroups, and whether the estimated exponent is sufficiently precise for the intended use. A narrow range of observed weights may provide little information about the slope on the log-log scale. Predictions far beyond the weights represented in the data are extrapolations and can be unreliable even when the fitted curve looks reasonable within the data.

The main alternative is an additive linear relationship, such as \(P(W)=P_{\mathrm{ref}}+\beta(W-W_{\mathrm{ref}})\), where \(\beta\) represents an absolute change per unit of weight. That form and a power function answer different questions: the additive form assigns the same absolute change for each kilogram, whereas the power form assigns a change tied to proportional weight. A categorical weight-group model is another possible approach, but it does not impose a smooth relationship across weights. Choice among these forms should reflect the scientific question, data support, and predictive behavior—not merely which equation is easiest to write.

A power function has limitations. It is positive for positive \(P_{\mathrm{ref}}\) and \(W\), which may be unsuitable for a parameter that can legitimately be zero or negative. A single exponent may not represent different weight ranges equally well. The model also does not explain why the relationship occurs, nor does it imply that weight is the only relevant covariate. Finally, a reference value and exponent are estimated quantities in an analysis; their uncertainty should be considered when drawing conclusions from predictions.

Common pitfall: Saying “\(\theta=0.75\) means the parameter increases by 0.75 units per kilogram” confuses a dimensionless exponent with an additive slope. The correction is that \(\theta\) describes proportional sensitivity: a weight multiplier \(c\) corresponds to a parameter multiplier \(c^\theta\). For a particular finite change, calculate that multiplier and report the resulting prediction in the parameter's units.

Using the Model Responsibly

1
Specify the parameter and reference.
Name the parameter, its units, the reference weight, and the units used for weight.
2
State the exponent and its direction.
A positive exponent implies an increasing prediction with weight; a negative exponent implies a decreasing prediction; zero implies no modeled weight effect.
3
Calculate with the ratio.
Use \(W/W_{\mathrm{ref}}\) and preserve enough intermediate precision to avoid rounding errors.
4
Check scope and interpretation.
Compare the prediction with the data-supported weight range, report units, and describe it as a model-based relationship rather than a causal effect.

This order helps prevent two common reporting problems: presenting an exponent without stating the weight scale and reference, and interpreting a relative effect as an absolute change. A clear report gives \(P_{\mathrm{ref}}\), \(W_{\mathrm{ref}}\), \(\theta\), the predicted parameter with units, and the relevant weight comparison.

Check Your Understanding

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

  1. A model has \(\theta=0.50\). What parameter multiplier does it predict when weight doubles?
  2. For \(P_{\mathrm{ref}}=10.0\ \text{L/h}\), \(W_{\mathrm{ref}}=75\ \text{kg}\), and \(\theta=0.50\), calculate \(P(90)\).
  3. Parameter values of \(6.00\ \text{L/h}\) at \(50\ \text{kg}\) and \(9.00\ \text{L/h}\) at \(75\ \text{kg}\) are connected by a power function. Calculate the exponent.
  4. In one sentence, distinguish the interpretation of a power exponent from that of an additive slope measured in parameter units per kilogram.

Solutions. (1) The multiplier is \(2^{0.50}=\sqrt{2}\approx1.414\), so the predicted parameter is about \(41.4\%\) higher. (2) \(P(90)=10.0(90/75)^{0.50}=10.0\sqrt{1.20}\approx10.95\ \text{L/h}\). (3) \(\theta=\ln(9.00/6.00)/\ln(75/50)=\ln(1.50)/\ln(1.50)=1.000\). (4) A power exponent describes proportional sensitivity, whereas an additive slope describes the absolute parameter change associated with a one-kilogram increase.

Clinical Biostats