## You can do this with the qr functions. tmp <- rnorm(4)
tmp.qr <- qr(tmp) tmp.complete <- qr.Q(tmp.qr, complete=TRUE) tmp.complete[,1] / tmp ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.