我一直在尝试为单变量分布拟合高斯混合模型(GMM),其中每个观察值都有一个分配的权重(从 0 到 1)。
为基本的“密度”功能(来自“统计”)分配权重可以很容易地完成,但我无法为 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