From Exposure to Time-to-Event Efficacy
A binary efficacy endpoint records whether an event occurred by a chosen visit. A time-to-event endpoint also uses when it occurred, while retaining information from participants whose event has not been observed by the end of follow-up. This distinction matters when follow-up differs across participants: an event-free participant observed for two months does not provide the same information as one observed for twelve months.
An exposure-dependent hazard model links a participant’s exposure to the instantaneous rate of the efficacy event, conditional on remaining event-free up to that time. The model can then predict event-free survival at a specified follow-up time. Those predictions are conditional on the model and its assumptions; an exposure association alone does not establish that changing exposure would cause the predicted change in efficacy.
For a concrete model, let \(x\) denote a participant’s exposure, measured in a defined unit, and let \(x_{\mathrm{ref}}\) be a reference exposure. Use a Weibull baseline hazard with an exposure effect proportional on the hazard scale:
Here \(\eta>0\) is the Weibull scale parameter, in units of time, \(k>0\) is the shape parameter, and \(\beta\) is the log hazard ratio for a 10-unit increase in exposure. The exposure unit must be stated in context; for example, if \(x\) is measured in mg·h/L, the increment is 10 mg·h/L. At \(x=x_{\mathrm{ref}}\), the exposure multiplier is 1, so the hazard is the baseline hazard. The model implies
The baseline event risk by time \(t\), at the reference exposure, is \(1-S(t\mid x_{\mathrm{ref}})\). This is an absolute probability over a stated period. By contrast, \(\exp(\beta)\) is a hazard ratio: at any time when the model applies, it compares the instantaneous hazards for exposures differing by 10 units. It is not a risk ratio, and it does not by itself specify the difference in event probabilities.
Assumptions and Right-Censoring
The model is appropriate only when its structure and data collection are defensible. The Weibull baseline allows the hazard to increase over time when \(k>1\), decrease when \(k<1\), or remain constant when \(k=1\). The exposure effect is proportional: the hazard ratio for a given exposure contrast does not change with time. Exposure must be defined consistently, and the model treats the specified exposure value as the relevant predictor throughout the modeled risk period. If exposure changes over time, this fixed-exposure model does not represent that change.
The analysis also assumes that event times are defined and observed consistently, participants are independent conditional on included predictors, and right-censoring is noninformative conditional on those predictors. In practical terms, after accounting for modeled information, a participant’s censoring time should not provide additional information about their unobserved event time. Loss to follow-up related to worsening efficacy can violate this assumption. Administrative study end can often be reasonable censoring, but that judgment depends on study conduct and the event process.
For participant \(i\), let \(Y_i\) be the observed time, the smaller of the event time and censoring time, and let \(\delta_i=1\) if the event was observed and \(0\) if follow-up was right-censored. The likelihood contribution is
An observed event contributes its event-time density \(f=hS\); a censored participant contributes the probability of remaining event-free through the observed follow-up, \(S\). This follows from \(f(t)=-dS(t)/dt\) and \(h(t)=f(t)/S(t)\). Under the stated independent-censoring assumption, the censoring distribution does not need to be modeled to estimate the event-time parameters using these contributions. Treating a censored participant as though an event occurred at the censoring time would instead use the wrong contribution.
Worked Example: Baseline Risk and Alternative Exposures
Suppose a hypothetical study models time to first efficacy failure in months. Exposure \(x\) is average exposure in mg·h/L, the reference is \(x_{\mathrm{ref}}=30\) mg·h/L, and fitted parameters are \(\eta=12\) months, \(k=1.4\), and \(\beta=-0.40\) per 10 mg·h/L. Assume the event definition and follow-up origin are common across participants, the model’s exposure effect is proportional, and censoring is independent conditional on exposure. These assumptions are modeling conditions, not facts established by the parameter values.
At the reference exposure and 12 months, \(H(12\mid30)=(12/12)^{1.4}=1\). Therefore,
The model predicts 36.79% event-free survival and 63.21% event risk by 12 months at the reference exposure. These are model-based probabilities for the stated population and time point, not estimates of an individual’s certainty.
At \(x=40\) mg·h/L, the exposure contrast is 10 mg·h/L. The hazard ratio is \(\exp(-0.40)=0.6703\), and the 12-month cumulative hazard is \(1(0.6703)=0.6703\). Thus \(S(12\mid40)=\exp(-0.6703)=0.5114\), giving event risk \(1-0.5114=0.4886\). The model predicts higher event-free survival at this exposure. The hazard is 32.97% lower at each time under the model, but the 12-month event risk is 14.35 percentage points lower than at the reference exposure; these are different comparisons.
At \(x=50\) mg·h/L, the contrast is 20 units and the hazard ratio is \(\exp(2[-0.40])=0.4493\). The 12-month cumulative hazard is 0.4493, so survival is \(\exp(-0.4493)=0.6380\) and event risk is \(1-0.6380=0.3620\). The absolute risk difference from the reference is 0.3620 minus 0.6321, or -0.2701: 27.01 percentage points lower by 12 months.
Working Through Censored Follow-Up
Censoring preserves partial information. A participant censored at month 6 is known to have remained event-free through month 6, even though the later event time is unknown. The Weibull model assigns a survival contribution through that time rather than inventing an event time.
Worked Example: Event and Censoring Likelihood Contributions
Use the same hypothetical parameters as above. At exposure \(x=40\) mg·h/L and time \(t=6\) months, the cumulative hazard is
The survival probability is \(S(6\mid40)=\exp(-0.2540)=0.7757\). The hazard at month 6 is \(h(6\mid40)=(k/t)H(t\mid40)=(1.4/6)(0.2540)=0.05927\) per month.
For a participant who experiences the event at month 6, the likelihood contribution is the event density \(hS=(0.05927)(0.7757)=0.04598\) per month. For a participant censored at month 6, the contribution is \(S=0.7757\). The first observation says an event occurred at that time; the second says only that the participant remained event-free through that time. The calculation uses the same conditional model for both participants and assumes their censoring is noninformative given exposure.
For a data set, multiply the individual likelihood contributions, or equivalently sum their log contributions, and estimate \(\eta\), \(k\), and \(\beta\) by maximizing that likelihood. A reproducible analysis should record the event definition, time origin, exposure summary and units, censoring rule, model parameterization, and estimation settings. The fitted model should be checked before its predictions are used.
Interpreting Time-Based Predictions
Worked Example: Comparing Median Event Times
The median event time is the time \(t_{0.5}\) at which \(S(t_{0.5}\mid x)=0.5\). For the specified Weibull model, solving the survival equation gives
At the reference exposure of 30 mg·h/L, the median is \(12[\log(2)]^{1/1.4}=9.24\) months. At 40 mg·h/L, the exposure multiplier is 0.6703, giving \(12[0.6931/0.6703]^{1/1.4}=12.29\) months. Under this hypothetical model, median event time is longer at 40 mg·h/L. This is a model-based comparison; it should not be interpreted as a causal benefit unless the study design and analysis justify a causal interpretation.
A median is not the only useful prediction. Event-free survival at a clinically meaningful time, such as 6 or 12 months, may be more informative when the median is beyond the observed follow-up. Report the time point and exposure contrast, and distinguish a predicted survival probability from a hazard ratio.
The main alternatives answer related but different questions. A Kaplan–Meier estimate describes event-free survival over time without imposing this parametric Weibull shape; it is useful for descriptive comparisons but does not by itself provide the same exposure-adjusted parametric predictions. A Cox proportional-hazards model leaves the baseline hazard unspecified while estimating hazard ratios, so it is less restrictive about baseline shape but does not directly supply a Weibull-based absolute survival curve without additional baseline estimation. The Weibull model gives a compact baseline and direct survival predictions, at the cost of assuming that baseline form. These methods do not remove confounding in observational exposure comparisons.
Diagnostics, Failure Modes, and Limitations
Diagnostics should address both the event-time form and how exposure and censoring enter the analysis. Compare model-predicted survival with descriptive survival patterns where data support that comparison. For a Weibull baseline, a plot of log cumulative hazard against log time should be approximately linear if the Weibull shape is adequate; systematic curvature is a warning, not an automatic proof of failure. Assess whether exposure contrasts appear compatible with a constant hazard ratio over time. Marked time-dependent differences challenge proportional hazards and may require a different model.
Also inspect event counts and follow-up across the exposure range. Sparse events, little exposure overlap, or few participants remaining under observation late in follow-up can make estimates and long-term predictions unstable. Compare censoring patterns by exposure and relevant covariates; concentration of early loss to follow-up in a group may signal a problem with the conditional independent-censoring assumption. A model cannot recover event-time information that was never observed without assumptions.
- Confusing hazard and risk. A hazard ratio of 0.67 is not necessarily a 33% reduction in event probability by a fixed time. Calculate \(1-S(t\mid x)\) to report an absolute risk.
- Calling censored follow-up an event. A participant event-free through month 6 contributes \(S(6)\), not the event density at month 6.
- Ignoring time and units. A prediction at 12 months is not interchangeable with one at 12 weeks. State exposure units, the exposure contrast, and the time origin.
- Over-interpreting exposure association. Differences in exposure may reflect prognosis or other covariates. A fitted association does not by itself show that increasing dose will improve efficacy.
The model concerns time to a first event and right-censoring. It does not automatically handle competing event types, recurrent efficacy events, exposure changing over time, or informative dropout. Its predictions also depend on the population and exposure range represented in the data; extrapolating to unsupported exposure levels or follow-up times can be misleading.
Specify the efficacy event, time origin, exposure measure and units, and rules for event observation and censoring.
Estimate the Weibull baseline and exposure coefficient using event and censoring likelihood contributions; assess baseline shape, proportional hazards, follow-up, and censoring patterns.
Give event-free survival or event risk at stated times and exposures, alongside hazard-ratio interpretations and the assumptions that qualify them.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- In a Weibull exposure model, what does \(\exp(\beta)\) compare, and why is it not automatically a risk ratio?
- At the reference exposure, suppose \((t/\eta)^k=0.5\). Calculate event-free survival and event risk at time \(t\).
- A participant is event-free when last observed at month 8. What is their event-time likelihood contribution under independent right-censoring?
- With \(\beta=-0.30\) per 10 exposure units, calculate the hazard ratio for a 20-unit increase.
- Name one diagnostic concern for the Weibull baseline and one for the exposure effect.
Solutions. (1) It is the instantaneous hazard ratio for a 10-unit exposure increase, assuming proportional hazards. A risk is accumulated over time through the survival function, so its ratio depends on baseline risk and follow-up time. (2) Survival is \(\exp(-0.5)=0.6065\); event risk is \(1-0.6065=0.3935\). (3) The contribution is \(S(8\mid x)\), because the event is known not to have occurred through month 8. (4) The ratio is \(\exp(2[-0.30])=\exp(-0.60)=0.5488\). (5) Curvature in the log cumulative hazard versus log time can challenge the Weibull baseline; evidence that exposure hazard ratios change over time can challenge proportional hazards.