%###############histogramy##########################

>> ht = 0:10:60;

>> pocet = [6 12 20 13 3 1 0];

%nespravny histogram

>> bar(ht, pocet, 'histc')%;shg

%spravny histogram

>> bar([0 10 20 30 60],[6 12 20 ((13+3+1)/4) 0],'histc')

%###################################################

%#######empiricka distribucna funkcia###############

>> x = [3 1 2 5 6 9];

>> xx = sort(x);

>> n = length(x);

>> yy = ((1/(2*n)) : (1/n) : (1 - (1/(2*n))));

>> plot(xx,yy),axis(-1 10 -1 1)

>> shg

% - 25 percentil

>> z = interp1(yy,xx,0.25)


%###########Korelacny koeficient####################

>> jano = [10 20 30 50 15 20];

>> eva = [30 20 15 40 22 56];

>> kor = corrcoef(jano,eva);

>> korkoef = kor(1,2) 

%###################################################
 
% ##########regresna priamka########################
 
>> clear
>> x = [179 172 183 165 193 186 169 181 176 177];

>> y = [180 175 183 170 189 187 174 183 180 179];

>> covmat = cov(x,y);

>> ssxy = covmat(1,2);

>> b0 = ssxy/var(x);
>> a0 = mean(y) - b0*mean(x);
>> plot([169 186], [(b0*169 + a0) (b0*186 + a0)])

%####################################################


%###########vykreslenie distribuc funkcie############

>> xi = [-5 0 6];
>> pxi = [0.2 0.3 0.5];
>> plot([-20 -5],[0 0]);hold on;plot([-5 0],[0.2 0.2]);hold on;plot([0 6],[0.5 0.5]);plot([6 10],[1 1]),axis([-7 8 -1 1.5]);

%#####################################################

%###########normalne rozdelenie#######################
%
% nech X~N(3,4) Urc P(0<X<5) 
%
% P(((0-3)/2)<((X-3)/2)<((5-3)/2)) = P((-3/2)<X0<1) = Fn(1) - Fn(-1.5) = Fn(1) + Fn(1.5) - 1
%
>> normcdf(5,3,2) - normcdf(0,3,2) 

%######################################################

%############exponencialne rozdelenie##################

>>expcdf(predpokladany cas, parameter mi)

%############erlangovo rozdel###########################
% 13 zakazka do 50 min

% Xi~exp(lambda = 0.2)

% X~erl(12,0.2)

>>gamcdf(50,12,0.5) % 1/0.2=0.5

%########################################################