2017-01-30 3 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 

또한 Zapf를 등에서 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

아니면 GitHub의에서 사용할 수있는 kripp.boot 기능을 사용할 수 MikeGruz/kripp.boot

+0

왜 1000 개 복제가있다? 왜 그다지 많지 않습니까? – 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