> lm(y ~ x-1) > solve(crossprod(x), t(x))%*%y # probably this can be done more > efficiently
You could do crossprod(x,y) instead of t(x))%*%y ______________________________________________ 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.
