我一直试图将高斯混合模型(GMM)拟合为单变量分布,其中每个观测值都有一个指定的权重(从0到1)。
可以轻松地为基本密度'分配权重。功能(来自' stats'),但我无法对GMM包执行此操作,因为' mixtools'或者' mclust'。
后者甚至具有这样的功能 - ' me.weighted' - 在链接" http://cran.open-source-solution.org/web/packages/mclust/vignettes/mclust-weights.pdf"中举例说明,但我无法使其发挥作用。 在他们的示例和我的数据中都出现错误:
> fitnew <- do.call("me.weighted", c(list(data = X, weights = w), fit))
Error in me.weighted(data = list(Sepal.Length = c(5.1, 4.9, 4.7, 4.6, : formal argument "data" matched by multiple actual arguments