问题标签 [kruskal-wallis]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
2 回答
3033 浏览

r - R - stat_compare_means 从 Kruskal-Wallis 测试返回不同的值

我想使用packageggplot中的 R 函数绘制 Kruskal-Wallis 测试的 p 值。stat_compare_meansggpubr

但是,如果我只是运行该函数,则绘制的值与该值不同:

kruskal.test(value ~ type, data = Profile_melt)

我绘制 p 值的代码是:

这是我的数据样本数据

0 投票
1 回答
129 浏览

python - 我的 Kruskal-Wallis 课程出了什么问题?

我试图建立一个可以执行 Kruskal-Wallis 测试的类。该类使用以下公式计算 H:

计算 H 的公式

但是,它产生的 H 值与kruskalscipy 的函数不同。有谁知道为什么会这样?

0 投票
1 回答
66 浏览

r - Kruskal.wallis 给出相等的 p 值

朋友们,

我在 r 中的 Kruskal wallis 测试存在问题,使用 Kruskal-wallis 测试测试稳定的季节性。为每个变量测试的 p 值是相同的。Kruskal.test(formula, data = mydata)library(stats)包装中使用。我很难相信 pvalues 会是相同的。

我的数据集是一个月度数据集,包含 163 个 obs、模型中的 3 个宏观经济变量和两个季节性虚拟变量。

我正在以下列方式用因变量测试每个独立的宏观经济变量Kruskal.test(y~x, data = mydata)。因此,对于下面的数据示例,它将是Kruskal.test(pr~mev06_mp_lag2, data = mydata). 并为数据集中的每个 mev 重复。用于测试 3 mev ( mev06_mp_lag2, mev29_lag2, mev108_lag1) 的所有 pvalues 都是这个输出: data: pr by mev29_lag2 Kruskal-Wallis chi-squared = 162, df = 162, p-value = 0.4852

这是数据:

我错过了这个测试的细微差别吗?想法?

0 投票
1 回答
936 浏览

r - 如何在不更改 pvalue 的情况下为 kruskal wallis 测试执行循环,比较多列上的 3 个组?

我对 R 相当陌生,我正在尝试进行 kruskal wallis 测试,以查看在查看不同基因时三组之间是否存在差异。我有 3 组和 127 种蛋白质。我已经能够创建一个可以做到这一点的代码,

样本数据”

不幸的是,获得的 pvalue 与我一次只对一个基因进行 kruskal wallis 测试的结果不同。例如:

对于 Gene1,从循环中获得的 pvalue 是 0.0389,但是当我运行 kruskal.test(Gene1,group, data=df) 时,我得到的 pvalue 是 0.84。

我遇到了这个问题,因为在进行 kruskal wallist 测试后,我进行了成对的 Mann Whitney 测试,并注意到 Kruskal wallis 的“显着”pvalues 与 Mann Whitney 的“显着”pvalues 无关。

此外,我继续使用 VassarStats 和 minitab,得到的 p 值为 0.84(平局调整)。我想知道如何在不影响 p 值的情况下循环运行这个 Kruskal wallis 测试。有什么我没有看到我做错了吗?

此外,我使用了我在上一篇文章中看到的 getAnywhere(kruskal.test.default),但我找不到在一遍又一遍地执行测试时会导致这种情况发生的原因。

0 投票
0 回答
210 浏览

r - 排名上的方差分析 VS kruskal wallis,有多大不同

我不确定这是解决此类问题的理想场所,但也许您可以帮助我。

我想检查 3 次治疗之间定量变量的差异,即执行方差分析。
不幸的是,我的模型的残差不是正态分布的。

我通常有两种解决方案:转换我的数据或使用我的测试的非参数等效项(这里是 kruskal wallis 等级测试)。

我尝试过的所有转换都无法满足常态(log、1/x、平方根、tukey 和 boxcox 幂),所以我想使用 kruskal 并继续前进。

然而,我的项目经理坚持只使用 ANOVA,并谈到 ANOVA 作为一个神奇的解决方案。

在 RI 上工作查找了一些示例,并artARTool包中找到了一个在等级上执行 anova 的函数。

基本上它需要你的变量并用它的等级替换它(通过平均等级来处理平局):

给出完全相同的输出。

我不是专业的统计学家,我想知道这个解决方案/转换有多有效。对我来说,这似乎很残酷,并且与 kruskal-wallis 测试的逻辑相距不远,即使统计数据不是直接根据等级计算的。

我发现有一个与 kruskal-wallis 不同的“等级方差分析”测试非常令人困惑(也称为等级单向方差分析),我不知道如何在这两个测试之间进行选择。

我不知道我是否说得很清楚,是否有人可以帮助我,但无论如何,感谢您的关注和评论!

PS:这是一个关于虚拟数据的例子

0 投票
1 回答
377 浏览

r - 如何在 R 中获得 Kruskal-Wallis 检验的精确 p 值?

如何在 R 中获得 Kruskal-Wallis(例如 3 组)测试的精确 p 值?

数据示例:

我尝试了coin使用该功能的包kruskal_test

虽然产生了错误:

如果我改变它的"exact"运行"approximate",但它不是确切的分布......

有什么想法吗?

0 投票
1 回答
181 浏览

boxplot - 箱线图坏了,只显示一行

所以我的数据集中在不同的处理方法以及它们如何影响发芽日。 狡猾的箱线图数据的图像

不久前,在 R 中制作小提琴图以显示根据处理发生发芽时间的分布时,我试图添加一个箱线图作为描述性统计数据,但只遇到了一条线。

我联系了许多根本不知道问题出在哪里的人,我在另一个小提琴图中使用了相同的数据,作为更大数据收集的一部分,包括这个数据在内的更多处理。

我继续前进,发现它很奇怪,现在当我开始在 SPSS 中执行统计测试时,我遇到了与下图相同的问题。当我尝试 Mann Whitney U 检验时,我被告知“无法计算”,因为没有两个变量超过 3 个测试字段(即 2 个)。

我完全不知所措,我尝试过重写数据,使用“分层”“strat”“s”等复制数据标签,但我不知道问题出在哪里,如果有人能给我任何指导的话非常感谢!

谢谢

0 投票
1 回答
77 浏览

r - 在函数中指定变量的准则

考虑dat1在这里创建:

dat1有 5 个变量的测量值,观察值 ( IDs) 可以根据 3 个分组变量进行 分组Loc:并保持我的分析整洁。我正在使用该包执行多项操作。以下函数将对我指定的数据进行 Kruskal Wallis 测试,计算效果大小并在单个数据框中返回结果:StateRegionrstatixefszres

现在我可以使用该函数循环遍历每个响应变量并Region在单个数据框中获取分组变量的结果(示例显示它),这就是我需要的:

这对 Kruskal Wallis 测试很有用,现在我想做一个非常相似的函数来做一个 duns 测试,但是看看会发生什么:

dunn_test()在用户编写的函数之外,您可以kruskal_test()以完全相同的方式指定数据。那么在这两个函数中指定变量有什么区别,为什么第一个起作用而第二个不起作用?

0 投票
0 回答
40 浏览

r - 问题:两组的 Kruskal-Wallis 检验

我使用聚类工具将样本分成两组。接下来,我想看看它们在临床特征方面的统计差异:淋巴结(连续变量)和肿瘤分期(分类变量)

现在,我对上述两个临床特征使用了 Kruskal Wallis 检验,而不是用于淋巴结的 Wilcoson 和秩检验和用于肿瘤分期的 Pearson 卡方检验。理论上能接受吗?

0 投票
1 回答
229 浏览

r - 在我的函数中应用巴特利特测试时出错

我是 R 的新手,刚开始在 R 中编写函数。我的目的是创建一个函数来改进检查变异的过程并进一步进行 ANOVA 或 Kruskal wallis 检验,还计算显着项目的平均值或中位数。

将变量放入函数后,出现错误

我已经查了一段时间了,complete.case() 在我的变量中都是 TRUE 的。我无法弄清楚我的功能有什么问题,也不确定我的功能的其他部分是否可以解决......希望你们能帮助我,谢谢!