数据data中Y为无序四分类因变量，已经有真实标签、预测标签，怎样绘制ROC曲线？使横纵坐标为TPR、FPR？请输出R代码。

请见以下R代码：

library(pROC)

#这里假设真实标签存储在y_true中，预测标签存储在y_pred中
#分别为0、1、2、3四个类别，将其转换为二分类标签
roc_data <- roc(as.factor(y_true==0), as.numeric(y_pred==0))
roc_data <- roc(as.factor(y_true==1), as.numeric(y_pred==1), add=TRUE)
roc_data <- roc(as.factor(y_true==2), as.numeric(y_pred==2), add=TRUE)
roc_data <- roc(as.factor(y_true==3), as.numeric(y_pred==3), add=TRUE)

#绘制ROC曲线
plot(roc_data, main="ROC Curve")

这段代码将四分类问题转换为四个二分类问题，使用pROC包中的roc函数计算每个分类器的ROC曲线，最后绘制在同一张图中。横轴为FPR，纵轴为TPR。注意，这里ROC曲线的斜率并不能反映分类器的性能优劣，只有曲线下面积（AUC）才能做出判断。