An Inverse Solution

Share:

Description

QR decomposition solution to Ax=b

Usage

1
Ainv(GAB, x, tol = 1e-12)

Arguments

GAB

design matrix

x

right hand side

tol

tolerance for singularity

Details

I needed something to make up for the lame-o matlab code that does this h = G\x to get the inverse

Value

Inverse Solution

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

Examples

1
2
3
4
5
6
7
8
9
set.seed(2015)
GAB = matrix(runif(36), ncol=6)
truex =rnorm(ncol(GAB))
rhs = GAB %*% truex

rhs = as.vector(rhs )

tout = Ainv(GAB, rhs, tol = 1e-12)
tout - truex 

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.