matlab - 高斯混合模型概率matlab

Question

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

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

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

我需要的是当我给出一个新数据时，YI 应该能够得到似然概率 $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的值都将是那个数量级，或者更小。