不好意思,本例答案已经解决。
思路和命令如下和大家分享:
因为我只想看某指标与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%)的灵敏度,筛查出阳性患者。所以也提示:进一步寻求联合的分子靶标,有利于提高实验阳性诊断的。
希望对大家有启示。