Caret中的附加指标 – PPV、敏感性、特异性

我在R中使用caret进行逻辑回归:

  ctrl <- trainControl(method = "repeatedcv", number = 10, repeats = 10,                        savePredictions = TRUE)  mod_fit <- train(Y ~ .,  data=df, method="glm", family="binomial",                   trControl = ctrl)  print(mod_fit)

默认输出的指标是准确率和Cohen kappa。我想提取匹配的指标,如敏感性、特异性、阳性预测值等,但找不到简单的方法来实现。最终模型是提供的,但它是用所有数据训练的(据我从文档中了解),所以我无法用它来进行新的预测。

混淆矩阵可以计算所有所需的参数,但将其作为摘要函数传递时不起作用:

  ctrl <- trainControl(method = "repeatedcv", number = 10, repeats = 10,                        savePredictions = TRUE, summaryFunction = confusionMatrix)  mod_fit <- train(Y ~ .,  data=df, method="glm", family="binomial",                   trControl = ctrl)Error: `data` and `reference` should be factors with the same levels. 13.stop("`data` and `reference` should be factors with the same levels.",     call. = FALSE) 12.confusionMatrix.default(testOutput, lev, method) 11.ctrl$summaryFunction(testOutput, lev, method) 

有没有办法除了准确率和kappa之外,还能提取这些信息,或者在caret train返回的train_object中找到这些信息?

提前感谢!


回答:

Caret已经有摘要函数可以输出您提到的所有指标:

defaultSummary 输出准确率和Kappa
twoClassSummary 输出AUC(ROC曲线下面积 – 见答案的最后一行)、敏感性和特异性
prSummary 输出精确度和召回率

为了获得组合指标,您可以编写自己的摘要函数,将这三个函数的输出结合起来:

library(caret)MySummary  <- function(data, lev = NULL, model = NULL){  a1 <- defaultSummary(data, lev, model)  b1 <- twoClassSummary(data, lev, model)  c1 <- prSummary(data, lev, model)  out <- c(a1, b1, c1)  out}

让我们在Sonar数据集上试试:

library(mlbench)data("Sonar")

在定义训练控制时,设置classProbs = TRUE 很重要,因为有些指标(ROC和prAUC)不能基于预测类别计算,而是基于预测概率计算。

ctrl <- trainControl(method = "repeatedcv",                     number = 10,                     savePredictions = TRUE,                     summaryFunction = MySummary,                     classProbs = TRUE)

现在拟合您选择的模型:

mod_fit <- train(Class ~.,                 data = Sonar,                 method = "rf",                 trControl = ctrl)mod_fit$results#output  mtry  Accuracy     Kappa       ROC      Sens      Spec       AUC Precision    Recall         F AccuracySD   KappaSD1    2 0.8364069 0.6666364 0.9454798 0.9280303 0.7333333 0.8683726 0.8121087 0.9280303 0.8621526 0.10570484 0.21620772   31 0.8179870 0.6307880 0.9208081 0.8840909 0.7411111 0.8450612 0.8074942 0.8840909 0.8374326 0.06076222 0.12218443   60 0.8034632 0.6017979 0.9049242 0.8659091 0.7311111 0.8332068 0.7966889 0.8659091 0.8229330 0.06795824 0.1369086       ROCSD     SensSD    SpecSD      AUCSD PrecisionSD   RecallSD        FSD1 0.04393947 0.05727927 0.1948585 0.03410854  0.12717667 0.05727927 0.084829632 0.04995650 0.11053858 0.1398657 0.04694993  0.09075782 0.11053858 0.057723883 0.04965178 0.12047598 0.1387580 0.04820979  0.08951728 0.12047598 0.06715206

在这个输出中,ROC实际上是ROC曲线下面积 – 通常称为AUC,而AUC是精确度-召回率曲线在所有截止点下的面积。

Related Posts

使用LSTM在Python中预测未来值

这段代码可以预测指定股票的当前日期之前的值,但不能预测…

如何在gensim的word2vec模型中查找双词组的相似性

我有一个word2vec模型,假设我使用的是googl…

dask_xgboost.predict 可以工作但无法显示 – 数据必须是一维的

我试图使用 XGBoost 创建模型。 看起来我成功地…

ML Tuning – Cross Validation in Spark

我在https://spark.apache.org/…

如何在React JS中使用fetch从REST API获取预测

我正在开发一个应用程序,其中Flask REST AP…

如何分析ML.NET中多类分类预测得分数组?

我在ML.NET中创建了一个多类分类项目。该项目可以对…

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注