
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
% Badamy dwie nowe karmy dla swin; nazwijmy je A i B. 
% Mamy dwie grupy po 12 zwierzat. 
% Uzyskane przyrosty masy sa nastepujace: 
%  A: 31 34 29 26 32 35 38 34 31 29 32 31 
%  B: 26 24 28 29 30 29 31 29 32 26 28 32 
%  Czy która? z karm daje istotnie wi?ksze przyrosty masy? 
A=[ 31 34 29 26 32 35 38 34 31 29 32 31 ];
B=[ 26 24 28 29 30 29 31 29 32 26 28 32 ];
ranksum(A,B)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
disp(' Czy stosunek do marihuany sie zmienil? ')
% 
% Rozwazmy dwie ankiety przeprowadzone w USA, pytano 
% 1500 respondentów o stosunek do legalizacji marihuany. 
% Pierwsza ankiete przeprowadzono w 1980, wówczas za legalizacja
% opowiadalo sie 52% a druga w 1985 i za legalizacja bylo 46%. 
% Czy wyniki tych dwóch ankiet sa istotnie rózne? 
% 
% Z jaka proporcja powinnismy porównywac te wyniki? 
% Jaka jest hipoteza zerowa? 

ank1=zeros(1500,1);
ank1(1:0.52*1500)=1; %za
ank2=zeros(1500,1);
ank2(1:0.46*1500)=1; %za
ranksum(ank1,ank2)


% ank=zeros(1500,1);
% ank(1:0.49*1500)=1; %za
% Nboot=10000;
% for i=1:Nboot
%     a1=ank(ceil(1500*rand(1500,1)));
%     a2=ank(ceil(1500*rand(1500,1)));
%     w(i)=abs(sum(a1)-sum(a2))/1500;
% end
% p=length(find(w>=(0.52-0.46)))/Nboot


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
disp('Firmy')
% Ponizsze dane przedstawiaja wyniki finansowe uzyskiwane przez 33 zestawy firm.
% Pierwszy typ to firmy po polaczeniu, drugi to firmy o rozmiarze zblizonym do rozmiaru firm, 
% które ulegly polaczeniu, trzeci to firmy o rozmiarze zblizonym do rozmiaru firmy powstalej 
% po polaczeniu. 



A=[    1   -0.2000         0.0256          0
    2   -0.3483         -0.1250         0.0805
    3   0.0751          0.0632          -0.0231
    4   0.1261          -0.0420         0.1647
    5   -0.1017         0.0800          0.2778
    6   0.0378          0.1491          0.4302
    7   0.1162          0.1518          0.1429
    8   -0.0984         0.0377          0.0400
    9   0.0214          0.0766          0.0111
   10   -0.0171         0.2843          0.1891
   11   -0.3648         0.1391           0.0389
   12   0.0881          0.0387          0.0948
   13   -0.2632         0.0564          0.0451
   14   -0.0494         0.0537          0.0083
   15   0.0115          0.0481          0.0948
   16   0.0097          0.1982          0.0609
   17   0.0714          0.4208          -0.0248
   18   0.0018          0.0743          0.0532
   19   0.0048          -0.0071         0.0501
   20   -0.0540         0.1715          0.1095
   21   0.0227          0.0279          -0.0225
   22   0.0598          0.0486          0.1671
   23   -0.0599         0.0264          0.0207
   24   -0.0886         -0.0593         0.0771
   25   -0.0248         -0.0184         0.0596
   26   0.0764          0.0126          0.0346
   27   -0.0017         -0.0455         0.0536
   28   -0.2198         0.3431          0.0428
   29   0.3824          0.2210          0.1158
   30   -0.0068         0.2549          0.2370
   31   -0.1630         0.0112          0.1905
   32   0.1918          0.1505          0.1520
   33   0.0612          0.1704          0.0935]

laczone=A(:,2);

male=A(:,3);
duze=A(:,4);
disp('laczone male duze ')
disp([ median(laczone) median(male) median(duze)])
m=signtest(laczone,male) 
m2=2*m
d=signtest(laczone,duze)
d2=2*d
disp('poniewaz robimy dwa testy powinnismy przynajmniej pomnozyc m i d przez 2')
