R/supplbasis.R

supplbasis <-
function(A,B,tol=sqrt(.Machine$double.eps)) {
  if (qr(A,tol=tol)$rank < qr(cbind(A,B),tol=tol)$rank) return(basis(B-proj(A,tol=tol)%*%B))
  else return(NULL)
}
marutter/lmreg documentation built on May 8, 2019, 1:51 p.m.