R：列表中線性模型的殘差

Question

我已經運行了一系列線性模型，我希望能夠輸出特定年份的殘差並將這些殘差匹配到公司。

數據類似於下圖所示

Year    Company         pctEarn
1990    3M Company      0.295918367
1991    AT&T            0.251497006
1992    Ford            0.293233083
1993    Microsoft       0.264705882

運行LM的循環如下所示

fits <- list()

for(i in 2000:2012){
  dtSub <- subset(dt, Year <= i)
  fit <- lm(pctEarn ~ Year + Company, data = dtSub)
  fits[[i - 1999]] <- fit
}

我遇到的問題是從輸出中檢索特定公司的殘差。

理想情況下，我將能夠為如下所示的每個回歸創建輸出

Year Company    pctEarn    Predicted  Residual
1990 3M Company 0.2959     0.4523     0.16000
Etc......

任何提示將不勝感激。

Answer 1

因此，我發現對於大型回歸而言，不理想的方法是簡單地建立一個累積回歸結果的數據框，以后再將其剪切。

out <- data.frame()
for(i in 2000:2012){
  dtSub <- subset(dt, Year <= i)
  fit <- lm(pctEarn ~ y2 + Company, data = dtSub)
  fits[[i - 1999]] <- fit
  dtSub$fitted <- fitted(fit)
  dtSub$resid <- residuals(fit)
  dtSub$reg <- i
  out <- rbind(out, dtSub)
}