我正在尝试将混合模型拟合到时间序列中以进行预测（平均一个arima、一个ets和一个stlf模型）。

有没有可以处理这个问题的 R 包？

我设法找到了许多可以处理混合模型的包，但它们不包括预测。

如果没有这样的包，是否有一种计算预测区间的好方法（预测本身只是平均值）。

谢谢！

我想使用 GMM（高斯混合模型对二值图像进行聚类，并且还想在二值图像本身上绘制聚类质心。

我以此为参考：http: //in.mathworks.com/help/stats/gaussian-mixture-models.html

这是我的初始代码

但问题是我无法在主二进制图像中绘制从 GMM 接收到的簇质心。我该怎么做？

**现在假设我有 10 个这样的图像序列并且我想将它们的平均位置信息存储在两个单元格数组中，那么我该怎么做。这是我的新问题的代码 **

我能够分割 10 个图像，但没有存储在单元格数组 xc 和 xb 中的均值值。它们仅存储 [] 代替均值值

我有一个大型数据集（我无法将整个数据放入内存）。我想在这个数据集上拟合 GMM。

我可以在小批量数据上重复使用GMM.fit()（ ）吗？sklearn.mixture.GMM

我试图在 R 中拟合高斯分布和均匀分布的混合物，但成功有限。

我依赖于flexmix提供一个 EM 框架，并使用 package.json 中匹配矩估计器的略微修改版本来拟合fitdistrplus组件。（适应只是用加权版本替换样本均值和方差。）

因为我flexmix对代码（如下）不是很熟悉，所以这是一个 hack。

看起来均匀分布的权重太小了。

非常感谢任何有关如何改善合身性的帮助。

首先是例子

生成玩具数据

理想的合身

合身

功能

使用下面的代码，我将数据集分成两个类，然后ggplot是点，将每个类标记为 1 类或 2 类。

问题：如何改用渐变着色来显示data_f.k2$posterior属于第 2 类的每个点的后验概率 ？在问题的底部，我分享了我在 using 方面的尝试scale_colour_gradient，但这是行不通的。

这将产生下图：

这是我尝试获得后渐变色的尝试，但它不起作用。我希望颜色褪色，比如说从红色到蓝色 - 中间的点（可能属于两种混合物）为紫色。

我想加载存储在 html 文件的各个子目录中的图像，例如在我的 index.html 文件所在的 images/ 文件夹中。

当我尝试显示带有此标签的图像时

如果我将项目创建为 html 但不适用于本地视图，则此方法有效。

我犯了什么愚蠢的错误吗？

我试图估计特威迪（或复合泊松伽马）分布的有限混合。我搜索了我能想到的任何资源，但没有找到任何关于如何做到这一点的资源。

我目前正在尝试使用 R 中的flexmix包编写不同的 M-step 驱动程序，如第 12-14 页的flexmix小插图中所述。这是我的代码，它依赖于cplm包：

但是，这会导致以下错误：

dtweedie(y, xi = p, mu = mu, phi = phi) 中的错误：对不一致数组的二元运算

有没有人做过或见过 tweedie 分布的有限混合？您能指出我使用flexmix还是其他方式来完成此任务的正确方向吗？

当我从中删除一些字母时，尝试使用此文件名和混合物给出页面未找到错误，然后它显示页面。我认为这不是文件名的限制。至少可以使用 255 个字母。这是混合物的问题吗？但我的服务器也没有显示文件

我试图将高斯混合模型的轮廓线可视化到我的散点图，但它们没有正确显示，首先我认为它们被点重叠，但即使我先画线然后绘制散点，问题仍然存在。

我生成图的代码：

这就是它的外观：注意轮廓线的切割

我目前使用混合模型包EMMIXskew，并且在我的数据（一些数值向量）上安装了一个 skew-t 分布。这个包有一些密度函数ddmst，但我没有在这个包中看到概率密度函数，我需要一些！

我以为我能做的是

1. 使用其他一些为 skew-t 分发提供 pdf 的包，例如 package snwith pst，但问题是这个分发对 skew-t 分发有不同的定义，或者
2. 我可以使用integrateon ddmst，但到目前为止它不起作用。

我尝试了类似的东西

并且如果我假设我的数据具有参数的 3 模态

我收到这个错误

Error in ddmix(dat, n, p, 1, "mst", mean, cov, nu, del) : dat does not match n and p.

我真的不知道我做错了什么！