ginv.1: Generalized inverse - test version

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

ginv.1(x, tol = .Machine$double.eps * d[1] * max(nrow(x), ncol(x)))

Arguments

x
tol

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, tol = .Machine$double.eps * d[1] * max(nrow(x), 
    ncol(x))) 
{
    r = svd(x)
    d = r$d
    dinv = ifelse(d < tol, 0, 1/d)
    ret = r$v %*% (dinv * t(r$u))
    dimnames(ret) = dimnames(x)[2:1]
    ret
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.