楼主: kanzongxuan
4038 3

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

  • 0关注
  • 0粉丝

大专生

6%

还不是VIP/贵宾

-

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

楼主
kanzongxuan 发表于 2019-4-16 17:29:19 |AI写论文
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
拉您进交流群
GMT+8, 2026-2-7 23:52