title: How to simulate bimodal distribution?
body: <p>I have the following code to generate bimodal distribution but when I graph the histogram. I don't see the 2 modes.  I am wondering if there's something wrong with my code.</p>
  
  <pre><code>mu1 &lt;- log(1)   
  mu2 &lt;- log(10)
  sig1 &lt;- log(3)
  sig2 &lt;- log(3)
  cpct &lt;- 0.4   
  
  bimodalDistFunc &lt;- function (n,cpct, mu1, mu2, sig1, sig2) {
    y0 &lt;- rlnorm(n,mean=mu1, sd = sig1)
    y1 &lt;- rlnorm(n,mean=mu2, sd = sig2)
  
    flag &lt;- rbinom(n,size=1,prob=cpct)
    y &lt;- y0*(1 - flag) + y1*flag 
  }
  
  bimodalData &lt;- bimodalDistFunc(n=100,cpct,mu1,mu2, sig1,sig2)
  hist(log(bimodalData))
  </code></pre>