powered by
Moore-Penrose generalized inverse
geninv(X, tol = .Machine$double.eps^(2/3))
matrix
tolerance for determining bad entries
A matrix of the same dimension as X is returned, the Moore-Penrose generalized inverse.
X
# NOT RUN { # create a positive definite 5x5 matrix x <- crossprod(matrix(rnorm(25),5)) # make it singular x[,2] <- x[,3]+x[,5] geninv(x) # }
Run the code above in your browser using DataLab