Stata:预测模型中的临床决策曲线-dca

发布时间:2021-11-17 阅读 3452

Stata连享会   主页 || 视频 || 推文 || 知乎 || Bilibili 站

温馨提示: 定期 清理浏览器缓存,可以获得最佳浏览体验。

New! lianxh 命令发布了:
随时搜索推文、Stata 资源。安装:
. ssc install lianxh
详情参见帮助文件 (有惊喜):
. help lianxh
连享会新命令:cnssc, ihelp, rdbalance, gitee, installpkg

课程详情 https://gitee.com/lianxh/Course

课程主页 https://gitee.com/lianxh/Course

⛳ Stata 系列推文:

PDF下载 - 推文合集

作者:Bethune0507
邮箱bethune0507@163.com


目录


1. 临床决策曲线分析

诊断试验中根据灵敏度和特异度可以做出 ROC 曲线,得到 AUC,来评价诊断试验的准确性。但数理统计上的准确并不代表患者一定会受益。比如下图所示,通过某观测值预测患者是否患了某病,无论选取哪个值为临界值,都会遇到假阳性和假阴性的可能。

但临床上有时候避免假阳性受益更大,有时候避免假阴性受益更大。所以希望找到一个净受益最大的方法。DCA 就应运而生了。纪念斯隆凯特琳癌症研究所的官方网站上设立了专门的网页,来介绍 DCA 的相关理论、发表文章以及软件实现方法。「网站」提供了 Stata、R 和 SAS 三种软件实现 DCA 曲线分析,并配有详细的教程。

2. 怎么解读 DCA 的含义

横坐标为阈概率 (ThresholdProbability):患者 i 诊断为某种疾病的概率记为 Pi;当 Pi 达某个阈值 (记为 Pt),就界定为阳性,采取治疗措施。此时会有病人治疗的获益 (利),也会有非病人治疗的伤害 (误诊) 以及病人未治疗的损失 (弊)。而纵坐标就是利减去弊之后的净获益 (Net Benefit,NB)。

灰色横线与黑色斜线代表两种极端情况。灰色横线代表所有患者全部没病,都没有治疗。这种情况下,没有治疗受益。黑色斜线代表所有患者全都诊断有病,因此全部接受治疗。黑色虚线代表模型的情况,一般来讲,曲线离两种极端情况越远越好。

3. Stata 实操

3.1 命令介绍

*命令安装
cnssc install lxhget, replace
lxhget dca.pkg, install replace
*命令语法
dca depvar indepvars [if] [in] [, options]
  • depvar 定义结局事件,有无糖尿病;
  • indepvars 定义预测概率。

3.2 具体示例

. webuse lbw.dta, clear
. logit low age lwt

Logistic regression                                     Number of obs =    189
                                                        LR chi2(2)    =   7.53
                                                        Prob > chi2   = 0.0231
Log likelihood = -113.56932                             Pseudo R2     = 0.0321
------------------------------------------------------------------------------
         low | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
         age |     -0.040      0.032    -1.23   0.218       -0.103       0.023
         lwt |     -0.013      0.006    -2.05   0.040       -0.025      -0.001
       _cons |      1.747      0.997     1.75   0.080       -0.207       3.701
------------------------------------------------------------------------------

. predict predlow, pr

. * Decision Curve Analysis
. dca low predlow, smooth xstop(0.5) lcolor(black gs8 black) ///
>     lpattern(solid solid dash) title("Decision Curve Analysis Example", ///
>     size(4) color(red)) scheme(s1mono)

4. 相关推文

Note:产生如下推文列表的 Stata 命令为:
lianxh probit logit, m
安装最新版 lianxh 命令:
ssc install lianxh, replace

相关课程

免费公开课

最新课程-直播课

专题 嘉宾 直播/回看视频
最新专题 文本分析、机器学习、效率专题、生存分析等
研究设计 连玉君 我的特斯拉-实证研究设计-幻灯片-
面板模型 连玉君 动态面板模型-幻灯片-
面板模型 连玉君 直击面板数据模型 [免费公开课,2小时]
  • Note: 部分课程的资料,PPT 等可以前往 连享会-直播课 主页查看,下载。

课程主页

课程主页

关于我们

  • Stata连享会 由中山大学连玉君老师团队创办,定期分享实证分析经验。
  • 连享会-主页知乎专栏,700+ 推文,实证分析不再抓狂。直播间 有很多视频课程,可以随时观看。
  • 公众号关键词搜索/回复 功能已经上线。大家可以在公众号左下角点击键盘图标,输入简要关键词,以便快速呈现历史推文,获取工具软件和数据下载。常见关键词:课程, 直播, 视频, 客服, 模型设定, 研究设计, stata, plus, 绘图, 编程, 面板, 论文重现, 可视化, RDD, DID, PSM, 合成控制法

连享会小程序:扫一扫,看推文,看视频……

扫码加入连享会微信群,提问交流更方便

✏ 连享会-常见问题解答:
https://gitee.com/lianxh/Course/wikis

New! lianxhsongbl 命令发布了:
随时搜索连享会推文、Stata 资源,安装命令如下:
. ssc install lianxh
使用详情参见帮助文件 (有惊喜):
. help lianxh