In [1]:
# rm(list=ls())
options(OutDec = ",") 
#==============================================================================
# Exemplo para gerar de uma mistura discreta de duas normais 
#==============================================================================
set.seed(54345)  # Semente
n      <- 100000 # Tamanho da amostra
omega  <- 0.3
mu1    <- 10
sigma1 <- 1
mu2    <- 15
sigma2 <- 1.5
y      <- rbinom(n,1,omega)
mu     <- y*mu1    + (1-y)*mu2
sigma  <- y*sigma1 + (1-y)*sigma2
x      <- rnorm(n,mu,sigma)
In [2]:
#==============================================================================
# Histograma
#==============================================================================
par(mfrow=c(1,1),lwd=2,cex.lab=1.5,cex.axis=1.5,lab=c(14,6,0),
    mar=c(4.5,5,2,1),bty="n")
hist(x,nclass=50,prob=T,main="",ylim=c(0,0.20),xlim=c(5,23),col="darkgreen",
     ylab=expression(f(x)),xlab=expression(x))
xseq <- seq(5,23,length=1000)
yseq <- ( omega*dnorm(xseq,mu1,sigma1)
        + (1-omega)*dnorm(xseq,mu2,sigma2) )
lines(xseq,yseq,col="red",lwd=3)
No description has been provided for this image
In [3]:
# rm(list=ls()) 
#==============================================================================
# graphics.off()
#==============================================================================
# Fim
#==============================================================================