Table of Contents
Fetching ...

Semi-supervised Clustering Through Representation Learning of Large-scale EHR Data

Linshanshan Wang, Mengyan Li, Zongqi Xia, Molei Liu, Tianxi Cai

TL;DR

This work tackles learning informative, transferable patient representations from high-dimensional, sparse EHR data with limited labels. It introduces SCORE, a semi-supervised framework that combines a Poisson-adapted latent factor mixture (PALM) with Gaussian variational approximation (GVA) and leverages pre-trained EHR embeddings to learn compact, transferable patient embeddings. The authors establish convergence and error results for a hybrid EM–GVA algorithm in a setting with diverging latent dimensionality, showing that unlabeled data improve accuracy and robustness to label scarcity. Through simulations and an MS disability-phenotyping application, SCORE outperforms existing methods in parameter estimation, embedding quality, and predictive power, and the learned embeddings generalize across related disability measures. Overall, SCORE offers a scalable, theory-backed approach for decision support and stratification in real-world healthcare, enabling better disease subtyping and outcome prediction with partially labeled, multi-domain data.

Abstract

Electronic Health Records (EHR) offer rich real-world data for personalized medicine, providing insights into disease progression, treatment responses, and patient outcomes. However, their sparsity, heterogeneity, and high dimensionality make them difficult to model, while the lack of standardized ground truth further complicates predictive modeling. To address these challenges, we propose SCORE, a semi-supervised representation learning framework that captures multi-domain disease profiles through patient embeddings. SCORE employs a Poisson-Adapted Latent factor Mixture (PALM) Model with pre-trained code embeddings to characterize codified features and extract meaningful patient phenotypes and embeddings. To handle the computational challenges of large-scale data, it introduces a hybrid Expectation-Maximization (EM) and Gaussian Variational Approximation (GVA) algorithm, leveraging limited labeled data to refine estimates on a vast pool of unlabeled samples. We theoretically establish the convergence of this hybrid approach, quantify GVA errors, and derive SCORE's error rate under diverging embedding dimensions. Our analysis shows that incorporating unlabeled data enhances accuracy and reduces sensitivity to label scarcity. Extensive simulations confirm SCORE's superior finite-sample performance over existing methods. Finally, we apply SCORE to predict disability status for patients with multiple sclerosis (MS) using partially labeled EHR data, demonstrating that it produces more informative and predictive patient embeddings for multiple MS-related conditions compared to existing approaches.

Semi-supervised Clustering Through Representation Learning of Large-scale EHR Data

TL;DR

This work tackles learning informative, transferable patient representations from high-dimensional, sparse EHR data with limited labels. It introduces SCORE, a semi-supervised framework that combines a Poisson-adapted latent factor mixture (PALM) with Gaussian variational approximation (GVA) and leverages pre-trained EHR embeddings to learn compact, transferable patient embeddings. The authors establish convergence and error results for a hybrid EM–GVA algorithm in a setting with diverging latent dimensionality, showing that unlabeled data improve accuracy and robustness to label scarcity. Through simulations and an MS disability-phenotyping application, SCORE outperforms existing methods in parameter estimation, embedding quality, and predictive power, and the learned embeddings generalize across related disability measures. Overall, SCORE offers a scalable, theory-backed approach for decision support and stratification in real-world healthcare, enabling better disease subtyping and outcome prediction with partially labeled, multi-domain data.

Abstract

Electronic Health Records (EHR) offer rich real-world data for personalized medicine, providing insights into disease progression, treatment responses, and patient outcomes. However, their sparsity, heterogeneity, and high dimensionality make them difficult to model, while the lack of standardized ground truth further complicates predictive modeling. To address these challenges, we propose SCORE, a semi-supervised representation learning framework that captures multi-domain disease profiles through patient embeddings. SCORE employs a Poisson-Adapted Latent factor Mixture (PALM) Model with pre-trained code embeddings to characterize codified features and extract meaningful patient phenotypes and embeddings. To handle the computational challenges of large-scale data, it introduces a hybrid Expectation-Maximization (EM) and Gaussian Variational Approximation (GVA) algorithm, leveraging limited labeled data to refine estimates on a vast pool of unlabeled samples. We theoretically establish the convergence of this hybrid approach, quantify GVA errors, and derive SCORE's error rate under diverging embedding dimensions. Our analysis shows that incorporating unlabeled data enhances accuracy and reduces sensitivity to label scarcity. Extensive simulations confirm SCORE's superior finite-sample performance over existing methods. Finally, we apply SCORE to predict disability status for patients with multiple sclerosis (MS) using partially labeled EHR data, demonstrating that it produces more informative and predictive patient embeddings for multiple MS-related conditions compared to existing approaches.

Paper Structure

This paper contains 25 sections, 5 theorems, 29 equations, 5 figures, 1 table.

Key Result

Proposition 1

Under the PALM (eq:model),

Figures (5)

  • Figure 1: Relative F-norm errors in estimating $\bm B$ and $\bm \Lambda$ across different methods. Shown are the relative errors $\|\widehat{\bm B} - \bm B_0\|_F/\|\bm B_0\|_F$ and $\|\widehat{\bm \Lambda} - \bm \Lambda_0\|_F/\|\bm \Lambda_0\|_F$ under varying values of $N$ and $n$, with $p = 400$ and $q = 20$ held fixed.
  • Figure 2: Relative F-norm errors in estimating $\bm B$ and $\bm \Lambda$ across different methods. Shown are the relative errors $\|\widehat{\bm B} - \bm B_0\|_F/\|\bm B_0\|_F$ and $\|\widehat{\bm \Lambda} - \bm \Lambda_0\|_F/\|\bm \Lambda_0\|_F$ under varying values of $p$ and $q$, with $N=5000$ and $n=100$ held fixed.
  • Figure 3: Mean cosine similarity between true and predicted patient embeddings across different methods. Results are shown for varying values of $n$ and $N$, with $p = 400$ and $q = 20$ held fixed.
  • Figure 4: Classification performance of different methods under various generative settings. When one generative parameter is varied, the others are fixed at $n=100$, $N=5000$, $p=400$, and $q=20$ with correct model specifications.
  • Figure 5: Classification performance with varying $n$ for phenotyping of disability for MS patients at MGB. PDDS: trained and validated using disability status based on PDDS, PDDS-EDSS: trained using disability status based on PDDS and validated on EDSS-based labels.

Theorems & Definitions (9)

  • Proposition 1
  • Remark 1
  • Remark 2
  • Remark 3
  • Remark 4
  • Theorem 1: Preliminary supervised estimation
  • Lemma 1: Linear contraction of EM-GVA
  • Theorem 2: Convergence of SCORE
  • Corollary 1: Quality of phenotyping and embedding for some new subject