[R] about the Choleski factorization
93354504
93354504 at nccu.edu.tw
Fri Mar 27 16:46:19 CET 2009
Hi there,
Given a positive definite symmetric matrix, I can use chol(x) to obtain U where U is upper triangular
and x=U'U. For example,
x=matrix(c(5,1,2,1,3,1,2,1,4),3,3)
U=chol(x)
U
# [,1] [,2] [,3]
#[1,] 2.236068 0.4472136 0.8944272
#[2,] 0.000000 1.6733201 0.3585686
#[3,] 0.000000 0.0000000 1.7525492
t(U)%*%U # this is exactly x
Does anyone know how to obtain L such that L is lower triangular and x=L'L? Thank you.
Alex
More information about the R-help
mailing list