From a Concentration-QTc Model to an Exposure-Specific Estimate
A fitted concentration-QTc relationship is useful for safety assessment only when its implications are evaluated at relevant exposures. The central calculation is the model-based mean QTc effect at a specified concentration, together with an interval that reflects uncertainty in that estimate. Results are often examined at therapeutic exposure and at a higher, supratherapeutic exposure because the latter can show whether the exposure-response relationship remains reassuring beyond the intended clinical range.
Here, the effect is the model-estimated mean, time-matched, baseline-adjusted, placebo-corrected QTc change, written \(\delta(C)\), at concentration \(C\). This is a population-average contrast under the model, not a prediction of one particular patient's ECG. The earlier material on baseline-adjusted and time-matched QTc change and on concentration-QTc models establishes how this contrast and its exposure-response relationship are formed.
This calculation assumes that the fitted model is an appropriate representation of the exposure-response relationship over the concentration range being evaluated; the contrast and concentration are defined on the intended scales; and the estimated coefficient covariance matrix reflects the model fit. For a mixed-effects analysis, the target here is the fixed-effect population mean contrast, with the model's appropriate uncertainty calculation. If concentration itself is uncertain, treating \(C^*\) as fixed does not account for that additional uncertainty.
Estimating the Effect and Its Uncertainty
Write the coefficient estimate as the vector \(\hat{\boldsymbol{\beta}}=(\hat{\beta}_0,\hat{\beta}_1)^\mathsf{T}\), and let \(\widehat{\mathrm{Var}}(\hat{\boldsymbol{\beta}})\) be its estimated covariance matrix. For a fixed concentration \(C^*\), the estimate is a linear combination of the coefficients. Its estimated variance is therefore:
The covariance term matters. An intercept and slope estimated from the same data are generally correlated, so adding their variances alone is not correct. The standard error describes uncertainty in the estimated mean effect; it is not the standard deviation of individual QTc responses.
With approximately normal estimation uncertainty, a two-sided 90% confidence interval is the estimate plus or minus 1.645 standard errors. With finite degrees of freedom, use the appropriate \(t\) critical value \(t_{0.95,df}\) instead. The upper endpoint of this two-sided 90% interval is also a one-sided 95% upper confidence bound under the same reference distribution.
The 1.645 multiplier is the large-sample normal approximation. The interpretation of a 90% confidence interval is about the long-run coverage of the interval procedure under its assumptions; it is not a 90% probability that the fixed true mean effect lies in this particular interval. An upper bound below a prespecified threshold supports excluding a mean effect at or above that threshold under the model and analysis assumptions. It does not establish that the drug has no QTc effect or that every individual effect is below the threshold.
Worked Example: Checking an Upper Bound Near 10 ms
Suppose a hypothetical fitted model gives an estimated mean QTc effect of 4.8 ms at a specified exposure, with estimated variance 9.4196 ms². Use the large-sample 90% two-sided interval to calculate the upper endpoint, and compare it with 10 ms.
First, the standard error is \(SE=\sqrt{9.4196}=3.0691\) ms. The 90% interval is:
The upper endpoint is approximately 9.849 ms. It is just below 10 ms, so this calculation places the upper confidence bound below that threshold, narrowly. The result does not mean that the estimated effect is zero: the point estimate is 4.8 ms. Nor should the small margin be presented as strong evidence of a large safety margin. Rounding, model adequacy, and uncertainty in exposure assignment deserve attention when an endpoint is this close to the threshold.
Comparing Therapeutic and Supratherapeutic Exposures
Therapeutic and supratherapeutic concentrations should be selected to represent scientifically relevant exposure levels, not chosen simply because they make the fitted curve look reassuring. A supratherapeutic exposure may be informed by clinical pharmacokinetics, including conditions associated with increased exposure. The model-based contrast is evaluated at each concentration, and its standard error is recalculated there because the uncertainty in a slope contributes increasingly as concentration moves away from zero.
Worked Example: Effects at Two Prespecified Exposures
Consider a hypothetical linear model with \(\hat{\beta}_0=1.20\) ms and \(\hat{\beta}_1=0.0350\) ms per ng/mL. Suppose the fitted coefficient covariance information is \(\widehat{\mathrm{Var}}(\hat{\beta}_0)=1.00\) ms², \(\widehat{\mathrm{Var}}(\hat{\beta}_1)=0.000100\) (ms per ng/mL)², and \(\widehat{\mathrm{Cov}}(\hat{\beta}_0,\hat{\beta}_1)=-0.0060\) ms² per ng/mL. Use a large-sample 90% interval. Let 80 ng/mL represent the hypothetical therapeutic exposure and 200 ng/mL the hypothetical supratherapeutic exposure.
At 80 ng/mL, the estimated mean effect is:
Its estimated variance is \(1.00+2(80)(-0.0060)+(80)^2(0.000100)=1.00-0.960+0.640=0.680\) ms². Thus \(SE=\sqrt{0.680}=0.8246\) ms, and the interval is \(4.00\pm1.645\sqrt{0.680}=[2.643,\ 5.357]\) ms, to three decimal places.
At 200 ng/mL, the estimated mean effect is \(1.20+0.0350(200)=8.20\) ms. Its variance is \(1.00+2(200)(-0.0060)+(200)^2(0.000100)=1.00-2.400+4.000=2.600\) ms². Therefore \(SE=\sqrt{2.600}=1.6125\) ms, giving \(8.20\pm1.645\sqrt{2.600}=[5.548,\ 10.852]\) ms.
The therapeutic-exposure upper bound is below 10 ms; the supratherapeutic-exposure upper bound is above 10 ms. Under this hypothetical model, the specified calculation therefore does not exclude a mean effect of 10 ms or greater at 200 ng/mL. This is not proof that the true effect reaches 10 ms: the point estimate is 8.20 ms, and the interval reflects substantial uncertainty. The negative coefficient covariance reduces the variance at 80 ng/mL, while the slope variance contributes more at 200 ng/mL.
Why the Upper Bound Matters, and What It Does Not Say
A threshold-based assessment asks whether the data and model estimate are sufficiently precise to place the upper confidence bound below a prespecified level of concern. In concentration-QTc assessment, 10 ms is a commonly used reference point for evaluating the upper bound of the mean effect. It is a decision criterion in a defined assessment framework, not a biological boundary that separates harmless from harmful outcomes in every setting. The size of the point estimate, the uncertainty, the exposure range, and the broader safety context all matter.
The following short justification explains why the covariance formula is necessary. Let \(a=(1,C^*)^\mathsf{T}\). Then \(\hat{\delta}(C^*)=a^\mathsf{T}\hat{\boldsymbol{\beta}}\). For any fixed vector \(a\), the variance of a linear combination is \(a^\mathsf{T}\mathrm{Var}(\hat{\boldsymbol{\beta}})a\). Expanding that matrix product gives:
The quadratic form also makes clear why the covariance matrix must be coherent: the resulting variance cannot be negative. A negative result signals an implementation, transcription, or covariance-matrix problem rather than a meaningful negative standard error.
An upper confidence bound for the mean must not be confused with a prediction interval for one future observation. A prediction interval adds variability in individual observations (and, when relevant, between-person variability) to uncertainty in the estimated mean. That interval is usually wider. Which quantity is needed depends on the question: a threshold assessment of the modelled mean effect is not answered by an interval for an individual ECG.
Worked Example: Mean-Effect Interval Versus Individual Prediction
At 160 ng/mL, use the hypothetical model and covariance values from the preceding example. The mean estimate is \(\hat{\delta}(160)=1.20+0.0350(160)=6.80\) ms. Its variance is \(1.00+2(160)(-0.0060)+(160)^2(0.000100)=1.00-1.92+2.56=1.64\) ms², so its standard error is 1.2806 ms. The upper endpoint of its 90% interval is \(6.80+1.645(1.2806)=8.907\) ms.
Now, only for illustration, suppose the conditional standard deviation of a future individual observation around the modelled mean is 5.0 ms, and assume normal, independent residual variation. A simple large-sample upper prediction endpoint is then:
The 8.907 ms endpoint describes uncertainty in the population mean effect; the illustrative 15.29 ms endpoint also includes individual residual variation. They answer different questions. This simplified prediction calculation is not a substitute for a fitted model's subject-level prediction method, especially when repeated measurements, random effects, or non-independent residuals are present.
Diagnostics, Alternatives, and Limitations
Before interpreting an exposure-specific upper bound, check that the target concentration is supported by the data and that the model's assumptions remain credible there. The earlier material on diagnostics for concentration-QTc model misspecification is relevant: curvature, influential observations, time-related patterns, or systematic residual trends can make a precise-looking linear-model bound misleading. A therapeutic exposure inside a well-sampled range and a supratherapeutic exposure far beyond the observed concentrations do not have equal evidential support.
- Check the concentration range. Mark the exposure values used for predictions against the observed concentration distribution. Predictions far outside that range are extrapolations and depend strongly on the assumed model form.
- Check the uncertainty calculation. Use the covariance matrix from the fitted model, including the intercept-slope covariance, and use the appropriate critical value and degrees of freedom where required.
- Check the target and estimand. Confirm that the predicted contrast is the intended time-matched, baseline-adjusted, placebo-corrected mean effect and that the exposure represents the clinical scenario being assessed.
- Check sensitivity to model choices. If plausible nonlinear forms or covariate specifications yield materially different upper bounds, the conclusion depends on model assumptions and should not be reported as robust.
The linear model is not the only possible approach. A nonlinear exposure-response model can represent a plateau or other curvature, but its prediction and uncertainty depend on the nonlinear parameter estimates and may be less stable when data are sparse. A categorical exposure analysis avoids imposing a straight-line relationship but can lose precision and does not naturally interpolate between exposure groups. Either alternative still requires a clearly defined mean contrast and an uncertainty calculation appropriate to its fitted model.
If exposure estimates are uncertain, a fixed-concentration calculation may understate total uncertainty. A model that propagates PK uncertainty, or a sensitivity analysis over plausible concentrations, may be needed; the earlier material on concentration assignment and PK uncertainty addresses that issue. Similarly, an upper confidence bound is only as reliable as the study design, ECG measurement quality, model specification, and uncertainty method. It cannot repair inadequate sampling or a poorly supported exposure range.
Common Mistake
A related mistake is to interpret an upper bound below 10 ms as proof of no effect, or as a guarantee that no individual will exceed 10 ms. It supports a narrower conclusion: under the stated model and assumptions, the upper confidence bound for the specified population mean contrast is below the assessment threshold. Other safety findings and the clinical context remain relevant.
Check Your Understanding
Try each question, then compare your reasoning with the solutions below.
- A linear model estimates \(\hat{\beta}_0=0.5\) ms and \(\hat{\beta}_1=0.020\) ms per ng/mL. What is the estimated mean effect at 100 ng/mL?
- At that exposure, suppose the coefficient variances are 0.64 ms² and 0.000064 (ms per ng/mL)², and their covariance is \(-0.0040\) ms² per ng/mL. Calculate the variance and standard error of the estimated mean effect.
- If the estimated mean effect is 2.5 ms with \(SE=2.0\) ms, calculate the large-sample 90% interval and state whether its upper bound is below 10 ms.
- Why is it incorrect to calculate the variance of \(\hat{\beta}_0+C^*\hat{\beta}_1\) by adding only the intercept and slope variance contributions?
- What additional variability is represented by an individual prediction interval that is not represented by a confidence interval for the modelled mean?
Solutions. (1) \(\hat{\delta}(100)=0.5+0.020(100)=2.5\) ms. (2) The variance is \(0.64+2(100)(-0.0040)+(100)^2(0.000064)=0.64-0.80+0.64=0.480\) ms², so \(SE=\sqrt{0.480}=0.6928\) ms. (3) \(2.5\pm1.645(2.0)=2.5\pm3.290\), giving \([-0.790,\ 5.790]\) ms; the upper bound is below 10 ms. (4) The intercept and slope estimates may be correlated, so the variance must include \(2C^*\widehat{\mathrm{Cov}}(\hat{\beta}_0,\hat{\beta}_1)\). (5) An individual prediction interval also includes variation among future individual observations, such as residual and potentially between-person variation, whereas the mean-effect confidence interval quantifies uncertainty in the estimated population mean.