提取指标并填充来自多个confusionMatrix的数据框

时间:2019-04-26 19:23:32

标签: r machine-learning confusion-matrix

对于许多不同的分类模型,我有许多confusionMatrix输出。我想从每个混淆矩阵中提取模型指标(即灵敏度->平衡精度,或者从cf_1 $ byClass中获取所有指标)以填充数据框。我希望数据框包含所有这些指标,以及它来自哪个混淆矩阵的标签。谁能帮我弄清楚该怎么做?请参见下面的可复制示例:

library(caret)
library(tidyverse)

df_1 <- data.frame(x = sample(LETTERS[1:2], 20, replace = T),
                      y = sample(LETTERS[1:2], 20, replace = T))

df_2 <- data.frame(x = sample(LETTERS[1:2], 20, replace = T),
                      y = sample(LETTERS[1:2], 20, replace = T))

cf_1 <- confusionMatrix(df_1$y, df_1$x)
cf_2 <- confusionMatrix(df_2$y, df_2$x)

# I would like a dataframe with each of these, as well as column for the associated confusion matrix name
cf_1$byClass
cf_2$byClass

除了所有相关指标外,我期望的输出如下所示:

x = c(0.2, 0.4)
y = c(0.5, 0.6)
z = c("cf_1", "cf_2"
      )

output <- data.frame(Sensitivity = x, Specificity = y, Model = z)

更新

这是我想出的解决方案,虽然不是很漂亮,但是可以用。话虽如此,仍然有待观察是否有人能提出一种更清洁,更有效的方式来做到这一点!

model_names <- c("cf_1", "cf_2")
model_metrics <- list()
for (i in 1:length(model_names)) {
     for (j in model_names) {
          model <- eval(as.name(j))
          results <- model$byClass %>% round(digits = 3)
          results <- c(model = j, results)
          model_metrics[[j]] <- results
     }
}

model_metrics <- do.call(rbind, model_metrics) %>% 
     data.frame %>% 
     `rownames<-`(NULL)

model_metrics

2 个答案:

答案 0 :(得分:1)

那又怎么样:

library(caret)
library(tidyverse)

df_1 <- data.frame(x = sample(LETTERS[1:2], 20, replace = T),
                      y = sample(LETTERS[1:2], 20, replace = T))

df_2 <- data.frame(x = sample(LETTERS[1:2], 20, replace = T),
                      y = sample(LETTERS[1:2], 20, replace = T))

cf_1 <- confusionMatrix(df_1$y, df_1$x, mode = "sens_spec")
cf_2 <- confusionMatrix(df_2$y, df_2$x, mode = "sens_spec")

bind_rows(stack(cf_1$byClass), stack(cf_2$byClass), .id = "id") %>% 
  filter(ind %in% c("Sensitivity", "Specificity")) %>% 
  spread(id, values) %>% 
  rename(metric = ind, mod1 = "1", mod2 = "2")


#>        metric      mod1 mod2
#> 1 Sensitivity 0.4285714  0.5
#> 2 Specificity 0.6153846  0.4

编辑 要获取数据框,只需堆叠cf_1列表:

data.frame(stack(cf_1$byClass))

答案 1 :(得分:1)

两种方法(顺便说一句,我认为您的示例中有一个类型,df_1使用了两次)

# Base R
output <- data.frame(rbind(cf_1$byClass, cf_2$byClass))
output$model = c('cf_1', 'cf_2')

# dplyr version
output <- bind_rows(cf_1$byClass, cf_2$byClass, .id = 'model')