数据库代码公开:二分类ROC计算
作者:半步博导 时间:2024年6月29日 21:40
阅读量:66
#以MKI67为例
geneName = "MKI67"
#以LUAD为例
CancerName = "LUAD"
#引用R包
library(pROC)
#读取基因表达数据
data = read.table(paste0("Expression.",geneName,".zscore.txt"),
header=T,
sep="\t",
check.names=F,
row.names=1)
data = data[(data[,"CancerType"] == CancerType),]
y = ifelse(data$Type == "Normal", 0, 1)
roc1 = roc(y, as.numeric(data[,geneName]))
ci1 = ci.auc(roc1, method="bootstrap")
ciVec = as.numeric(ci1)
ciVec[2]
paste0("95% CI: ",
sprintf("%.03f",ciVec[1]),
"-",
sprintf("%.03f",ciVec[3]))