2013-01-13 41 views
0

我想绘制拟合值与观察值之间的拟合值,并希望将直线显示出拟合的好处。但是,我不想使用abline(),因为我没有使用lm命令计算拟合值,因为我使用了R未涵盖的模型。我计算了系数并用它们来计算拟合值。那么,我能做些什么来在R或winbugs中获得这样的情节?在R或winbugs中绘制拟合值与观察值的拟合值

这里是我想要的

enter image description here

+0

我找不到任何重要的东西 –

+0

你对这条线有什么了解?拦截和坡度会要求'abline'。两点要求计算截距和斜率,然后使用'abline'。还有什么可以呢? – Thilo

+0

我的线性模型是单向误差分量随机效应模型,其规格为y_i = 1 -xbeta + delta_i + e_i其中x是设计矩阵,beta是回归系数的向量,delta_i是第i个随机效应, e是第i个常见错误。我获得了该模型的可能性,然后通过使用R中的优化值和拟合值获得了β。另外请注意,没有拦截 –

回答

1

仍然没有提供的数据,但也许使用curve功能这个简单的例子将通知过程:

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

enter image description here

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

+0

明确的想法Dwin。谢谢 –