2017-01-30 11 views
0

これはBootstrapping Krippendorff's Alphaに翻訳されていると確信しています。しかし、私はその質問や答えを理解していませんでした。そして、たとえ答えとコメントでさえ互いに矛盾しているように見えます。RのKrippendorf's Alphaの信頼区間を計算するには?

set.seed(0) 
df <- data.frame(a = rep(sample(1:4),10), b = rep(sample(1:4),10)) 
kripp.alpha(t(df)) 

これは出力です。

Krippendorff's alpha 

Subjects = 40 
    Raters = 2 
    alpha = 0.342 

ここでは信頼区間をどのように計算できますか?

答えて

2

あなたは正しくブートストラップに接続されています。あなたは、信頼区間を次のように計算することができます:

Bootstrap Statistics : 
     original  bias std. error 
t1* 0.3416667 -0.01376158 0.1058123 

BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS 
Based on 1000 bootstrap replicates 

CALL : 
boot.ci(boot.out = b, type = "perc") 

Intervals : 
Level  Percentile  
95% (0.1116, 0.5240) 
Calculations and Intervals on Original Scale 

もツァップらからRスクリプトがあります:

library(irr) 
library(boot) 

alpha.boot <- function(d,w) { 
     data <- t(d[w,]) 
     kripp.alpha(data)$value 
} 

b <- boot(data = df, statistic = alpha.boot, R = 1000) 
b 
plot(b) 
boot.ci(b, type = "perc") 

これが出力されます。 2016 look for Additional file 3 at the bottom of the page just before the references

または1000回の反復があるなぜあなたはgithubのMikeGruz/kripp.boot

+0

にkripp.boot機能が利用可能に使用することができますか?なぜそれ以上のものはありませんか? – buhtz

+1

[Zapf et al 2016](https://bmcmedresmethodol.biomedcentral.com/articles/10.1186/s12874-016-0200-9)は[Efron](https://statistics.stanford.edu/sites/)のため1000を使用しました。デフォルト/ファイル/ BIO%20139.pdf)。 – Kev

+0

私は紙を正しく理解すると、あなたの答えはFleiss Kappaにも使えますか? – buhtz

関連する問題