r - 在 R 或 winbugs 中绘制拟合值与观察值

Question

我想绘制拟合值与观察值，并希望用直线显示拟合优度。但是，我不想使用abline()，因为我没有使用lm命令计算拟合值，因为我使用的模型R不包括在内。我计算了系数并用它们来计算拟合值。那么，我该怎么做才能在 R 或 winbugs 中获得这样的情节？

这就是我想要的

Answer 1

仍然没有提供数据，但也许这个使用该curve函数的简单示例将通知该过程：

x <- 1:10
y <- 2+ 3*(1:10) + rnorm(10)
plot(1:10, y)
curve( 2+3*x, 0, 10,  add=TRUE)

新 R 用户请注意。表达式y_i = 1 - xbeta + delta_i + e_i在 R 中会失败，部分原因是 x 和 beta 没有被运算符分隔。但是，如果您确实了解 R 的矩阵语法，即使“X”是多维的，它也可能是一个非常紧凑的表达式。所有 htis 都取决于我们迄今为止缺乏的细节。