
% Czy przedstawione wyniki potwierdzaj±, ¿e w grupie leczonej 
% wagi narodzeniowe s± wy¿sze ni¿ w grupie kontrolnej?
% Skonstruowaæ test permutacyjny (mieszamy losowo przynale¿no¶æ do grup) i bootstrapowy (losujemy z powtórzeniami ze wspólnego \emph{universum}).
X=[6.9	 6.4
	7.6 	6.7
 	7.3 	5.4 
	7.6 	8.2 
	6.8 	5.3 
	7.2 	6.6 
	8.0 	5.8 
	5.5 	5.7 
	5.8 	6.2 
	7.3 	7.1 
	8.2 	7.0 
	6.9 	6.9 
	6.8 	5.6 
	5.7 	4.2 
	8.6	6.8];

Lek=X(:,1);
NLek=X(:,2);
w=[Lek' NLek'];
N=length(w);
S=mean(Lek)-mean(NLek);

% permutacyjnie

Lek_ix=1:N/2;
NLek_ix=N/2+1:N;

N_rep=1e5;
s_dist=zeros(1,N_rep);
tic
for i=1:N_rep
    ww=w(randperm(N));
    s_dist(i)=mean(ww(Lek_ix))-mean(ww(NLek_ix));
end
p=length(find(s_dist>=S))/N_rep
 toc
 tic

% bootstrapowo

for i=1:N_rep
    ww=w(ceil(N*rand(1,N)));%randsample(w,N,true);
    s_dist(i)=mean(ww(Lek_ix))-mean(ww(NLek_ix));
end
p=length(find(s_dist>=S))/N_rep
toc