零膨胀高斯混合模型

问题描述

mclust 或其他 R 包是否有任何扩展，可以对零膨胀的连续数据进行建模？ Mclust 无法处理膨胀的 0，请参阅下面的代码。

library(tidyverse)
library(mclust)

n <- 250
y0 <- rep(0,n)
y1 <- rnorm(n/2,4,0.5)
y2 <- rnorm(n/2,5,2)
y3 <- rnorm(n/2,12,1)
y4 <- rnorm(n/2,10,2)
y <- c(y0,y1,y2,y3,y4)
hist(y)

#including zeros
mc <- densityMclust(data=y)
summary(mc)
plot(mc,what = "BIC")

#excluding zeros
y2 <- c(y1,y4)
mc <- densityMclust(data=y2)
summary(mc)
plot(mc,what = "BIC")

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）