# rm(list=ls())
options(OutDec = ",") 
#==============================================================================
# Exemplo para gerar da Bernoulli a partir de U(0,1)
#==============================================================================
set.seed(54345) # Semente
p  <- 0.9       # Probabilidade de X=1
n  <- 1000      # Tamanho da amostra

#==============================================================================
# Primeira maneira 
#==============================================================================
x  <- rep(0,n)
for(i in 1:n){
	u <- runif(1)
	if(u <= 0.9) x[i] <- 1
}
print(mean(x))

[1] 0,903

#==============================================================================
# Histograma
#==============================================================================
par(mfrow=c(1,1),lwd=2,cex.lab=1.5,cex.axis=1.5,lab=c(2,6,0),
    mar=c(4.5,5,1,1),bty="n")
barplot(table(x)/n,ylim=c(0,1),col="darkgreen",xlab=expression(x),
       ylab="Frequência relativa")

#==============================================================================
# Segunda maneira
#==============================================================================
# Tudo em bloco!
w  <- as.integer( runif(n) <= p )
print(mean(w))

[1] 0,894

#==============================================================================
# Histograma
#==============================================================================
par(mfrow=c(1,1),lwd=2,cex.lab=1.5,cex.axis=1.5,lab=c(2,6,0),
    mar=c(4.5,5,1,1),bty="n")
barplot(table(w)/n,ylim=c(0,1),col="darkblue",xlab=expression(x),
       ylab="Frequência relativa")

# rm(list=ls()) 
#==============================================================================
# graphics.off()
#==============================================================================
# Fim
#==============================================================================