p_fisher.m
p_fisher.m — Objective-C source code, 674 bytes
Contenuto del file
p1=p_cut1; %inizializzo la probabilit� (1) p2=p_cut2; %inizializzo la probabilit� (1) %ciclo per tutte le tabelle di contingenza equivalenti for n=1:size(M,3) a=M(1,1,n); %calcolo elemeti (solo per maggior chiarezza ei passaggi successivi) b=M(1,2,n); c=M(2,1,n); d=M(2,2,n); test1=nchoosek(a+b,a)*nchoosek(c+d,c)/nchoosek(a+b+c+d,a+c) %calcolo la prob per la tabella selezionata test2=pdf('Hypergeometric',a,a+b+c+d,a+b,a+c) if test1<=p_cut1 %verifico se la prob della tabella � pi� estrema o ugualmente probabile p1=p1+test1; end if test2<=p_cut2 p2=p2+test2; end end p1 p2