用R中的顏色和頻率繪制混淆矩陣

Question

我想繪制一個混淆矩陣，但是，我不想只使用熱圖，因為我認為它們的數值分辨率很差。 相反，我還想在正方形的中間繪制頻率。 例如，我喜歡這個的輸出：

library(mlearning);
data("Glass", package = "mlbench")
Glass$Type <- as.factor(paste("Glass", Glass$Type))

summary(glassLvq <- mlLvq(Type ~ ., data = Glass));
(glassConf <- confusion(predict(glassLvq, Glass, type = "class"), Glass$Type))

plot(glassConf) # Image by default

但是，1。）我不明白“01,02等”是指每個軸。 我們怎樣才能擺脫這種局面？ 2.）我希望'Predicted'作為'y'維度的標簽，'Actual'作為'x'維度的標簽3.）我想用頻率/概率替換絕對計數。

或者，是否有其他包可以做到這一點？

從本質上講，我希望在R中：

http://www.mathworks.com/help/releases/R2013b/nnet/gs/gettingstarted_nprtool_07.gif

要么：

http://c431376.r76.cf2.rackcdn.com/8805/fnhum-05-00189-HTML/image_m/fnhum-05-00189-g009.jpg

Answer 1

通過繪制混淆矩陣， mlearning包似乎非常不靈活。

從您的glassConf對象開始，您可能希望執行以下操作：

prior(glassConf) <- 100 
# The above rescales the confusion matrix such that columns sum to 100.
opar <- par(mar=c(5.1, 6.1, 2, 2))
x <- x.orig <- unclass(glassConf)
x <- log(x + 0.5) * 2.33
x[x < 0] <- NA
x[x > 10] <- 10
diag(x) <- -diag(x)
image(1:ncol(x), 1:ncol(x),
      -(x[, nrow(x):1]), xlab='Actual', ylab='',
      col=colorRampPalette(c(hsv(h = 0, s = 0.9, v = 0.9, alpha = 1), 
                             hsv(h = 0, s = 0, v = 0.9, alpha = 1), 
                             hsv(h = 2/6, s = 0.9, v = 0.9, alpha = 1)))(41), 
      xaxt='n', yaxt='n', zlim=c(-10, 10))
axis(1, at=1:ncol(x), labels=colnames(x), cex.axis=0.8)
axis(2, at=ncol(x):1, labels=colnames(x), las=1, cex.axis=0.8)
title(ylab='Predicted', line=4.5)
abline(h = 0:ncol(x) + 0.5, col = 'gray')
abline(v = 0:ncol(x) + 0.5, col = 'gray')
text(1:6, rep(6:1, each=6), 
     labels = sub('^0$', '', round(c(x.orig), 0)))
box(lwd=2)
par(opar) # reset par

上面的代碼使用了plot.confusion調用的confusionImage函數的plot.confusion 。

Answer 2

這是一個繪制混淆矩陣的函數，我從jbaums開發出色的答案。
它是相似的，但看起來更好（IMO），並沒有轉置你喂它的混亂矩陣，這可能會有所幫助。

### Function for plotting confusion matrices
confMatPlot = function(confMat, titleMy, shouldPlot = T) {
  #' Function for plotting confusion matrice
  #' 
  #' @param confMat: confusion matrix with counts, ie integers. 
  #' Fractions won't work
  #' @param titleMy: String containing plot title
  #' @return Nothing: It only plots

  ## Prepare data
  x.orig = confMat; rm(confMat)  # Lazy conversion to function internal variable name
  n = nrow(x.orig)  # conf mat is square by definition, so nrow(x) == ncol(x)
  opar <- par(mar = c(5.1, 8, 3, 2))
  x <- x.orig
  x <- log(x + 0.5)  # x<1 -> x<0 ,  x>=1 -> x>0
  x[x < 0] <- NA
  diag(x) <- -diag(x)  # change sign to give diagonal different color
  ## Plot confusion matrix
  image(1:n, 1:n,  # grid of coloured boxes
        # matrix giving color values for the boxes
        # t() and [,ncol(x):1] since image puts [1,1] in bottom left by default
        -t(x)[, n:1],  
        # ylab added later to avoid overlap with tick labels
        xlab = 'Actual', ylab = '',
        col = colorRampPalette(c("darkorange3", "white", "steelblue"), 
                               bias = 1.65)(100),
        xaxt = 'n', yaxt = 'n'
        )
  # Plot counts
  text(rep(1:n, each = n), rep(n:1, times = n), 
       labels = sub('^0$', '', round(c(x.orig), 0)))
  # Axis ticks but no lables
  axis(1, at = 1:n, labels = rep("", n), cex.axis = 0.8)
  axis(2, at = n:1, labels = rep("", n), cex.axis = 0.8)
  # Tilted axis lables
  text(cex = 0.8, x = (1:n), y = -0.1, colnames(x), xpd = T, srt = 30, adj = 1)
  text(cex = 0.8, y = (n:1), x = +0.1, colnames(x), xpd = T, srt = 30, adj = 1)
  title(main = titleMy)
  title(ylab = 'Predicted', line = 6)
  # Grid and box
  abline(h = 0:n + 0.5, col = 'gray')
  abline(v = 0:n + 0.5, col = 'gray')
  box(lwd = 1, col = 'gray')
  par(opar)
}

輸出示例：