如何找到R中的真实混合物密度？

Question

我不知道如何找到真实混合物密度的图。在下面，我提供了为普通混合物创建函数的代码。而且我了解如何创建每个单独的密度图，但是我不知道如何找到 true 密度。

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
}

Answer 1

我和我的教授写了这篇论文：D.S. Young，X.Chen，D.C。Hewage和R.N. Poyanco（2018）。 “有限的伽马混合准备分布：估计，推断和基于模型的聚类。”

他还有一个R库，称为mixtools：https://cran.r-project.org/web/packages/mixtools/vignettes/mixtools.pdf

例如，可以在第6页上使用wait1 <- normalmixEM(waiting, lambda = .5, mu = c(55, 80), sigma = 5)来估算参数。

此库不需要您提供均值和标准差，它将根据您提供的范围为您进行估算。而且，您可以使用本机R绘图功能来生成漂亮的绘图。

希望对您有帮助。