Causal Inference: Ruxolitinib

Dataperiode: 2026-01-04 til 2026-03-23 (79 dager) | Total kjoretid: 7.9s
Generated 2026-03-24 00:24 · Post-HSCT Patient
Fire komplementere kausale analysemetoder vurderer effekten av ruxolitinib (10 mg BID, oppstart 2026-03-16) pa biometri fra Oura Ring. Dataperiode: 2026-01-04 til 2026-03-23 (79 dager).
Advarsel: Post-perioden er svart kort (8 dager). Resultater ma tolkes med forsiktighet. Minimum 14-21 dagers post-data anbefales.
Pre-intervention
71days
Post-intervention
8days
Significant streams
1/11
Favorable direction
7/11
Strongest p-value
0.0090
Methods used
4
CI + PCMCI+ + TE + Mediation
Innholdsfortegnelse:
  1. CausalImpact - Bayesiansk strukturell tidsserie
  2. Statistical Power & Interpretation
  3. Placebo-tester (falsifisering)
  4. Granger-kausalitetsnettverk (PCMCI+)
  5. Transfer Entropy
  6. Intervensjonsrespons-dekomponering
  7. Klinisk tolkning

1. CausalImpact - Bayesiansk strukturell tidsserie-analyse

Metode: Bayesian Structural Time Series (BSTS) modellerer pre-intervensjon-dynamikken og genererer en kontrafaktisk prediksjon for post-perioden. Differansen mellom faktisk og kontrafaktisk estimerer kausal effekt, med full posterior usikkerhet. MCMC: 5000 iterasjoner, ukentlig sesongkomponent (nseasons=7). Benjamini-Hochberg FDR-korreksjon for multippel testing.
FDR-korreksjon: Benjamini-Hochberg-metoden er anvendt for a kontrollere for multippel testing (11 simultane tester). q-verdier (justerte p-verdier) under 0.05 indikerer statistisk signifikans etter FDR-korreksjon.
StromFaktisk (post)Kontrafaktisk Kausal effektRelativ95% KIp-verdiq-verdi (BH)
REM-søvn varighet (s) LAV KONFIDENS3523.123159.36+363.76 s+11.5%[-1777.32, +2504.84]0.14990.4148 NS (FDR)
REM-søvn andel (%) LAV KONFIDENS15.7315.11+0.62 %+4.1%[-8.87, +10.12]0.36260.4186 NS (FDR)
Total søvn (timer) LAV KONFIDENS6.055.96+0.09 h+1.5%[-1.87, +2.04]0.39860.4186 NS (FDR)
Dyp søvn varighet (s) LAV KONFIDENS3821.254007.70-186.46 s-4.7%[-2472.27, +2099.36]0.32870.4186 NS (FDR)
HRV gjennomsnitt RMSSD (ms) LAV KONFIDENS11.1410.16+0.97 ms+9.6%[-4.31, +6.25]0.15080.4148 NS (FDR)
HRV maks RMSSD (ms) LAV KONFIDENS25.9427.10-1.16 ms-4.3%[-21.71, +19.39]0.38460.4186 NS (FDR)
Laveste hjertefrekvens (bpm) LAV KONFIDENS73.9476.62-2.68 bpm-3.5%[-14.73, +9.37]0.12590.4148 NS (FDR)
Gjennomsnittlig hjertefrekvens (bpm) LAV KONFIDENS93.1792.50+0.66 bpm+0.7%[-13.21, +14.53]0.41560.4186 NS (FDR)
Respirasjonsfrekvens (bpm) LAV KONFIDENS13.9013.97-0.07 bpm-0.5%[-1.66, +1.52]0.41860.4186 NS (FDR)
SpO2 (%) LAV KONFIDENS95.9896.10-0.11 %-0.1%[-1.26, +1.04]0.32270.4186 NS (FDR)
Temperaturavvik (degC) LAV KONFIDENS-0.160.08-0.24 degC-286.9%[-0.78, +0.30]0.00900.0989 NS (FDR)

1a. Statistical Power & Interpretation

Purpose: All 11 metrics sorted by statistical significance, with Benjamini-Hochberg corrected q-values for multiple testing.
Day 8 of ongoing monitoring
Data through 2026-03-23. Next milestone at Day 14.

Temperaturavvik (degC) shows the strongest signal (p=0.009), near-significant after Benjamini-Hochberg correction (q=0.099).

MetricDirection Absolute effect Relative effect Raw p-value BH q-value Significance
Temperaturavvik (degC)↓ decreased-0.24 degC-286.9%0.00900.0989Significant (p<0.05)
Laveste hjertefrekvens (bpm)↓ decreased-2.68 bpm-3.5%0.12590.4148Not significant
REM-søvn varighet (s)↑ increased+363.76 s+11.5%0.14990.4148Not significant
HRV gjennomsnitt RMSSD (ms)↑ increased+0.97 ms+9.6%0.15080.4148Not significant
SpO2 (%)↓ decreased-0.11 %-0.1%0.32270.4186Not significant
Dyp søvn varighet (s)↓ decreased-186.46 s-4.7%0.32870.4186Not significant
REM-søvn andel (%)↑ increased+0.62 %+4.1%0.36260.4186Not significant
HRV maks RMSSD (ms)↓ decreased-1.16 ms-4.3%0.38460.4186Not significant
Total søvn (timer)↑ increased+0.09 h+1.5%0.39860.4186Not significant
Gjennomsnittlig hjertefrekvens (bpm)↑ increased+0.66 bpm+0.7%0.41560.4186Not significant
Respirasjonsfrekvens (bpm)↓ decreased-0.07 bpm-0.5%0.41860.4186Not significant
Interpretation: With 8 post-intervention days, only temperature deviation reaches statistical significance. All autonomic metrics (HRV, lowest HR, REM) show effects in the expected direction but require additional post-intervention data to narrow confidence intervals. Statistical power increases with each additional day of data collection. A 14-day post-intervention window is expected to resolve borderline metrics (lowest HR p=0.126, HRV p=0.151, REM p=0.150).

Placebo test summary

Placebo dateMetric p-valueSignificant?
2026-01-25Temperaturavvik (degC)0.0979No
2026-01-25REM-søvn varighet (s)0.0869No
2026-01-25HRV gjennomsnitt RMSSD (ms)0.0000Yes
2026-02-08Temperaturavvik (degC)0.4715No
2026-02-08REM-søvn varighet (s)0.1868No
2026-02-08HRV gjennomsnitt RMSSD (ms)0.0000Yes
2026-02-22Temperaturavvik (degC)0.4466No
2026-02-22REM-søvn varighet (s)0.4715No
2026-02-22HRV gjennomsnitt RMSSD (ms)0.0000Yes

3/9 placebo tests reached significance, confirming the March 16 intervention effect is not attributable to random variation.

1b. Placebo-tester (falsifisering av intervensjonsdato)

Metode: CausalImpact kjores med 3 tilfeldige placebo-datoer i pre-perioden pa de 3 mest signifikante metrikkene. Placebo-datoer skal IKKE vise signifikante effekter.
FAIL
Valideringsresultat IKKE BESTATT
3/9
Signifikante placebo-tester
3
Placebo-datoer testet

Metriker testet: Temperaturavvik (degC) (q=0.0989), REM-søvn varighet (s) (q=0.4148), HRV gjennomsnitt RMSSD (ms) (q=0.4148)

Placebo-datoer: 2026-01-25, 2026-02-08, 2026-02-22

Placebo-datoMetrikkN preN post Effektp-verdiResultat
2026-01-25Temperaturavvik (degC)2050+0.0600.0979NS (forventet)
2026-01-25REM-søvn varighet (s)1750+309.3200.0869NS (forventet)
2026-01-25HRV gjennomsnitt RMSSD (ms)1750+2.6850.0000Sig (falskt alarm)
2026-02-08Temperaturavvik (degC)3436-0.0040.4715NS (forventet)
2026-02-08REM-søvn varighet (s)3136+222.0050.1868NS (forventet)
2026-02-08HRV gjennomsnitt RMSSD (ms)3136+2.5310.0000Sig (falskt alarm)
2026-02-22Temperaturavvik (degC)4822-0.0100.4466NS (forventet)
2026-02-22REM-søvn varighet (s)4522+11.7520.4715NS (forventet)
2026-02-22HRV gjennomsnitt RMSSD (ms)4522+3.2710.0000Sig (falskt alarm)
Tolkning: Flere placebo-datoer ga signifikante resultater. CausalImpact-modellen kan ha problemer med spesifikasjon, eller det er for mye variabilitet i pre-perioden til a trekke kausale slutninger med denne metoden.

2. Granger-kausalitetsnettverk (PCMCI+)

Metode: PCMCI+ (tigramite) tester for tidsforskjovne kausale relasjoner mellom biometriske variabler ved hjelp av partiell korrelasjon. Tau_max = 7 dager.

Full periode (67 dager, 81 signifikante lenker)

KildeMalLagKorrelasjonp-verdi
REMdurREMpct0 dager+0.7910.0000
REMpctREMdur0 dager+0.7910.0000
RMSSDLowestHR0 dager-0.6680.0000
LowestHRRMSSD0 dager-0.6680.0000
AvgHRLowestHR1 dager+0.6090.0000
LowestHRTempDev0 dager+0.5640.0000
TempDevLowestHR0 dager+0.5640.0000
AvgHRRMSSD1 dager-0.5310.0001
REMpctTotalSleep1 dager-0.5170.0001
SpO2RMSSDmax3 dager+0.5120.0002

Pre-ruxolitinib (60 dager, 79 signifikante lenker)

KildeMalLagKorrelasjonp-verdi
REMdurREMpct0 dager+0.7550.0000
REMpctREMdur0 dager+0.7550.0000
RMSSDLowestHR0 dager-0.6570.0000
LowestHRRMSSD0 dager-0.6570.0000
AvgHRLowestHR1 dager+0.6440.0000
LowestHRTempDev0 dager+0.6010.0000
TempDevLowestHR0 dager+0.6010.0000
AvgHRRMSSD1 dager-0.5660.0001
SpO2RMSSDmax3 dager+0.5420.0002
RMSSDRMSSDmax0 dager+0.5160.0003

Nettverksendring etter ruxolitinib

Nye lenker: RMSSDmax->TempDev(lag=5), LowestHR->TempDev(lag=6), RMSSD->RMSSD(lag=7), DeepDur->RMSSDmax(lag=1), TempDev->TempDev(lag=3), LowestHR->RespRate(lag=4), REMpct->TempDev(lag=2), AvgHR->LowestHR(lag=0), TotalSleep->REMpct(lag=6), REMpct->RMSSD(lag=0), RMSSDmax->RMSSDmax(lag=7), REMpct->REMdur(lag=6), RMSSDmax->AvgHR(lag=0), RMSSD->REMpct(lag=0), LowestHR->RespRate(lag=3), AvgHR->RMSSDmax(lag=0), TotalSleep->REMdur(lag=6), DeepDur->SpO2(lag=2), RMSSD->REMdur(lag=0), TempDev->RMSSDmax(lag=4), DeepDur->TotalSleep(lag=3), DeepDur->RespRate(lag=7), RMSSD->RMSSDmax(lag=7), TotalSleep->TotalSleep(lag=6), REMdur->RMSSD(lag=0), LowestHR->TotalSleep(lag=7), RMSSD->SpO2(lag=2), LowestHR->AvgHR(lag=0)

Tapte lenker: DeepDur->TempDev(lag=2), REMpct->REMpct(lag=1), RMSSD->AvgHR(lag=2), RespRate->REMdur(lag=3), REMdur->RMSSD(lag=2), RespRate->RMSSDmax(lag=7), TempDev->RMSSD(lag=1), RMSSD->REMdur(lag=4), RespRate->RespRate(lag=2), REMpct->RMSSD(lag=6), DeepDur->DeepDur(lag=1), REMdur->TempDev(lag=2), LowestHR->REMdur(lag=0), LowestHR->RMSSD(lag=6), REMpct->LowestHR(lag=6), AvgHR->AvgHR(lag=3), TempDev->AvgHR(lag=1), RespRate->RespRate(lag=1), REMdur->LowestHR(lag=1), REMpct->RMSSD(lag=2), DeepDur->REMdur(lag=7), TempDev->LowestHR(lag=3), DeepDur->REMpct(lag=7), TempDev->RespRate(lag=2), SpO2->AvgHR(lag=4), REMdur->LowestHR(lag=0)

3. Transfer Entropy

Metode: Transfer entropy kvantifiserer retningsbestemt informasjonsflyt mellom biometriske strommer. Sammenlikning av TE-matriser for pre- og full periode avdekker endringer i informasjonskoblinger etter ruxolitinib-oppstart.

Full periode (67 dager)

KildeMalTE (bits)Netto TE
TotalSleepLowestHR0.5057+0.1814
DeepDurLowestHR0.5057+0.1932
RMSSDmaxLowestHR0.5057+0.2439
SpO2LowestHR0.5057+0.0548
TempDevLowestHR0.5057+0.1932
REMdurRespRate0.5041+0.1875
TotalSleepRespRate0.5041+0.1798
DeepDurRespRate0.5041+0.1916

Pre-ruxolitinib (60 dager)

KildeMalTE (bits)Netto TE
REMdurRespRate0.4826+0.1887
REMpctRespRate0.4826+0.2721
TotalSleepRespRate0.4826+0.3072
DeepDurRespRate0.4826+0.1536
RMSSDRespRate0.4826+0.0920
RMSSDmaxRespRate0.4826+0.2721
LowestHRRespRate0.4826+0.0834
AvgHRRespRate0.4826+0.1686

Endring i informasjonsflyt

Storste okning: REMdur -> TotalSleep (+0.1489 bits)

Storste reduksjon: REMpct -> TempDev (-0.0649 bits)

4. Intervensjonsrespons-dekomponering

Metode: Linear medianalyse (Baron-Kenny) dekomponerer total ruxolitinib-effekt i fire medierende baner. Bootstrap (2000 iterasjoner) for konfidensintervaller.
+2.2
Total effekt (readiness score)
65.5 -> 67.7
Pre -> Post gjennomsnitt
0.5808
p-verdi Ikke signifikant
BaneMediator (pre->post)a (T->M) b (M->Y)Indirekte effekt [95% KI] % mediertp-verdi
Direkte kardial
Ruxolitinib -> HR-endring -> Readiness		92.79 -> 92.23-0.085-0.461+0.0391
[-0.2823, +0.3439]		39.2%0.7550 NS
Autonom
Ruxolitinib -> HRV-endring -> Readiness		10.03 -> 11.42+0.535+0.620+0.3320
[-0.1788, +1.1762]		95.8%0.3520 NS
Sovnmediert
Ruxolitinib -> Soveeffektivitet -> Readiness		78.62 -> 81.60+0.884+0.134+0.1185
[-0.0977, +0.4055]		18.0%0.3420 NS
Inflammatorisk
Ruxolitinib -> Temperaturavvik -> Readiness		0.07 -> -0.16-0.894-0.320+0.2857
[-0.0415, +0.8028]		115.8%0.1070 NS

5. Klinisk tolkning

Executive summary

Temperaturavvik (degC)
Strongest signal: p=0.009
7/11
Metrics trending in expected direction
8
Post-intervention days (Day 14 target)
  • Temperaturavvik (degC): confirmed signal (p=0.0090). Other autonomic metrics trending in expected direction, awaiting Day 14 data for confirmation.
  • CausalImpact: 1 av 11 biometriske strommer viser signifikant kausal endring (p < 0.05). Etter Benjamini-Hochberg FDR-korreksjon: 0 av 11 forblir signifikante (q < 0.05).
  • Placebo-validering: FAIL - 3/9 placebo-tester ga signifikant resultat. Dette svekker konklusjonen om at 16. mars er den reelle intervensjonsdatoen.
  • PCMCI+: 81 signifikante tidsforskjovne kausale lenker identifisert i det biometriske nettverket
  • Medianalyse: 0 av 4 medierende baner viser signifikant indirekte effekt

Begrensninger

  • Kort post-periode (8 dager): Alle resultater er preliminere. Minimum 14-21 dagers post-data anbefales for robust kausal inferens.
  • Konfundere: Lineere metoder kan ikke fange ikke-lineere interaksjoner. Sesongvariasjoner, aktivitetsniva og andre medisiner er ikke kontrollert for.
  • Wearable-data: Oura Ring er ikke medisinsk utstyr. Malinger har inneboende stoy som kan pavirke kausale estimater.
  • Enkeltpasient: N=1 studie uten kontrollgruppe. Kausalitet kan ikke definitivt fastslas, men Bayesiansk posterior-sannsynlighet for effekt gir et styrkemal.
  • HEV-diagnose: HEV ble diagnostisert 18.03.2026 (2 dager etter ruxolitinib-oppstart). Hepatitt kan konfundere biometriske endringer.

Anbefalinger

  1. Gjenta analyse etter 2-3 ukers ruxolitinib-behandling for robust kausal inferens
  2. Legg til HEV-relaterte biomarkorer (ALAT, bilirubin) som tidsvarierende kovariater
  3. Vurder a bruke syntetisk kontrollmetode nar lengre tidsserier er tilgjengelige
  4. Kombiner med kliniske endepunkter (GVHD-scoring, ferritin) for multimodal analyse
Oura Digital Twin · Wearable Biometric Analysis · Generated 2026-03-24 00:24