全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4497 1
2014-11-18
悬赏 1 个论坛币 未解决
通过LOGIT 命令可以得到回归,然后拟合ROC曲线,通过LSENS可以画出敏感度和特异度交会的曲线图,这是可以理解的。

寻找的CUTOFF通常是医学诊断上一个重要目的,CUTOFF即为切点最大值的约登指数。

LSENS显示了交会图,但是却不显示CUTOFF,有没有选项能够显示呢?

如果不行是否需要编程?

请大家百忙之中不吝赐教!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2014-11-18 10:34:51
不好意思,本例答案已经解决。

思路和命令如下和大家分享:
因为我只想看某指标与CASE的关系,所以只需要用ROCTAB即可。

格式如下和范例:

  . webuse hanley
        . roctab disease rating
        . roctab disease rating, graph
        . roctab disease rating, graph summary
        . roctab disease rating [fw=pop]
        . roctab disease rating, table detail
        . roctab disease rating, lorenz
        . roctab disease rating, lorenz graph

在这里我们可以用DETAIL这个选项。他会给出如下的结果:Detailed report of Sensitivity and Specificity
------------------------------------------------------------------------------
                                                        Correctly
Cutpoint      Sensitivity   Specificity     Classified          LR+          LR-
------------------------------------------------------------------------------
( >= .56 )        100.00%         0.00%       46.56%       1.0000     
( >= .78 )         99.45%         0.00%       46.31%       0.9945     
( >= .81 )         98.91%         0.00%       46.06%       0.9891     
( >= .82 )         98.91%         0.48%       46.31%       0.9938       2.2951
( >= .83 )         98.91%         1.43%       46.82%       1.0034       0.7650
( >= .9 )          98.91%         1.90%       47.07%       1.0083       0.5738
( >= .91 )         98.36%         2.38%       47.07%       1.0076       0.6885
( >= .94 )         97.81%         2.38%       46.82%       1.0020       0.9180
( >= .95 )         97.81%         3.33%       47.33%       1.0119       0.6557
( >= .96 )         97.27%         3.33%       47.07%       1.0062       0.8197

。。。。。

在CORRECTLY CLASSIFIED里面会LIST出最大的VALUE,这时候就是CUTOFF点。
大家不信可以用EXCEL进行验证,用函数:YOUDENINDEX=SENSIT+SPECIF-1,然后SORT。找出那个点。
我验证过,结果一致。

但是有个问题,这时候的灵敏度和特异度均不高,而会弱化诊断的价值。这是个矛盾。
这个现象的出现可能出本身的ROC AERA不高有关,点估计只有0.58。

                      ROC                    -Asymptotic Normal--
           Obs       Area     Std. Err.      [95% Conf. Interval]
         --------------------------------------------------------
           XXX     0.5843       0.0288        0.52790     0.64067



对于临床研究,我希望能够有较好(例如大于85%)的灵敏度,筛查出阳性患者。所以也提示:进一步寻求联合的分子靶标,有利于提高实验阳性诊断的。

希望对大家有启示。

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群