R 中的 QR 分解 - 强制正对角线

Question

我对 R 中的 qr 函数有一个问题。我的输入矩阵是正定的，所以 R 应该给 r 函数一个对角线均为正的三角矩阵。然而，我发现对角线中有一些负值。我该如何解决这个问题？

假设我们有一个y如下所示的矩阵：

[1,] 0.07018171 -0.07249188 -0.01952050 
[2,] -0.09617788 0.52664014 -0.02930578 
[3,] -0.01962719 -0.09521439 0.81718699 

它是正定的：

> eigen(y)$values
[1] 0.82631283 0.53350907 0.05418694

我在 R 中应用 qr() ，它给了我 Q =

          [,1]      [,2]        [,3] 
[1,] -0.5816076 -0.6157887 0.5315420 
[2,] 0.7970423 -0.5620336 0.2210021 
[3,] 0.1626538 0.5521980 0.8176926

且 R =

[1,] -0.1206685 0.4464293 0.1209139    
[2,] 0.0000000 -0.3039269 0.4797403    
[3,] 0.0000000 0.0000000 0.6513551 

其中对角线不为正。

非常感谢。

这是矩阵：

structure(c(0.07018171, -0.09617788, -0.01962719, -0.07249188, 
0.52664014, -0.09521439, -0.0195205, -0.02930578, 0.81718699), .Dim = c(3L, 
3L))

Answer 1

我可以简单地将对角矩阵与符号（R）相乘，以强制对角线条目为正，然后调整 Q 的相应值。Q 仍然是一个正交矩阵。

示例代码

qr.decom <- qr(A)  
Q <- qr.Q(qr.decom)
R <- qr.R(qr.decom)
sgn <- sign(diag(R))
R.new <- diag(sgn) %*% R
Q.new <- Q %*% diag(sgn)

那么R.new有一个正对角线元素。

我们可以使用问题部分中的示例在 R 中进行尝试。