我已经成功创建了2x2列联表。列联表是通过以下过程创建的:
plus <- ua.mask + vm.c.m.s;
minus <-ua.mask - vm.c.m.s;
这里,ua.mask和vm.cms是栅格图层(可以视为矩阵),其中单元格的值为0(非城市区域)和1(城市区域).ua.mask是参考图和vm.cms是一个精确度是我要检查的地图。
m <- matrix(0, nrow=2, ncol=2)
m[1,1] <- length(which(plus[]==2))
m[1,2] <- length(which(minus[]==-1))
m[2,1] <- length(which(minus[]==1))
m[2,2] <- length(which(plus[]==0))
因此,当我添加和减去这些栅格图层时,单元格的值将为1,-1,0和2.这意味着,我将有4种可能的结果,其中2 = TP,-1 = FP,1 = TP,0 = TN。 在最后,我得到一个2x2列联表,如下所示:
[,1] [,2]
[1,] 2 -1
[2,] 1 0
然后,我计算了灵敏度(TPR)和1-特异性(FPR)。
现在,我想绘制所有灵敏度和1-特异性对,为所有城市创建ROC曲线(灵敏度与1-特异性),但我无法弄清楚如何做到这一点。
我在R中浏览了ROCR包,但无法理解函数(预测和执行)在我的情况下是如何工作的。
或者,我的输入数据是否适用于ROC-Curves?
非常感谢您提供的任何帮助!
谢谢!