凌晨两点,如此艰难而又幸运地完成了一张ROC曲线图:
使用的是B站德民的代码。由于发现了列名不匹配这一问题,经过一番折腾后问题被巧妙解决了,上代码:
for(i in 1:length(hubgenes)){
x_color = c('blue','red','green','black','yellow','orange')
plot.roc(group,as.numeric(hubgenes_expression[i,2:67]),main = hubgenes[i],
col = x_color[i],print.auc = T,percent = T,cex.lab = 1.5,print.auc.cex=1.5)
}
感谢B站爱学习又乐于分享的热心网友。