porting R's lm over