gsi.svdsolve {compositions} | R Documentation |
Internal function: Solves singular and non square equation systems
Description
Based on the singular value decomposition, a singular equation system ax=b is solved.
Usage
gsi.svdsolve(a,b,...,cond=1E-10)
Arguments
a |
the matrix of ax=b (a.k.a. left-hand side matrix) |
b |
the vector or matrix b of ax=b (a.k.a right-hand side, independent element) |
cond |
the smallest-acceptable condition of the matrix. Smaller singular values are truncate |
... |
additional arguments to svd |
Value
The "smallest" vector or matrix solving this system with minimal joint error among all vectors.
Note
Do not use gsi.* functions directly since they are internal functions of the package
Author(s)
K.Gerald v.d. Boogaart http://www.stat.boogaart.de
Examples
#A <- matrix(c(0,1,0,0,0,0),ncol=2)
#b <- diag(3)
#erg <- gsi.svdsolve(A,b)
#erg
#A %*% erg
#diag(c(0,1,0)) # richtig
[Package compositions version 2.0-8 Index]