Я пытаюсь адаптировать код, который я использовал для создания кривых ROC в прошлом, для очень похожих данных, но создается новое сообщение об ошибке, которое не происходит с другим набором данных, даже если они отформатированы одинаково. Вот код:
control<-subset(thesis_rep_dat,view==0 & nso==0)
gvnso<-subset(thesis_rep_dat,view==0 & nso==1)
pvcont<-subset(thesis_rep_dat,view==1 & nso==0)
pvnso<-subset(thesis_rep_dat,view==1 & nso==1)
##Now let's create ROCs for each condition
rocobj1<roc(controls=control$conf3[control$ta_tp==0],cases=control$conf3[control$ta_tp==1],
levels=c(0,1),direction="<",partial.auc=FALSE)
rocobj2 <-roc(controls=gvnso$conf3[gvnso$ta_tp==0],cases=gvnso$conf3[gvnso$ta_tp==1],
levels=c(0,1),direction="<",partial.auc=FALSE)
rocobj3 <- roc(controls=pvcont$conf3[pvcont$ta_tp==0],cases=pvcont$conf3[pvcont$ta_tp==1],
levels=c(0,1),direction="<",partial.auc=FALSE)
rocobj4 <- roc(controls=pvnso$conf3[pvnso$ta_tp==0],cases=pvnso$conf3[pvnso$ta_tp==1],
levels=c(0,1),direction="<",partial.auc=FALSE)
R возвращает следующее сообщение об ошибке:
Ошибка в roc.default (controls = control $ conf3 [gvnso $ ta_tp == 0],case = control $ conf3 [gvnso $ ta_tp ==: кейсы и элементы управления должны быть числовыми или упорядоченными.
Я пытался использовать аргумент:
control$ta_tp<-as.numeric(control$ta_tp)
, предшествующийконструирование объектов ROC, но оно возвращает ту же ошибку.
Извините за наивность с этим!