r - Solving a system of unknowns in terms of unknowns -


i trying solve 5x5 cholesky decomposition (for variance-covariance matrix) in terms of unknowns (no constants).

a simplified version, sake of giving example, 2x2 decomposition:

[[a,0],[b,c]]*[[a,b],[0,c]]=[[u1,u2],[u2,u3]] 

is there software (i'm proficient in r, if r can great) solve above yield answer of left-hand variables in terms of right-hand variables? i.e. final answer:

a = sqrt(u1) b = u2/sqrt(u1) c = sqrt(u3+u2/u1) 

take @ this wikipedia section.

the symbolic definition of (i,j)th entry of decomposition defined recursively in terms of entries above , left. implement these recursions using matlab's symbolic math toolbox , apply them (symbolically) obtain formulas 5x5 case. warned you'll end extremely complicated formulas of unknowns, , - excepting unusual circumstances - fine implement decomposition iteratively fixed size 5x5 matrix.


Comments

Popular posts from this blog

how to proxy from https to http with lighttpd -

android - Automated my builds -

python - Flask migration error -