function [TWald,TPomVer,kvantil,p1,p2,zamitnuti1,zamitnuti2]=test_shody1(nj,npj,alfa)
% funkce test_shody provede test dobre shody 
% pro pozorovane a ocekavane cetnosti
% synatxe: 
% [TWald,TPomVer,kvantil,p1,p2,zamitnuti1,zamitnuti2]=test_shody1(nj,npj,alfa)
% vstupni parametry:
% nj ... sloupcovy vektor pozorovanych cetnosti
% npj ... sloupcovy vektor ocekavanych cetnosti
% alfa ... hladina vyznamnosti testu
% vystupni parametry:
% TWald ... testova Waldova statitika
% TPomVer ... testova statistika pomerem verohodnosti
% kvantil ... prislusny kvantil chi-kvadrat rozlozeni
% p1 ... p-hodnota v pripadu Wald. testu
% p2 ... p-hodnota v pripadu testu pomerem verohodnosti
% zamitnuti{1/2} ...=0, kdyz H0 nezamitame
%                   =1, kdyz H0 zamitame{ve Wald.testu/testu pom.verohod.}
% overeni podminek dobre aproximace....hodnota 1 pro poruseni
if sum(npj<5)>0
    poruchy_podminek=(npj<5)'
    error('Nejsou splneny podminky dobre aproximace.')
end;
TWald=sum((nj-npj).^2./npj);
TPomVer=2*sum(nj.*log((nj/sum(nj))./(npj/sum(nj))));
kvantil=chi2inv(1-alfa,size(nj,1)-1);
p1=1-chi2cdf(TWald,size(nj,1)-1);
p2=1-chi2cdf(TPomVer,size(nj,1)-1);
zamitnuti1=(p1<alfa);
zamitnuti2=(p2<alfa);