R/qlr.R

Defines functions qlr

Documented in qlr

qlr <-
function(X,thres){
  X.svd = svd(X);
  d = X.svd$d - thres;
  k = sum(d > 1e-8);
  if (k==1){
    D = as.matrix(d[1]);
  }else{
    D = diag(d[1:k]);
  }
  X.lr = X.svd$u[ ,1:k] %*% D %*% t(X.svd$v[ ,1:k]);
  rownames(X.lr) = rownames(X);
  colnames(X.lr) = colnames(X);
  return(list("lr"=X.lr,"rank"=k));
}

Try the QRFCCA package in your browser

Any scripts or data that you put into this service are public.

QRFCCA documentation built on Nov. 17, 2017, 7:20 a.m.