Advanced Search

Question: Cholesky Decomposition

Question:Cholesky Decomposition

MrMarc 3158 Maple
cholesky
October 21 2011
0

Ok If R is the Cholesky decomposition of the Covariance matrix then I need:

i) Cov=R*R'
ii) R'=R

I can show that i) indeed holds but ii) does not seem to work. Why?

restart:
with(LinearAlgebra):
Cov := Matrix([[.1, .2], [.2, 1.3]]);
Ap := LUDecomposition(Cov, 'method' = 'Cholesky');

Cov = Ap.Transpose(Ap);
Ap = Transpose(Ap)
Please Wait...