r - Rまたはwinbugsでの適合値と観測値のプロット

Question

適合値と観測値をプロットし、適合度を示す直線を配置したいと考えています。ただし、カバーしないモデルを使用したため、コマンドをabline()使用して適合値を計算しなかったため、使用したくありません。係数を計算し、それらを使用して適合値を計算しました。では、R または winbugs でそのようなプロットを取得するにはどうすればよいでしょうか?lmR

これが私が欲しいものです

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_ix とベータが演算子で区切られていないため、式は R では部分的に失敗します。しかし、R の行列構文を理解していれば、"X" が多次元であっても、非常にコンパクトな式になる可能性があります。すべての htis は、これまでのところ欠けている詳細に依存します。