请问ranksum检验的p值怎么记录？难道stata不返回这个东西？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问ranksum检验的p值怎么记录？
Eg：
     webuse fuel2
     ranksum mpg, by(treat)

Two-sample Wilcoxon rank-sum (Mann-Whitney) test


       treat |      obs    rank sum    expected
-------------+---------------------------------
           0 |       12         128         150
           1 |       12         172         150
-------------+---------------------------------
    combined |       24         300         300


unadjusted variance      300.00
adjustment for ties       -4.04
                     ----------
adjusted variance        295.96


Ho: mpg(treat==0) = mpg(treat==1)
             z =  -1.279
    Prob > |z| =   0.2010  怎么能让stata记录下这个值好在后面的程序里调用？


看帮助文件发现返回值列表里只有Z值，在r(z)里，但是怎么没有p值呢？难道不返回p值啊？怎么能让stata把这个p值也记录吗？
因为我要搞好多个变量的秩和检验，要是不返回p值的话，我就得一个个手动复制p值，就悲催了……


谢谢哈~@

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：ranksum Ranks Stata tata Rank 记录 检验

有了z
可以自己用函数求p值

蓝色 发表于 2012-7-30 17:06
有了z
可以自己用函数求p值

真是惭愧 我不知道公式是什么。。。麻烦斑竹指点一下啦    谢谢！

. sysuse auto
(1978 Automobile Data)

.
. ranksum price,by(foreign)

Two-sample Wilcoxon rank-sum (Mann-Whitney) test

     foreign |      obs    rank sum    expected
-------------+---------------------------------
    Domestic |       52        1862        1950
     Foreign |       22         913         825
-------------+---------------------------------
    combined |       74        2775        2775

unadjusted variance     7150.00
adjustment for ties        0.00
                     ----------
adjusted variance       7150.00

Ho: price(foreign==Domestic) = price(foreign==Foreign)
             z =  -1.041
    Prob > |z| =   0.2980

.
. scalar P_value=2*normal(r(z))

. scalar list P_value
   P_value =  .29801025

sysuse auto

ranksum price,by(foreign)

scalar P_value=2*normal(r(z))
scalar list P_value

    normal(z)
       Domain:       -8e+307 to 8e+307
       Range:        0 to 1
       Description:  returns the cumulative standard normal distribution.

难道是NORMSDIST（）？

蓝色 发表于 2012-7-30 18:57
. sysuse auto
(1978 Automobile Data)

谢谢斑竹~！！！

一个通用的公式：
loc z = abs(r(z))
loc q = 2*(1-normal(`z'))
暂元q里的就是P值啦

后来人可以按需要套用

http://www.stata.com/statalist/archive/2004-12/msg00622.html