% 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;
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
