设置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中根据glm标记为'culture_positive'的列，但是无论如何，我不断收到此消息，提示未找到'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，所以它们可以放在一起。