- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
您好,我是 R 编程的新手,我在不使用任何包的情况下绘制 ROC 曲线时遇到了麻烦。
我使用以下方法生成我的数据:
d=rpearsonIII(100,0.5,360,20)
nd=rnorm(100,450,25)
我需要一个 d 值 <400 和 nd 值 >400 的向量,所以我做到了:
spec = (cumsum(nd[nd>400])/sum(nd))*100
sens = (cumsum(d[d<400])/sum(nd))*100
我是这样绘制的:
plot(1-spec,sens)
但是剧情和我想象的完全不一样
编辑:感谢给出的建议,我的代码现在看起来像这样:
sc2 = c(rnorm(50,450,25),rpearsonIII(50,0.5,360,20))
scF = sc2 < 395
thresholds <- sort(sc2)
pos <- sum(scF);pos
neg <- sum(!scF);neg
tn <- cumsum(!scF);tn
spec <- tn/neg;spec
tp <- pos - cumsum(scF);tp
sens <- tp/pos;sens
plot(1 - spec, sens, type = "l", col = "red",
ylab = "Sensitivity", xlab = "1 - Specificity")
abline(c(0,0),c(1,1))
绘制的 roc 曲线如下所示: roc curve
我现在的问题是,如果改变生成数据的顺序(rnorm 和 rpearsonIII),曲线就会反转。
最佳答案
我不知道 rpearsonIII
是什么,所以我将制作一个样本随机数据,其中包含实际类别 actuals
以及预测分数 分数
。
set.seed(100)
actuals <- sample(c(TRUE,FALSE), 100, replace = TRUE)
scores <- runif(100,-1,1)
如果您的数据中的实际值是字符串或因子而不是逻辑值,您需要使用以下方法将它们转换为逻辑值:
actuals <- actuals == "postiveClass"
接下来我们要根据分数对实例进行排序。我们可以使用:
actuals <- actuals[order(scores)]
如果您想跟踪灵敏度和特异性的阈值,可以使用以下方法使它们保持一致:
thresholds <- sort(scores)
现在我们需要了解我们的敏感性和特异性。灵敏度为 TP/P
,特异性为 TN/N
。获取正数 P
的总数很容易,因为我们的实际值是合乎逻辑的,我们可以只使用 sum(actuals)
。相似性,我们可以使用 sum(!actuals)
得到负数 N
。
pos <- sum(actuals)
neg <- sum(!actuals)
首先让我们在每个阈值处得到我们的真底片。这很简单,它只是等于或低于每个阈值的 FALSE
值的数量。由于我们的数据按阈值排序,我们可以使用以下方法计算(和特异性):
tn <- cumsum(!actuals)
spec <- tn/neg
真正的阳性数量稍微难一些,因为我们正在寻找大于阈值的阳性数量,所以单独使用 cumsum
是行不通的。然而,由于高于阈值的数量等于低于或等于阈值的总负数,我们可以使用以下方法获得真阳性:
tp <- pos - cumsum(actuals)
sens <- tp/pos
现在我们需要做的就是绘制两者。
plot(1 - spec, sens, type = "l", col = "red",
ylab = "Sensitivity", xlab = "1 - Specificity")
abline(c(0,0),c(1,1))
要获得曲线的 AUC,我们只需计算每个实际值的曲线高度(灵敏度)乘以宽度(1 的差异 - 特异性)。我们已经有了敏感性,我们只需要特异性。 diff
函数将为我们提供相邻值的特异性差异,但是,我们需要在开头放置一个 0
值以获得第一列的宽度。
width <- diff(c(0, 1 - sens))
auc <- sum(spec*width)
actuals <- actuals[order(scores)]
sens <- (sum(actuals) - cumsum(actuals))/sum(actuals)
spec <- cumsum(!actuals)/sum(!actuals)
plot(1 - spec, sens, type = "l", col = "red",
ylab = "Sensitivity", xlab = "1 - Specificity")
abline(c(0,0),c(1,1))
(auc <- sum(spec*diff(c(0, 1 - sens))))
关于r - 在不使用任何包的情况下在 R 中绘制 ROC 曲线,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/41754220/
我是 Java 新手,这是我的代码, if( a.name == b.name && a.displayname == b.displayname && a.linknam
在下面的场景中,我有一个 bool 值。根据结果,我调用完全相同的函数,唯一的区别是参数的数量。 var myBoolean = ... if (myBoolean) { retrieve
我是一名研究 C++ 的 C 开发人员: 我是否正确理解如果我抛出异常然后堆栈将展开直到找到第一个异常处理程序?是否可以在不展开的情况下在任何 throw 上打开调试器(即不离开声明它的范围或任何更高
在修复庞大代码库中的错误时,我观察到一个奇怪的情况,其中引用的动态类型从原始 Derived 类型更改为 Base 类型!我提供了最少的代码来解释问题: struct Base { // some
我正在尝试用 C# 扩展给定的代码,但由于缺乏编程经验,我有点陷入困境。 使用 Visual Studio 社区,我尝试通过控制台读出 CPU 核心温度。该代码使用开关/外壳来查找传感器的特定名称(即
这可能是一个哲学问题。 假设您正在向页面发出 AJAX 请求(这是使用 Prototype): new Ajax.Request('target.asp', { method:"post", pa
我有以下 HTML 代码,我无法在所有浏览器中正常工作: 我试图在移动到
我对 Swift 很陌生。我如何从 addPin 函数中检索注释并能够在我的 addLocation 操作 (buttonPressed) 中使用它。我正在尝试使用压力触摸在 map 上添加图钉,在两
我设置了一个详细 View ,我是否有几个 Nib 文件根据在 Root View Controller 的表中选择的项目来加载。 我发现,对于 Nibs 的类,永远不会调用 viewDidUnloa
我需要动态访问 json 文件并使用以下代码。在本例中,“bpicsel”和“temp”是变量。最终结果类似于“data[0].extit1” var title="data["+bpicsel+"]
我需要使用第三方 WCF 服务。我已经在我的证书存储中配置了所需的证书,但是在调用 WCF 服务时出现以下异常。 向 https://XXXX.com/AHSharedServices/Custome
在几个 SO 答案(1、2)中,建议如果存在冲突则不应触发 INSERT 触发器,ON CONFLICT DO NOTHING 在触发语句中。也许我理解错了,但在我的实验中似乎并非如此。 这是我的 S
如果进行修改,则会给出org.hibernate.NonUniqueObjectException。在我的 BidderBO 类(class)中 @Override @Transactional(pr
我使用 indexOf() 方法来精细地查找数组中的对象。 直到此刻我查了一些资料,发现代码应该无法正常工作。 我在reducer中尝试了上面的代码,它成功了 let tmp = state.find
假设我有以下表格: CREATE TABLE Game ( GameID INT UNSIGNED NOT NULL, GameType TINYINT UNSIGNED NOT NU
代码: Alamofire.request(URL(string: imageUrl)!).downloadProgress(closure: { (progress) in
我是一名优秀的程序员,十分优秀!