В биостатере по касательной затронули тему вычмата и вопрос того, что неплохо бы знать, а что там происходит внутри lm() и как подгоняются модели. Вот тут пост, который рассматривает это с более широкого ракурса, можно с него начать ощупывание темы.