Как создать ROC-кривую для модели, в которой в качестве ответа используется журнал шансов? - PullRequest
0 голосов
/ 24 апреля 2019

У меня есть вопрос о графике ROC curve для моей модели, в которой log of odds имеет значение response.Например:

model<-lm((ln(y/1-y)~Temp+RH+DmaxT, data=fit) #'y' is a proportion

Предсказанный ответ был получен для нового набора данных:

Predicted_model<-predict(model, newdata, type = 'response')

Предсказанные значения были преобразованы обратно для получения значений в пропорции

Iу меня есть новые наблюдения в пропорции, и я использовал 0.05 предельное значение для представления контроля (<0,05) и случаев (> 0,05)

newdata$observed<-ifelse(newdata$observed > 0.05, "cases", "controls")

Я построил кривую ROC, используя следующую формулу

roc(newdata$observed, predicted_model_backtrans, legacy.axes = TRUE, plot = TRUE, print.auc = TRUE)

С этой формулой я получил AUC значение 1, и график отличается от ожидаемого.Я не мог понять, что было бы лучшим способом создать ROC curve для моей модели.Буду признателен за любую помощь.

Я также попытался создать кривую ROC, где я изменил наблюдаемую и предсказанную пропорцию на бинарные характеристики (контроль (<0.05) и случаи (>0.05)), что дало мне прямую кривуюа не гладкий.

...