我有一个“长”数据框定义为:
q <- data.frame(Indicator.Code=factor(),Year=numeric(),Value=numeric())
我试图在单个 xyplot 中绘制每个不同年份的值作为年份的函数Indicator.Code
, 如下
xyplot( Value~Year,data=q,group=Indicator.Code)
到目前为止,一切都很好。现在我尝试添加与线性回归相对应的线
rlm(q$Value[q$Indicator.Code==a]~q$Year[q$Indicator.Code==a])
对于所有的值Indicator.Code
.
我不知道怎么做。添加回归线的常用方法,即
xyplot( Value~Year,data=q,group=Indicator.Code),
panel = function(x, y) {
panel.xyplot(x, y)
panel.abline(rlm(y ~ x))
}))
无法正常工作(它计算单个回归,并为整个数据集添加单个回归线)。此外,我已经计算了回归(除了图形之外,我也需要它们),并且讨厌必须重新计算它们的想法。
新手可以遵循的任何提示吗?