Why the Distinction Matters
A biomarker measured during follow-up is not usually observed without error. A laboratory value can vary because of assay imprecision, short-term biological fluctuation, or specimen handling, even when the subject’s underlying marker level has changed little. If an event hazard is meant to reflect the subject’s evolving marker status, treating each recorded value as that status can make the event association respond to noise rather than to the trajectory of interest.
In the basic shared-random-effects joint model, the longitudinal measurements inform an underlying trajectory, and the event submodel relates risk to that trajectory. The distinction is important: a noisy measurement is evidence about the latent marker level, not necessarily the level itself. This tutorial focuses on that distinction, not on the later likelihood calculations.
A standard simple measurement-error specification is \(\epsilon_{ij}\sim N(0,\sigma_\epsilon^2)\), independent across measurements and independent of the event process conditional on the latent trajectory and shared random effects. The mean-zero assumption says that measurements are centered on the latent value; the variance \(\sigma_\epsilon^2\) quantifies residual measurement variability. These are modeling assumptions, not facts guaranteed by collecting repeated measurements.
The latent trajectory is commonly represented by a longitudinal model with subject-specific random effects, for example a subject-specific intercept and slope. Those random effects induce the trajectory’s variation across subjects and are shared with the event submodel. Under the current-value link, for two subjects with the same baseline hazard and other covariates, the hazard ratio at time \(t\) corresponding to latent marker values differing by \(\Delta m\) is \(\exp(\alpha\Delta m)\). Here \(\alpha\) is the log-hazard association per one-unit increase in the latent marker.
How Measurement Noise Changes the Link
Suppose two subjects have the same latent marker value at a given time but their recorded measurements differ because of error. A model that substitutes the recorded values directly into the hazard link would assign different risks to them. The latent-current-value model does not make that change solely because one assay happened to read high and the other low. Their longitudinal records may still provide different evidence about their latent trajectories, but a single residual deviation is not automatically interpreted as a true marker difference.
Worked Example: A Noisy Reading Is Not a One-Unit Trajectory Change
Assume that the current-value association is \(\alpha=0.40\) per marker unit. At a particular time, subject A’s latent marker value is 6 units, while subject B’s is 5 units. For this comparison, suppose their other covariates and baseline hazard are the same. Subject A’s observed value happens to be 9 units, because the measurement error is \(+3\).
Latent-trajectory comparison. The marker difference for the current-value link is \(6-5=1\) unit. The hazard ratio is
Thus, under the model, A’s hazard is about 1.49 times B’s hazard at that time, for the stated one-unit latent difference and with the other comparison conditions held fixed.
Direct substitution of the observed reading. If the observed value 9 were substituted for A’s latent value, the apparent difference would be \(9-5=4\) units and the resulting ratio would be
That larger ratio treats a three-unit measurement error as though it were a true trajectory difference. The calculation does not show that the observed-value model is always inappropriate for every scientific question; it shows that it answers a different question from the latent-current-value model.
There is also a timing problem with directly using observations. Measurements are generally taken at discrete visits, while a hazard is defined over follow-up time. Between visits there may be no current observed value to insert. Carrying the last value forward or interpolating between readings creates an additional rule that can affect the apparent risk relationship. A latent trajectory provides a modeled marker value across time, including between measurement occasions.
A Short Derivation: Averaging a Noisy Plug-In Hazard
The consequences of substituting a noisy value can be made explicit under a Gaussian error model. Consider a hazard multiplier \(\exp(\alpha m_i(t))\) based on the latent value. If instead one uses \(Y_{ij}=m_i(t)+\epsilon_{ij}\), the multiplier becomes \(\exp\{\alpha(m_i(t)+\epsilon_{ij})\}\). Averaging this plug-in multiplier over a new measurement error gives
The second equality uses the moment-generating function of a mean-zero normal random variable with variance \(\sigma_\epsilon^2\). Unless \(\alpha=0\) or \(\sigma_\epsilon^2=0\), the average of the exponentiated noisy value differs from the exponentiated latent value. This is a consequence of the nonlinear exponential link, not evidence that the latent-value model should include this extra factor. It describes what happens when one averages a different, noisy plug-in multiplier over measurement error.
Worked Example: Quantifying the Noise Effect Under a Gaussian Model
Suppose the marker measurement error has standard deviation \(0.50\) units, so \(\sigma_\epsilon^2=0.25\), and the association parameter is \(\alpha=0.40\) per marker unit. The multiplier introduced by averaging the exponentiated noisy plug-in value is
Under these assumptions, the average plug-in multiplier is about 2.02% larger than the multiplier based on the latent value. This is a mathematical comparison of the two expressions, not a general estimate of bias in a fitted joint model. In practice, the impact on an estimated association depends on the measurement schedule, the longitudinal model, and other features of the data.
Repeated Measurements and the Latent Trajectory
Repeated readings help distinguish persistent subject-specific levels and trends from occasion-specific deviations. They do not eliminate measurement error, and an individual trajectory cannot generally be known exactly from a finite set of readings. Instead, the longitudinal submodel uses the observed record to estimate or infer the trajectory under its stated assumptions. The event submodel then uses that latent trajectory as the marker history relevant to the current-value association.
Worked Example: A Trajectory Between Visits
Consider hypothetical measurements for one subject at times 0, 1, and 2 years: \(4.8\), \(6.3\), and \(6.9\) marker units. For illustration, fit a straight line by least squares to these three values. The mean time is \(\bar{t}=1\), and the mean observed marker is \(\bar{Y}=6.0\). The fitted slope is
marker units per year. The fitted intercept is \(\hat{\beta}_0=\bar{Y}-\hat{\beta}_1\bar{t}=6.0-1.05(1)=4.95\) units, so the illustrative fitted trajectory is \(\hat{m}_i(t)=4.95+1.05t\). At year 1, it gives \(\hat{m}_i(1)=6.00\) units, whereas the observed value is 6.3 units. The 0.3-unit difference is a residual from this fitted line, not automatically a persistent increase in the subject’s marker.
For an association of \(\alpha=0.25\) per marker unit, the hazard ratio comparing a latent marker value of 6 units with a value of 5 units is \(\exp\{0.25(6-5)\}=\exp(0.25)=1.284\). Substituting the single observed value 6.3 instead would yield \(\exp\{0.25(6.3-5)\}=\exp(0.325)=1.384\). The simple fitted line is only an illustration; a joint model typically represents uncertainty in trajectories through its longitudinal model and shared random effects rather than treating this hand-fitted line as known without error.
Assumptions, Diagnostics, and Alternatives
The latent-versus-observed distinction is useful only when the measurement and event models are credible. The standard formulation relies on assumptions such as the following:
- The longitudinal model represents the systematic marker pattern adequately, including relevant time trends and subject-to-subject variation.
- Conditional on the latent trajectory and modeled random effects, the measurement errors have the specified mean and variability and do not contain unmodeled systematic patterns.
- The current-value association is an appropriate description of how the latent marker relates to the hazard; it is not automatically a causal effect.
- The shared random-effects structure and conditional independence assumptions are adequate for connecting the longitudinal and event processes.
Diagnostics should target the parts of the model that make the latent trajectory meaningful. Examine longitudinal residuals for changing variance, trends over time, or outliers that suggest an inadequate measurement model. Compare observed marker profiles with fitted trajectories, and assess whether the trajectory form is plausible rather than relying on a single average trend. Consider whether assay precision differs across visits or subjects. A pattern of residuals related to time, risk factors, or impending events may indicate that the assumed measurement error or longitudinal structure is inadequate.
Also inspect the visit schedule and its relation to the event process. The simple argument that repeated readings inform a latent history can fail if visits occur more often when a subject is deteriorating and this visit process is not adequately addressed. Sparse measurements can leave trajectories weakly determined; substantial measurement error can make the latent association imprecise. Neither a sophisticated link nor repeated measurements can recover information that the data do not contain without relying on model assumptions.
| Approach | What enters the event association | Main distinction |
|---|---|---|
| Latent current-value link | The modeled underlying marker value \(m_i(t)\) | Separates systematic trajectory from occasion-specific measurement error; requires a credible longitudinal model. |
| Observed-value plug-in | A recorded marker value, or a value created by carrying forward or interpolating observations | Can make risk depend on measurement noise and on the rule used between visits; may suit a different question about available readings. |
| Measurement-error-corrected analysis | A corrected or calibrated marker representation | Can be useful when error properties are known or estimable, but depends on the correction model and need not represent a full joint model. |
These approaches are not interchangeable merely because each uses a marker in an event analysis. The latent link targets association with the modeled underlying trajectory. An observed-value approach may instead target association with the value available to a clinician at a visit, which can be scientifically relevant but must be specified and interpreted accordingly. Measurement-error correction is another option when its assumptions and data support it. The key is to state which marker quantity the risk association concerns.
Common Pitfall and Practical Interpretation
A common mistake is to say, “The assay was 9, so the subject’s current marker level was 9; therefore the hazard should use 9.” That reasoning equates a noisy observation with the latent quantity and ignores the measurement model. The correction is to treat the recorded value as an observation centered on the latent trajectory, then use the trajectory in the current-value association if that is the model being fitted.
A second mistake is to assume that replacing observations with a latent trajectory removes all uncertainty. It does not: the trajectory is estimated from finite, noisy data, and its adequacy depends on the model and available measurements. The association describes a conditional statistical relationship under that model. It does not establish that changing the biomarker would cause a change in event risk.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- If two subjects have identical latent marker values at time \(t\) but different observed readings, what does the latent-current-value link imply about the part of their hazard comparison attributable to the marker?
- For \(\alpha=0.30\) per marker unit, calculate the hazard ratio associated with a 2-unit difference in latent marker value, holding other comparison conditions fixed.
- With Gaussian measurement error of variance \(0.16\) and \(\alpha=0.50\), calculate the multiplicative factor from averaging the exponentiated noisy plug-in value over measurement error.
- Name one diagnostic that could reveal a problem with the longitudinal measurement or trajectory model, and explain why it matters for the hazard association.
Solutions. (1) The latent-current-value component assigns no marker-based difference because the latent values are equal; other covariates or random effects could still distinguish their hazards. (2) The ratio is \(\exp\{0.30(2)\}=\exp(0.60)=1.822\), so the higher latent value corresponds to a hazard 1.822 times as large under the stated comparison. (3) The factor is \(\exp\{(0.50)^2(0.16)/2\}=\exp(0.0200)=1.020\). (4) For example, a systematic trend in longitudinal residuals over time suggests the fitted trajectory may miss time structure. Since the hazard association uses that trajectory, misspecification can distort the marker history that enters the risk relationship.