R glm многомерный синтаксис короче

Question

У меня есть пример данных.

DF<- data.frame(x=c(1,1,2,2,3,3,3), y1=c(1,2,3,4,8,8,4), y2=c(12,11,5,7,4,7,8), y3=c(2,3,4,5,3,2,1))

DF1 <- glm(DF$y1 ~ -1 + factor(DF$x)) 
DF2 <- glm(DF$y2 ~ -1 + factor(DF$x)) 
DF3 <- glm(DF$y3 ~ -1 + factor(DF$x)) 
result <- cbind(DF1$res, DF2$res, DF3$res) 

Можете ли вы помочь, как я могу GLM синтаксис модели взять короче?

Answer 1

Поскольку GLM с family = gaussian просто обычная линейная регрессия:

fit <- lm(cbind(y1,y2,y3) ~ -1 + factor(x), data=DF) 
result <- fit$res