如何在 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 绘图函数生成漂亮的绘图。

希望对您有所帮助。

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

How do I find a true mixture density in R?

r

density-plot

mixture