function [K,kvantil,p,zamitnuti]=test_shody(nj,npj,alfa)
% funkce test_shody provede test dobre shody 
% pro pozorovane a ocekavane cetnosti
% synatxe: [K,kvantil,p,zamitnuti]=test_shody(nj,npj,alfa)
% vstupni parametry:
% nj ... sloupcovy vektor pozorovanych cetnosti
% npj ... sloupcovy vektor ocekavanych cetnosti
% alfa ... hladina vyznamnosti testu
% vystupni parametry:
% K ... testova statitika
% kvantil ... prislusny kvantil chi-kvadrat rozlozeni
% p ... p-hodnota
% zamitnuti ... =0, kdyz H0 nezamitame
%               =1, kdyz H0 zamitame
% overeni podminek dobre aproximace....hodnota 1 pro poruseni
if sum(npj<5)>0
    poruchy_podminek=(npj<5)'
    error('Nejsou splneny podminky dobre aproximace.')
end;
K=sum((nj-npj).^2./npj);
kvantil=chi2inv(1-alfa,size(nj,1)-1);
p=1-chi2cdf(K,size(nj,1)-1);
zamitnuti=(p<alfa);