问题标签 [rstatix]

所以，我正在计算两个值之间的统计数据，一个是随机值，另一个是“真实”分数。因为它是成对的，所以我绘制了不同条件和状态的差异分布，facet_grid但我也对边缘（全局分布，即 margin = T）感兴趣。现在我为 t 检验添加置信水平（顺便说一句，这有意义吗？），因为它们对于每个条件都是正态分布的。尽管如此，我不希望水平条挤在边缘分布(all)中。可能吗？最好只添加一个水平条来指示全局分布的置信限。

产生问题的代码以及输出图如下：

这真让我抓狂。我正在使用 rstatix 的 anova_test ，它告诉我我的列不存在，而它们显然存在。

这是我的数据框的样子：

ID、Form、Pen 是因子，Time 是数字。所以每个受试者用红色、蓝色和绿色的笔完成了表格 A 和 B，我测量了每个人完成表格的时间。

这是我故意提出来提出这个问题的一个假数据集。实际上，此数据框源自具有更多变量的更大数据集。每个变量都有更多的观察值（因此，在本例中，主题 1 和表格 A 和红笔不仅仅是一个数据点，而是多个数据点），所以我总结了平均时间。

我想测试主要和交互效果，所以我正在做一个 2 x 3 重复测量 ANOVA（双向 ANOVA，因为 Form 和 Pen 是两个自变量）。

我不断收到这个错误：

为什么？！任何帮助将不胜感激。我一直在寻找 HOURS 的解决方案，但我感到非常沮丧。

考虑以下数据框：

按样本分组，我希望对每种细菌进​​行逐行双边 Fisher 精确检验。（例如 HM_001 如下所示）。

依此类推，本质上为数据帧中的 1350 行中的每一行生成一个 p 值。

以下是我的尝试：

但它不起作用，输出以下错误消息：

任何指针将不胜感激！

我正在尝试从 4 路混合方差分析中计算组均值。代码运行良好——我可以将结果输出到我的控制台中——但是如何将此表保存到文件中（即导出为 .csv）？我尝试使用该capture.output()命令，但找不到对象分数。

我有一个调查，在学期开始和结束时（前和后）在 3 个不同的课程（数学、物理、生物）中进行。在调查中，有3组问题（A、B、C）和李克特式量表。我将所有答案转换为数字分数

我想针对每种问题类型测试每门课程的分数是否存在差异。我还想在此处添加 Bonferroni 校正以考虑多重比较：

但是，p.adjust.method = "bonferroni"不产生任何输出。我不确定为什么。

当我尝试执行 leven_test() 时，我以前没有问题地执行过，现在我收到以下错误：

我尝试更改公式中的术语，并调用库（统计），但我不知道可能发生什么。有人可以帮忙吗？

所以我正在寻求创造这样的东西：

示例输出

使用我自己的数据，特别是使用我在这里找到的 p 值：

KS 检验 p 值

现在，我能够用不正确的方法制作类似的东西。具体来说，我能够使用 T 检验产生类似的东西： T 检验 p 值

我通过编写以下代码产生了这个：

换句话说，我使用 stat_compare_means() 使用 ggplot2/tidyverse/ggpubr/rstatix。

但是，如果我修改代码中的方法，它似乎可以正确显示 Wilcoxon 和 T 检验，但不能正确显示 anova 和 kruskal wallis 检验。此外，似乎 stat_compare_means() 仅支持这四个而不支持 KS，但我特别感兴趣的是从我的 KS 测试输出中将平均值比较绘制到我的小提琴图上。还有其他我可以使用的包吗？

另请注意：对于 KS 测试，“UpScorePos”“DownScorePos”等是按规则比较 ARE 分数（就像我对 T 测试中的图表所做的那样）。

我试图在 R 中的 3x12 ANOVA 之后可视化 Tukey HSD 成对比较的主要影响和交互。我已经能够使用类似于下面的代码添加第二个（12 级）因子的成对比较。（下面是一个小得多的数据子集，让你们的生活更轻松。）

第二个因素的主效应示例

但是，当尝试以类似方式添加交互时，它们不会正确映射到图上。如果有某种方式我应该改变测试和/或情节以使它们彼此“交谈”，或者如果有更好的方法来解决这个问题，我会很高兴听到它。

交互不会根据需要映射到绘图上

注意：我知道建议tukeyfull %>% add_xy_position(y="ftwo")在添加测试的 p 值之前执行类似的操作，但是当我尝试运行该示例时出现此错误：Error: Must group by variables found in '.data'. * Column 'fone:ftwo' is not found.

我正在尝试使用 RStatix's 执行一系列 T 检验t_test()，其中因变量在每个测试中都是相同的，并且分组变量会发生变化。我在一个循环中进行这些测试，所以我想用列号而不是列名来选择分组变量。我尝试使用 执行此操作colnames(dataframe)[[columnnumber]]，但出现以下错误：“无法提取不存在的列”。如何使用列号而不是列名选择分组变量？

下面是一个带有虚构数据框的最小可复制示例；当指示分组变量的名称（性别）时，测试正常工作，但在指示列号时不能正常工作。

我有一个不平衡的重复测量设计，我想为每个时间段（即曲线）运行单独的方差分析，然后 Bonferroni 纠正结果。

这是数据，其中曲线是重复测量：

我遇到的这个问题与 'rstatix' 包和 anova_test() 函数有关。它对 Topt 变量运行得很好。

这给出了：

但是，相同的代码对 Aopt 变量给出了一个奇怪的结果，其中 DFn 和 DFd 对于 Curve = "Fifth" 不正确。

这给出了：

有任何想法吗？谢谢