--- a
+++ b/MATLAB/analysis/findDuplicateSubjects.m
@@ -0,0 +1,61 @@
+%FINDDUPLICATESUBJECTS finds duplicate subjects in the database.
+%
+% AUTHOR: Maximilian C. M. Fischer
+% COPYRIGHT (C) 2023 Maximilian C. M. Fischer
+% LICENSE: EUPL v1.2
+%
+
+clearvars; close all
+
+addpath(genpath('..\src'))
+VSD_addPathes('..\..\..\..\..\')
+
+% Load subjects & meta data
+subjectXLSX = '..\res\VSD_Subjects.xlsx';
+Subjects = readtable(subjectXLSX);
+
+%% Import
+NoS = size(Subjects, 1);
+patchProps.FaceAlpha = 0.5;
+TFM2 = repmat(struct('APP',[]),NoS,1);
+Sacrum = repmat(struct('vertices',[],'faces',[]),NoS,1);
+for s=1:NoS
+    % Import the bones
+    load(['..\..\Bones\' Subjects.ID{s} '.mat'], 'B')
+    % Import the transformation to the APP coordinate system
+    load(['D:\Biomechanics\Hip\Code\AcetabularMorphologyAnalysis\data\' Subjects.ID{s} '.mat'], 'TFM2APP')
+    TFM2(s).APP = TFM2APP;
+    % Transform the sacrum to the APP coordinate system
+    Sacrum(s) = transformPoint3d(splitMesh(...
+        B(ismember({B.name},'Sacrum')).mesh,'maxBoundingBox'), TFM2(s).APP);
+end
+clearvars B TFM2APP
+
+% Compare the sacrum of each subject using an ICP registration
+patchProps.FaceAlpha = 0.5;
+SS = 1:NoS;
+rigidRegRMSE = nan(NoS);
+for s=1:NoS
+    SS(SS==s)=[];
+    for ss=SS
+        rigidReg = icp(Sacrum(s),Sacrum(ss),'Plot',false,'ChangeRate',1e-2,'Delay',0);
+        rigidRegRMSE(s,ss) = rigidReg.rmse;
+    end
+end
+
+save('rigidRegRMSE','rigidRegRMSE')
+
+%% Evaluate results
+load('rigidRegRMSE','rigidRegRMSE')
+
+medianRMSE = median(rigidRegRMSE(:), 'omitnan');
+outlierCutOff = prctile(rigidRegRMSE(:),25)-1.5*iqr(rigidRegRMSE(:));
+[row,col] = ind2sub(size(rigidRegRMSE),find(rigidRegRMSE<outlierCutOff));
+for d=1:length(row)
+    disp([...
+        'Subject ' Subjects.ID{row(d)} ' and ' Subjects.ID{col(d)} ...
+        ' are dupes based on a comparison of the sacrum!'])
+end
+
+% [List.f, List.p] = matlab.codetools.requiredFilesAndProducts([mfilename '.m']);
+% List.f = List.f'; List.p = List.p';
\ No newline at end of file