В R функция lm позволяет одновременно устанавливать несколько отдельных линейных моделей - из документации для lm:
Если ответМатрица. Линейная модель подбирается отдельно по методу наименьших квадратов к каждому столбцу матрицы.
Но как вы можете подгонять несколько отдельных моделей одновременно для GLM, например, для логистической регрессии?
glm.fit требует, чтобы ответ y был вектором.Если вы зададите ему матрицу в аргументе y, то он просто использует первый столбец матрицы y и создаст одну модель GLM.