%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Native samples %%%%%%%%%%%%%%%%%%%%%% >> Y = pdist(X,'spearman'); squareform(Y) Z = linkage(Y,'average'); dendrogram(Z, 'Labels',{'Lame397';'Lame3982';'iTunes';'FhG';'l3enc';'Helix'}) Coph = cophenet(Z,Y) trineqtest(Y) ans = 0 0.0428 0.2790 0.1730 0.2526 0.1360 0.0428 0 0.2230 0.1519 0.2209 0.0903 0.2790 0.2230 0 0.1784 0.3661 0.1448 0.1730 0.1519 0.1784 0 0.3074 0.0847 0.2526 0.2209 0.3661 0.3074 0 0.2849 0.1360 0.0903 0.1448 0.0847 0.2849 0 Coph = 0.8814 Triangle inequality holds for 116 cases from 120 possible (96.67%). >> % Helix Lame3982 FhG Lame397 XX = [X(6,:); X(2,:); X(4,:); X(1,:)]; Y = pdist(XX,'spearman'); AvrgD = mean(Y) MaxD = max(Y) trineqtest(Y) AvrgD = 0.1131 MaxD = 0.1730 Triangle inequality holds for 22 cases from 24 possible (91.67%). % ----------------------------------------------------------------- >> Y = pdist(X,@distancecorr); squareform(Y) Z = elinkage(Y); % e-link dendrogram(Z, 'Labels',{'Lame397';'Lame3982';'iTunes';'FhG';'l3enc';'Helix'}) Coph = cophenet(Z,Y) trineqtest(Y) ans = 0 0.0611 0.3051 0.1918 0.3152 0.1762 0.0611 0 0.2549 0.1675 0.2593 0.1130 0.3051 0.2549 0 0.1873 0.4179 0.1559 0.1918 0.1675 0.1873 0 0.3347 0.0954 0.3152 0.2593 0.4179 0.3347 0 0.3302 0.1762 0.1130 0.1559 0.0954 0.3302 0 Coph = 0.8197 Triangle inequality holds for 118 cases from 120 possible (98.33%). >> % Helix Lame3982 FhG Lame397 XX = [X(6,:); X(2,:); X(4,:); X(1,:)]; Y = pdist(XX,@distancecorr); AvrgD = mean(Y) MaxD = max(Y) trineqtest(Y) AvrgD = 0.1342 MaxD = 0.1918 Triangle inequality holds for 22 cases from 24 possible (91.67%). % Helix Lame3982 FhG Lame397 df = [-24.8237 -22.1923 -23.8205 -20.7318] sq = [4.59 4.51 4.44 4.28] sqa = sepsy3(df, sq, df) ((sqa-sq) ./ sq) * 100 RMSE = sqrt(mean((sqa-sq).^2)) % 4.5900 4.5100 4.4400 4.2800 % 4.5712 4.4129 4.5108 4.3251 % -0.4100 -2.1527 1.5955 1.0530 % 0.0649 Helix Lame3982 FhG Lame397 RMSEr True QS 4.5900 4.5100 4.4400 4.2800 Computed QS 4.5712 4.4129 4.5108 4.3251 Error -0.41% -2.15% +1.60% +1.05% 0.0649 %%%%%%%%%%%%%%%%%%%%%%%%%% The random mix %%%%%%%%%%%%%%%%%%%%%%% >> X = [df_lame97, df_lame982, df_itunes, df_fhg, df_l3enc, df_helix]; >> X = X'; >> Y = pdist(X,'spearman'); squareform(Y) Z = linkage(Y,'average'); dendrogram(Z, 'Labels',{'Lame397';'Lame3982';'iTunes';'FhG';'l3enc';'Helix'}) Coph = cophenet(Z,Y) trineqtest(Y) ans = 0 0.0371 0.4628 0.1196 0.1632 0.0674 0.0371 0 0.4863 0.1168 0.1513 0.0490 0.4628 0.4863 0 0.5655 0.5064 0.4663 0.1196 0.1168 0.5655 0 0.2530 0.0900 0.1632 0.1513 0.5064 0.2530 0 0.1670 0.0674 0.0490 0.4663 0.0900 0.1670 0 Coph = 0.9861 Triangle inequality holds for 118 cases from 120 possible (98.33%). >> % Helix Lame3982 FhG Lame397 XX = [X(6,:); X(2,:); X(4,:); X(1,:)]; Y = pdist(XX,'spearman'); AvrgD = mean(Y) MaxD = max(Y) trineqtest(Y) AvrgD = 0.0800 MaxD = 0.1196 Triangle inequality holds for 24 cases from 24 possible (100.00%). % ----------------------------------------------------------------- >> Y = pdist(X,@distancecorr); squareform(Y) Z = elinkage(Y); % e-link dendrogram(Z, 'Labels',{'Lame397';'Lame3982';'iTunes';'FhG';'l3enc';'Helix'}) Coph = cophenet(Z,Y) trineqtest(Y) ans = 0 0.0459 0.4765 0.1494 0.1955 0.0883 0.0459 0 0.4916 0.1424 0.1770 0.0658 0.4765 0.4916 0 0.5773 0.5116 0.4772 0.1494 0.1424 0.5773 0 0.2942 0.1103 0.1955 0.1770 0.5116 0.2942 0 0.1993 0.0883 0.0658 0.4772 0.1103 0.1993 0 Coph = 0.9827 Triangle inequality holds for 120 cases from 120 possible (100.00%). >> % Helix Lame3982 FhG Lame397 XX = [X(6,:); X(2,:); X(4,:); X(1,:)]; Y = pdist(XX,@distancecorr); AvrgD = mean(Y) MaxD = max(Y) trineqtest(Y) AvrgD = 0.1004 MaxD = 0.1494 Triangle inequality holds for 24 cases from 24 possible (100.00%). % Helix Lame3982 FhG Lame397 df = [-30.7113 -27.1183 -29.4716 -25.3550] sq = [4.59 4.51 4.44 4.28] sqa = sepsy3(df, sq, df) ((sqa-sq) ./ sq) * 100 RMSE = sqrt(mean((sqa-sq).^2)) % 4.5900 4.5100 4.4400 4.2800 % 4.5674 4.4088 4.5127 4.3310 % -0.4919 -2.2429 1.6376 1.1922 % 0.0682 Helix Lame3982 FhG Lame397 RMSEr True QS 4.2530 4.0390 3.6980 1.5450 Computed QS 4.5674 4.4088 4.5127 4.3310 Error -0.49% -2.24% +1.64% +1.19% 0.0682 Native samples The random mix Average / Max Spearman's correlation distance 0.1131 / 0.1730 0.0800 / 0.1196 Average / Max Distance correlation distance 0.1342 / 0.1918 0.1004 / 0.1494