当我运行下面的代码时，我可以计算每个类别c的回归系数。现在我想知道如何应用这些估计系数来计算所有观测值的残差。例如，这里只有25个观测值属于c = 1，但我需要根据该类别的估计系数计算所有50个观测值的拟合值/残差。
A<-cars$speed
B<-cars$dist
c<-rep(1:2,25)
S<-data.frame(A,B,c)
library(plyr)
lmodel <- dlply(S,"c", function(d) lm(B~A, data = d))
答案 0 :(得分：0)
我不是100％确定我理解你的意思，但下面的代码会给你一个残差列表。列表的第一个元素包含使用c=1的系数和c=2的第二个系数的所有50个观测值的残差。
residuals<- lapply(lmodel, function(x) B - coef(x)[1] - coef(x)[2]*A)