*DISTRIBUIÇÃO DE POISSON

*Verdadeira probabilidade
di 50 / 30000

*Probabilidade utilizada: p=0.002

*Número de estudantes que você passa em um dia típico (n)
*Caminhada de 30 minutos
*5 estudantes por minuto
di 5 * 30

*Média = mu = np
di 150 * 0.002

*P(x) = (mu^x * e^-mu) / x!

*P(0)
di ((0.3^0)*(2.71828^(-0.3)))/1
di ((0.3^0)*(exp(-0.3)))/1

*P(1)
di ((0.3^1)*(2.71828^(-0.3)))/1
di ((0.3^1)*(exp(-0.3)))/1

*P(2)
di ((0.3^2)*(2.71828^(-0.3)))/(2*1)
di ((0.3^2)*(exp(-0.3)))/(2*1)

*P(3)
di ((0.3^3)*(2.71828^(-0.3)))/(3*2*1)
di ((0.3^3)*(exp(-0.3)))/(3*2*1)

*P(4)
di ((0.3^4)*(2.71828^(-0.3)))/(4*3*2*1)
di ((0.3^4)*(exp(-0.3)))/(4*3*2*1)

*P(5)
di ((0.3^5)*(2.71828^(-0.3)))/(5*4*3*2*1)
di ((0.3^5)*(exp(-0.3)))/(5*4*3*2*1)

*P(>5) = 1 - [P(1)+P(2)+P(3)+P(4)+P(5)]
di 1- (0.74081837 + 0.22224551 + 0.03333683 + 0.00333368 + 0.00025003 + 0.000015)

*Criando banco de dados
clear
input x p
0 0.74081837
1 0.22224551
2 0.03333683
3 0.00333368
4 0.00025003
5 0.000015
6 0.00000058
end

*Verificando se somatório é igual a uma unidade
egen sump = total(p)
tab sump

*Distribuição de probabilidade em tabela
list x p

*Distribuição de probabilidade em histograma
graph bar p, over(x, gap(0)) bar(1, fcolor(gray) lcolor(black))