testing_speed.R

microbenchmark::microbenchmark({lm(ys~xs, data = testdf)$residuals},
                               {fastLm(ys~xs, data = testdf)$residuals},
                               {homemade_residuals(indeps, deps)}, times = 500)

identical(homemade_residuals(indeps, deps), homemade_residuals2(indeps, deps))
# all.equal()

testdf <- data_frame(xs = runif(200, 0, 15),
                     ys = xs * 5 + 15 + rnorm(200, sd = 3))
profvis::profvis({homemade_residuals(testdf)})

	pro.residuals1.1<-function(Y,X){
		X<-sweep(X,2,colMeans(X,na.rm=TRUE),check.margin = FALSE)
		Y<-sweep(Y,2,colMeans(Y,na.rm=TRUE),check.margin = FALSE)
		Yfit<-X%*%solve(t(X)%*%X)%*%(t(X))%*%Y
		Yres<-Y-Yfit
	return(Yres)
	}

http://happygitwithr.com/
http://r-pkgs.had.co.nz/r.html

https://www.dropbox.com/s/bd2le8lc2d6uvb5/traits_with_NA_multiv.txt?dl=0