Priment Statistics and Methodology seminar
Department of Statistical Science, UCL
Resources
Slides and code here: github.com/n8thangreen/HCM-ICD-cost-eff-model-talk
HCM is the most common heritable cardiac condition and is the most common cause of sudden cardiac arrest in the young, affecting about 1 in 500 people. Familial hypertrophic cardiomyopathy is a heart condition characterized by thickening (hypertrophy) of the heart muscle, more specifically the ventricle. The thickened heart muscle, can make it challenging to keep up with the oxygenation demands of the body and the heart muscle itself. Many people with HCM have few, if any, symptoms and can lead normal lives without significant symptoms. However, this condition can also have serious consequences. Life threatening arrhythmias resulting in cardiac arrest can sometimes be the first symptom.
Symptoms include: Shortness of breath, Chest pain, Fainting, Palpitations
An ICD is a small battery-powered device placed under the skin—usually near the collarbone. It constantly monitors your heart rhythm and delivers life-saving electrical shocks or pacing to correct dangerously fast or irregular heartbeats (arrhythmias), preventing sudden cardiac arrest
O’Mahony et al. (2014)
Hypertrophic cardiomyopathy (HCM) is a leading cause of sudden cardiac death (SCD) in young adults. Current risk algorithms provide only a crude estimate of risk and fail to account for the different effect size of individual risk factors.
The aim of this study was to develop and validate a new SCD risk prediction model that provides individualized risk estimates.
The prognostic model was derived from a retrospective, multi-centre longitudinal cohort study. The model was developed from the entire data set using the Cox proportional hazards model and internally validated using bootstrapping.
The cohort consisted of 3675 patients from six centres. During a follow-up period of 24,313 patient-years (median 5.7 years), 198 patients (5%) died suddenly or had an appropriate ICD shock.
Of eight pre-specified predictors, associated with SCD/appropriate ICD shock at the 15% significance level age, maximal left ventricular wall thickness, left atrial diameter, left ventricular outflow tract gradient, family history of SCD, non-sustained ventricular tachycardia, and unexplained syncope
final model to estimate individual probabilities of SCD at 5 years.
Define the Prognostic Index (PI) using the coefficients from the standard Cox model
\[ \begin{align*} \text{PI}_{\text{Naive}} &= 0.159 \times \text{Maximal wall thickness} \\ &\quad - 0.003 \times \text{Maximal wall thickness}^2 \\ &\quad + 0.026 \times \text{Left atrial diameter} \\ &\quad + 0.004 \times \text{Maximal LVOT gradient} \\ &\quad + 0.458 \times \text{Family history SCD} \\ &\quad + 0.826 \times \text{NSVT} \\ &\quad + 0.717 \times \text{Unexplained syncope} \\ &\quad - 0.018 \times \text{Age} \end{align*} \]
Risk Calculator
\[ \begin{equation*} \hat{P}_{\text{SCD at 5 years}} = 1 - 0.998^{\exp(\text{PI})} \end{equation*} \]
Green et al. (2024)
A discrete-time Markov model was used to determine the cost-effectiveness of different ICD decision-making rules for implantation.
Several scenarios were investigated, including the reference scenario of implantation rates according to observed real-world practice.
A 12-year time horizon with an annual cycle length was used.
Transition probabilities used in the model were obtained using Bayesian analysis.
\[ \begin{pmatrix} p_{11}^s & p_{12}^s & 0 & 0 & p_{15}^s \\ 1 - p_{15}^s & 0 & 0 & 0 & p_{15}^s \\ 0 & 0 & p_{33} & p_{34}^s & p_{35}^s \\ 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 & 1 \end{pmatrix} \]
Denote \(x\) as the observed number of transitions, \(p\) the probability of a transition and \(n\) as the total number of transitions from a given state. The hyperparameters \(\alpha\) characterise the prior knowledge on \(p\). Superscripts indicate the decision rule used.
Likelihood
\[ \begin{align*} x_{i.}^{(1)} &\sim \text{Multinomial}\left(p_{i.}^{(1)}, n_{i}^{(1)}\right), \quad i = 1, 3 \\ x_{i.}^{(2)} &\sim \text{Multinomial}\left(p_{i.}^{(2)}, n_{i}^{(2)}\right), \quad i = 1, 3 \\[1em] \end{align*} \]
Priors
\[ \begin{align*} p_{i.}^{(1)} &\sim \text{Dirichlet}(\alpha^{(1)}), \quad i = 1, 3 \\ p_{i.}^{(2)} &\sim \text{Dirichlet}(\alpha^{(2)}), \quad i = 1, 3 \end{align*} \]
For all final sink states,
\[ p_{ij}^{(s)} = \begin{cases} 1 & \text{if } i = j \\ 0 & \text{if } i \neq j \end{cases} \]
(main effect / first-order Sobol index)
\[ S_i = \frac{V_{X_i}(E_{\mathbf{X}_{\sim i}}(Y \mid X_i))}{V(Y)} \]
Structural Model Extension: Adapt the existing Markov framework to differentiate between Transvenous (TV-ICD) and Subcutaneous (S-ICD) devices.
Differential Cost-Effectiveness: Evaluate the economic trade-off between the higher upfront capital costs of S-ICDs against the potential long-term savings.
Subgroup Stratification: Younger cohorts?
From Model Uncertainty to Decision Uncertainty
For a specific parameter of interest (\(X_i\)) is calculated as:
EVPPI
\[ \text{EVPPI}_{X_i} = E_{X_i} \left[ \max_{d \in \{0, 1\}} E_{\boldsymbol{X}_{\sim i}} \left[ \text{NMB}(d, X_i, \boldsymbol{X}_{\sim i}) \right] \right] - \max_{d \in \{0, 1\}} E_{\boldsymbol{\theta}} \left[ \text{NMB}(d, \boldsymbol{\theta}) \right] \]
Fine-Gray Model (Subdistribution Hazards)
\[ \begin{equation*} \text{CIF}_{\text{SCD}}^{\text{FG}}(5) = 1 - \exp\left(-\hat{H}_{\text{SCD}}(5) \exp(\text{PI}_{\text{FG}})\right) \end{equation*} \]
Cause-Specific Model
\[ % \hat{h}_{SCD}(t) and \hat{h}_{Other}(t) are the baseline hazards at time t. % PI_{cs,SCD} and PI_{cs,Other} are the linear predictors from their respective Cox models. \begin{equation*} \text{CIF}_{\text{SCD}}^{\text{cs}}(5) = \int_{0}^{5} \hat{h}_{\text{SCD}}(u) \exp(\text{PI}_{\text{cs,SCD}}) \exp\left(-\int_{0}^{u} \left[ \hat{h}_{\text{SCD}}(s) \exp(\text{PI}_{\text{cs,SCD}}) + \hat{h}_{\text{Other}}(s) \exp(\text{PI}_{\text{cs,Other}}) \right] ds\right) du \end{equation*} \]
Individualized Risk Stratification: The novel risk prediction model provides individualized 5-year probabilities for Sudden Cardiac Death, offering a more precise tool for clinical decision-making than previous crude guidelines.
Robust Economic Evaluation: By embedding this prediction model within a Bayesian discrete-time Markov framework, we can formally quantify the cost-effectiveness of different ICD implantation thresholds under real-world uncertainty.
Future Work:
Nathan Green | UCL | n.green@ucl.ac.uk