13-Obteniendo una muestra para distribuciones discretas- caso Poisson

Muestra aleatoria Poisson

David J. Santana

17 de septiembre de 2016

Generando una muestra Poisson y después la comparamos con una muestra generada con rpois.

set.seed(123)
a=dpois(0:50,25)
b=100000*a
c=round(b)
d=25
for(i in 1:51){
  if(c[i]>0){d=c(d,rep(i-1,c[i]))}
}
ee=sample(d,10000,replace = T)
aa=rpois(10000,25)
brbr=seq(-0.5,50.5,1)
par(bg='lightblue')
hist(aa,breaks = brbr,freq=F, col='yellow',
     main = 'Poisson Random Sample',xlab = 'X')
hist(ee,breaks = brbr,freq = F,add = T,
     col="green")
lines(0:50,dpois(0:50,25),col=2,lw=2)
text(40,0.08,'La curva roja',col=4)
text(40,0.074,' une los puntos de la',col=4)
text(40,0.068,'densidad verdadera.',col=4)