R/compbasis.R

Defines functions compbasis

Documented in compbasis

compbasis <-
function(M,tol=sqrt(.Machine$double.eps)) {
  if (qr(M,tol=tol)$rank < dim(M)[1]) {
    return(as.matrix(svd(diag(1,dim(M)[1]) - projector(M,tol=tol))$
                       u[,1:(dim(M)[1] - qr(M,tol=tol)$rank)]))}
  else return(NULL)
}

Try the lmreg package in your browser

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

lmreg documentation built on May 2, 2019, 9:29 a.m.