我想绘制拟合值与观察值之间的拟合值,并希望将直线显示出拟合的好处。但是,我不想使用abline()
,因为我没有使用lm
命令计算拟合值,因为我使用了R
未涵盖的模型。我计算了系数并用它们来计算拟合值。那么,我能做些什么来在R或winbugs中获得这样的情节?在R或winbugs中绘制拟合值与观察值的拟合值
这里是我想要的
我想绘制拟合值与观察值之间的拟合值,并希望将直线显示出拟合的好处。但是,我不想使用abline()
,因为我没有使用lm
命令计算拟合值,因为我使用了R
未涵盖的模型。我计算了系数并用它们来计算拟合值。那么,我能做些什么来在R或winbugs中获得这样的情节?在R或winbugs中绘制拟合值与观察值的拟合值
这里是我想要的
仍然没有提供的数据,但也许使用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”是多维的,它也可能是一个非常紧凑的表达式。所有这些都取决于我们迄今为止缺乏的细节。
明确的想法Dwin。谢谢 –
我找不到任何重要的东西 –
你对这条线有什么了解?拦截和坡度会要求'abline'。两点要求计算截距和斜率,然后使用'abline'。还有什么可以呢? – Thilo
我的线性模型是单向误差分量随机效应模型,其规格为y_i = 1 -xbeta + delta_i + e_i其中x是设计矩阵,beta是回归系数的向量,delta_i是第i个随机效应, e是第i个常见错误。我获得了该模型的可能性,然后通过使用R中的优化值和拟合值获得了β。另外请注意,没有拦截 –