我有一个nnet:multinom
由5个类构成的用于26个变量的多项式模型:
mirna_multinom_0 = multinom(formula_0, data= clase_training, maxit=10000 )
然后,我使用以下命令创建ROCS:
multiclass.roc(clase_training$clase, mirna_multinom_0$fitted.values,plot=TRUE)
我情节。
pred_test_inter_multinom_5 = predict(interaction_multinom_model_5, newdata = clase_test, "probs")
multiclass.roc(clase_test$clase, pred_test_inter_multinom_5,plot=TRUE)
为了理解它们,我将其存储为一个对象e
,并将其称为对比,names(e$roc)
以查看对比。
e = multiclass.roc(clase_training$clase, mirna_multinom_0$fitted.values)
names(e$rocs)
[1] "Control/Idiop_grave" "Control/Idiop_leve" "Control/Isquem_grave"
[4] "Control/Isquem_leve" "Idiop_grave/Idiop_leve" "Idiop_grave/Isquem_grave"
[7] "Idiop_grave/Isquem_leve" "Idiop_leve/Isquem_grave" "Idiop_leve/Isquem_leve"
[10] "Isquem_grave/Isquem_leve"
这给了我每个图2个图,一个在>方向,另一个在<方向。
现在。我可以以某种方式在绘图中绘制每个对比的标题吗?
而且,是否有一种方法可以获取每个ROC对比度的曲线AUC下的面积?我仅在有关多项式的消息中获得它。Wich没有PROC图。我能否获得多项式ROC,或者仅获得其无图形表示的构造?
我可以以某种方式在绘图中绘制每个对比的标题吗?
您将需要自己遍历曲线,但是可以轻松完成,例如:
for (contrast in names(e$rocs)) {
plot(e$rocs[[contrast]][[1]], col = "green", main = contrast)
plot(e$rocs[[contrast]][[2]], col = "blue", add = TRUE)
}
有没有一种方法可以针对每种ROC对比度获取曲线AUC下的面积?
您可以使用该auc
功能执行类似的操作:
for (contrast in names(e$rocs)) {
print(contrast)
print(auc(e$rocs[[contrast]][[1]]))
print(auc(e$rocs[[contrast]][[2]]))
}
我能否获得多项式ROC,或者仅获得其无图形表示的构造?
这是Hand&Till在doi:10.1023 / A:1010920819831中描述的AUC平均值。没有相应的曲线可表示。
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句