Как найти истинную плотность смеси в 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

Я и мой профессор написали эту статью: Д. С. Янг, Х. Чен, Д. К. Хьюэйдж и Р. Н. Поянко (2018). «Конечные распределения подготовки смеси гаммы: оценка, вывод и кластеризация на основе моделей».

И у него также есть библиотека 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 plot для создания хороших графиков.

Надеюсь, это поможет.