This function calculates matrix inverses for symmetric matrices.
Usage
Rmatsolve(m)
Arguments
m
a (symmetric) matrix.
Value
inv
the matrix inverse of m.
Details
This function uses the eigenvalue decomposition of a matrix m to work out its inverse. The function is used here since standard matrix inverse algorithms do not cope well with matrices which are near singular (this often happens in the regression stages of the forward transforms.