3

vrtest パッケージ (および Lo.Mac コマンド) を使用して分散比テストを実行しました。私の問題は、Rが計算されたテスト統計を提供するだけで、テストのフロンティア値、p値などを提供しないことです。誰かがこれを行う方法を知っていますか? これが私がやったことです：

> Lo.Mac(rcorr,kvec)
$Stats
             M1         M2
k=2   2.2525920  2.1072370
k=5   0.8903658  0.8383312
k=10 -0.5530420 -0.5388423

> str(Lo.Mac(rcorr,kvec))
List of 1
 $ Stats: num [1:3, 1:2] 2.253 0.89 -0.553 2.107 0.838 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:3] "k=2" "k=5" "k=10"
  .. ..$ : chr [1:2] "M1" "M2"

どうもありがとうございました ！

ジェラルディン

4

1 に答える 1

1

見る?Boot.test

この関数は、Lo-MacKilay (1988) および Chow-Denning (1993) 検定のブートストラップ p 値を返します。

Boot.test(rcorr,kvec,nboot=500,wild="Normal")
$Holding.Period
[1]  2  5 10

$LM.pval
[1] 0.012 0.284 0.810

$CD.pval
[1] 0.028

$CI
          2.5%    97.5%
k=2  -2.013676 1.733034
k=5  -1.725987 1.745693
k=10 -1.562838 1.815892
于 2013-01-07T07:49:37.557 に答える