% wymaga zainstalowanego Statistics Toolbox (funkcje ttest2 i ranksum)
x=[12 9 21 14 7 17];
y=[5  9 3  11 8 19];
[h p ci stats]= ttest2(x,y,.05,1);
fprintf('pierwsze dane/test t: p=%.4f\n', p);
[p, h, stats]=ranksum(x,y);
fprintf('pierwsze dane/test Wilcoxona-Manna-Whitneya: p=%.4f\n', p);
plot([x y], '+');
line([length(x)+.5, length(x)+.5], [min([x y])-1, max([x y])+1]);
line([0, length(x)+.5], [mean(x) mean(x)], 'linestyle','--');
line([length(x)+.5 length(x)+length(y)+1], [mean(y) mean(y)], 'linestyle','--');
ylabel('liczba braków'); xlabel('dni');text(length(x),min(y),'remont')
% test permutacyjny
[roznice, p]= perm_test(x,y); 
%;disp([ mean(x) mean(y) mean(x)-mean(y) p]);
fprintf('pierwsze dane/test permutacyjny: p=%.4f\n', p);
figure(2); hist(roznice);
figure(3); hist(roznice,54);
%dorysowujemy obszar krytyczny za 5.2:
p=patch([5.2 10 10 5.2], [0 0 40 40], 'c', 'erasemode', 'xor');
figure(4); subplot 211; hist(roznice,54);
p=patch([5.2 10 10 5.2], [0 0 40 40], 'c', 'erasemode', 'xor');
subplot 212; hist(roznice,54);
p=patch([6.3 10 10 6.3], [0 0 40 40], 'c', 'erasemode', 'xor');
p=patch([-6.3 -10 -10 -6.3], [0 0 40 40], 'c', 'erasemode', 'xor');

% po zebraniu dodatkowych danych:
y=[5  9 3  11 8 19 7 5 9 12 5 11 9 6 8 17 9 12];
[h p ci stats]= ttest2(x,y,.05,1);
fprintf('wiêcej danych/test t: p=%.4f\n', p);
[p, h, stats]=ranksum(x,y,.05);
fprintf('wiêcej danych/test Wilcoxona-Manna-Whitneya: p=%.4f\n', p);
figure(5); plot([x y], '+');
line([length(x)+.5, length(x)+.5], [min([x y])-1, max([x y])+1]);
line([0, length(x)+.5], [mean(x) mean(x)], 'linestyle','--');
line([length(x)+.5 length(x)+length(y)+1], [mean(y) mean(y)], 'linestyle','--');
ylabel('liczba braków'); xlabel('dni');text(length(x),min(y),'remont')
set(gca,'xlim', [0 length(x)+length(y)+1], 'ylim', [min([x y])-1 max([x y])+1])
% test permutacyjny
clear roznice; [roznice, p]=perm_test(x,y);
fprintf('wiêcej danych/test permutacyjny: p=%.4f\n', p);
figure(6); hist(roznice,20);
