似然分解

机器算法验证 最大似然 协方差 可能性 高斯过程
2022-03-22 16:08:30

优化全似然函数有时很耗时,并且包含许多数值问题和不稳定性,尤其是在需要矩阵求逆时。如果我们有 3 个输入向量基于此,我想建模其中 K 是所有三个输入向量之间的联合协方差矩阵。输入连接如下图所示:是独立的,而取决于y1,y2,y3y=[y1,y2,y3]y

yN(0,K)
y1y2y1y2y3换句话说,协方差矩阵表示为 其中是任何半正定协方差函数。我有什么特定的方法可以利用这种独立性并分解完整的似然函数,例如我可以写并分别优化每个部分?的表达式是什么K
K=(k110k130k22k23k13k23k33)
k
f(y1,y2,y3)=f(y3|y1,y2)f(y1)f(y2)
f(y3|y1,y2)给定上面的协方差矩阵。还有一种特定的方法可以使用复合似然方法,例如成对似然法。

非常感谢有关此类正常似然分解的任何良好参考

在此处输入图像描述

1个回答

我不知道这最终是否是一件好事,但是您可以使用https://en中显示的标准 Schur 补码方法 条件,这是一维法线。 wikipedia.org/wiki/Multivariate_normal_distribution#Conditional_distributionsy3y1,y2

的均值。μ1,μ2,μ3y1,y2,y3

将K_表示为具有条目的对角矩阵,即的协方差矩阵。它的逆是通过反转对角元素来获得的。K12k11k22[y1,y2]T

y3为条件的 y_3是正常的并且具有 [y1,y2]T=[x1,x2]T

mean=μ3+[k13,k23]K121[x1μ1,x2μ2]T
variance=k33[k13,k23]K121[k13,k23]T