/****************************************************
  ECO375 Class1 Example: Distribution  
*****************************************************/

/* Set Number of observations */
clear all

set obs 100000

/* Binomial Distribution (number of tries, prob(x=1) for each try ) */
gen x = rbinomial(1, 0.5)
gen x2 = rbinomial(2, 0.5)
gen x5 = rbinomial(5, 0.5)
gen x1k = rbinomial(1000, 0.5)
label var x "Once coin"
label var x2 "two coins"
label var x "a coin"
label var x5 "five coins"
label var x1k "1000 coins"
hist x , discrete name(fig1, replace)
hist x2 , discrete name(fig2, replace)
hist x5 , discrete name(fig3, replace)
hist x1k , discrete name(fig4, replace)
graph combine fig1 fig2 fig3 fig4

/* Poisson Distribution (Lambda) Lambda = mean = variance */
gen p1 = rpoisson(1)
gen p2 = rpoisson(2)
gen p5 = rpoisson(5)
gen p100 = rpoisson(100)
label var p1 "poisson with lambda = 1"
label var p2 "poisson with lambda = 2"
label var p5 "poisson with lambda = 5"
label var p100 "poisson with lambda = 100"
hist p1 , discrete name(figp1, replace)
hist p2 , discrete name(figp2, replace)
hist p5 , discrete name(figp5, replace)
hist p100 , discrete name(figp100, replace)
graph combine figp1 figp2 figp5 figp100
global white_back graphregion(color(white))

/* Normal Distribution (Mean, Standard Deviation)  */
gen n1 = rnormal(10, 5)
gen n2 = rnormal(0,1)
label var n1 "Normal(10,5)"
label var n2 "Normal(0,1)"
hist n2, $white_back name(fign2) normal
hist n1, $white_back name(fign1, replace) normal
graph combine fign1 fign2, col(1)


 /* End of Code */