请求计算ROC曲线面积以及两两比较的stata程序及其菜单操作

我有如下数据，请高手写出 可求出 x  y  z  对转归的ROC曲线 面积以及 x    y    z 三条ROC曲线面积两两比较的stata程序（三条ROC曲线不交叉和曲线交叉的各自情况），以及满足前述要求的STATA菜单操作（三条ROC曲线不交叉和曲线交叉的各自情况）。谢谢高手指点！

      x            y            z       转归
  2.5800   110.20      .00        1
15.2600   133.79      .00        1
  2.2300   143.26     4.50        1
  5.6000   203.20     5.30        1
   .5130   206.03     3.10        1
  8.3000   209.40     5.80        1
  1.3100   215.50     4.40        1
   .5900   223.65     6.90        1
  4.5600   225.10     5.30        1
  2.5900   236.57     4.00        1
  7.1300   237.02     3.30        1
  3.4100   254.43     1.50        1
  3.5700   259.60     3.80        1
  3.8500   262.64     3.50        1
10.5600   279.56     2.80        0
  3.6800   282.00     4.80        1
  3.1900   286.72      .00        1
  5.8000   287.60      .00        1
   .2800   290.45     6.00        0
  1.6500   292.27     4.60        1
  3.5000   296.07     2.90        1
  6.5900   307.01     5.00        1
  1.4200   314.19     3.70        1
  1.0800   316.75     7.50        0
  1.3900   320.04     4.90        1
  1.0300   326.63     1.70        1
13.5500   343.30     4.70        1
  8.6300   353.52     6.50        1
  2.7000   363.97     4.90        1
  3.7800   375.38     4.30        1
  3.9000   385.56      .00        1
  5.2700   405.35     5.50        1
  1.0300   420.56      .00        1
  3.0200   421.91     6.00        1
  2.0800   430.27     7.60        0
  6.0200   434.94     6.40        1
  2.9000   435.83     3.20        1
  3.7000   436.70      .00        1
  8.9800   447.89     3.00        1
  8.9800   447.89     4.70        1
  8.0200   454.84     2.20        1
  9.0200   482.91     3.40        1
  1.3800   493.71      .00        1
  2.4000   510.92     4.10        1
  1.0400   512.36     2.70        1
24.8500   516.81     6.20        1
  2.7000   523.19     4.20        1
  6.7800   529.00     2.00        1
  7.5100   531.08     3.60        1
  8.5000   535.20     5.10        1
  3.6000   536.81     3.20        1
  4.5000   544.21     4.90        1
  1.2400   549.24      .00        1
  7.3000   556.00     6.30        1
  5.3000   586.96     2.90        1
11.8300   598.00     4.30        1
   .3900   619.67     3.90        1
  4.1700   619.90     2.00        1
  7.1500   619.90     3.70        1
  4.7900   624.27     6.30        1
  9.2000   635.30     5.10        1
  8.0200   636.14     3.40        1
  1.7500   643.16     4.20        1
  1.2000   650.87     4.80        1
  4.2200   666.42     2.60        0
  3.4600   687.34     7.10        1
  1.8700   693.26     4.30        1
  8.8300   697.08     3.20        1
38.4300   711.40     3.40        1
16.8900   716.53     4.00        1
  1.7000   731.20     5.90        1
  2.5000   756.40     4.80        0
  5.7000   779.51     4.90        1
  1.0800   783.39     2.80        1
  7.8000   792.54     5.80        0
  6.0900   792.88     6.20        1
  7.5000   800.61     1.80        1
  3.3800   809.14     4.00        1
  6.7000   810.60     4.50        1
  2.5200   813.35     1.20        1
12.6200   821.40     3.00        1
10.0000   823.10     6.60        1
  3.3700   829.62     7.00        1
  3.9500   853.61     2.70        1
11.5300   860.83     4.10        1
  4.1500   871.25     5.90        0
  9.4500   895.55     5.10        1
  5.6000   906.63     6.30        1
10.6400   908.98     1.80        1
  1.8200   910.70      .00        1
  2.6300   919.12     3.00        1
  5.4500   926.40     4.80        1
15.6000   931.20     3.10        0
  1.2700   931.91     4.30        1
  1.3000   975.38      .00        1
  6.9600   995.63     5.30        1
21.9600  1040.99     4.10        0
  6.8700  1074.00     5.00        1
  6.7900  1084.00     2.00        1
  7.5000  1084.00     6.80        1
  8.3000  1104.89     7.40        1
  5.3000  1106.59     5.60        1
  4.1700  1124.31     4.20        1
  2.9200  1124.75     5.80        1
  6.8900  1161.30     4.40        0
  5.3200  1170.00     5.80        1
  4.8700  1195.50     3.70        1
  3.7200  1220.50     3.60        1
  6.2500  1239.00     2.00        1
  2.0500  1262.00     5.60        1
  4.8400  1274.00     3.90        1
  9.7200  1278.92     4.00        1
  2.4700  1343.00     5.20        1
  3.1000  1411.30     5.30        1
  9.2500  1457.27     5.40        0
  8.4500  1491.68     4.90        0
  6.7800  1497.65     5.40        1
  3.6700  1514.31     5.00        1
  7.5600  1528.00     3.90        0
  2.8600  1542.65     2.90        1
19.3000  1553.95     5.80        1
  6.0200  1582.32     2.90        1
20.9200  1618.00     2.90        1
  1.5700  1754.00     3.00        1
   .7500  1799.64     3.00        1
  8.5000  1916.37     6.40        0
10.1100  1929.06     7.80        0
  6.2600  1931.66     5.20        1
  4.5600  1939.00     4.70        0
  2.2300  1954.43      .00        1
  1.5700  2016.57     4.40        1
  8.2300  2061.54     2.20        1
  6.7800  2165.52     5.50        0
  2.2800  2207.15     2.50        1
11.2500  2218.21     6.20        1
  7.9500  2265.61     4.50        0
  3.0500  2284.00     5.70        1
  8.8000  2321.98     2.90        1
20.3900  2329.60     6.30        1
  9.2000  2341.66     6.50        1
  2.1800  2367.19     6.90        0
  2.9500  2376.34     3.40        1
  3.9300  2377.23     6.60        1
  3.5600  2456.76     2.20        1
  8.0600  2616.63     7.20        0
  8.0600  2625.57     2.00        1
  3.9300  2691.76     2.00        1
31.3200  2731.50     4.10        1
  2.4800  2867.06     8.70        1
12.2500  3024.00     4.60        0
12.2500  3041.59     6.70        1
   .2400  3059.00     3.00        1
18.6700  3258.00     4.90        0
17.0000  3429.52     3.90        0
16.1800  3438.00     3.90        0
  6.9200  3805.26     6.50        1
18.4500  3841.21     3.90        0
18.6700  3923.00     2.30        0
  5.1100  3929.92     4.30        1
  9.1100  3932.21     2.60        1
  3.9900  4149.90     5.80        0
20.0800  4400.35     6.80        0
  7.5800  4552.60     3.90        1
  1.0400  4617.81     4.30        1
18.0400  4751.00     4.30        0
  3.2000  4905.00     5.20        1
  9.0500  5547.39     6.00        1
11.8400  5706.72     9.00        0
24.8300  5749.03     3.80        0
  4.2900  5914.44     3.10        0
22.7100  6421.22     4.80        0

在对ROC曲线面积比较时，不需要分不交叉和交叉的，因为比较的是面积大小。

stata 本身有计算器功能，你可以 在上面 的程序 后面加上 一个 p值乘以3 的公式，这样就 自动输出了调整后的P值。

或者你把 显著性 水平 0.05 除以 3, 得到的 p 值 直接 和 调整后 的 显著 性 水平 比， 这样也可以。

sas里面的两两比较也是要对p-value进行校正的， 不过它的矫正公式在程序里面，不用你手工校正了。

但是道理 都是一样的 。

我今天刚刚下载stata，并且刚刚学会用stata的菜单操作进行几条ROC曲线（不相交）面积的计算以及两两比较，如果有高手，可以用stata程序或者菜单操作进行几条相交的ROC曲线面积的比较，我就把他的答案设为最佳答案（按照回答的顺序）！谢谢高手们！

1 先把数据输入
2任一个自变量为 var1 ,给结局变量起名为var2

.logit var2 var1

Iteration 0:
log likelihood =
-117.336

Iteration 1:
log likelihood =
-116.0535

Iteration 2:
log likelihood = -116.05053

Iteration 3:
log likelihood = -116.05053


Logistic regression
Number of obs
=
189

LR chi2(1)
=
2.57

Prob > chi2
=
0.1088
Log likelihood = -116.05053
Pseudo R2
=
0.0110

------------------------------------------------------------------------------

var2 |
Coef.
Std. Err.
z
P>|z|
[95% Conf. Interval]
-------------+----------------------------------------------------------------

var1 |
-.0542063
.0340834
-1.59
0.112
-.1210085
.0125959

_cons |
.7062236
.9465123
0.75
0.456
-1.148906
2.561354
------------------------------------------------------------------------------
. lsens


. lroc


Logistic model for var2
number of observations =
189
area under ROC curve
=
0.5698

有些乱，我把结果整理成了PDF，见附件。

需要说明一下，PDF中蓝色标注的是命令，一共用了三个命令

.logit Y X
在你的题目中 转归是Y， X是自变量

.lroc
会绘制ROC曲线，给出曲线下面积

.lsens
提供敏感度，特异度 分别对概率分割点的制图

PDF中其它的都是stata自己输出的结果。

你像这样做三次，每一次放入一个自变量（也就是你题目中的X，Y，Z 分别放入自变量就好了）

另外强调一下，我用的数据和你不同，但是做法是一样的。

你好wangxinbio26！首先谢谢你，我发现你专门在我的疑难贴上回答我的难题！希望你以后多帮我解决难题！如果你愿意，悬赏价格都可以商量——当然对所有愿意帮助我的高手，悬赏价格都可以商量！！第二，你给我的程序只是画出ROC曲线和算出ROC曲线的面积，并没有对两条ROC曲线进行比较，当然也没有对有交叉的ROC曲线进行两两比较—这是我最需要的，也可以是我的唯一要求（对交叉的ROC曲线面积进行两两比较，据胡良平说，可以用rockit0.9β解决）！

roccomp y xb1 xb2, graph summary
                              ROC                    -Asymptotic Normal--
                   Obs       Area     Std. Err.      [95% Conf. Interval]
-------------------------------------------------------------------------
xb1                200     0.5785       0.0254        0.52870     0.62828
xb2                200     0.8330       0.0305        0.77329     0.89274
-------------------------------------------------------------------------
Ho: area(xb1) = area(xb2)
    chi2(1) =    50.18       Prob>chi2 =   0.0000

说明：面积的比较用的程序如下：

logit  y  x1
predict xb1
logit  y  x2
predict xb2
roccomp  y  xb1  xb2, graph summary



y 是两分类的结局变量， x1和x2 分别是两个logisitic 回归的自变量，xb1 和 xb2是两个logisitic 回归的预测值

然后用 roccomp 命令得出ROC面积值，及面积的 标准误 ， 置信区间 及 两个面积比较的 P-value .  那个graph summary

命令也可以不用，用的话会多一张图(见附件）

用的是上面粗体的命令， 今天看书时正好看到，就直接回复个帖子过来。再强调一下，数据和上面不同。

另外应该没有胡良平说的那么复杂。

最后要说的是：三组的两两比较，最后还应当对P-value 进行校正。

你好wangxinbio26！你所说的是两条ROC曲线比较，如多条比较，则最后应当对P值进行校正，我想请问你，可否编一个stata程序，可直接对几条ROC曲线进行两两比较，就像SAS程序一下，之后我们不需要对P值进行校正！最后还问一下，在对ROC曲线面积比较时，真的需要向胡良平所说的那样，需要分不交叉和交叉的情况吗？谢谢！