pinv: Pseudo-inverse of a matrix

In abarbour/pborepo: Access to data-products and metadata from the Plate Boundary Observatory geodetic network

Pseudo-inverse of a matrix

Description

Uses singular value decomposition to calculate the pseudo-inverse of a matrix.

Usage

pinv(Sij, ...)

Arguments

Sij	matrix, or an object to be coerced into a matrix
...	additional arguments

Details

This can be useful when predictions of tensor strains are used to make predictions on a gauge-by-gauge basis with strain data

The attribute "dimnames" is recycled; any other attribute will be lost.

Author(s)

A. Barbour

See Also

pseudoinverse for details of the method

calibration_matrix

Examples

# get a calibration table
pbotbl <- get_caltbl("pbo")

# and the first available calibration matrix...
m <- any_calibration(pbotbl, "B082")

# Calculate the pseudoinverse
print(mi <- pinv(m))

# we can get back the original, less attributes of course
mo <- pinv(mi)
try( all.equal(m, mo, check.attributes=FALSE) )

abarbour/pborepo documentation built on April 27, 2022, 8:47 a.m.