This version continues optimizing by calling faster bcoef function rather than lm().