Tutorials › Joint Modeling of Longitudinal and Time-to-Event Data › Interpreting the Current-Value Association Parameter

The Basic Joint Model · Tutorial 131 of 697

Interpreting the Current-Value Association Parameter

Interpret the current-value coefficient on the log-hazard scale, express it as a hazard ratio for meaningful marker differences, and recognize the assumptions and limits of that interpretation.

Advanced 9 min read

What You'll Learn

  • State the model conditions needed to interpret the current-value association coefficient.
  • Derive the hazard ratio for a specified difference in latent marker values.
  • Interpret positive, negative, and zero coefficients in context and with units.
  • Transform a coefficient confidence interval to the hazard-ratio scale.
  • Distinguish conditional hazard comparisons from risk, causal, and observed-marker comparisons.
  • Identify common diagnostics and limitations of the current-value interpretation.

What the Association Parameter Compares

In a joint model, the longitudinal submodel represents a subject’s underlying marker trajectory, while the event submodel relates that latent trajectory to the event hazard. The current-value association parameter describes how the hazard changes with the marker’s latent value at the time being considered. Its interpretation depends on the scale of the marker and on which other hazard predictors are held fixed.

Let \(m_i(t)\) denote the latent current value of the marker for subject \(i\) at time \(t\). In the current-value link, the relevant part of the conditional hazard model has the form

$$ h_i(t\mid b_i)=h_0(t)\exp\{\gamma^\mathsf{T}w_i+\alpha m_i(t)\}. $$

Here \(h_i(t\mid b_i)\) is the event hazard conditional on the subject’s random effects \(b_i\), \(h_0(t)\) is the baseline hazard, \(w_i\) collects other hazard predictors, and \(\gamma\) contains their coefficients. The parameter \(\alpha\) is the current-value association coefficient. This equation specifies a linear contribution of the latent marker to the log hazard; it is that specification that gives \(\alpha\) its simple per-unit interpretation.

Definition: In a current-value link, \(\alpha\) is the change in the conditional log hazard associated with a one-unit increase in the latent marker value \(m_i(t)\), with time and the other hazard predictors held fixed. The corresponding hazard ratio is \(\exp(\alpha)\) for a one-unit increase.

The comparison is conditional: it is made at the same time \(t\), with the same values of the other hazard predictors and the same baseline hazard. It describes a relative instantaneous event rate, not a difference in event probabilities. If two subjects have different latent marker values, their hazards may also differ for reasons represented by other predictors; the interpretation of \(\alpha\) isolates the marker contrast while holding those predictors fixed.

From the Log Hazard to a Hazard Ratio

The association parameter appears inside an exponential. To interpret it, compare two otherwise identical conditional hazards at time \(t\), one with latent marker value \(m\) and the other with value \(m+\Delta\). Their ratio is

$$ \frac{h(t\mid m+\Delta,\text{ other predictors})} {h(t\mid m,\text{ other predictors})} = \frac{h_0(t)\exp\{\gamma^\mathsf{T}w+\alpha(m+\Delta)\}} {h_0(t)\exp\{\gamma^\mathsf{T}w+\alpha m\}} = \exp(\alpha\Delta). $$

The baseline hazard and the contributions of the other predictors cancel. Taking logarithms gives a log-hazard difference of \(\alpha\Delta\). For a one-unit increase, \(\Delta=1\), so the hazard ratio is \(\exp(\alpha)\). For a difference of \(d\) marker units, it is \(\exp(\alpha d)\). Thus \(\alpha\) itself is not a hazard ratio; it is a log-hazard change per marker unit.

A positive \(\alpha\) means that higher latent marker values are associated with a higher conditional hazard under this model. A negative \(\alpha\) means that higher values are associated with a lower conditional hazard. If \(\alpha=0\), the specified model has no current-value contribution to the conditional log hazard. This does not establish that the marker and event time are unrelated in every sense: other associations, omitted predictors, or misspecification may still matter.

Worked Example: Interpreting a Positive Coefficient

Suppose a hypothetical joint model uses a biomarker measured in milligrams per litre and estimates \(\alpha=0.35\) per mg/L. Interpret the association for a one-unit increase and for a 2 mg/L difference.

For a 1 mg/L increase, the log-hazard difference is \(0.35\times1=0.35\). The hazard ratio is

$$ \exp(0.35)=1.419. $$

Conditional on the same time and other hazard predictors, a subject whose latent biomarker value is 1 mg/L higher has an estimated hazard about 1.42 times as large, or about 41.9% higher. This is a relative instantaneous hazard comparison, not a 41.9 percentage-point increase in event probability.

For a 2 mg/L difference, the log-hazard difference is \(0.35\times2=0.70\), so the hazard ratio is

$$ \exp(0.70)=2.014. $$

Under the model, that 2 mg/L latent-value difference corresponds to an estimated conditional hazard about 2.01 times as large, with the other predictors and time held fixed. The two-unit comparison is not obtained by doubling the one-unit hazard ratio; it is obtained by exponentiating twice the log-hazard change.

Units, Scale, and Meaningful Contrasts

The numerical value of \(\alpha\) depends on the units used to express the marker. If the marker is recorded in mg/L, \(\alpha\) is a log-hazard change per mg/L. If the same values are instead recorded in µg/L, the coefficient per µg/L is one-thousandth as large, because 1 mg/L equals 1000 µg/L. The corresponding hazard ratio for the same physical marker difference does not change when the conversion is done consistently.

This makes it useful to report hazard ratios for interpretable differences, such as 5 units or one standard deviation, as well as for a one-unit increase when that unit is meaningful. For any chosen difference \(d\), calculate \(\exp(\alpha d)\), and state what \(d\) represents. Do not describe \(\exp(\alpha)\) as the effect of a “typical” increase unless one marker unit actually has that meaning.

Worked Example: Rescaling the Marker

Suppose the hypothetical estimate is \(\alpha=0.004\) per µg/L. A difference of 250 µg/L is equivalent to 0.25 mg/L. Find the hazard ratio for this difference.

Using the coefficient in its stated units, the log-hazard difference is

$$ 0.004\ \text{per µg/L}\times250\ \text{µg/L}=1.00. $$

Therefore the hazard ratio is \(\exp(1.00)=2.718\). Equivalently, the coefficient per mg/L is \(0.004\times1000=4.00\), and the same contrast is 0.25 mg/L, giving \(4.00\times0.25=1.00\) on the log-hazard scale. Both calculations yield the same hazard ratio. The units cancel in the product \(\alpha d\), as they must.

A change in the marker’s origin has a different role from a change in its units. Replacing \(m_i(t)\) by \(m_i(t)-c\) shifts the marker’s reference value and changes the corresponding intercept contribution, but it does not change the slope \(\alpha\) or the hazard ratio for a fixed marker difference. It is the difference between values, not the chosen zero point, that determines \(\exp(\alpha\Delta)\).

Uncertainty and Interpretation in Practice

An estimated association coefficient is subject to sampling uncertainty. If an estimated coefficient is approximately normally distributed with standard error \(SE(\hat{\alpha})\), a Wald confidence interval on the log-hazard scale is \(\hat{\alpha}\pm z_{1-\alpha_{\mathrm{level}}/2}SE(\hat{\alpha})\), where \(\alpha_{\mathrm{level}}\) denotes the chosen significance level. Exponentiating both endpoints gives a confidence interval for the one-unit hazard ratio. For a \(d\)-unit contrast, multiply both coefficient endpoints by \(d\) before exponentiating when \(d>0\).

Worked Example: Coefficient and Hazard-Ratio Uncertainty

In a hypothetical analysis, \(\hat{\alpha}=0.28\) per mg/L and \(SE(\hat{\alpha})=0.10\) per mg/L. Use a 95% Wald confidence interval and interpret the estimate.

The estimated coefficient is 2.8 standard errors above zero, since \(0.28/0.10=2.80\). Under the usual large-sample normal approximation, the two-sided Wald test has \(z=2.80\) and \(p\approx0.0051\). The 95% coefficient interval is

$$ 0.28\pm1.96(0.10) = 0.28\pm0.196 = (0.084,\ 0.476) \quad\text{per mg/L}. $$

For a one mg/L increase, exponentiate the estimate and both interval endpoints:

$$ \exp(0.28)=1.323,\qquad \bigl(\exp(0.084),\exp(0.476)\bigr)=(1.088,\ 1.610). $$

The estimated conditional hazard is 1.32 times as large per 1 mg/L higher latent marker value; the 95% Wald confidence interval for that hazard ratio is 1.09 to 1.61. Under the model and the normal approximation, the data are inconsistent with a zero association at the 5% two-sided level. The interval is for the model parameter transformed to the hazard-ratio scale; it is not a range containing 95% of individual hazard ratios.

The interpretation relies on more than the arithmetic. The current-value term must be an adequate representation of how the latent marker relates to log hazard: a one-unit change has the same log-hazard increment across the marker range, and that increment is represented by the same \(\alpha\) over time. The other parts of the event model must also be appropriate for the comparison, and the longitudinal model must provide a meaningful estimate of the latent trajectory. If these conditions are doubtful, the coefficient may be a poor summary even if its standard error is small.

What the Coefficient Does Not Say

A hazard is an instantaneous event rate among subjects still at risk at time \(t\); it is not the probability of experiencing the event by \(t\). Consequently, a hazard ratio of 1.42 does not mean that the probability of an event is 42% higher. Absolute event probabilities also depend on the baseline hazard, follow-up time, and the distribution of predictors. Likewise, a positive current-value association does not by itself imply a fixed difference in median survival.

The coefficient is a conditional model association, not a causal effect. A causal interpretation would require additional design and identification assumptions that are not supplied by the current-value link itself. Further, \(m_i(t)\) is the latent marker value in the model, not simply a single recorded measurement substituted into the hazard equation. Observed measurements contain measurement error and may not equal the underlying value. The distinction matters when explaining what contrast the fitted coefficient represents.

The current-value link is also not interchangeable with every other way of associating longitudinal and event processes. A term based on the marker’s rate of change, an accumulated marker history, or shared random effects answers a different conditional comparison. Their coefficients need not have the same units or interpretation as \(\alpha\). The current-value coefficient specifically describes the log-hazard change per unit difference in the marker value at time \(t\), under the stated model.

Diagnostics and Common Pitfalls

The interpretation is most credible when the assumed relationship is plausible over the observed marker range and follow-up period. Useful checks include examining whether the longitudinal model describes the observed marker patterns adequately, assessing whether the fitted event model shows systematic lack of fit, and considering whether a constant linear current-value contribution is reasonable. A pattern of poor fit at high marker values, for example, could indicate that one log-hazard increment per marker unit is too restrictive. The association estimate should also be treated cautiously when data provide little information about parts of the marker range or time period to which an interpretation is being extended.

  • Wrong: “\(\alpha=0.35\) means a 0.35% increase in risk.” Correction: It is a 0.35-unit increase in log hazard per marker unit; the one-unit hazard ratio is \(\exp(0.35)\).
  • Wrong: “A hazard ratio of 1.42 means 42% more patients will have the event.” Correction: It is a conditional relative instantaneous hazard, not a probability or event-count ratio.
  • Wrong: “The association proves that raising the marker causes the hazard to rise.” Correction: The coefficient is a conditional model association; causal claims require additional assumptions.
  • Wrong: “A one-unit increase always means the same physical change.” Correction: State the marker units and the contrast. Rescaling changes the numerical coefficient.
  • Wrong: “The coefficient applies to an observed reading at any time.” Correction: The current-value link concerns the model’s latent marker value at that time, not an error-free interpretation of a single noisy observation.
Key takeaway: For a latent marker difference \(\Delta\), the current-value model implies a conditional hazard ratio of \(\exp(\alpha\Delta)\). Interpret \(\alpha\) on the log-hazard scale, report the marker units, and keep the comparison distinct from absolute risk and causal effect.

Check Your Understanding

For each question, give the calculation and a contextual interpretation. The solutions are included below.

  1. A hypothetical model estimates \(\alpha=-0.20\) per unit of a latent marker. What is the one-unit hazard ratio, and how should its direction be described?
  2. For \(\alpha=0.15\) per unit, calculate the hazard ratio for a 3-unit difference.
  3. A coefficient is reported as 0.006 per µg/L. What is the hazard ratio for a difference of 100 µg/L?
  4. In one sentence, explain why a hazard ratio of 1.5 is not an event probability 50% higher.

Solutions. (1) \(\exp(-0.20)=0.819\). At the same time and with other hazard predictors held fixed, a one-unit higher latent marker value corresponds to an estimated conditional hazard about 0.82 times as large, or 18.1% lower. (2) The log-hazard difference is \(0.15\times3=0.45\), so the hazard ratio is \(\exp(0.45)=1.568\); the 3-unit higher latent value corresponds to a conditional hazard about 1.57 times as large. (3) The log-hazard difference is \(0.006\times100=0.60\), so the hazard ratio is \(\exp(0.60)=1.822\) for the 100 µg/L contrast. (4) A hazard ratio compares instantaneous hazards among those at risk at a given time, whereas an event probability depends on the baseline hazard and follow-up time as well.

Clinical Biostats