% Chcemy zbadaæ czy cztery metody resocjalizacji 
% w istotny sposób daj± ró¿ne wyniki. Badamy ka¿d± na
% grupie 20 m³odocianych przestêpców. Jako wynik negatywny  
% traktujemy ilo¶æ osób, które po terapii ponownie pope³ni³y 
% przestêpstwa.
% Otrzymane wyniki s± nastêpuj±ce: 
 x(1)=3 ;
x(2)=10;
x(3)=10;
x(4)= 13;
%Pytanie: czy która¶ z terapii jest lepsza? \\
% Wskazówka: 
% H0: Wszystkie wyniki pochodz± z tej samej populacji
%  --- nasza najlepsza wiedza o tej populacji to proporcja 
%      ponownych przestêpców (3+10+10+13)/80.  
%      Mo¿emy badaæ najwiêksz± ró¿nicê pomiêdzy czterema próbami 
%      lub odstêpstwo od ¶redniej ka¿dej z prób osobno. 

w=zeros(1,80);
w(1:36)=1;
S=max(x)-min(x);

N_rep=1e4;
N=20;
S_dist=zeros(1,N_rep);
for i=1:N_rep
    for j=1:4
        x_rep(j)=sum(randsample(w,N,true));
    end
    S_dist(i)=max(x_rep)-min(x_rep);
end

p=length(find(S_dist>=S))/N_rep