Вы можете использовать такой пакет, как pROC, что может быть проще, чем создавать его самостоятельно. Вы используете функцию roc () для создания объекта roc, а затем plot () для построения объекта и создания кривой ROC.
Кроме того, не похоже, что вы делаете это правильно. Кривая ROC рассчитывается из значений чувствительности и специфичности, взятых в разных точках разреза в диапазоне возможных вероятностей вашей модели, т.е. в точках разреза между 0 и 1. Вам не нужно создавать 100 наборов предсказаний модели, только одна модель будет достаточно.
Попробуйте что-то вроде этого, где у - переменная ответа, а р - вектор значений вероятности, выведенных вашей моделью:
plot(roc(y, p)), print.auc = TRUE)