R/gmi.R

In RKUM: Robust Kernel Unsupervised Methods

Defines functions gmi

Documented in gmi

gmi <-
function(X, tol = sqrt(.Machine$double.eps))
 {
  dnx <- dimnames(X)
  if(is.null(dnx)) dnx <- vector("list", 2)
  s <- svd(X)
  nz <- s$d > tol * s$d[1]
  structure(
    if(any(nz)) s$v[, nz] %*% (t(s$u[, nz])/s$d[nz]) else X,
    dimnames = dnx[2:1])
 }

Try the RKUM package in your browser