R/normX.R

normX <- function (x,w)
	{
          #qq<-La.svd(sqrt(w)*x); list(q=(1/sqrt(w))*(qq$u),r=qq$d)
          qq <- qr((1/sqrt(w)) * x)
	  list(q = (1/sqrt(w)) * qr.Q(qq), r=abs(diag(qr.R(qq))))
	}

Try the homals package in your browser

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

homals documentation built on May 2, 2019, 6:51 p.m.