如何計算Krippendorf在R中的Alpha的置信度？

Question

我確信這是引用Krippendorff的阿爾法 。 但我不明白這個問題，也沒有答案。 看起來連答案和評論都是相互矛盾的。

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 

我如何計算置信區間？

Answer 1

你是對的，它連接到bootstrapping。 您可以通過以下方式計算置信區間：

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

這是輸出：

 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

還有來自Zapf等人的R腳本。 2016 在參考文獻之前的頁面底部查找附加文件3

或者您可以使用github MikeGruz / kripp.boot上提供的kripp.boot函數