function [] = Estatisticas (password, Delta, Fourrier,fator,dirIn, listaIn, dirOut, listaOut,tempo)

%Recebendo a lista de arquivos de audio
dirRelat = 'C:\warter\Relatorio\';

%delta = strcat('Delta: ', num2str(Delta));
%fat = strcat('Fator: ', num2str(fator));
%disp (delta);
%disp (fat);



segmento = 1;
%Inicializando variaveis para verificaçao de ruido
 maxSNR = 0;
 maxk = 0;
 
for i=3:length (listaIn)
    
    %Lendo os dados do arquivo do arquivo NAO Marcado
    nomeIn = strcat(dirIn,listaIn(i).name);
    [tam Fs bits] = wavread(nomeIn,'size');
    %disp (nomeIn);

        
    %Nome do arquivo Marcado
    nomeOut = strcat (dirOut,listaOut(i).name);
    %disp (nomeOut);


    %Dividindo o arquivo original em pedaços
    resto = rem(tam(1),tempo); 
    
    %Gerando os dados estatisticos
    for j=1:tempo:(tam(1)-resto)

   
        y = wavread(nomeIn,[j (j+tempo-1)]);
        
        % Gerando os conjuntos A e B que sao indices do vetor y
        [A,B] = GeraConjuntos(length(y), (password+j), fator);
        
        S1 = y(A) - y(B);
        meanS1(segmento) = mean(S1);
        stdS1(segmento) = std (S1);

        yMarcado = wavread(nomeOut,[j (j+tempo-1)]);
        
        S2 = yMarcado(A) - yMarcado(B); 
        meanS2(segmento) = mean(S2);
        stdS2(segmento) = std (S2);
  
        %[votar(segmento),pValue(segmento)] = ttest(S2,0);  

        %Caluculando o SNR(signal to noise ratio)
        % x2 - Nao Marcado (Referencial)
        % y2 - Marcado
        x2 = y .* y;
        x2 = sum(x2);
        x2 = sqrt(x2);
        
        
        y2 = yMarcado .* yMarcado;
        y2 = sum(y2);
        y2 = sqrt(y2);
        
        if y2 == 0
            SNR = 0;
        elseif x2 == 0
            SNR = 0;
        else
            % Calculando o valor do SNR em decibel
            SNR = 10 * log10 (y2/x2);
        end
        
        if (SNR > maxSNR)
                maxSNR = SNR;
                maxk = segmento;
        end
        segmento = segmento + 1;
    end
    %Posicionado j na primeira posiçao ainda nao lida do audio original
    j = j + tempo;
    
    %Gravando o restante do audio, caso exista. 
    if (j < tam(1))
        
        y = wavread(nomeIn, [j tam(1)]);
        
        % Gerando os conjuntos A e B que sao indices do vetor y
        [A,B] = GeraConjuntos(length(y), password, fator);
        
        S1 = y(A) - y(B);
        meanS1(segmento) = mean(S1);
        stdS1(segmento) = std (S1);

        yMarcado = wavread(nomeOut,[j tam(1)]);
        
        S2 = yMarcado(A) - yMarcado(B); 
        meanS2(segmento) = mean(S2);
        stdS2(segmento) = std (S2);

        %[votar(segmento),pValue(segmento)] = ttest(S2,0);  

        %Caluculando o SNR(signal to noise ratio)
        % x2 - Nao Marcado (Referencial)
        % y2 - Marcado
        x2 = y .* y;
        x2 = sum(x2);
        x2 = sqrt(x2);
        
        y2 = yMarcado .* yMarcado;
        y2 = sum(y2);
        y2 = sqrt(y2);
        
        if y2 == 0
            SNR = 0;
        elseif x2 == 0
            SNR = 0;
        else
            % Calculando o valor do SNR em decibel
            SNR = 10 * log10 (y2/x2);
        end

        if (SNR > maxSNR)
            maxSNR = SNR;
            maxk = segmento;
        end        
        segmento = segmento + 1;
    end
end

%Arquivo de Saida
nomeRelat = strcat (dirRelat,'ResultadosTeste.xls');
fid = fopen(nomeRelat,'a');

fprintf (fid,'%.3f\t%d\t%d\t%.5f\t%.5f\t', Delta, ((floor (tempo/2)*fator)-1),segmento, mean(meanS1), mean(meanS2));
fprintf (fid,'%.5f\t%.5f\t',mean(stdS1), mean(stdS2));

%Testando se distribuiçao e normal
[JmeanS1, LmeanS1] = TestaNormal(meanS1);
[JmeanS2, LmeanS2] = TestaNormal(meanS2);

fprintf (fid,'%d\t%d\t%d\t%d\t',JmeanS1, LmeanS1, JmeanS2, LmeanS2);


%Teste Anova
%NomesGrafico = {'Media Alt.' 'Media Nao Alt.' };
%aux = [meanS1' meanS2'];
%ValorAnova = anova1(aux, NomesGrafico, 'off');
%if (ValorAnova < 0.05)
%marcado = 1;
%else 
%    marcado = 0;
%end

%fprintf (fid,'%.5f\t%d\t',ValorAnova,marcado);

%Valor dos SNR
fprintf (fid,'%.5f\t',maxSNR);
%Realizando o teste de hipotese
% 0 - Nao Marcado (Nao posso rejeitar a hipostes)
% 1 - Marcado (Posso rejeita a hipotese)
[marcado,pValue] = ttest(meanS2,0);  
%marcado = sum(votar)/length (votar);
%disp (marcado);
%if (marcado > 0.9)
    %marcado = 1;
    %else
    %marcado = 0;
    %end
%pValue = mean (pValue);

fprintf (fid,'%d\t%.5f\t\t%d\n',marcado, pValue,password);

%testando se Media alterada e a media nao alterada sao iguais
%Assumimos que os std sao iguais apesar de desconhecidos
%1 - pode rejeitar a hipotese
%0 - nao pode rejeitar a hipotese
%Resultado esperado para marcado e 1
%[marcado,pValue] = ttest2(meanS1,meanS2);  
%fprintf (fid,'%d\t%.5f\t',marcado, pValue);

%testando se o Media nao alterada e zero
% Resultado esperado para marcado e 0
%[marcado,pValue] = ttest(meanS1,0);  
%fprintf (fid,'%d\t%.5f\t\t%d\n',marcado, pValue,password);



fclose(fid);

