%myszy
clear;
myszy_spec=[94 38 23 197 99 16 141];
myszy_reg=[52 10 40 104 51 27 146 30 46];
Nboot=100000;
A=zeros(1,Nboot);
%B=zeros(1,Nboot);
for i=1:Nboot
%wartosci1=randsample(myszy_spec,7,true);
%wartosci2=randsample(myszy_reg,9,true);
%A(i)=mean(wartosci1)-mean(wartosci2);
A(i)=mean(randsample([myszy_spec myszy_reg],7,true))-mean(randsample([myszy_spec myszy_reg],9,true));
%B(i)=mean(wartosci2);

end 
%roznica=mean(A)
lo=prctile(A,2.5);
hi=prctile(A,97.5);
rozn=mean(myszy_spec)-mean(myszy_reg);
%disp(sprintf('na populacji: %.3f\n',rozn))
disp(sprintf('srednia i 95 procentowy przedzial %.3f %.3f %.3f\n',lo,rozn,hi))
%-mean(B)
nowa_delta=sqrt(1/length(myszy_spec)*var(myszy_spec)+1/length(myszy_reg)*var(myszy_reg));
f=length(myszy_spec)+length(myszy_reg)-2;
lo2=tinv(0.025,f);
hi2=tinv(0.975,f);
disp(sprintf('srednia i 95 procentowy przedzial z zalozeniem normalnosci %.3f %.3f %.3f\n',lo2*nowa_delta,rozn,hi2*nowa_delta))