我不明白如何找到真实混合密度的图。在下文中，我提供了为正常混合创建函数的代码。而且我了解如何创建每个单独的密度图，但我不知道如何找到真实的密度。
NormalMix <- function(n,omega,mu1,mu2,sigma1,sigma2) {
z <- sample(2,n,replace=T,prob=c(omega,1-omega))
n1 <- sum(z==1)
n2 <- sum(z==2)
z[z==1] <- rnorm(n1,mu1,sigma1)
z[z==2] <- rnorm(n2,mu2,sigma2)
z
}
我和我的教授写了这篇论文：DS Young、X. Chen、DC Hewage 和 RN Poyanco (2018)。“Gamma 准备分布的有限混合：估计、推理和基于模型的聚类。”</p>
而且他还有一个 R 库，叫做 mixtools：https ://cran.r-project.org/web/packages/mixtools/vignettes/mixtools.pdf
例如，您可以使用wait1 <- normalmixEM(waiting, lambda = .5, mu = c(55, 80), sigma = 5)第 6 页上的示例来估算参数。
这个库不需要你提供平均值和标准差，它会根据你提供的范围为你估计。您可以使用原生 R 绘图功能生成漂亮的绘图。
我希望它有所帮助。