我有一个混淆矩阵,我想从中计算生产者和用户的准确性,并获得整体 Kappa 和 AIC 值。我知道 R 有几个包可以计算 Cohen 的 kappa 统计量,但是我找不到任何处理混淆矩阵的包(或者我可能没有进行足够广泛的搜索)。我的混淆矩阵是 50 x 50 矩阵,因此手动计算容易出错。
所以对于矩阵的样本:
V.1 V.2 V.3 V.4 V.5
V.1 10 0 2 1 0
V.2 1 12 0 0 0
V.3 0 0 9 4 0
V.4 0 2 0 11 3
V.5 4 1 0 0 11
因此,用户对 V.1 的准确度将为(10/13)*100%
. 同样,生产者的 V.1 (10/15)*100%
。我必须cbind
为用户的准确性添加一个新列,并rbind
为生产者的准确性添加一行。我假设我可以使用melt
包reshape2
中的数据来格式化可用于计算 kappa 的表中的数据。但理想情况下,我想以混淆矩阵格式保存数据。
理想的最终结果如下所示:
V.1 V.2 V.3 V.4 V.5 User Acc
V.1 10 0 2 1 0 76.92
V.2 1 12 0 0 0 92.30
V.3 0 0 9 4 0 69.23
V.4 0 2 0 11 3 68.75
V.5 4 1 0 0 11
Prod 66.66
kappa = 价值,AIC = 价值
任何指针将不胜感激