设置ROC的y值

Question

为一个非常基本的问题道歉。我很难让R识别ROC的y值

我正在尝试做一个基本的ROC，但似乎无法为y设置向量。

fullmodel= glm(culture_positive ~  No_symptoms + sex + art_status_v1 +current_cd4 + 
               bmi_v1 +nurse_tb_diagnosis_crp_v1 + temperature_v1,
               family="binomial",  data= Data1)
roc(y , fullmodel$fitted.values, plot=TRUE)

roc中的错误（y，fullmodel $ fitted.values，plot = TRUE）：找不到对象'y'

所以'y'是我的数据集Data1中的一列，标记为'culture_positive'，按照glm，但无论我尝试什么，我都会收到这条消息'y'找不到。

再一次，为一个基本问题道歉，但它真的让我感到高兴。

Answer 1

由于y不在您的全球环境中，您需要指定在哪里找到y。您可以使用您用于拟合模型的值：

roc(culture_positive , fullmodel$fitted.values, plot=TRUE)

或存储在glm对象中的响应

roc(fullmodel$y , fullmodel$fitted.values, plot=TRUE)

我会推荐第二种选择，它有点安全，因为你从同一个物体中取出y和fitted.values，所以它们会合在一起。