pinv: Moore-Penrose pseudo-inverse of a matrix

In rchemo: Dimension Reduction, Regression and Discrimination for Chemometrics

Moore-Penrose pseudo-inverse of a matrix

Description

Calculation of the Moore-Penrose (MP) pseudo-inverse of a matrix X.

Usage

pinv(X, tol = sqrt(.Machine$double.eps))

Arguments

X	X-data (n, p).
tol	A relative tolerance to detect zero singular values.

Value

Xplus	The MP pseudo-inverse.
sv	singular values.

Examples

n <- 7 ; p <- 4
X <- matrix(rnorm(n * p), ncol = p)
y <- rnorm(n)

pinv(X)

tcrossprod(pinv(X)$Xplus, t(y))
lm(y ~ X - 1)

rchemo documentation built on Sept. 11, 2024, 8:05 p.m.