GvHD Prediction Model

Oura Ring, 2026-01-08 to 2026-03-23 (74 days) | Ruxolitinib started 2026-03-16 | HEV diagnosed 2026-03-18
Generated 2026-03-24 00:24 · Post-HSCT Patient
Peak GVHD Composite
72.9
on 2026-03-02
rSLDS Event Detection
PARTIAL
Feb 9 classified as: Recovery
Alert Lead Time
30d
RED alerts before Feb 9 event
Combined GVHD + BOS
36.6
LOW
Clinical Context: Chronic GVHD affecting skin, liver, mouth (NIH 2014: moderate). Known acute decompensation on Feb 9, 2026 (validation target). Temperature deviation from Oura is the primary inflammatory signal. This analysis tests whether wearable data can provide early warning of GVHD flares.

1. Temperature Fluctuation Analysis

Temperature deviation from Oura measures deviations from personal baseline body temperature. GVHD causes micro-inflammatory signatures visible as increased night-to-night variability. CUSUM change-point detection identifies regime shifts in temperature dynamics.

Regime changes detected: 2026-01-31, 2026-02-11, 2026-03-02

2. Multi-Stream GVHD Composite Score

Six biometric streams weighted by clinical relevance: temperature (25%), HRV (20%), SpO2 (15%), sleep fragmentation (15%), resting HR (15%), activity (10%). Each component z-scored against first 14 days baseline, normalized to 0-100 (higher = more GVHD-like).

3. Recurrent Switching Linear Dynamical System

4-state rSLDS (Remission, Pre-flare, Active Flare, Recovery) fitted via Laplace-EM (Linderman et al. 2017). Each discrete state governs linear dynamics in a continuous latent space, with recurrent transitions that depend on the latent state. Observations: composite score + temperature deviation + HRV (3 features). Falls back to Gaussian HMM (hmmlearn) if ssm library is unavailable.

rSLDS State Distribution

StateDaysPercentage
Remission34.0%
Pre-flare2736.0%
Active Flare2330.7%
Recovery2229.3%

4. Early Warning Alerts

YELLOW alert: P(pre-flare) > 0.3 for 3+ consecutive days.
RED alert: P(pre-flare) > 0.5 OR P(active flare) > 0.2.
Retrospective validation against Feb 9, 2026 acute decompensation.

Alert History (50 total)

DateLevelP(Pre-flare)P(Flare)Composite
2026-01-10RED0.01.059.4
2026-01-11RED0.01.058.8
2026-01-12RED0.01.042.1
2026-01-13RED0.01.066.7
2026-01-14RED0.01.040.0
2026-01-15RED1.00.052.1
2026-01-16RED1.00.044.5
2026-01-17RED1.00.044.3
2026-01-18RED1.00.055.2
2026-01-19RED1.00.022.1
2026-01-20RED1.00.056.3
2026-01-21RED0.01.054.3
2026-01-22RED0.01.052.2
2026-01-23RED1.00.049.7
2026-01-24RED1.00.039.6
2026-01-25RED1.00.040.4
2026-01-26RED1.00.028.7
2026-01-27RED1.00.036.6
2026-01-28RED1.00.034.5
2026-01-29RED1.00.062.5

5. Predictive Feature Importance

Features ranked by combined score: point-biserial correlation (30%), mutual information (30%), composite correlation (20%), Cohen's d effect size (20%). Binary target: composite score > 65th percentile.

Top Predictive Features for GVHD State

FeatureImportanceCorrelationMutual InfoCohen's d
HRV Median (RMSSD)0.520-0.5480.266-1.45
Sleep Heart Rate0.503+0.5940.124+1.46
Lowest HR0.502+0.5710.170+1.41
Readiness Score0.353-0.4100.079-0.89
Temperature Deviation0.282+0.3260.127+0.65
Temperature Gradient0.209+0.2550.064+0.53
HRV Variability0.200-0.2440.044-0.57

6. BOS Risk Integration

Bronchiolitis obliterans syndrome (BOS) is the pulmonary manifestation of chronic GVHD. Existing BOS risk score (24.9) integrated with systemic GVHD composite. SpO2 trend analysis provides supplementary pulmonary assessment.
Disclaimer: This analysis is for research purposes only and should not be used as a sole basis for clinical decisions. N=1 retrospective case study - all detection metrics are descriptive, not inferential. Sensitivity/specificity cannot be computed from a single patient. Validation requires an external multi-patient cohort. Temperature deviation from consumer wearables has limited precision compared to clinical thermometry. All clinical decisions should be made in consultation with the treating hematologist.