% Calculo del umbral de decision T que iguala las dos areas por debajo de
% las funciones densidad de probabilidad de la distribucion normal y de la
% distribucion de Weibull

%Fecha de creación: 25 de Febrero de 2007

n=load('Resultados\prueba_eleccion_umbral_pdf_12anuncios.mat'); %variable de partida que contiene el resultado de la prueba
                                            %prueba_eleccion_umbral_pdf
long=length(n.n);
corte=700;           % indice a partir del cual se discrimina entre valores para el calculo de Pf y Pn

%Seleccionamos la parte correspondiente a la probabilidad Pn                                            
n1=n.n;
n1(corte:long)=zeros(1,long-corte+1);

%Seleccionamos la parte correspondiente a la probabilidad Pf
n2=n.n;
n2(1:corte)=zeros(1,corte);

%Transformamos n1 y n2 en vectores que me permitan el calculo de los
%parametros de las dos distribuciones consideradas (normal y weibull)

vect1=transforma_formato(n1);
vect2=transforma_formato(n2);

%Calculo de los parametros que mejor aproximan las distribuciones
[a,b]=wblfit(vect1);
[c,d]=normfit(vect2);

i=1;
%Calculo de las probabilidades Pf y Pn en funcion de alpha
for alpha=0.03:0.01:0.45
    prob1(i)=1-wblcdf(alpha,a(1,1),a(1,2));
    prob2(i)=normcdf(alpha,c,d);
    i=i+1;
end

plot([0.03:0.01:0.45],prob1);
hold on;
plot([0.03:0.01:0.45],prob2);
hold off;

figure
plot_dist_weibull(a(1,1),a(1,2),long);
hold on;
plot_dist_normal(c,d,long);
hold off;

