計算同類相關系數（ICC）R

計算同類相關系數（ICC）R語言

需要的包：readr，irr

library(readr)
library(irr)


fpath <- choose.files()

feature_1 <- read_csv(fpath)

fpath <- choose.files()

feature_3 <- read_csv(fpath)
len <- N# N是名額的數量
icc_val<-vector(length=Len)
thr <- 0.75
selected <- feature_1[feature_1$ID %in% feature_3$ID,]#擷取id相同的樣本
for (i in 2:len){#len是特征的數量，第一列是ID，從第二列開始進行名額的比較
  ratings <- cbind(selected[,i],feature_3[,i])
  icc <- icc(ratings, model = "twoway", 
             type = "agreement", 
             unit = "single", r0 = 0, conf.level = 0.95)
  icc_val[i] <- icc$value
}
Index <- which(icc_val > thr)
dim(icc_val)=c(1,len)
write.csv( icc_val,file = "output.csv",row.names = F)
           

icc處的參數：type：agreement 是不同評分者的評價是否一緻 consistency不同評分者的評分是否存在相關性

                      unit：single 是單個評分 average是取平均的評分