Computes a stable pseudoinverse of a symmetric positive semi-definite matrix using singular value decomposition (SVD) and ridge regularization. This is useful when the matrix is ill-conditioned or rank-deficient.
Arguments
- m
A symmetric numeric matrix (e.g., from
crossprod()).- alpha
Non-negative scalar specifying the ridge penalty. A small positive value stabilizes the inversion by shrinking large singular values.