带有调整 p 值的 R 的 Kruskal-Wallis 检验
Posted
技术标签:
【中文标题】带有调整 p 值的 R 的 Kruskal-Wallis 检验【英文标题】:Kruskal-Wallis Test on R with ajusted p-value 【发布时间】:2019-06-14 23:52:11 【问题描述】:我正在对我的数据集执行 Kruskal-Wallis 测试,我正在尝试调整 p 值,但它似乎不起作用,这是我的代码:
> kruskal.test(df$Folate_biosynthesis, df$Group, p.adj="holm")
Kruskal-Wallis rank sum test
data: df$Folate_biosynthesis and df$Group Kruskal-Wallis chi-squared
= 8.5144, df = 5, p-value = 0.1301
> kruskal.test(df$Folate_biosynthesis, df$Group, p.adj="none")
Kruskal-Wallis rank sum test
data: df$Folate_biosynthesis and df$Group Kruskal-Wallis chi-squared
= 8.5144, df = 5, p-value = 0.1301
如您所见,如果我输入p.adjust = "none"
,我得到的结果完全相同。怎么可能?
在此先感谢所有愿意提供帮助的人。
安德烈亚
【问题讨论】:
P 值调整通常在您有多个 p 值时进行。您只有一个 p 值,所以我想知道您希望调整在这里做什么。 【参考方案1】:当您有多个 p 值时,通常会进行 P 值调整。您只有一个 p 值,所以我想知道您希望在这里进行什么调整。
话虽如此,kruskal.test
也没有出现p.adj
参数。该函数有一个点参数,但据我所知,它不会使用或将这些参数传递给任何其他函数,因此任何不是命名参数的输入基本上都将被忽略。
如果您想调整来自 kruskal.test
的多个输出的 p 值,您可以将 p 值收集到一个向量中,然后使用适当的方法将它们直接传递给 p.adjust
。
但话虽如此,尚不清楚您希望实现什么 - 但很明显,尝试在 kruskal.test
中使用 p.adj
参数并不是实现目标的方法。
【讨论】:
以上是关于带有调整 p 值的 R 的 Kruskal-Wallis 检验的主要内容,如果未能解决你的问题,请参考以下文章
R语言mgcv包中的gam函数拟合广义加性模型(Generalized Additive Model)GAM(对非线性变量进行样条处理计算RMSER方调整R方可视化模型预测值与真实值的曲线)