6

Rで回帰を実行するとき、返される係数の順序は? 例えば：

 coef(lm(y ~ x + z, data=data.frame(x=1:10, y=10:1, z=1:5)))

xに関連付けられた係数が常に に関連付けられた係数の前に返されることが保証されていzますか? 順序とは、返された係数のベクトル内の順序を意味します。これが私にとって重要な理由は、モデルの係数に関する線形仮説をテストしたいため、によって返される分散共分散行列の係数の順序と、問題vcovによって返される実際の推定値coefです。

4

1 に答える 1

1

位置ではなく、名前でインデックスを作成します。そうすれば、必ず正しい答えが得られます。

coef(lm(y ~ x+z, data=data.frame(x=1:10, y=10:1, z=1:5)))['x']
##  x 
## -1 
coef(lm(y ~ x+z, data=data.frame(x=1:10, y=10:1, z=1:5)))['z']
##             z 
## -1.855301e-16

そして、それらの両方を、望ましい順序で：

coef(lm(y ~ x+z, data=data.frame(x=1:10, y=10:1, z=1:5)))[c('x', 'z')]
##             x             z 
## -1.000000e+00 -1.855301e-16
于 2013-05-19T03:38:57.743 に答える