cách mô phỏng phân phối lưỡng cực trong R?

Question

Tôi có mã sau để tạo phân phối hai mặt nhưng khi tôi vẽ đồ thị biểu đồ. Tôi không thấy 2 chế độ. Tôi tự hỏi nếu có gì sai với mã của tôi.

mu1 <- log(1) 
mu2 <- log(10) 
sig1 <- log(3) 
sig2 <- log(3) 
cpct <- 0.4 

bimodalDistFunc <- function (n,cpct, mu1, mu2, sig1, sig2) { 
    y0 <- rlnorm(n,mean=mu1, sd = sig1) 
    y1 <- rlnorm(n,mean=mu2, sd = sig2) 

    flag <- rbinom(n,size=1,prob=cpct) 
    y <- y0*(1 - flag) + y1*flag 
} 

bimodalData <- bimodalDistFunc(n=100,cpct,mu1,mu2, sig1,sig2) 
hist(log(bimodalData)) 

Answer 1

Vấn đề có vẻ là quá nhỏ n và chênh lệch quá nhỏ giữa mu1 và mu2, lấy mu1=log(1), mu2=log(50) và n=10000 cho này: