Objetivo: Construir funções para gerar amostras das distribuições Binomial e Bernoulli.
# Gerando uma Bernoulli
rbernoulli<- function(nc,p)
{
u<-runif(nc)
x<- (u>1-p)*1
}
amostra_bernou <- rbernoulli(30,0.9)
barplot(table(amostra_bernou), col = rainbow(3))# Gerando uma binomial
rbinomial<-function(nc,n,p)
{
X<-c(1:nc)
for (i in 1:nc)
{
X[i]<-sum(rbernoulli(nc,p))
}
X
}
amostra_binom <- rbinomial(30, 10, 0.9)
barplot(table(amostra_binom), col = rainbow(10))