taqa-interface/boxplotCarlos.m at master · fofoni/taqa-interface · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
%Rodando script do Fonini
resultados;

notas = NOTAS;
sinais = SINAIS;

notas(notas == -Inf) = NaN;
[~, I] = sort(nanmean(notas));
sinais = sinais(I);
notas = notas(:,I);
boxplot(notas)
axis([0,45, 0.9, 5.1]);
xlabel({ 'Número do Sinal', '(olhar tabela de correspondência)'});
ylabel('Nota na escala MOS');
title('Notas para cada sinal de teste');


nomeSinal = {};
parametros = {};
tipoTeste = {};
label = {};
notasColuna = [];

for i = 1:length(sinais)
    notasTeste = notas(:,i);
    notasTeste = notasTeste(~isnan(notasTeste));
    a = cellstr(repmat(sinais(i).sinal, length(notasTeste),1));
    b = cellstr(repmat(sinais(i).teste, length(notasTeste),1));
    c = cellstr(repmat(sinais(i).parametro, length(notasTeste),1));

    nomeSinal = vertcat(nomeSinal , a);
    tipoTeste = vertcat(tipoTeste, b);
    parametros = vertcat(parametros, c);

    notasColuna = [notasColuna; notasTeste];
    label = vertcat(label, strcat(sinais(i).sinal,'_', sinais(i).teste,'_', sinais(i).parametro));
end

%Sinal = strcat(nomeSinal,'_', tipoTeste,'_', parametros);
%figure;
%boxplot(notasColuna, Sinal);
%figure;
%boxplot(notasColuna, parametros);
figure;
for i= 1: length(tipoTeste)
    teste_temp = tipoTeste{i};
    if strcmp(teste_temp, 'CODEC4KHZ')
        parametros{i} = 'CODEC 4K';
    end
end
boxplot(notasColuna, parametros);
xlabel('Parâmetro testado');
ylabel('Nota na escala MOS');
title('Notas por parâmetro testado');
figure;
boxplot(notasColuna, parametros, 'GroupOrder', {'20dB', '35dB', '50dB', '64k', '128k', '320k', 'CODEC 4K'});
xlabel('Parâmetro testado');
ylabel('Nota na escala MOS');
title('Notas por parâmetro testado');