如何在 R 中找到真正的混合密度？答案

【问题标题】：How do I find a true mixture density in R?如何在 R 中找到真正的混合密度？
【发布时间】：2020-02-24 00:21:06
【问题描述】：

我不明白如何找到真实混合密度的图。在下文中，我提供了为正常混合创建函数的代码。而且我了解如何创建每个单独的密度图，但我不知道如何找到 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
}

【问题讨论】：

这是一个统计问题吗？就像你试图找到组合 pdf 的功能形式？
也许这会有所帮助：stackoverflow.com/questions/23480529/… 或此stackoverflow.com/questions/39945473/… 或此stackoverflow.com/questions/20452650/…

标签： r density-plot mixture

【解决方案1】：

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

而且他还有一个 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 绘图功能来生成漂亮的绘图。

希望对你有帮助。

【讨论】：