2013-01-06 31 views
3

我使用vrtest軟件包(和Lo.Mac命令)運行了Variance Ratio測試。我的問題是,R只是給我提供了計算的測試統計數據,而不是測試的邊界值,p值等。有人知道如何做到這一點?這是我所做的:vrtest包:擴展文本答案?

> 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" 

非常感謝!

傑拉爾丁

回答

1

看到?Boot.test

該函數返回羅MacKilay(1988)和 周-丹寧(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 
+0

好了,然後你用你從Lo.Mac函數得到的結果進行比較呢?這對我來說有點困惑。 – Geraldine