powered by
Inverts a symmetric positive-definite matrix without requiring the Matrix package.
sinv(obj)
A matrix of the same size as the input object
The symmetric positive-definite matrix
Claus Ekstrom, claus@rprimer.dk.
This function does no error checking and it is up to the user to ensure that the input is indeed symmetric, positive-definite, and a matrix.
m <- matrix(c(1, 0, .5, .5, 0, 1, .5, .5, .5, .5, 1, .5, .5, .5, .5, 1), 4) sinv(m)
Run the code above in your browser using DataLab