请选择 进入手机版 | 继续访问电脑版
楼主: kanzongxuan
3383 3

[问答] 求助:如何提取置换检验中的P值 [推广有奖]

  • 0关注
  • 0粉丝

大专生

5%

还不是VIP/贵宾

-

威望
0
论坛币
283 个
通用积分
2.3358
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
412 点
帖子
22
精华
0
在线时间
42 小时
注册时间
2015-6-8
最后登录
2023-6-6

kanzongxuan 发表于 2019-4-16 17:29:19 |显示全部楼层 |坛友微信交流群
20论坛币
请教大神:
在《R语言实战》中提到可用coin包中的oneway_test函数做置换检验,但怎么提取出里面的P值来?


library(coin)
score <- c(40, 57, 45, 55, 58, 57, 64, 55, 62, 65)
treatment <- factor(c(rep("A", 5), rep("B", 5)))
mydata <- data.frame(treatment, score)
t.test(score ~ treatment, data = mydata, var.equal = TRUE)
fit<-oneway_test(score ~ treatment, data = mydata, distribution = "exact")
fit


        Exact Two-Sample Fisher-Pitman Permutation Test

data:  score by treatment (A, B)
Z = -1.9147, p-value = 0.07143
alternative hypothesis: true mu is not equal to 0



另外,这里的Z值为-1.9147,但根据Z值来求P值的话:pnorm(-1.9147)*2=0.0555,与上面的P值也不同,但与下面的置换检验(参数 distribution取“asymptotic”)中P值相同,想问一下:为何两个不同的参数,得到同一个Z值,但得到不同的P值,谢谢!
> oneway_test(score ~ treatment, data = mydata, distribution = "asymptotic")

        Asymptotic Two-Sample Fisher-Pitman Permutation Test

data:  score by treatment (A, B)
Z = -1.9147, p-value = 0.05553
alternative hypothesis: true mu is not equal to 0

关键词:置换检验
kanzongxuan 发表于 2019-4-17 08:54:55 |显示全部楼层 |坛友微信交流群
自己查了下coin的文档,已经解决了。使用pvalue函数即可。
pvalue(fit)

使用道具

zzz1999x 发表于 2019-4-23 14:58:26 |显示全部楼层 |坛友微信交流群
我也有同样的问题,我没有找到coin这个包啊?

使用道具

卧雨禅 发表于 2020-2-24 16:01:25 |显示全部楼层 |坛友微信交流群
你好,,那个关于置换检验里的p值如何提取请问你解决了吗

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 04:49