如何用R画Logistic回归模型的ROC曲线

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

老师们好！

我用LASSO回归做出了Logistic回归模型，下一步需要该模型的ROC曲线，请问应该如何实现？

还要画该模型对测试的ROC曲线，咋画呢？  谢谢老师啦

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：回归模型 如何实现 归模型 下一步

帮你顶顶帖子让更多人看到哈~~

啊啊啊啊啊吖 发表于 2018-12-30 15:48
帮你顶顶帖子让更多人看到哈~~

谢谢！！

可以用pROC package

wudinbg 发表于 2019-1-20 22:26
可以用pROC package

请问是这个代码吗？> install.packages("pROC")

> library("pROC")

> roc

> plot(roc, print.auc=TRUE, auc.polygon=TRUE, grid=c(0.1,0.2),grid.col=c("green","red"),

有点看不明白，我是做了随进森林，但是不会计算AUC和ROC曲线

老师们，就像这样的logistic模型，怎么画他的ROC曲线呢？？ 求指点

你可以用AUC package来画，我下面这个代码是做了一个logistic regression，然后画出了roc并标出了auc的值，你可以参考下运用到你的数据中

LR_model <- glm(as.factor(Diagnosis) ~ Season+Age+Chd_d+Acci_tra+High_fever+
                  Ach_con+Smoke,data=training.set,
                family=binomial(link = "logit"))
glm.probs <- predict(LR_model,test.set[,-10],type="response")
sum(ifelse(glm.probs>0.5,1,0)==test.set$Diagnosis)/30 #accuracy

library(AUC)
plot(roc(glm.probs,as.factor(test.set$Diagnosis)), col=rainbow(10),
     main="ROC curve for Logistic Regression model")
lr.auc <- auc(roc(glm.probs,as.factor(test.set$Diagnosis)))
text(0.8, 0.2, labels=paste('AUC=', round(lr.auc, 3), sep=''), cex=1.2)

你可以用AUC包来画，下面的代码是先建立了一个logistic regression，然后用AUC包里的roc,auc函数来画roc图并标注auc的值，你可以参考下用到你的数据中。（发完贴发现粘贴的代码排版被网站搞乱了，可以看附件图片）

##Logistic regression
LR_model <- glm(as.factor(Diagnosis) ~ Season+Age+Chd_d+Acci_tra+High_fever+
                  Ach_con+Smoke,data=training.set,
                family=binomial(link = "logit"))
glm.probs <- predict(LR_model,test.set[,-10],type="response")
sum(ifelse(glm.probs>0.5,1,0)==test.set$Diagnosis)/30 #accuracy
##plot roc and indicate auc
library(AUC)
plot(roc(glm.probs,as.factor(test.set$Diagnosis)), col=rainbow(10),
     main="ROC curve for Logistic Regression model")
lr.auc <- auc(roc(glm.probs,as.factor(test.set$Diagnosis)))
text(0.8, 0.2, labels=paste('AUC=', round(lr.auc, 3), sep=''), cex=1.2)

thanks for sharing

谢谢老师指点！