Boa noite
Estou precisando de ajuda para fazer a curva ROC no R.
Estou utilizando o seguinte comando:
library(ROCR) pred <- prediction(predictions = banco$ic, labels = banco$SM) perf <- performance(pred, measure = "tpr", x.measure = "fpr") plot(perf, main = "ic", col = "blue", lwd = 3) abline(a = 0, b = 1, lwd = 2, lty = 2) perf.auc <- performance(pred, measure = "auc") unlist([email protected])
library(ROCR)
pred <- prediction(predictions = banco$ic, labels = banco$SM)
perf <- performance(pred, measure = "tpr", x.measure = "fpr")
plot(perf, main = "ic", col = "blue", lwd = 3)
abline(a = 0, b = 1, lwd = 2, lty = 2)
perf.auc <- performance(pred, measure = "auc")
unlist([email protected])
A variável "ic" é contínua, e significa índice de conicidade. A variável SM é categórica (binária), e significa síndrome metabólica (0 - sem SM; 1 - com SM).
Fiz a curva com essas variáveis no SPSS e ela fica acima da linha diagonal com AUC = 0,7519.
Fiz também no R, com as mesmas variáveis, mas a curva fica abaixo da linha diagonal com AUC = 0,2481. Esse valor de AUC que deu no R é complementar ao que deu no SPSS.
No SPSS sei que é necessário identificar o estado da variável (foi colocado 1, pois 1 representa os doentes) e em opções é necessário identificar a direção (a opção marcada foi "um resultado do teste menor indica um teste mais positivo").
No entanto, no R não sei qual o comando que deve ser utilizado para identificar a direção, como se deve fazer no SPSS. Testei a curva no R também com 0 = com SM e 1 = sem SM, e deu a mesma coisa.
Alguém poderia me informar qual comando devo usar no R para que a curva fique igual ao SPSS?, ou seja, para que a curva AUC no R fique invertida!!