Invert of covariance matrices
Tries to invert a matrix by solve. If this fails because of singularity, an eigenvector decomposition is computed, and eigenvalues below 1/cmax are replaced by 1/cmax, i.e., cmax will be the corresponding eigenvalue of the inverted matrix.
solvecov(m, cmax = 1e+10)
m |
a numeric symmetric matrix. |
cmax |
a positive value, see above. |
A list with the following components: inv the inverted matrix, coll TRUE if solve failed because of singularity.
solvecov code was taken from package fpc: Christian Hennig
Christian Hennig
solve, eigen
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.