関数を使用して、R で次の行列の下三角コレスキー分解を取得しようとしていchol()
ます。ただし、上三角分解を返し続け、ドキュメントを調べた後でも、下三角分解を取得する方法を見つけることができないようです。以下は私が使用しているコードです -
x <- matrix(c(4,2,-2, 2,10,2, -2,2,5), ncol = 3, nrow = 3)
Q <- chol(x)
Q
# [,1] [,2] [,3]
# [1,] 2 1 -1.000000
# [2,] 0 3 1.000000
# [3,] 0 0 1.732051
Q
基本的に、 のような行列を見つける必要がありますQQ' = X
。ご協力いただきありがとうございます!