У меня есть образец 1м записей, полученных из моих исходных данных. (Для справки вы можете использовать эти фиктивные данные, которые могут генерировать примерно аналогичное распределение

b <- data.frame(matrix(rnorm(2000000, mean=c(8,17), sd=2)))
c <- b[sample(nrow(b), 1000000), ]

Я полагал, что гистограмма представляет собой смесь двух логнормальных распределений, и я попытался согласовать суммированные распределения, используя алгоритм EM, используя следующий код:

install.packages("mixtools")
lib(mixtools)
#line below returns EM output of type mixEM[] for mixture of normal distributions
c1 <- normalmixEM(c, lambda=NULL, mu=NULL, sigma=NULL) 
plot(c1, density=TRUE)

Первый график представляет собой график логарифмического правдоподобия, а второй (если вы нажмете клавишу возврата еще раз) дает аналогичные следующие кривые плотности:

Как я уже упоминал, c1 имеет тип mixEM [], и функция plot () может это учитывать. Я хочу заполнить кривые плотности цветами. Это легко сделать с помощью ggplot2 (), но ggplot2 () не поддерживает данные типа mixEM [] и выдает это сообщение:

«ggplot не знает, как обращаться с данными класса mixEM». Могу ли я использовать какой-либо другой подход для решения этой проблемы? Любые предложения с благодарностью!

Спасибо!