我已经实现了exp。mu
最大化算法,它会收敛并正确返回 和 的值sigma
,我已经检查了各种示例。
我试图绘制对数似然,但我不知道它会以正确的形式出现吗?这是等式: 我的情节,y是对数似然值,x是迭代次数。
负值很奇怪,也许我应该标准化可能性?对数似然在 Exp 中意味着什么。最大化?
logLikelihood = 0;
for i = 1 : n
logTemp = 0;
for j = 1 : k
logTemp = logTemp + p(j) * mvnpdf(x(i,:), mu(j,:), sigma(:,:,j));
end
logLikelihood = logLikelihood + log(logTemp);
end
plot(iteration, logLikelihood,'r*');
hold on;