% Badamy skuteczno¶æ leku na raka. 
% Mamy grupê 12 chorych: 
% 6 osobom podajemy lek --- poprawa wyst±pi³a u 5 osób,  
% pozosta³ym sze¶ciu osobom podajemy placebo --- poprawa wyst±pi³a u 2 osób.
% Czy te wyniki upowa¿niaj± do stwierdzenia, 
% ¿e lek istotnie zwiêksza szansê poprawy?

% Zakladamy hipotezy 
% H0: lek nie daje poprawy 
% H1: lek daje poprawe

% zgodnie z H0 obie proby pochodza ze swiata:
w=[ones(1,7) zeros(1,5)]; % jedynki = wystapila poprawa

% reprobkujemy
N_rep=1e5;

n_l=5; % ilosc popraw w grupie leku
n_p=2; % ilosc popraw w grupie placebo
% jako statystyke testowa przyjmiemy roznice w poprawach miedzy grupami
% w tym problemie istotne jest zwiekszenie ilosci popraw wiec stosujemy test
% jednostronny

st_0= n_l - n_p;
st_rep=zeros(1,N_rep);
for i=1:N_rep
    n_l_rep   = sum(randsample(w,6,true));
    n_p_rep   = sum(randsample(w,6,true));
    st_rep(i) = n_l_rep - n_p_rep; % wartosc statystyki w i-tym reprobkowaniu
end

% jaka proporcja statystyk reprobkowanych daje wynik taki jak oryginalnie
% lub bardziej ekstremalny?
p=length(find(st_rep>=st_0))/N_rep


