r语言 - 设置 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">

因此，根据glm，"y"是我的数据集Data1中的一列，标记为"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，所以它们会放在一起。