高斯混合模型概率matlab

Question

我有一个尺寸为50x100000的数据。 （100000个功能，每个功能的维度为50）。

我想使用这些数据拟合高斯混合模型。我使用了以下代码。

               obj = gmdistribution.fit(X',3);

我需要的是当我提供新数据时我应该能够得到概率$ p（Y | \ theta）$，其中$ \ theta $是高斯混合模型参数。

我使用以下代码来获取概率值。

               P = pdf(obj,X');

但我得到的价值非常低，所有人都是0.这真是幸福吗？如何获得适当的概率值？

Answer 1

在一个维度中，高斯分布的pdf的最大值是1 / sqrt（2 * PI）。所以在50个维度中，最大值将是1 /（sqrt（2 * PI）^ 50），大约是1E-20。因此，pdf的值都将达到该数量级或更小。