I believe another workaround to this is to use glm() with the option family="gaussian" instead of lm()