% matlab code

% Scores and analysis for listening test included with
% "Decoding and Compression of Channel and Scene Objects for Spatial Audio"

% 8 subjects, 3 conditions

similarity_score = flip ([ 9 5 3; 10 7 4; 8 3 1; 6 6 5; 10 5 1; 10 7 9; 10 7 1; 10 5 2 ], 2 ); 
    
stability_score = flip( [ 9 3 2; 10 5 2; 8 8 2; 7 6 6; 10 5 1; 10 7 5; 10 6 1; 10 4 2 ], 2 );


figure( 'Position', [0, 0, 450, 200] );

subplot(1,2,1);
boxplot(similarity_score, 'Labels', {'VL','RQ','VBAP'}, 'Whisker',1000);  %  'Notch','on'
title('Localisation')

%figure( 'Position', [0, 0, 200, 200]);
subplot(1,2,2);
boxplot(stability_score, 'Labels', {'VL','RQ','VBAP'}, 'Whisker',1000 );
title('Stability')


[h,p] = ttest(similarity_score(:,1), similarity_score(:,2))

[h,p] = ttest(stability_score(:,1), stability_score(:,2))


N = size(similarity_score, 1);
M = sum( similarity_score(:,2) > similarity_score(:,1) );

p_similarity2GT1 = 1 - binocdf(M-1, N, 0.5)


N = size(stability_score, 1);
M = sum( stability_score(:,2) > stability_score(:,1) );

p_stability2GT1 = 1 - binocdf(M-1, N, 0.5)



[h,p] = ttest(similarity_score(:,2), similarity_score(:,3))

[h,p] = ttest(stability_score(:,2), stability_score(:,3))


N = size(similarity_score, 1);
M = sum( similarity_score(:,3) > similarity_score(:,2) );

p_similarity2GT1 = 1 - binocdf(M-1, N, 0.5)


N = size(stability_score, 1);
M = sum( stability_score(:,3) > stability_score(:,2) );

p_stability2GT1 = 1 - binocdf(M-1, N, 0.5)