1

R では多重回帰に対して返されることは知っていますが、βi=0 の仮説検定を返しますが、βi=1 のような検定をテストしたい場合はどうすればよいでしょうか。これには簡単なコマンドがありますか、そうでない場合は、係数の標準誤差、係数の値、回帰の自由度をどのように呼び出すので、t 分布 cdf を使用して p 値を計算できますか。一般的なプログラムで複数のデータを実行するためにこれを行いたい

4

2 に答える 2

3

R には、係数が 0 以外の値と異なるかどうかをテストできるパッケージがいくつかあります。たとえば、https://www.rforge.net/doc/packages/FSA/hoCoef.htmlです。この場合、hoCoef 関数で bo = 1 を指定します。

于 2015-03-04T23:44:26.080 に答える
2

これを行うパッケージはありますが、実行するのは非常に簡単で、小さな関数を書くことができます。

slこれは、 の値に対する同等性に対する両側検定の回帰要約の p 値を返しますb0。

testb0=function(sl,b0) {
  slm=sl$coefficients                               #$
  t0=(slm[,1]-b0)/slm[,2]
  pt(abs(t0),sl$df[2],lower.tail=FALSE)
}

その機能のテスト:

testb0( summary(lm(dist~speed+I(speed^2),cars)), b0=c(0,1,0) )

3 つの p 値を返します

(Intercept)       speed  I(speed^2) 
 0.43415754  0.48308979  0.06820122
于 2015-03-04T23:53:14.363 に答える