计算同类相关系数（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是取平均的评分