R Rrank 求上三角矩陣的秩

說明

通過 rank 塊估計上 rank 的條件數，並減少 rank 直到該值達到可接受的低值，從而找到上三角矩陣 R 的秩。假設 R 是通過使用旋轉的方法計算的，通常是旋轉 QR 或 Choleski。

用法

Rrank(R,tol=.Machine$double.eps^.9)

參數

R	上三角矩陣，通過旋轉 QR 或旋轉 Choleski 獲得。
tol	用於判斷排名的容差。

細節

該方法基於 Cline 等人。 (1979) 如 Golub 和 van Loan (1996) 中所述。

例子

  set.seed(0)
  n <- 10;p <- 5
  x <- runif(n*(p-1))
  X <- matrix(c(x,x[1:n]),n,p)
  qrx <- qr(X,LAPACK=TRUE)
  Rrank(qr.R(qrx))

作者

Simon N. Wood simon.wood@r-project.org

參考

Cline, A.K., C.B. Moler, G.W. Stewart and J.H. Wilkinson (1979) An estimate for the condition number of a matrix. SIAM J. Num. Anal. 16, 368-375

Golub, G.H, and C.F. van Loan (1996) Matrix Computations 3rd ed. Johns Hopkins University Press, Baltimore.