在R中使用后期时获得未格式化的p值

Question

我目前正在分析一个数据集，其中的组是异方差的。因此，我使用韦尔奇变异数来检验不同收入水平之间的显著差异。这些结果是重要的，因此，我希望执行一个后特别游戏-豪厄尔测试，以确定哪些组之间的差异。我使用的是posthocTGH函数，它是对用户友好的科学包的一部分。这根据APA标准格式化p值，但我想知道未格式化的值是什么。我尝试使用formatPvalue = FALSE参数来修正这个问题，但是这不会改变显示的p值。

目前的代码：

posthoc<-posthocTGH(y=dataset$climb.distance, 
           x=dataset$genotype, 
           method=c("games-howell"),
           formatPvalue = FALSE)

Answer 1

如果我正确理解了您的问题，您可以得到如下未格式化的p值：

phoc <- posthocTGH(y=ChickWeight$weight, x=ChickWeight$Diet, method="games-howell")
(outmtx <- phoc$output$games.howell)

         diff       ci.lo    ci.hi         t       df            p   p.adjusted
2-1 19.971212   0.3575276 39.58490 2.6378339 201.3839 4.420348e-02 0.1768139320
3-1 40.304545  17.5439718 63.06512 4.5930381 175.9190 4.883308e-05 0.0002929985
4-1 32.617257  13.4519674 51.78255 4.4086175 203.1602 9.873989e-05 0.0004936994
3-2 20.333333  -6.2029177 46.86958 1.9829842 229.9425 1.974058e-01 0.5922172720
4-2 12.646045 -10.9092650 36.20136 1.3891079 235.8785 5.074188e-01 1.0000000000
4-3 -7.687288 -33.8987923 18.52422 0.7590752 226.1567 8.726938e-01 1.0000000000

可以使用以下方法提取p-值的向量：

outmtx$p.adjusted

[1] 0.1768139320 0.0002929985 0.0004936994 0.5922172720 1.0000000000 1.0000000000

我希望这能帮到你。