geninv: Moore-Penrose generalized inverse

Description Usage Arguments Value Examples

View source: R/utils.R

Description

Moore-Penrose generalized inverse

Usage

1
geninv(X, tol = .Machine$double.eps^(2/3))

Arguments

X

matrix

tol

tolerance for determining bad entries

Value

A matrix of the same dimension as X is returned, the Moore-Penrose generalized inverse.

Examples

1
2
3
4
5
# create a positive definite 5x5 matrix
x <- crossprod(matrix(rnorm(25),5))
# make it singular
x[,2] <- x[,3]+x[,5]
geninv(x)

sgaure/durmod documentation built on Dec. 19, 2019, 4:57 a.m.