Как я могу добавить заголовок к графикам multiclass.ro c в pRO C? - PullRequest
Новая
       30

Как я могу добавить заголовок к графикам multiclass.ro c в pRO C?

1 голос
/ 24 марта 2020

У меня есть полиномиальная модель, построенная из 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)

Plot of first contrast Ctrl vs Group1

Чтобы понять их, я сохраняю их как объект 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 графика для каждого из них, 1 в> направлении и другой в <направлении. </p>

Сейчас. Могу ли я каким-то образом построить названия каждого контраста на графиках?

А также, есть ли способ получить области под кривой AU C для каждого один из RO C контраст? Я получаю его только в сообщении для полинома. У которого нет PRO C сюжета. Могу ли я получить полиномиальный RO C, или это просто конструкция без графического представления?

1 Ответ

0 голосов
/ 25 марта 2020

Могу ли я каким-то образом нанести названия каждого контраста на графиках?

Вам нужно будет l oop по кривым самостоятельно, но это можно легко сделать, например: 

for (contrast in names(e$rocs)) {
    plot(e$rocs[[contrast]][[1]], col = "green", main = contrast)
    plot(e$rocs[[contrast]][[2]], col = "blue", add = TRUE)
}

есть ли способ получить области под кривой AU C для каждого контраста RO C?

Можно сделать что-то похожее с функцией auc: 

for (contrast in names(e$rocs)) {
    print(contrast)
    print(auc(e$rocs[[contrast]][[1]]))
    print(auc(e$rocs[[contrast]][[2]]))
}

Можно ли получить полиномиальный RO C, или это просто конструкция без графического представления?

Это среднее значение AU C, описанное Hand & Till в doi: 10.1023 / A: 1010920819831 . Нет соответствующей кривой для представления.

...