我想pairwise.wilcox.test一次申请多个自变量，然后想要长格式的输出。对于特定的波长，我可以使用以下代码来完成

我想要的最终输出是

现在一次将它应用于所有波长，我使用了dplyr包（最新版本 1.0.0），例如

返回我

现在如何以长格式输出，例如

数据

如果我使用插入符号的两种方法（NN 和 KNN）然后我想提供显着性检验，我该如何进行 wilcoxon 检验。

我提供了我的数据样本如下

如何在wilcox.test()这里表演。

我有一个记录在 2 个位置的 11 种植物的变量数据框。对于每个物种，我正在尝试使用 t.test（或 wilcoxon 测试）比较两个不同位置之间变量的平均值。

这是我的数据的前几行

例如，对于每个物种，我想比较（测试显着差异）豪登省和普马兰加省孢子密度的平均值。请问有什么帮助吗？

我想对跨多个组的两种治疗使用 Wilcoxon 2 面检验，即对几个样本位点中的每一个都有一个治疗前后 (Conc)。我想按站点将数据集拆分为一个列表，然后应用测试，这样我就可以分别为每个站点提供一个输出，但是，我无法将其设置为可以重复的函数。

我有多个站点 (Site) 和两个级别的治疗 (Scenario)，结果得分 (Conc)：

每个站点/场景组合中有多个 Conc 数据点 (~60)。我选择 Wilcoxon 测试的原因主要是因为我在每个站点的处理（场景）之间的样本数量略有不均匀。

当我将此代码用于整个数据集时，我得到了一个合理的结果：

但是，此代码不会单独为每个站点应用测试。

我查看了我能找到的所有类似示例（在 SO 和其他地方），这是我能想到的最好的代码：

这段代码为每个站点提供了一个输出，但所有测试输出都是相同的，即使是 p 值：

谁能看到我做错了什么？感谢您的帮助

为了找到 iris 数据集的两列之间的关系，我正在执行 kruskal.test 并且 p.value 显示了这两列之间的有意义的关系。

结果如下：

散点图也显示了某种关系。 plot(iris$Petal.Length, iris$Petal.Width)

为了找到这两个变量的有意义的边界，我运行了pairwise.wilcox.test测试，但是为了使这个测试起作用，其中一个变量需要是分类的。如果我将两个连续变量都传递给它，那么结果将不符合预期。

作为输出，我需要一个明确的切入点，这两个变量有某种关系以及这种关系在哪里结束（如上图红线所示）

我不确定是否有任何统计测试或其他编程技术来找到这些边界。

例如手动我可以做这样的事情来标记边界 -

但是，R 中有没有一种编程技术或库来找到这样的边界？

重要的是要注意我的实际数据是有偏差的。

谢谢， 索拉布

在 R 中，我试图创建一个具有各种 p 值的向量，如下面的代码所示，但wilcox.test无法确定四个 p 值之一，因为只有一个类别。如何使代码工作，以便在无法执行NA时生成带有 an 的向量？wilcox.test

我正在尝试在 R 中进行一些分析。我需要对我的数据框执行非配对 Wilcox 测试。首先我为分析做了一个脚本：

但是，我想通过将列 Obs1 更改为 Obs2 然后更改为 Obs3 来执行此分析的 for 循环，但不知道如何开始。

输入（df）

结构（列表（样本 = c（A，A，B，B，C，C，D，D，E，E，F，F，G，G，H，H，I，I），Var1 = c（" 1”、“2”、“1”、“2”、“1”、“2”、“1”、“2”、“1”、“2”、“1”、“2”、“1” , "2", "1", "2", "1", "2"), Var2 = c("2", "2", "3", "3", "1", "1", “1”、“1”、“2”、“2”、“2”、“2”、“1”、“1”、“3”、“3”、“3”、“3”）、Obs1 = c(3.12, 6.54, 4.14, 5.40, 8.56, 4.92, 3.36, 8.57, 7.56, 7.89. 2.34, 5.43, 6.21, 7.43, 5.83, 9.28, 10.30, 3.42), Obs2 = c(4.13, 2.6) 9.42、4.95、3.21、9.99、1.25、9.78、9.54、7.23、6.24、1.23、2.35、4.73、4.12, 5.30, 0.12), Obs3 = c(5.35, 1.11, 8.98, 6.72, 2.62, 3.97, 2.62, 6.74, 9.41, 5.37, 7.43, 9.62, 4.69, 4.27, 2.74, 3.53, 5.4 = 7) c(8.52, 9.59, 3.45, 9.40, 8.77, 4.26, 9.78, 5.55, 2.31, 5.12. 2.35, 4.33, 7.61, 5.37, 6.84, 9.98, 8.65, 0.43), row.names = c(NA, -6L) , 类 = c("tbl_df", "tbl", "data.frame"))

有人可以帮我解决循环吗？提前致谢。

我想使用基本 R 箱线图显示非参数成对比较的结果。我有多个组（物种）和多个变量（比率）。

我将提供一些虚拟数据，因为我的数据集非常大：

我的想法是使用基础 R 创建具有 3 行（速率 1-3）和 3 列（处理 1-3）的面板箱线图。然后，我想显示我的非参数统计检验的结果：Kruskal Wallace 处理之间的差异和事后 Wilcoxon 成对比较（字母指数）以显示不同处理中物种之间的差异。我还没有找到如何有效地显示这些字母。
我想显示显着比较的字母。

在此先感谢您的帮助！

我有一个数据集，可以测量来自多个样本站点的大型无脊椎动物的丰度。我希望将最近几年的抽样结果与同一地点所有前几年的抽样结果进行比较。

我的数据如下所示：

s_abundance1 代表最近站点的站点丰度，s_abundance2 代表先前采样站点的站点丰度

整个数据集大约有 4000 行，由许多不同流域的样本数据组成。

我想执行 mann-whitney u 测试，比较 s_abundance1 和 s_abundance2，但在单个输出中按盆地分组

我一直在使用的代码是：

它似乎有效，只是所有的 p 值都完全相同。这是输出：

我需要更改哪些内容才能为每个盆地获得不同的结果？

对于我的理学硕士 项目我尝试按中位数对（由用户）给定数据帧的所有列进行排序，通过特定模式（稍后提到）在列上应用 wilcox.test，然后在箱须图中绘制每列的值。

排序和绘图工作得很好，但我很难找到一种方法将 wilcox.test 应用于以下模式中的数据框：

而i=1andj=2和 都递增直到j=ncol(dataframe)。所以我想用参数第 1 列和第 2 列运行函数，然后用第 2 列和第 3 列等等，直到j是数据框的最后一列。

我也想将所有 p 值存储在一个数据框中，其中包含一行（包含 p 值），并且每一行都具有作为 wilcox.test 中参数的两列的名称，因为我不想只绘制所有列（每列都代表一个“解决方案”），但我也想在控制台中打印每个测试的 p 值（例如：“Wilcoxon-test with 'Solution1' and 'Solution2'导致 p 值：'来自解决方案 1 和解决方案 2 的 wilcox.test 的 p 值'，这意味着解决方案是/不是显着不同的“）。

我试图调整其他帖子中有关此问题的一些代码，但没有任何结果。不幸的是，我在 R 方面也非常缺乏经验，所以我希望我上面写的也不是完全胡说八道。我也尝试以 Java 方式使用 for 循环和增量迭代数据帧的列，因为这是我所学的唯一编程语言，但这根本不起作用（真令人惊讶）。

我的代码基于具有非常不同值的数据框创建的图：

谢谢大家能给我的任何建议，非常感谢！