在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都取决于我们迄今为止缺乏的具体细节。