私が実行した一連の線形モデルがあり、特定の年からの残差を出力し、それらの残差を企業に一致させたいと考えています。
データは次のようになります。
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......
ヒントをいただければ幸いです。