Siamese neural network-enhanced electrocardiography can re-identify anonymized healthcare data

Figures

• 

  Figure 1. The SNN architecture and triplet loss function: the SNN receives three ECG inputs: a reference ECG from a given subject A (anchor), another ECG from the same subject A (positive), and an ECG from a different subject B (negative) and encodes embeddings for these using three identical CNNs. The Euclidean distances between anchor and positive embeddings and between anchor and negative embeddings are calculated, with the triplet loss function updating the weights of the CNN to encode anchor-positive embeddings closer together than anchor-negative embeddings. AN, anchor-negative; AP, anchor-positive; CNN, convolutional neural network; SNN, Siamese neural network.

• 

  Figure 2. Re-identification of a subject from an anonymized dataset: (A) A visualization of the re-identification process during which Euclidean distances are calculated between a subject’s reference ECG and all other ECGs in an anonymized dataset (where only one also belongs to that subject). The Euclidean distance for the only correct ECG pair (i.e. that belonging to the same subject) is shown in green, and the Euclidean distances for all other incorrect ECG pairs is shown in red. (B) SNN-ECG is only successful at re-identification when the Euclidean distance between the two ECGs belonging to the reference subject is the shortest. (C) In a real-world setting, where the correct ECG from the anonymized dataset is not known, a certainty score is output alongside the model’s top choice to aid interpretation. As the discrepancy between top two ECGs from the anonymized dataset that are closest to the reference ECG increases (i.e. their Euclidean distances are most different), the model is more confident of its top choice. SNN, Siamese neural network.

• 

  Figure 3. Distributions of AP and AN pair normalized Euclidean distances: these are shown for (A) all hold-out test subjects (n subjects = 21 737) and subsets of (B) outpatient (n = 17 085), (C) normal (n = 3438), (D) LBBB (n = 1438), (E) RBBB (n = 2139), and (F) AF (n = 4206) subjects. The binary thresholds, overall accuracies, and confusion matrices are shown. If there was perfect discrimination between AP and AN pairs, there would be no overlap between the two histograms. AF, atrial fibrillation; AN, anchor-negative; AP, anchor-positive; LBBB, left bundle branch block; RBBB, right bundle branch block.

• 

  Figure 4. Evolution of ECG similarity over time: data for two subjects with multiple ECGs, from the hold-out test set, are shown. (A) and (B) are from one subject, and (C) and (D) are from a second subject. (A) and (C) show the normalized Euclidean distances for all ECG pairs. A point represents a unique ECG pair and is coloured by normalized Euclidean distance as output by SNN-ECG with smaller distances in green and larger distances in brown. The chronologically successive ECG pairs are plotted along the bottom-left to top-right diagonal. The normalized Euclidean distance generally increases towards the top-left corner with increasing time between the acquisition of ECGs in a given AP pair. (B) and (D) show four ECG median beat traces for each subject. In (B), these correspond to the first recorded ECG, subsequent ECGs from dates i and ii, and the last recorded ECG. In (D), these correspond to the first recorded ECG, a subsequent ECG from date iii, a further ECG also from date iii, and the last recorded ECG. The normalized Euclidean distance and temporal difference (in days) between the four ECGs is shown. A normalized Euclidean distance < 1 indicates that the model correctly predicts…