abc Matlab - elektronski priručnik
III dio Matematika u Matlabu - 3.4. Vjerovatnoća i statistika

3.4.5. Vrste raspodjela
3.4.5.1. Binomna raspodjela
3.4.5.1.1. Binomna raspodjela - Primjer 1

Pretpostavimo da postoji 20% šansi da odrasli Amerikanac pati od psihijatrijskog poremećaja. Uzima se slučajni uzorak od 25 odraslih Amerikanaca. Ako X predstavlja broj ljudi koji imaju psihijatrijski poremećaj, tada je X binomna slučajna varijabla sa parametrima (25, 0.20). Interesuje nas vjerovatnoća da barem 3 od odabranih ljudi imaju takav poremećaj. Možemo upotrijebiti MATLAB Statistics Toolbox funkciju binocdf da odredimo , na sljedeći način:

                prob = binocdf(3,25,0.2);

Možemo takođe sabrati pojedinačne vrijednosti funkcije mase vjerovatnoće od X=0 do X=3:

                prob2 = sum(binopdf(0:3,25,0.2));
Obe ove komande vraćaju vjerovatnoću od 0.234. Sada ćemo pokazati kako se generiše binomna raspodjela prikazana na Slici 1.

% Uzmi vrijednosti za domen, x.
x = 0:6;
% Uzmi vrijednosti funkcije mase vjerovatnoće.
% Prvo za n = 6, p = 0.3:
pdf1 = binopdf(x,6,0.3);
% Sada za n = 6, p = 0.7:
pdf2 = binopdf(x,6,0.7);

Sada imamo vrijednosti za funkciju mase vjerovatnoće (ili visinu kolona (engl. bar)). Crteži se dobijaju pomoću sljedećeg koda.

                % Crtanje.
subplot(1,2,1),bar(x,pdf1,1,'w')
title(' n = 6, p = 0.3')
xlabel('X'),ylabel('f(X)')
axis square
subplot(1,2,2),bar(x,pdf2,1,'w')
title(' n = 6, p = 0.7')
xlabel('X'),ylabel('f(X)')
axis square

Vrste raspodjela - Binomna raspodjela    <    Index    >    Vrste raspodjela - Poasonova raspodjela