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

12917

收藏 2012-07-30

请问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值，就悲催了……

谢谢哈~@

扫码加我拉你入群

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

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

全部回复

蓝色

2012-7-30 17:06:18

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

扫码加我拉你入群

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

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

carweed

2012-7-30 18:52:33

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

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

谢谢！

扫码加我拉你入群

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

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

蓝色

2012-7-30 18:57:54

. 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

扫码加我拉你入群

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

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

蓝色

2012-7-30 18:58:06

sysuse auto

ranksum price,by(foreign)

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

扫码加我拉你入群

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

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

蓝色

2012-7-30 18:59:22

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

扫码加我拉你入群

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

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

点击查看更多内容…

carweed

2012-7-30 18:59:27

难道是NORMSDIST（）？

扫码加我拉你入群

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

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

carweed

2012-7-30 19:02:08

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

谢谢斑竹~！！！

扫码加我拉你入群

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

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

carweed

2012-7-30 19:23:04

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

后来人可以按需要套用

扫码加我拉你入群

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

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

sungmoo

2012-7-31 13:47:46

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

扫码加我拉你入群

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

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

carweed

2012-8-1 09:46:08

sungmoo 发表于 2012-7-31 13:47
http://www.stata.com/statalist/archive/2004-12/msg00622.html

O(∩_∩)O谢谢~！
请问一下你是怎么搜索到这个页面的？（http://www.stata.com/statalist/archive/2004-12/msg00622.htm）
谢谢！

扫码加我拉你入群

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

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

mxcz

2013-7-1 15:17:21

求助楼主！我用Stata做Wann-Whitney检验的时候，总出错，请问你输入的命令是什么啊？

扫码加我拉你入群

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

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

alasaa

2017-5-30 09:33:45

sungmoo 发表于 2012-7-31 13:47
http://www.stata.com/statalist/archive/2004-12/msg00622.html

老师您好，这个链接打不开了...请问您还有原版的资料吗？

扫码加我拉你入群

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

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

huasheng121

2018-5-23 18:08:24

膜拜一下

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群