Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Reverse Engineering Orthognathic Surgery and Orthodontics in Individuals with Cleft Lip and/or Palate: A Case Report
Previous Article in Journal
On the Fused Deposition Modelling of Personalised Bio-Scaffolds: Materials, Design, and Manufacturing Aspects
Previous Article in Special Issue
HGTMDA: A Hypergraph Learning Approach with Improved GCN-Transformer for miRNA–Disease Association Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Unsupervised Error Detection Methodology for Detecting Mislabels in Healthcare Analytics

1
Department of Systems Design Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
2
AI Engineering Team, SpassMed Inc., Toronto, ON M5H 2S6, Canada
*
Author to whom correspondence should be addressed.
Bioengineering 2024, 11(8), 770; https://doi.org/10.3390/bioengineering11080770
Submission received: 30 June 2024 / Revised: 19 July 2024 / Accepted: 25 July 2024 / Published: 31 July 2024

Abstract

:
Medical datasets may be imbalanced and contain errors due to subjective test results and clinical variability. The poor quality of original data affects classification accuracy and reliability. Hence, detecting abnormal samples in the dataset can help clinicians make better decisions. In this study, we propose an unsupervised error detection method using patterns discovered by the Pattern Discovery and Disentanglement (PDD) model, developed in our earlier work. Applied to the large data, the eICU Collaborative Research Database for sepsis risk assessment, the proposed algorithm can effectively discover statistically significant association patterns, generate an interpretable knowledge base for interpretability, cluster samples in an unsupervised learning manner, and detect abnormal samples from the dataset. As shown in the experimental result, our method outperformed K-Means by 38 % on the full dataset and 47 % on the reduced dataset for unsupervised clustering. Multiple supervised classifiers improve accuracy by an average of 4 % after removing abnormal samples by the proposed error detection approach. Therefore, the proposed algorithm provides a robust and practical solution for unsupervised clustering and error detection in healthcare data.

Graphical Abstract

1. Introduction

Machine Learning (ML) models are used to process big healthcare datasets, and various ML models have proven their effectiveness for predictive analytics and diagnostic applications [1]. For instance, Alvin et al. [2] demonstrate how deep learning models improve healthcare quality when applied to electronic health records (EHR) data. Usharani et al. [3] introduced a specialized Fuzzy C-mean segmentation model to improve brain tumor identification. Similarly, another deep learning model is used for rectal cancer detection [4]. In [5], classical and deep learning models were used to estimate the compensatory index for hypovolemia, enhancing emergency medicine assessment. In [6], an LSTM was developed for early diabetes detection. However, there are several challenges that machine learning may face in healthcare data analysis [7].
First, machine learning models always require accurate and complete datasets to ensure effective functionality. Thus, data quality significantly influences reliability and can impact clinical decision-making and patient care [8]. Specifically, the reliability of labels in healthcare data is crucial for the effectiveness of predictive models. When class labels are inaccurate or unreliable, even models with high predictive accuracy may not be reliable due to their dependence on incorrectly labelled data. Moreover, it is difficult for medical experts to manually detect these elusive errors due to lack of contextual information, limiting data privacy regulations, and the sheer scale of data to be reviewed [9].
Next, healthcare providers often need clear explanations for AI-driven decisions, which is challenging with complex models like deep neural networks. Many of these models are black box models, such as genetic algorithms solely focusing on performance [10], lacking transparency and accountability, which can lead to severe consequences [11]. For high-stake applications such as healthcare, the judicial system, and human resource recruitment, machine learning models must be transparent to enable the interpretation of their predictions and explanation of decision boundaries [11,12,13]. As suggested by Rudin [11], creating an interpretable model is a better solution than trying to explain black box models. This allows for synergistic improvement as the interpretable model both integrates and enhances expert knowledge.
Finally, many medical conditions are rare, resulting in imbalanced datasets where the number of cases for one outcome is much smaller than those for others. This class imbalance problem [14] has posed a challenge for years, with machine learning models still struggling to overcome it.
To address these challenges, we have developed a model based on statistical theory, named Pattern Discovery and Disentanglement (PDD) [15,16], to efficiently discover compact sets of significant association patterns linked to their primary sources. It generates a knowledge base interlinking patterns, entities, and sources for visualization, interpretation, and further exploration. In this paper, we propose an unsupervised error detection method using patterns discovered by the PDD to improve the quality of the data.
Some existing algorithms also focus on interpretable patterns. For example, patterns from four gene expressions were explored to classify patients effectively [17], but the number of explored patterns is relatively small, considering the initial input containing 23 features. However PDD’s explainability is more succinct, direct, and extensive, which can use the full dimensions of its input. In addition, the local rule-based explainability method described in Metta et al. [18] uses a genetic algorithm to create a synthetic neighborhood upon which an interpretable model such as a decision tree is applied. This method is limited by the scope of rules that can be extracted from models like decision trees, which are relatively simple and outperformed by PDD.
To identify errors, existing error correction strategies often rely on prediction results from a trained model. These strategies include those implemented through Cleanlab, isolation forest, and AdaBoost. Cleanlab’s system [19] can provide model-specific insights, which are highly dependent on the model’s quality and may inherit biases. Isolation Forest can be used for anomaly detection by determining which points require fewer partitions to separate from other data points [20]. AdaBoost calculates the weighted error of the classifier for each iteration, giving higher weights to misclassified samples to emphasize their importance in the next iteration [21]. Unlike the above strategies, our proposed error detection approach does not depend on a specific model and avoids the issue of model bias. Moreover, it provides an interpretable explanation of dataset and error properties, ensuring a higher-quality dataset before any model training occurs.
To demonstrate the effectiveness of the proposed algorithm, we applied it to the eICU Collaborative Research Database (eICU-CRD) [22] for assessing sepsis risk based on vital signs. Sepsis is a condition that requires prompt diagnosis [23]. Our analysis confirms that the algorithm successfully addresses imbalances in healthcare data, identifies errors and mislabels, and enhances sepsis detection.
Therefore, the contributions of this paper are significant both algorithmically and practically. Algorithmically, (1) we introduce an unsupervised algorithm that effectively discovers patterns or groups of patterns associated with distinct classes, which are functional and statistically independent from the patterns of other classes, as demonstrated by our experimental results; and (2) we propose an error detection method to detect mislabels in the dataset, the effectiveness of which is validated through comparing classification results before and after error correction. Practically, (1) the application of our algorithm to the disease sepsis reveals insightful characteristics of patients and patient groups, enhancing our understanding of the disease’s impact; and (2) the error detection process improves the quality of the dataset, helping to identify outliers or mislabelled patient records.

2. Materials and Methods

2.1. Dataset and Preprocessing

The eICU-CRD [22] includes data from critical care unit patients across the United States for 2014 and 2015. We used the same dataset as the study mentioned in [23], which was collected from key tables such as patient demographics, diagnoses, nursing assessments, and vital periodic records. This reference dataset consisted of 138 features related to patient information, baseline measurements, labels, and temporal vital signs. There are features related to patient data, baseline physical measurements, diagnosis, label, and temporal vital signs. Patient data includes parameters like gender, age, ethnicity, weight and height. Baseline physical measurements include Glasgow coma scale (GCS), systolic blood pressure, diastolic blood pressure, mean blood pressure, pulse pressure, heart rate, respiration, and SpO2 oxygen saturation. Diagnosis includes diagnosisoffset and diagnosispriority. Label includes ICD9 codes and categories. Temporal vital sign data includes feature-lagged physical measurements.
To clean up the dataset, we removed 12 features with excessive missing values, reducing the feature count from 138 to 126. Further preparing the dataset for analysis and preventing data leakage, we excluded an additional six label-related features (e.g., ICD9, categories 1–4, and binary label) and six features used solely for preprocessing (e.g., patientunitstayid, observationoffset, activeupondischarge, diagnosisoffset, diagnosispriority, hospitaladmitoffset). Labels were generated from clinician-evaluated ICD9 codes corresponding to specific diagnoses, and the records were categorized into two groups: sepsis and other. It was notably imbalanced, with 24.4 % of the samples being sepsis-positive [23]. The final dataset used to compare PDD’s error detection comprised of 10,743 records, 113 features, and one label feature. Notably, while the full 114-feature dataset was compared before and after error detection using supervised learning models in Section 3.4, only a 28-feature subset was used by PDD to interpret the dataset and compare unsupervised learning methods in Section 3.1 and Section 3.2. Detailed processing steps are included in the Appendix A with relevant features defined in Table A1.

2.2. Methodology

In this section, we presented the proposed algorithm applied to the dataset descriptive in Section 2.1. The algorithm completes three main tasks: (1) Pattern discovery via PDD for interpretability, as discussed in the [16], (2) Unsupervised learning applied to the dataset; (3) Detecting errors or abnormal samples (i.e., mislabelled records, outliers, and undecided records).

2.2.1. Interpreting Dataset

Our early work with PDD focused on interpretable pattern discovery. In this section, we presented how PDD works for interpretability. As detailed in Appendix A, the processed dataset of 114 features will be used to compare the supervised learning models. This dataset includes features related to patient information, physical assessments, and temporal features capturing two observations (i.e., time points 2 and 14) of six physical measurements. Statistics such as mean, standard deviation, minimum, maximum, kurtosis, and skewness are calculated for each observation. In the pattern discovery process, we aim to discover associations between feature values. Given that many of the statistical features convey the same information, to simplify the analysis and enhance interpretability, we only retained the standard deviation for each observation, as it best describes the variability of data. Therefore, to avoid the complexity introduced by an overwhelming number of features (attributes), we finally retained a subset of 27 features describing patient demographics and their physical assessments. This includes 15 features related to physical measurements and patient information (i.e., GCS, systolic blood pressure, diastolic blood pressure, mean blood pressure, pulse pressure, heart rate, respiration rate, SpO2, age, gender, ethnicity, discharge status, admission weight, discharge weight, and height); another two observation features for each of the six physical measurements, totaling 12 features; and one label feature. Consequently, the dimension of the dataset is 10,743 by 28.
Due to the unlimited degrees of freedom inherent in numerical features, correlating these features with the target variable and interpreting the associations present significant challenges. So, the first step of the PDD process is discretizing the numerical features into event-based or discrete categories according to clinical standards as Table 1 shows. Other numerical features without clear clinical significance such as age and admission weight, were discretized into intervals using the Equal Frequency [24]. This method ensures an equal distribution of data points (i.e., records) within each interval. It is achieved by first sorting the numerical data in ascending order, and then dividing the data into three intervals (bins) each containing an equal number of data points. In information theory, this method indirectly maximizes the data’s informational content. It also handles skewed distributions and simplifies analysis with interpretable categories that are more understandable than precise numerical values. Hence, by applying this discretization approach, all features are converted to categorical ones.
Then, on the discrete dataset, using PDD [16], we construct a statistical residual matrix (SR-Matrix) to account for the statistical strength of the associations among feature values. In pattern discovery, the term “attribute” is used instead of “feature”, so “attribute value” (AV) will be used subsequently. Since the meaning of the attribute values (AVs) and their class labels are implicit, the discovery of a statistically significant association of the AVs is unaffected by prior knowledge or confounding factors. To evaluate the association between each pair of AVs in the SR-Matrix, we calculated the statistical measure of adjusted standardized residual to represent the statistical weights of the association between distinct AV pairs. For instance, if A V 1 represents the value of attribute s y s t o l i c as H, and A V 2 represents the value of attribute d i a s t o l i c as L. Then the adjusted standard residual of the association between A V 1 and A V 2 is calculated in Equation (1), which is denoted as S R ( A V 1 , A V 2 ) .
S R ( A V 1 , A V 2 ) = O c c ( A V 1 , A V 2 ) E x p ( A V 1 , A V 2 ) E x p ( A V 1 , A V 2 ) · ( 1 O c c ( A V 1 ) N ) · ( 1 O c c ( A V 2 ) N )
where O c c ( A V 1 ) and O c c ( A V 2 ) represent the number of occurrences of each attribute value; O c c ( A V 1 , A V 2 ) is the total number of co-occurrences of the two attribute values; and E x p ( A V 1 , A V 2 ) refers to the expected frequency of co-occurrences of the two attribute values; N is the total number of entities. Further details on the calculation are provided in Appendix B.
The key idea of the pattern discovery and disentanglement are briefly described below. We assume that certain attribute values (AVs) or associations co-occur on samples due to certain primary causes (referred to as primary sources). Other significant or insignificant AV associations may also occur in the samples, entangled with those originating from the primary sources. Hence, the goal of PDD is to discover statistically significant association patterns from AV groups on a disentangled principal component (PC), then cluster each AV group into subgroups, denoting them as Disentangled Statistical Units (DSUs). These DSUs are functionally independent and statistically connected within the same group but not associated with other groups.
To achieve such a goal, PDD first applies a linear transformation, Principal Component Analysis (PCA), to decompose the SR-Matrix into Principal Components (PCs). Each PC is functionally independent, capturing associations uncorrelated to those obtained in other PCs. To look for AVs with statistically significant associations with other AVs on the same PC, PDD reprojects each PC back onto an SR-Matrix to generate a Reconstructed SR-Matrix (RSR-Matrix) for each distinct PC. If the maximum residual between a pair of attribute values (AVs) within an RSR-Matrix exceeds a statistical threshold, such as 1.96 for a 95% confidence interval, the association is considered statistically significant. Notably, the associations identified within each RSR-Matrix (or PCs) remain functionally independent from those in other RSR-Matrices (or PCs). Then, on the PCs with statistically significant associations, we obtain AV groups where the AVs within are statistically connected. From entities containing these AVs, we use a hierarchical clustering algorithm to obtain entity clusters containing AV subgroups with statistically connected AVs within, but not outside of that subgroup. They are the DSUs originate from distinct primary sources.
Figure 1 illustrates the AV association disengagement concept. After Principal Component Decomposition of the SR-Matrix, PCs and their corresponding RSR-Matrices are obtained. Only RSR-Matrices containing SR values exceed the statistical significance threshold (i.e., 1.96), and their corresponding PCs are retained. For each of those retaining PCs, three groups of projected AVs can be identified along its axis, each showing different degrees of statistical associations.
  • Those close to the origin are not statistically significant and thus do not associate with distinct groups/classes (marked as (a) in Figure 1), They are close to the origin because all their coordinates (SRs) are low and insignificant.
  • Those at one extremal end of the projections; and those at the opposite extremal end (both marked as (b)).
  • The AV groups or subgroups in (b), if their AVs within are statistically connected but disconnected from other groups, may associate with distinct sources (i.e., classes) (marked as (c)).
As a result, two AV groups at opposite extremes were discovered. Each AV within such a group is statistically linked to at least one other AV within, and none of them is statistically connected to AVs in other groups. Furthermore, it is possible that among the AVs in the AV groups, some of the AV Subgroups may only occur on subsets of the entities uncorrelated to other subgroups. Hence, to achieve a more detailed separation of groups, several subgroups are separated in each AV group based on their appearance in entity groups. This is done using a similarity measure defined by the overlapping of entities each AV can cover. Then AV subgroups found in different entity groups can be identified. We denote such an AV subgroup by a three-digit code [#PC, #Group, #SubGroup] and refer to it as a Disentangled Space Unit (DSU). We hypothesize that these DSUs originate from distinct functional sources.
Therefore, in each subgroup denoted by DSU, a set of AVs is included, which are referred to as pattern candidates. We then developed a pattern discovery algorithm to grow high-order patterns, called comprehensive patterns, from the pattern candidates. In the end, a set of high-order comprehensive patterns is generated within each DSU, and they are all associated with the same distinct source.
The interpretable output of the PDD is organized in a PDD Knowledge Base. This framework is divided into three parts: the Knowledge Space, the Pattern Space, and the Data Space. Firstly, the Knowledge Space lists the disentangled AV subgroups referred to as a Disentangled Space Unit (DSU) (denoted by a three-digit code, [#PC, #Group, #SubGroup] shown in the three columns of the knowledge space to indicate different levels of grouping) linking to the patterns discovered by PDD on the records. Secondly, the Pattern Space displays the discovered patterns, detailing their associations and their targets (the specified class or groups). Thirdly, Data Space shows the record IDs of each patient, linking to the knowledge source (DSU) and the associated patterns. Thus, this Knowledge Base effectively links knowledge, patterns, and data together. If an entity (i.e., a record) is labelled as a class, we can trace the “what” (i.e., the patterns it possesses), the “why” (the specific functional group it belongs to), and the “how” (by linking the patterns to the entity clusters containing the pattern(s)).
The novelty and uniqueness of PDD lie in its ability to discover the most fundamental, explainable, and displayable associations at the AV level from entities (i.e., records) associated with presumed distinct primary sources. This is based on robust statistics, unbiased by class labels, confounding factors, and imbalanced group sizes, yet its results are trackable and verifiable by other scientific methods.

2.2.2. Clustering Patient Records

Without specifying the number of clusters to direct the unsupervised process, PDD can cluster records based on the disentangled pattern groups and subgroups.
As described in Section 2.2.1, the output of PDD is organized into a Knowledge Base, where each pattern subgroup is represented by a D S U [ # P C , # G r o u p , # S u b G r o u p ] . As defined in Section 2.2.1, the set of AVs displayed in each DSU is a summarized pattern, representing the union of all the comprehensive patterns on entities discovered from that subgroup. We denote the number of comprehensive patterns discovered from the summarized pattern in DSU as # C P D S U . For example, in D S U [ 1 , 1 , 1 ] , if 10 comprehensive patterns are found, then # C P [ 1 , 1 , 1 ] = 10 . Each record may possess none, one, or multiple comprehensive patterns for each DSU. We denote the number of comprehensive patterns possessed by a record in a specific DSU as # I D D S U . For example, # 1 [ 1 , 1 , 1 ] = 5 and # 1 [ 2 , 1 , 1 ] = 6 represent the record with I D = 1 possess 5 comprehensive patterns in D S U [ 1 , 1 , 1 ] and 6 comprehensive patterns in D S U [ 2 , 1 , 1 ] .
Each DSU can represent a specific function or characteristic in the data, potentially associated with a particular class. For example, in this study, D S U [ 1 , 1 , 1 ] is associated with s e p s i s , while D S U [ 1 , 2 , 1 ] is associated with o t h e r . The fact that D S U [ 1 , 1 , 1 ] and D S U [ 1 , 2 , 1 ] appear as two opposite groups in P C 1 (Figure 1) indicates that their AV associations have significant differences as captured by P C 1 . Some DSUs might reveal rare patterns not associated with any class while the class label is not in the association.
Based on the definitions described above, we cluster the records by assigning each record to the class that matches the most comprehensive patterns compared to any other class. To provide a more detailed explanation of the clustering process, consider the following example. The DSU outputted by PDD are D S U [ 1 , 1 , 1 ] , D S U [ 1 , 2 , 1 ] , D S U [ 2 , 1 , 1 ] , and D S U [ 2 , 2 , 1 ] , which are associated with s e p s i s , o t h e r , s e p s i s , and o t h e r , respectively. The total number of comprehensive patterns in these DSUs are # C P [ 1 , 1 , 1 ] = 100 , # C P [ 1 , 2 , 1 ] = 100 , # C P [ 2 , 1 , 1 ] = 200 , and # C P [ 2 , 2 , 1 ] = 200 respectively. Consider a record ( I D = n ) with comprehensive patterns possessed by this record are # n [ 1 , 1 , 1 ] = 50 , # n [ 1 , 2 , 1 ] = 60 , # n [ 2 , 1 , 1 ] = 0 , and # n [ 2 , 2 , 1 ] = 150 .
Due to the variation in the number of comprehensive patterns across DSUs, we use a percentage rather than an absolute value to measure the association of the record with pattern groups. Hence, to determine how the record ( I D = n ) is associated with the s e p s i s patterns, we calculate the average percentage of the number of comprehensive patterns associated with a specific class possessed by the record, denoted as # I D C l a s s . Due to # n [ 2 , 1 , 1 ] = 0 indicates that the record is not covered by the DSU[2,1,1], it is excluded from the calculation to avoid the significant impact of a zero value on the final percentage. Hence, the association of the record ( I D = n ) with the s e p s i s patterns is calculated as # n s e p s i s = 50 / 100 = 0.5 . Similarly, # n o t h e r s = m e a n ( 60 / 100 , 150 / 200 ) = 0.675 . Since # n o t h e r s is greater than # n s e p s i s , the record is assigned as o t h e r . To evaluate the accuracy of this assignment for all records, we compare the assigned class label with the original implicit class label.

2.2.3. Detecting Abnormal Records

The evaluation of the classification or prediction involves comparing the predicted labels with the original class labels. However, this comparison is unreliable if mislabels exist in the original data. To address this issue, we proposed an error detection method to identify abnormal records using the patterns discovered by PDD. In our early work on PDD, we integrated both the supervised and unsupervised methods for error detection and class association. In this paper, we simplify the process by using only a novel unsupervised method on a dataset with implicit class labels as the ground truth, making the error detection process more succinct.
To determine whether a record is abnormal, the proposed algorithm compares the class assigned by PDD with its original labels, evaluating the consistency of discovered patterns with their respective explicit or implicit class labels. We define three statuses to an abnormal record: Mislabelled, Outlier, and Undecided, which are detailed below.
  • Mislabelled: If a record is categorized into one class but matches more patterns from a different class according to the PDD output, it suggests the record may be mislabelled. For example, consider the same record with I D = n described in Section 2.2.2 with the same setting of the pattern groups where # n s e p s i s = 0.5 and # n o t h e r s = 0.675 . If the record is originally labelled as s e p s i s in the dataset, but the relative difference ( | # n o t h e r s # n s e p s i s | | # n o t h e r s + # n s e p s i s | ) is greater than 0.1 , this suggests that the record ( I D = n ) is more associated with o t h e r than with s e p s i s . The relative difference is used instead of absolute difference because it provides a scale-independent comparison of the number of patterns associated with one class to another. A value greater than 0.1 indicates that the number of patterns associated with one class is statistically significantly greater than the number associated with another class. Hence, the record ( I D = n ) may be mislabelled.
  • Outlier: If a record possesses no patterns or very few patterns, it may indicate the record is an outlier. For example, a record with I D = m uses the previously described pattern group settings. The comprehensive patterns possessed by this record are: # m [ 1 , 1 , 1 ] = 1 , # m [ 1 , 2 , 1 ] = 0 , # m [ 2 , 1 , 1 ] = 0 , and # m [ 2 , 2 , 1 ] = 1 . Calculating the percentages, # m s e p s i s is 1 / 100 = 0.01 and # m o t h e r s is 1 / 200 = 0.005 . Both percentages are less than or equal to 0.01 , suggesting that record m possesses fewer than 1% of the patterns associated with either class, which may indicate it is an outlier.
  • Undecided: If the number of possessed patterns for a record is similar across different classes, the record should be classified as undecided. For example, a record with I D = k uses the previously described pattern group settings. The comprehensive patterns possessed by this record are: # k [ 1 , 1 , 1 ] = 50 , # k [ 1 , 2 , 1 ] = 60 , # k [ 2 , 1 , 1 ] = 110 , and # k [ 2 , 2 , 1 ] = 100 . Calculating the percentages, # k s e p s i s is the mean of 50 / 100 and 110 / 200 , which is 0.55 ; and # k o t h e r s is the mean of 60 / 100 and 100 / 200 , which is 0.55 . Since the difference between the two percentages is zero or less than 0.001 , record k may be associated with both classes, suggesting it is undecided
To avoid adding new incorrect information, mislabelled, undecided, and outliers are removed from the dataset. Hence, to validate the effectiveness of abnormal records detected by PDD, we compared classification results from the original dataset to those from a dataset without abnormal records when various classifiers were applied.

3. Results

The experimental results are presented in this section. First, Section 3.1 demonstrates the interpretable output derived from PDD. Next, Section 3.2 compares PDD’s clustering accuracy with the K-means algorithm which serves as our baseline. Then, Section 3.3 presents our error detection results, showing interpretable outcomes with discovered patterns after mislabeling detection. Finally, Section 3.4 demonstrates the efficacy of error detection by comparing the performance of supervised learning algorithms both before and after the application of error detection.

3.1. Interpretable Results

In this section we present interpretable results, via the PDD knowledge base, containing all summarized patterns from different DSUs in Figure 2. The knowledge base contains three parts: the knowledge space, the pattern space, and the data space. First, the Knowledge Space uses the triple code of the DSU to indicate the primary source that gives rise to the disentangled patterns from the AV Subgroup. In this specific study, as shown in Figure 2, the first row displays the group discovered as the first disentangled space denoted by the code [ 1 , 1 , 1 ] , indicating in reverse order the first Sub-Pattern group (SubGroup) of the first Pattern Group (Group) discovered in the first Principal Component (PC). Here, we identify two groups corresponding to two types of diseases: s e p s i s and o t h e r . The results show all DSUs discovered by PDD representing all the disentangled spaces with statistically significant patterns obtained from the first three PCs.
Second, the Pattern Space reveals the specific AV association patterns and their targets (the specified class or groups, if given) obtained from the DSU. We find that the disease sepsis exhibits different patterns than other. For example, the sepsis class was associated with lower blood pressure, higher heart rate, higher respiration, and moderate or severe GCS. Inversely, the other class was associated with higher blood pressure, normal heart rate, lower respiration, and mild GCS. Since the patterns outputted by PDD are disentangled, they clearly distinguish between the different classes of sepsis and other.
The patterns related to sepsis demonstrate clinical relevance as hypotension, tachycardia, elevated respiratory rate, and altered mental status are all considered key criteria defined in the Third International Consensus for Sepsis and Septic Shock (Sepsis-3) [25]. Singer et al. [25] discussed how septic shock is associated with underlying circulatory, cellular, and metabolic abnormalities that substantially increase mortality. These characteristic abnormalities helped to improve understanding of sepsis pathobiology and informed treatment such as a vasopressor requirement to constrict blood vessels and maintain a minimum mean arterial pressure in the absence of hypovolemia or adequate fluid resuscitation. This shows evidence of a positive correlation between low blood pressure and the severity of sepsis among patients with septic shock, which is consistent with patterns discovered by PDD. In addition to this, sepsis is impacted by various biomarkers and physiological parameters, such as C-reactive protein, white blood cell count, and central venous pressure [26,27]. Thus, with a more detailed dataset, PDD can further examine associations between sepsis and these factors.
Thirdly, the Data Space displays the records with their IDs from the original dataset and the discovered patterns each record possesses. As shown in Figure 2, the number in the column of each ID on the row linking to the pattern(s) and the DSU(s) represents the number of the comprehensive patterns in a DSU that the record possesses. For example, 180 indicates there are 180 comprehensive patterns in the D S U [ 1 , 1 , 1 ] covered by record 1. Due to space limitations, only the first nine records ( I D = 1 to 9 ) are displayed.
In the original data, records labelled s e p s i s are listed before those labelled o t h e r , so the first nine records belong to the s e p s i s class. It is also evident in the Data Space that the records possess more patterns associated with s e p s i s than with o t h e r .

3.2. Comparison of Unsupervised Learning

In addition, we compared the PDD’s clustering results with the classical clustering algorithm K-mean as the baseline. For a clear comparison, both the same subset ( 10 , 743 × 28 ) and the full dataset ( 10 , 743 × 114 ) were used to evaluate K-means performance. To run K-means, we used the s k l e a r n . c l u s t e r s package in Python 3.0 [28] with all default parameter settings and assigned the number of clusters as two, as the original data only contains the two classes of s e p s i s and o t h e r .
To evaluate an imbalanced dataset, balanced accuracy is more reliable than regular accuracy. Using regular accuracy, a model might frequently predict the majority class and still achieve a high accuracy rate, despite failing to accurately predict the minority class [29]. Balanced accuracy, however, considers the true positive rate of each class, offering a more honest evaluation of model performance, which is essential when one class significantly outnumbers another. Therefore, we utilize both balanced accuracy and accuracy for evaluating classification or clustering performance, in addition to precision, recall, and a weighted F1-score for comparison.
Based on the unsupervised clustering results in Table 2, PDD significantly outperformed K-Means on both the 28-feature subset and 114-feature dataset. On the 28-feature subset, the balanced accuracy of 0.89 achieved with PDD was significantly higher than the 0.42 achieved with K-means. On the 114-feature dataset, K-means exhibited improved performance, demonstrating the benefit of including more features. However, even when comparing PDD’s performance on the reduced feature subset to K-means performance on the full feature dataset, PDD still outperformed K-means by 38 % .
The two clustering methods were also compared in terms of their effectiveness in handling dataset imbalance, given that only 24.4% of samples were sepsis-positive. As shown in Table 2, PDD’s balanced accuracy on the 28-feature subset exceeds regular accuracy indicating a more effective prediction of the minority class s e p s i s . In contrast, K-means’ regular accuracy exceeds balanced accuracy on both the 28 and 114-feature datasets indicating a more effective prediction of the majority class o t h e r . Moroever, the difference between the accuracy and balanced accuracy of PDD results is a minimal 2% in comparison to K-means’ results with differences of 10% and 24%. This indicates that PDD is better suited for handling imbalanced data and is a valuable tool for datasets where accurate prediction of the minority class is crucial. Also, that PDD has learned meaningful patterns in the minority class despite its scarcity in the dataset while K-means struggles to accurately capture these characteristics and is therefore challenged by the imbalanced distribution. In this case, PDD not only predicts more accurate clustering results for the minority class but also provides patterns that reveal the characteristics of the cluster. PDD can display all the statistically significant patterns with their source (DSU) possessed by each cluster. It also corrects mislabels based on its error detection capability and the records’ implicit class labels.

3.3. Error Detection Results

This section presents the error detection results using the PDD knowledge base. As mentioned in Section 2.2.3 and Section 3.1, we used the comprehensive pattern count from the PDD data space to detect abnormal records in the original dataset. For example, I D = 5510 has # 5510 [ 1 , 1 , 1 ] = 30 , # 5510 [ 2 , 1 , 1 ] = 13 , # 5510 [ 3 , 1 , 1 ] = 2 , # 5510 [ 5 , 1 , 1 ] = 6 , and all DSU is associated with s e p s i s . Thus, # 5510 s e p s i s = m e a n ( 30 / 1357 , 13 / 2617 , 2 / 453 , 6 / 940 ) 0.0095 . Similarly, # 5510 [ 1 , 2 , 1 ] = 2 , # 5510 [ 2 , 2 , 1 ] = 51 , # 5510 [ 3 , 2 , 1 ] = 3 , # 5510 [ 5 , 2 , 1 ] = 8 , thus # 5510 o t h e r = m e a n ( 2 / 696 , 51 / 2514 , 3 / 561 , 8 / 1712 ) 0.0083 . Both # 5510 s e p s i s and # 5510 o t h e r are less than the threshold 0.01 specified in Section 2.2.3, I D = 5510 is categorized as an outlier.
In total, PDD identified 996 mislabelled records, 2 outliers, and 159 undecided records. To avoid adding new incorrect information, we finally removed the 1157 records, and 9586 records were retained. For simplicity, only a partial result is shown in Figure 3.
For better visualization, in Figure 3a, all the DSUs with l a b e l = s e p s i s are highlighted in red, and all the DSUs with l a b e l = o t h e r are highlighted in green. An important observation after such highlighting is that attributes are consistent within its class. For example, for all DSUs with l a b e l = s e p s i s , their blood pressure features (i.e., s y s t o l i c b p , d i a s t o l i c b p , m e a n b p ) are mostly low. And for all DSUs with l a b e l = o t h e r , their blood pressure attributes are mostly high. Due to this consistency, we assigned the same highlighting scheme to specific AV in Figure 3b. For example, when s y s t o l i c b p = l o w , we highlighted the AV in red, implying its positive correlation with sepsis through PDD’s Pattern Space. With such a highlighting scheme, one can easily observe the pattern entanglement among the original dataset and how well PDD can disentangle the source attributes with a statistically supported, interpretable knowledge base.

3.4. Error Detection Applied to Supervised Learning Models

To demonstrate the effectiveness of the error detection methods, we compared the classification results of various classifiers on the dataset before and after removing errors. The classifiers tested include Random Forest, SVM, Neural Networks, Logistic Regression, LightGBM, and XGBoost. 80 % of the data are used for training and 20 % for testing. To evaluate the results, the same metrics (i.e., Precision, Recall, Accuracy, Balanced Accuracy, and Weighted F1-Score) are used for comparison.
Table 3 shows the comparison results of existing ML models before and after error detection. This error detection process by PDD resulted in an accuracy improvement of approximately 1–9% and 4 % on average. The error detection results show how PDD can remove those abnormal records (outliers/undecided/mislabelled) to improve the accuracy of class association for building a classifier based on a cleaner training set since the balanced accuracy is increased for all classifiers.
In Table 3, it is also interesting to note how different classifiers handle overfitting and bias in imbalanced data with varying capabilities. First, for all classifiers, the regular accuracies are greater than the balanced accuracies, indicating that the predictions are more effective for the majority class, particularly for the first four traditional classifiers (RF, SVM, NN, and LR). However, for models with class balancing capability (LightGBM and XGBoost), the differences between regular accuracy and balanced accuracy are small. Second, the higher accuracy before error detection indicates that models are performing better on the majority class, which can be a sign of overfitting to the majority class. So the accuracy slightly decreases after removing abnormal cases due to reducing the overfitting for the majority class. Lastly, after error detection, the balanced accuracies increase. This indicates that outliers and undecided cases were biasedly absorbed into the majority classes before error detection, but after error detection, the classifiers show better prediction results for the minority class, leading to better results for the balanced prediction.
Additionally, to effectively highlight comparison results, Figure 4 displays radar charts that visualize each classifier’s overall performance profile. These charts distinctly show that the blue areas, representing improved metrics, are larger than the red areas for all models, with particularly obvious differences in LightGBM and XGBoost.

4. Discussion

In this section, we discuss the results of interpretability and error detection using PDD from the perspective of previous studies.

4.1. Interpretability

Decision Trees are widely regarded as interpretable models, and similarly, Random Forests (RF) can also provide interpretable insights despite being less straightforward due to their ensemble nature. Various methods have been developed to extract feature importance or feature contributions to enhance its interpretability [30]. In this section, we compare the interpretability between RF and PDD.
First, according to the RF model illustrated in Figure 5, RF assesses the importance of features associated with class labels. Similarly, as shown in Figure 2, PDD discovers patterns that reveal the associations between feature values rather than just feature associations.
Second, the importance analysis provided by RF shows the association between one feature and the target, while PDD goes further by revealing the significance of high-order comprehensive patterns linked to the target, which provides a more direct, succinct, precise, and interpretable understanding of the problem. This capability allows PDD to provide a more in-depth and precise interpretation, which is crucial for high-stakes applications requiring transparency and explainability.
Lastly, RF can also output prediction rules similar to the high-order patterns in PDD. However, according to the supervised learning comparison results provided in Section 3.4, before error detection, the accuracy of Random Forest is 86%, while the balanced accuracy is only around 73%. This indicates that, for imbalanced classes, the rules outputted by RF are easily biased towards the majority class. In contrast, as demonstrated in the study [16], PDD can discover disentangled patterns even for imbalanced data. Furthermore, PDD uses hierarchical clustering to highlight the independent functions of each pattern group, enhancing interpretability. It generates a knowledge base linking patterns, entities, and classes, which can trace the origins of the discovered patterns. This capability ensures the accuracy of error detection. This level of detail and traceability is a limitation of the RF model.

4.2. Error Detection

Real-world data may contain inconsistencies due to incorrect data collection, human errors, equipment inaccuracies, and missing values [31]. So removing errors and cleaning data is essential in developing a reliable ground truth for training and testing. In this study, the conditions of patients with sepsis may be various, leading to fluctuations that complicate establishing baseline values. The eICU-CRD dataset, collected in 2014 and 2015, predates the sepsis definition updated in 2016 [25], with the initial definition in 2001. Thus, the labels in the dataset may not reflect the latest medical insights. Furthermore, integrating data from various hospitals with different environments, procedures, equipment, and staff may introduce non-standardized variables.
PDD’s error detection results in Section 3.4 demonstrate removing abnormal data and creating a more reliable ground truth can improve model performance. The interpretability of PDD enables it to discover and disentangle patterns, which can be used to assign status—mislabelled, outlier, and undecided—to records. The traceability of PDD allows humans to trace patterns back to the data and clarify why a particular sample should be corrected. Furthermore, unlike other existing error detection methods that rely on model predictions, PDD is model-independent.
To better understand PDD’s error detection, sample scatter pair plots between the features of the eICU-CRD dataset are shown in Figure 6. The outlier, undecided, and mislabelled records are the 10 abnormal records presented in Section 3.3. In Figure 6a,b, regardless of whether the x-axis values are continuous or discrete, the decision boundaries for classification are not linear. Two classes are mixed, making classification and error detection based on the original features challenging due to entanglement.
However, the two classes can be visually disentangled in PDD’s Data Space as shown in Figure 7 and Figure 8 by treating each DSU as a distinct feature and the comprehensive pattern count as the feature value. The only difference between Figure 7 and Figure 8 is the proportion of the abnormal records that are highlighted. Figure 7 highlights the 10 abnormal records from Section 3.3. Figure 8 highlights all the 1157 abnormal records, which comprise 996 mislabelled, 2 outliers, and 159 undecided records. The DSU plot shows a more prominent classification decision boundary.
Specifically, in Figure 7a,b, the DSUs of the x and y axes are found at the two opposite ends of the same PC, as demonstrated in Figure 1. So, if one DSU captures the AV associations solely from class s e p s i s , the other DSU should capture the AV associations solely from class o t h e r . By observation, Figure 7b has a more prominent classification decision boundary than Figure 7a, suggesting that PC2’s clustering of the opposite groups of AV Associations is better aligned with the expected classification than PC1. The classification decision boundary can be curved as in Figure 7b, or approximated by a straight line as in Figure 7d. Both demonstrate the disentanglement capability of PDD.
The figures also revealed some insights about the abnormal records. First, some undecided records are close to the decision boundary, as shown in Figure 7a,b,d, since these records possess almost an equal number of patterns from both classes. This is more obvious in Figure 8b,d since the two classes are more disentangled and all the undecided records are highlighted. Second, most of the mislabelled records are consistently found in the s e p s i s cluster in Figure 7 and Figure 8. Specifically, 894 of the 996 ( 90 % ) mislabelled records were originally labelled as o t h e r but suggested as s e p s i s . This is more evident in Figure 8b,d since the two classes are more disentangled and all the mislabelled records are highlighted. Furthermore, all outlier records are situated close to the origin because they possess few patterns. Due to the unavoidable overlapping issue, this is more observable in Figure 7. Lastly, when the DSUs of both axes are associated with the same class but from different PCs, which is s e p s i s in Figure 7c and Figure 8c, the classification decision boundary is not as clear as in o t h e r figures when the DSU pairs are associated with different classes as in Figure 7a,b,d and Figure 8a,b,d.

5. Conclusions

In this study, we applied our Pattern Discovery and Disentanglement (PDD) methodology to healthcare data analysis, addressing challenges like mislabeling, class imbalance, and the need for explainable AI. The uniqueness and novelty of PDD lie in its ability to discover fundamental, explainable associations at the feature value level, unbiased by class labels, confounding factors, or imbalanced group sizes. PDD generates a comprehensive knowledge base that links patterns, entities, and sources/classes, providing deep insights and enhancing the interpretability of the results. After detecting the errors (i.e., outliers, undecided, and mislabelled), the original dataset can be cleaned. In this study, using the eICU-CRD data for sepsis risk assessment, PDD demonstrated significant improvements over traditional K-means clustering, achieving a 38 % improvement on the full dataset and a 47 % improvement on the reduced dataset. PDD’s error detection identified 996 mislabelled records, 2 outliers, and 159 undecided cases, resulting in the removal of 1157 abnormal records. This process improved the accuracy of multiple supervised ML models by an average of 4 % , with improvements ranging from 1 % to 9 % . The analysis highlighted PDD’s ability to handle patterns challenging for traditional models, particularly lower-order patterns with limited feature correlation.
In conclusion, PDD offers a robust and practical solution for unsupervised clustering and error detection in large-size and low-quality healthcare data, enhancing both the stand-alone performance and the results of other ML models. Its application to this large, complex dataset demonstrates its practicality and effectiveness. Future work will extend PDD to other medical conditions and datasets, further validating its utility and impact in diverse healthcare applications.

Author Contributions

Conceptualization, P.-Y.Z., A.K.C.W., A.B. and S.P.; methodology, P.-Y.Z., A.K.C.W. and T.J.W.; software, P.-Y.Z., F.L. and T.J.W.; validation, P.-Y.Z. and A.K.C.W.; formal analysis, P.-Y.Z.; investigation, P.-Y.Z., T.J.W. and F.L.; resources, A.B., S.P. and C.D.; data curation, A.B., S.P. and C.D.; writing—original draft preparation, P.-Y.Z., F.L., T.J.W. and A.K.C.W.; writing—editing and review, P.-Y.Z., F.L., T.J.W., A.B., S.P., C.D. and A.K.C.W.; visualization, P.-Y.Z. and T.J.W.; supervision, P.-Y.Z. and A.K.C.W.; project administration, P.-Y.Z.; funding acquisition, A.K.C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NSERC I2I Phase I Grant.

Institutional Review Board Statement

This study did not require ethical approval as it utilized The eICU Collaborative Research Database, a freely available multi-center database designed for critical care research.

Informed Consent Statement

Not applicable.

Data Availability Statement

Access to the eICU Collaborative Research Database can be requested at https://eicu-crd.mit.edu/ (accessed on 15 May 2023).

Conflicts of Interest

Author Anubhav Bhatti, Surajsinh Parmar, and Chen Dan were employed by the company SpassMed Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MLMachine Learning
EHRelectronic health records
PDDPattern Discovery and Disentanglement
eICU-CRDeICU Collaborative Research Database
GCSGlasgow coma scale
SPO2oxygen saturation
SRadjusted standard residual
SR-Matrixstatistical residual matrix
RSR-MatrixReprojected SR-Matrix
AVattribute value
PCAPrincipal Component Analysis
DSUDisentangled Space Unit
SubGroupSub-Pattern group
GroupPattern Group
PCPrincipal Component
RFRandom Forests

Appendix A. Dataset Description and Preprocessing Details

These steps outline the process by which the reference dataset, with 138 features, from Parmar et al. [23] was obtained, and how this dataset was further processed to obtain the final 114 feature dataset for this paper’s analysis.

Appendix A.1. Step 1: Generate Temporal Features Based on Vital Signs

First, the vital periodic table contains vital sign data generated from one minute averages and archived as five minute median values at timepoints that can be identified by their observationOffset. Notably, although the dataset claims to be standardized, many of these values are missing data and not actually aligned to five-minute intervals. Every “group” contains temporal vital sign observations over a selected six-hour period. The period begins at a specific observationOffset and contains values recorded at subsequent time points. It must also contain a diagnosisOffset. The diagnosisOffset is the time since admission before a diagnosis was recorded. Patients can have multiple diagnoses during their ICU stay. The objective was to select suitable six-hour periods and align time points to five-minute intervals. To ensure diagnosisOffset was always present within the period, the search began with the first diagnosisOffset after skipping the first hour after ICU admission, which is considered an unstable period. A six-hour total period was determined by retrieving a three-hour window before the diagnosisOffset followed by a three-hour backend length. The three-hour window would be used for prediction and the three-hour backend for labeling. This process was repeated for the next available diagnosisOffset after a minimum three-hour buffer time to avoid overlap. Time points were manually added using fill-forward and fill-backward to fill in missing values and standardize them to five-minute intervals. To prioritize reliance on observed rather than manually added data, groups were filtered to select those where the actual observation offsets constitute at least half of the duration covered by manually added time points. Moreover, groups with more than 20 missing time points, out of the total 72 five-minute intervals, were filtered. For each group, the mean blood pressure must be less than systolic and diastolic blood pressure. If not, these values were replaced with NaN and then filled with the same techniques for missing values. This data was then used to generate lagged features and lagged Fourier-transformed features to identify trends and cyclical patterns. Mean, standard deviation, minimum, maximum, kurt, and skew were also calculated and added as features without look-ahead bias. The final dataset used timepoints 2 and 14 which were one hour apart.

Appendix A.2. Step 2: Process Labels

Next, labels were first created using diagnosisString based on clinician evaluation at the time. They were relabelled afterward using the calculated qSOFA score for sepsis. The diagnosis string column contains a string made up of four different categorical diagnoses split by pipe. These were separated into four new columns named Categories 1–4. Only patients with diagnosisPriority as primary was considered. Groups labelled as sepsis included entries where the string ‘sepsis’ or ‘septic shock’ appeared under Category 3. Groups labelled as ‘other’ included entries with conditions within the top 20 diseases in the dataset under Category 3. The calculated qSOFA score was used to relabel all groups. The Third International Consensus Definition for Sepsis (Sepsis-3) was updated in 2016 with an acute increase in SOFA score of 2 or more, from baseline, in a patient with suspected infection associated with a minimum twofold increased risk of death. At least two of the following clinical criteria together would constitute a new bedside clinical score termed quick sofa (qSOFA): systolic blood pressure of 100 mmHg or less, respiratory rate of 22 breaths/min or higher, and Glasgow Coma scale (GCS) less than 15 [25]. The last three hours of each group were used to calculate qSOFA. Groups were removed if the qSOFA score did not match the doctor’s diagnosis for more than time.

Appendix A.3. Step 2: Filter Features

Further filtering was used to obtain the same dataset with 138 features used in Parmar et al.: dropping columns with all null values, caused by lagged feature generation and inputting ethnicity with a simple input using the ’most frequent’ strategy. Lastly, the final processed dataset with 10,743 records and 114 features, for this paper’s analysis, was obtained by removing encoded binary values, removing label-related features to avoid leaking information (ICD9, categories 1–4, and label), and removing unnecessary features only used for preprocessing (patientunitstayid, observationoffset, activeupondischarge, diagnosisoffset, diagnosispriority, hospitaladmitoffset).
Table A1. Preprocessing Features.
Table A1. Preprocessing Features.
FeatureTableDescription
patientUnitStayIDpatientA distinct ICU admission. The same patient can have multiple.
observationOffsetvital periodicThe number of minutes, since admission, before the periodic value was recorded.
diagnosisoffsetdiagnosisNumber of minutes, since admission, before diagnosis was recorded.
diagnosisPrioritydiagnosisDiagnosis was marked as Primary, Major, or Other
diagnosisStringdiagnosisThe full pathstring of diagnosis selected in Philips eCareManager software with sections separated with pipes (e.g., cardiovascular|shock/hypotension|sepsis )
ICD9CodediagnosisInternational Classification of Diseases code used to classify diagnoses (e.g., 518.81, 491.20, etc.)

Appendix B. Calculation of Adjusted Residual

In statistics, two events are dependent if their observed frequency of occurrences deviates significantly from the independence hypothesis. Such deviation is known as residual. The adjusted residual is the normalized residual and is defined as the standardized residual divided by an estimate of the standard error [32]. The standardized residual is the ratio of the difference between the observed count and the expected count to the standard deviation of the expected count [33].
Specifically, the adjusted standard residual of the association between A V 1 and A V 2 is calculated in Equation (A1), which is denoted as S R ( A V 1 , A V 2 ) . If the true value of the S R ( A V 1 , A V 2 ) exceeds 1.96, the A V 1 and A V 2 are associated and likely to occur together, with a confidence level of 95%. If S R ( A V 1 , A V 2 ) is less than 1.96, then A V 1 and A V 2 are unlikely to co-occur [34].
S R ( A V 1 , A V 2 ) = O c c ( A V 1 , A V 2 ) E x p ( A V 1 , A V 2 ) E x p ( A V 1 , A V 2 ) · ( 1 O c c ( A V 1 ) N ) · ( 1 O c c ( A V 2 ) N )
Consider the following example to illustrate the calculation of adjusted residuals. Suppose there are two features, systolic and diastolic blood pressure. Systolic blood pressure can be H i g h or L o w , while diastolic blood pressure can be H i g h , M e d i u m , or L o w . The Table A2 summarizes the observed frequencies.
Table A2. Example of an Observed Frequency Table for Two Attributes: Systolic and Diastolic.
Table A2. Example of an Observed Frequency Table for Two Attributes: Systolic and Diastolic.
Association TableSystolic—HighSystolic—LowTotal
Diastolic—High50 (39)40 (51)90
Diastolic—Medium30 (43.33)70 (56.67)100
Diastolic—Low50 (47.67)60 (62.33)110
Total130170300
Assume A V 1 represents Diastolic is High and A V 2 represents Systolic is High. To determine the association between these two attribute values, we can calculate the adjusted residual. According to Equation (1), O c c ( D i a s t o l i c = H i g h , S y s t o l i c = H i g h ) represents the co-occurrence of these two attribute values, which is given as 50 in the table. The Expected value can be calculated using Equation (A2). In this example E x p ( D i a s t o l i c = H i g h , S y s t o l i c = H i g h ) is 90 130 300 which results in 39.
E x p ( A V 1 , A V 2 ) = O c c ( A V 1 ) O c c ( A V 2 ) N
So,
S R ( D i a s t o l i c = H i g h , S y s t o l i c = H i g h ) = 50 39 39 · ( 1 90 300 ) · ( 1 130 300 ) = 11 39 · ( 0.7 ) · ( 0.57 ) = 2.7885
The adjusted residual result is 2.7885, which exceeds the threshold of 1.96, indicating that the association between High Systolic and High Diastolic is statistically significant.

References

  1. Weerasinghe, K.; Scahill, S.L.; Pauleen, D.J.; Taskin, N. Big data analytics for clinical decision-making: Understanding health sector perceptions of policy and practice. Technol. Forecast. Soc. Chang. 2022, 174, 121222. [Google Scholar] [CrossRef]
  2. Rajkomar, A.; Oren, E.; Chen, K.; Dai, A.M.; Hajaj, N.; Hardt, M.; Liu, P.J.; Liu, X.; Marcus, J.; Sun, M.; et al. Scalable and accurate deep learning with electronic health records. NPJ Digit. Med. 2018, 1, 1–10. [Google Scholar] [CrossRef] [PubMed]
  3. Bhimavarapu, U.; Chintalapudi, N.; Battineni, G. Brain tumor detection and categorization with segmentation of improved unsupervised clustering approach and machine learning classifier. Bioengineering 2024, 11, 266. [Google Scholar] [CrossRef]
  4. Lin, C.Y.; Wu, J.C.H.; Kuan, Y.M.; Liu, Y.C.; Chang, P.Y.; Chen, J.P.; Lu, H.H.S.; Lee, O.K.S. Precision Identification of Locally Advanced Rectal Cancer in Denoised CT Scans Using EfficientNet and Voting System Algorithms. Bioengineering 2024, 11, 399. [Google Scholar] [CrossRef] [PubMed]
  5. Bedolla, C.N.; Gonzalez, J.M.; Vega, S.J.; Convertino, V.A.; Snider, E.J. An explainable machine-learning model for compensatory reserve measurement: Methods for feature selection and the effects of subject variability. Bioengineering 2023, 10, 612. [Google Scholar] [CrossRef] [PubMed]
  6. Al Sadi, K.; Balachandran, W. Leveraging a 7-Layer Long Short-Term Memory Model for Early Detection and Prevention of Diabetes in Oman: An Innovative Approach. Bioengineering 2024, 11, 379. [Google Scholar] [CrossRef] [PubMed]
  7. Rudin, C.; Chen, C.; Chen, Z.; Huang, H.; Semenova, L.; Zhong, C. Interpretable machine learning: Fundamental principles and 10 grand challenges. Stat. Surv. 2022, 16, 1–85. [Google Scholar] [CrossRef]
  8. Raghupathi, W.; Raghupathi, V. Big data analytics in healthcare: Promise and potential. Health Inf. Sci. Syst. 2014, 2, 3. [Google Scholar] [CrossRef] [PubMed]
  9. Nguyen, T.; Diakiw, S.; VerMilyea, M.; Dinsmore, A.; Perugini, M.; Perugini, D.; Hall, J. Efficient automated error detection in medical data using deep-learning and label-clustering. Sci. Rep. 2023, 13, 19587. [Google Scholar] [CrossRef]
  10. Caponetto, R.; Fortuna, L.; Graziani, S.; Xibilia, M. Genetic algorithms and applications in system engineering: A survey. Trans. Inst. Meas. Control 1993, 15, 143–156. [Google Scholar] [CrossRef]
  11. Rudin, C. Stop explaining black box machine learning models for high stakes decisions and use interpretable models instead. Nat. Mach. Intell. 2019, 1, 206–215. [Google Scholar] [CrossRef]
  12. Mehrabi, N.; Morstatter, F.; Saxena, N.; Lerman, K.; Galstyan, A. A survey on bias and fairness in machine learning. ACM Comput. Surv. (CSUR) 2021, 54, 1–35. [Google Scholar] [CrossRef]
  13. Amann, J.; Blasimme, A.; Vayena, E.; Frey, D.; Madai, V.I.; Consortium, P. Explainability for artificial intelligence in healthcare: A multidisciplinary perspective. BMC Med. Inform. Decis. Mak. 2020, 20, 1–9. [Google Scholar] [CrossRef]
  14. Sun, Y.; Wong, A.K.; Kamel, M.S. Classification of imbalanced data: A review. Int. J. Pattern Recognit. Artif. Intell. 2009, 23, 687–719. [Google Scholar] [CrossRef]
  15. Wong, A.K.; Zhou, P.Y.; Butt, Z.A. Pattern discovery and disentanglement on relational datasets. Sci. Rep. 2021, 11, 5688. [Google Scholar] [CrossRef] [PubMed]
  16. Wong, A.K.; Zhou, P.Y.; Lee, A.E.S. Theory and rationale of interpretable all-in-one pattern discovery and disentanglement system. Npj Digit. Med. 2023, 6, 92. [Google Scholar] [CrossRef]
  17. Ghalwash, M.F.; Radosavljevic, V.; Obradovic, Z. Extraction of interpretable multivariate patterns for early diagnostics. In Proceedings of the 2013 IEEE 13th International Conference on Data Mining, Dallas, TX, USA, 7–10 December 2013; IEEE: Piscataway, NJ, USA, 2013; pp. 201–210. [Google Scholar]
  18. Metta, C.; Beretta, A.; Pellungrini, R.; Rinzivillo, S.; Giannotti, F. Towards Transparent Healthcare: Advancing Local Explanation Methods in Explainable Artificial Intelligence. Bioengineering 2024, 11, 369. [Google Scholar] [CrossRef]
  19. Northcutt, C.G.; Athalye, A.; Mueller, J. Pervasive label errors in test sets destabilize machine learning benchmarks. arXiv 2021, arXiv:2103.14749. [Google Scholar]
  20. Zhou, Z.H. Ensemble Methods: Foundations and Algorithms; CRC Press: Boca Raton, FL, USA, 2012. [Google Scholar]
  21. Khan, A.A.; Chaudhari, O.; Chandra, R. A review of ensemble learning and data augmentation models for class imbalanced problems: Combination, implementation and evaluation. Expert Syst. Appl. 2023, 244, 122778. [Google Scholar] [CrossRef]
  22. Pollard, T.J.; Johnson, A.E.; Raffa, J.D.; Celi, L.A.; Mark, R.G.; Badawi, O. The eICU Collaborative Research Database, a freely available multi-center database for critical care research. Sci. Data 2018, 5, 180178. [Google Scholar] [CrossRef]
  23. Parmar, S.; Shan, T.; Lee, S.; Kim, Y.; Kim, J.Y. Extending Machine Learning-Based Early Sepsis Detection to Different Demographics. In Proceedings of the 2024 IEEE First International Conference on Artificial Intelligence for Medicine, Health and Care (AIMHC), Laguna Hills, CA, USA, 5–7 February 2024; IEEE: Piscataway, NJ, USA, 2024; pp. 70–71. [Google Scholar]
  24. Zhou, P.Y.; Sze-To, A.; Wong, A.K. Discovery and disentanglement of aligned residue associations from aligned pattern clusters to reveal subgroup characteristics. BMC Med. Genom. 2018, 11, 35–49. [Google Scholar] [CrossRef]
  25. Singer, M.; Deutschman, C.S.; Seymour, C.W.; Shankar-Hari, M.; Annane, D.; Bauer, M.; Bellomo, R.; Bernard, G.R.; Chiche, J.D.; Coopersmith, C.M.; et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA 2016, 315, 801–810. [Google Scholar] [CrossRef] [PubMed]
  26. Hilarius, K.W.; Skippen, P.W.; Kissoon, N. Early recognition and emergency treatment of sepsis and septic shock in children. Pediatr. Emerg. Care 2020, 36, 101–106. [Google Scholar] [CrossRef] [PubMed]
  27. Dellinger, R.P.; Levy, M.M.; Rhodes, A.; Annane, D.; Gerlach, H.; Opal, S.M.; Sevransky, J.E.; Sprung, C.L.; Douglas, I.S.; Jaeschke, R.; et al. Surviving sepsis campaign: International guidelines for management of severe sepsis and septic shock: 2012. Crit. Care Med. 2013, 41, 580–637. [Google Scholar] [CrossRef] [PubMed]
  28. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  29. Haixiang, G.; Yijing, L.; Shang, J.; Mingyun, G.; Yuanyue, H.; Bing, G. Learning from class-imbalanced data: Review of methods and applications. Expert Syst. Appl. 2017, 73, 220–239. [Google Scholar] [CrossRef]
  30. Palczewska, A.; Palczewski, J.; Marchese Robinson, R.; Neagu, D. Interpreting random forest classification models using a feature contribution method. Integr. Reusable Syst. 2014, 263, 193–218. [Google Scholar]
  31. Marchenko, O.; Russek-Cohen, E.; Levenson, M.; Zink, R.C.; Krukas-Hampel, M.R.; Jiang, Q. Sources of safety data and statistical strategies for design and analysis: Real world insights. Ther. Innov. Regul. Sci. 2018, 52, 170–186. [Google Scholar] [CrossRef]
  32. Wong, A.K.; Au, W.H.; Chan, K.C. Discovering high-order patterns of gene expression levels. J. Comput. Biol. 2008, 15, 625–637. [Google Scholar] [CrossRef]
  33. Everitt, B.; Skkrondal, A. The Cambridge Dictionary of Statistics; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  34. Wong, A.K.; Wang, Y. High-order pattern discovery from discrete-valued data. IEEE Trans. Knowl. Data Eng. 1997, 9, 877–893. [Google Scholar] [CrossRef]
Figure 1. Principal Components ( P C 1 , P C 2 ) with AV Groups and Subgroups discovered by PDD. (a) On P C 1 , AVs with no significant statistical association with other AVs (SR less than the statistical significance threshold) are projected near the origin (marked as (a)). (b) AV Groups with statistically significant associations are projected further from the origin, indicating strong associations with other AVs (marked as (b)). (c) Subgroups, denoted by Disentangled Statistical Units (DSUs), with statistical associations among the AVs within but statistically uncorrelated to AVs of any other AV Subgroup (marked as (c)).
Figure 1. Principal Components ( P C 1 , P C 2 ) with AV Groups and Subgroups discovered by PDD. (a) On P C 1 , AVs with no significant statistical association with other AVs (SR less than the statistical significance threshold) are projected near the origin (marked as (a)). (b) AV Groups with statistically significant associations are projected further from the origin, indicating strong associations with other AVs (marked as (b)). (c) Subgroups, denoted by Disentangled Statistical Units (DSUs), with statistical associations among the AVs within but statistically uncorrelated to AVs of any other AV Subgroup (marked as (c)).
Bioengineering 11 00770 g001
Figure 2. Knowledge Base outputted by PDD. Rows highlighted in pink and green represent patterns associated with the classes sepsis and other, respectively.
Figure 2. Knowledge Base outputted by PDD. Rows highlighted in pink and green represent patterns associated with the classes sepsis and other, respectively.
Bioengineering 11 00770 g002
Figure 3. Error Detection results by PDD. (a) Knowledge Base of partial abnormal records. Rows highlighted in pink and green represent patterns associated with the classes sepsis and other, respectively. (b) Partial abnormal records of the original dataset. Blocks highlighted in pink and green represent feature values associated with the classes sepsis and other, respectively.
Figure 3. Error Detection results by PDD. (a) Knowledge Base of partial abnormal records. Rows highlighted in pink and green represent patterns associated with the classes sepsis and other, respectively. (b) Partial abnormal records of the original dataset. Blocks highlighted in pink and green represent feature values associated with the classes sepsis and other, respectively.
Bioengineering 11 00770 g003
Figure 4. Comparison of the Accuracy of Classifiers Before and After Error Detection.
Figure 4. Comparison of the Accuracy of Classifiers Before and After Error Detection.
Bioengineering 11 00770 g004
Figure 5. Top feature importances associated with class labels outputted by the Random Forest Model.
Figure 5. Top feature importances associated with class labels outputted by the Random Forest Model.
Bioengineering 11 00770 g005
Figure 6. Sample scatter pair plots of the eICU-CRD dataset. All records are plotted. The 10 abnormal records presented in Section 3.3 are highlighted, which comprise 2 outliers, 4 undecided, and 4 mislabelled records originally labelled as o t h e r , denoted by m i s l a b e l l e d _ o o . (a) d i a s t o l i c b p vs. p p . (b) d i a s t o l i c b p vs. g c s .
Figure 6. Sample scatter pair plots of the eICU-CRD dataset. All records are plotted. The 10 abnormal records presented in Section 3.3 are highlighted, which comprise 2 outliers, 4 undecided, and 4 mislabelled records originally labelled as o t h e r , denoted by m i s l a b e l l e d _ o o . (a) d i a s t o l i c b p vs. p p . (b) d i a s t o l i c b p vs. g c s .
Bioengineering 11 00770 g006
Figure 7. Sample scatter pair plots of the Data Space in the Pattern Discovery and Disentanglement (PDD)’s output. The Disentangled Space Units (DSU) in PDD are used as features, and the comprehensive pattern counts as the feature values. All records from the eICU-CRD dataset are plotted. The 10 abnormal records presented in Section 3.3 are highlighted, which comprise 2 outliers, 4 undecided, and 4 mislabelled records originally labelled as o t h e r , denoted by m i s l a b e l l e d _ o o . (a) D S U [ 1 , 1 , 1 ] vs. D S U [ 1 , 2 , 1 ] . (b) D S U [ 2 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] . (c) D S U [ 2 , 1 , 1 ] vs. D S U [ 3 , 1 , 1 ] . (d) D S U [ 1 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] .
Figure 7. Sample scatter pair plots of the Data Space in the Pattern Discovery and Disentanglement (PDD)’s output. The Disentangled Space Units (DSU) in PDD are used as features, and the comprehensive pattern counts as the feature values. All records from the eICU-CRD dataset are plotted. The 10 abnormal records presented in Section 3.3 are highlighted, which comprise 2 outliers, 4 undecided, and 4 mislabelled records originally labelled as o t h e r , denoted by m i s l a b e l l e d _ o o . (a) D S U [ 1 , 1 , 1 ] vs. D S U [ 1 , 2 , 1 ] . (b) D S U [ 2 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] . (c) D S U [ 2 , 1 , 1 ] vs. D S U [ 3 , 1 , 1 ] . (d) D S U [ 1 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] .
Bioengineering 11 00770 g007
Figure 8. Sample scatter pair plots of the Data Space in the Pattern Discovery and Disentanglement (PDD)’s output. The Disentangled Space Units (DSU) in PDD are used as features, and the comprehensive pattern counts as the feature value. All records from the eICU-CRD dataset are plotted. All the 1157 abnormal records are highlighted, which comprise 2 outliers, 159 undecided and 996 mislabelled records. m i s l a b e l l e d _ o o represents the 894 mislabelled records originally labelled as o t h e r , and m i s l a b e l l e d _ o s represents the 102 mislabelled records originally labelled as s e p s i s . (a) D S U [ 1 , 1 , 1 ] vs. D S U [ 1 , 2 , 1 ] . (b) D S U [ 2 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] . (c) D S U [ 2 , 1 , 1 ] vs. D S U [ 3 , 1 , 1 ] . (d) D S U [ 1 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] .
Figure 8. Sample scatter pair plots of the Data Space in the Pattern Discovery and Disentanglement (PDD)’s output. The Disentangled Space Units (DSU) in PDD are used as features, and the comprehensive pattern counts as the feature value. All records from the eICU-CRD dataset are plotted. All the 1157 abnormal records are highlighted, which comprise 2 outliers, 159 undecided and 996 mislabelled records. m i s l a b e l l e d _ o o represents the 894 mislabelled records originally labelled as o t h e r , and m i s l a b e l l e d _ o s represents the 102 mislabelled records originally labelled as s e p s i s . (a) D S U [ 1 , 1 , 1 ] vs. D S U [ 1 , 2 , 1 ] . (b) D S U [ 2 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] . (c) D S U [ 2 , 1 , 1 ] vs. D S U [ 3 , 1 , 1 ] . (d) D S U [ 1 , 1 , 1 ] vs. D S U [ 2 , 2 , 1 ] .
Bioengineering 11 00770 g008
Table 1. Clinician Feature Bins.
Table 1. Clinician Feature Bins.
FeaturesLow RangeNormal RangeHigh Range
GCS3–8 (Severe)9–12 (Moderate)13–15 (Mild)
Systolic BP 1<90 mm Hg90–120 mm Hg>120 mm Hg
Diastolic BP 1<60 mm Hg60–80 mm Hg>80 mm Hg
Mean BP 1<70 mm Hg70–93 mm Hg>93 mm Hg
Pulse Pressure<30 mm Hg30–40 mm Hg>40 mm Hg
Heart Rate<60 bpm60–100 bpm>100 bpm
1 Blood Pressure.
Table 2. Comparison between PDD and K-means for unsupervised learning.
Table 2. Comparison between PDD and K-means for unsupervised learning.
Unsupervised LearningPrecisionRecallAccuracyBalanced AccuracyWeighted F1-Score
PDD (28-features data)0.830.900.870.890.88
K-means (28-features data)0.430.420.520.420.54
K-means (114-features data)0.590.510.750.510.67
Table 3. Results of the Accuracy of Existing ML Models Before and After Error Detection Using PDD.
Table 3. Results of the Accuracy of Existing ML Models Before and After Error Detection Using PDD.
Classifiers (Before/After Error Detection)PrecisionRecallAccuracyBalanced AccuracyWeighted F1-Score
Random Forest0.85/0.870.73/0.760.86/0.860.73/0.760.84/0.85
SVM0.71/0.710.53/0.540.77/0.750.53/0.540.7/0.68
NN0.71/0.740.71/0.800.79/0.770.71/0.800.79/0.78
Logistic Regression0.84/0.830.79/0.800.87/0.860.79/0.800.87/0.86
LightGBM0.9/0.930.89/0.930.92/0.950.89/0.930.92/0.95
XGBoost0.89/0.940.88/0.930.92/0.950.88/0.930.92/0.95
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, P.-Y.; Lum, F.; Wang, T.J.; Bhatti, A.; Parmar, S.; Dan, C.; Wong, A.K.C. An Unsupervised Error Detection Methodology for Detecting Mislabels in Healthcare Analytics. Bioengineering 2024, 11, 770. https://doi.org/10.3390/bioengineering11080770

AMA Style

Zhou P-Y, Lum F, Wang TJ, Bhatti A, Parmar S, Dan C, Wong AKC. An Unsupervised Error Detection Methodology for Detecting Mislabels in Healthcare Analytics. Bioengineering. 2024; 11(8):770. https://doi.org/10.3390/bioengineering11080770

Chicago/Turabian Style

Zhou, Pei-Yuan, Faith Lum, Tony Jiecao Wang, Anubhav Bhatti, Surajsinh Parmar, Chen Dan, and Andrew K. C. Wong. 2024. "An Unsupervised Error Detection Methodology for Detecting Mislabels in Healthcare Analytics" Bioengineering 11, no. 8: 770. https://doi.org/10.3390/bioengineering11080770

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop