问题标签 [rstatix]

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 投票
1 回答
1266 浏览

r - R:从 conover 测试中绘制 p 值的问题(邓恩的作品很好)

我遇到了一些我无法解决的错误消息。

我使用以下代码在 R 中运行 Kruskal Wallis,然后使用 Dunn's 进行成对比较:

然后我使用这个结果来绘制出如下所示的重要性:

得到的情节很漂亮(万岁!)

在此处输入图像描述

但是封面测试显然比邓恩的更强大和更受欢迎......(也欢迎对此提出任何意见!)。运行以下代码会返回错误:

我无法弄清楚...我可以很好地运行 Conover 测试,但无法解决上述错误以使其与 ggboxplot 一起使用。

我实际上更喜欢使用更美观的 ggplot 和 geom_boxplot 但根本无法与 Dunn_test 交互......

任何解决方案都将受到欢迎!

谢谢

PS我在ggplot中使用它:http: //www.sthda.com/english/articles/24-ggpubr-publication-ready-plots/76-add-p-values-and-significance-levels-to-ggplots但是我不想要一个 wilcoxon 或 t.test 进行成对比较......

0 投票
3 回答
226 浏览

r - 如何获得 r 中多个组的平均 SD 和 Pvalue?

你好我所有的数据框看起来像

我想分别根据类和类别计算(test1和test2之间)的平均SD和Pvalue

我使用 dplyr 计算平均值和 SD 并努力计算 Pvalue,因为我的数据集包含 1000 行和 4 个不同的类别和 8 个类

这是使用 dplyr 表示 mean 和 sd 后得到的结果

而SD只是手动输入以显示它的错误值我想要的输出是

提前致谢

0 投票
1 回答
309 浏览

r - rstatix 包安装

我在安装“rstatix”包时遇到问题。

当安装完成并尝试运行库时,我收到以下消息:

错误:loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) 中“rstatix”的包或命名空间加载失败:命名空间“vctrs”正在加载 0.2.0,但需要 >= 0.3.2 另外:警告消息:包 'rstatix' 是在 R 版本 3.6.3 下构建的

你以前处理过这个问题吗?你是怎么解决的?谢谢!

0 投票
1 回答
38 浏览

anova - anova_test 在主题 ANOVA 中没有以三种方式返回 Mauchly

我正在使用一个名为 sleep 的数据集(在此处找到:https ://drive.google.com/file/d/15ZnsWtzbPpUBQN9qr-KZCnyX-0CYJHL5/view )在受试者 ANOVA 中运行三种方式,比较基于刺激、剥夺、和时间。在使用来自 rstatix 的 anova_test 之前,我已经成功地做到了这一点。我想查看球形输出,但它没有出现在输出中。我已经在主题数据集中想出了其他三种方式,所以我不确定为什么会这样。这是我的代码:

我还尝试将它保存到一个对象并使用 get_anova_table,但这看起来并没有什么不同。

这是我从互联网上提取的理想数据集,因此我开始认为数据的 W 为 1(完美的球形度),因此 rstatix 正在跳过此输出。这是 anova_test 做的吗?

这也是我使用返回 Mauchly 的数据集的代码:

0 投票
0 回答
165 浏览

r - 我可以将 anova_test (rstatix) 类的对象转换为 aov、anova 或 aovlist 类的对象吗?

重复测量方差分析的 DATANOVIA 指南中的一些样本数据:

我想使用report来自出版商 E​​asyStats ( https://github.com/easystats/report ) 的软件包。它接受 、 或 类的aovANOVAanova对象aovlist。ANOVA 中的对象rstatix属于类anova_test,因此与report.

有什么方法可以转换对象类,还是应该使用基 R 进行重复测量方差分析?

回复 1(已被删除):

class(res.aov) <- "aov"class(res.aov) <- c("aov", class(res.aov))。您可以添加或更改对象的类属性

所以我尝试了这个:

对象摘要

然后你的建议:

对象摘要 2 类更改似乎已经奏效,但将对象传递给report. 我该如何进行?

附带说明一下,如果有人可以澄清这些是否等效,我的问题可以通过将我的 ANOVA 传递给基础 R 来解决:

anova_test(response ~ A * B + Error(subject/(A * B)))rstatix

相当于

aov(response ~ A * B + Error(subject/(A * B)))stats

0 投票
1 回答
581 浏览

r - 使用 rstatix 在 ggplots 上的 p 值位置不正确

我无法使用 rstatix 将 p 值放置在 ggplot 的 y 轴上的正确位置。我可以让包作者博客上提供的示例正常工作,但是当我更改值时,位置不正确。工作版本:

但是当我更改值时,p 值的位置太高了。

我想导致问题的第二个数据帧存在差异,但我无法弄清楚。谢谢!

0 投票
2 回答
110 浏览

r - 如何映射数据框,这是一个 tidyeval 错误吗?

想要映射数据框中的列并针对固定列对每列执行 t 检验。所需的输出将是一个数据框,每行都是 t 检验结果 - 可以使用map_dfr一次映射过程

深入研究 tidy eval,不确定这是否是 tidy eval 错误 - 非常感谢任何帮助!

mtcars作为玩具数据集)

遵循 tidyeval 指导:https ://tidyeval.tidyverse.org/dplyr.html 尝试查看是否引用/取消引用是问题,但没有骰子

0 投票
1 回答
47 浏览

r - 循环变量以生成余额表

我正在使用rstatix. 我可以为每个变量生成我想要的结果,但无法遍历多个变量以一次性生成一个链接表。

这失败了:

0 投票
2 回答
144 浏览

r - 多个组和变量的卡方

我有一个包含多个组(2 到 6 个级别的因子)和二分变量(0、1)的庞大数据集。

示例数据

我想为所有变量的每个组做一个卡方。

我设法通过使用两个嵌套的 for 循环使其工作,但我确信有更好的解决方案

我想我需要某种应用功能,但我想不通

0 投票
1 回答
115 浏览

r - R中具有两路方差分析的分类协变量

这是我的数据集的随机样本:

看起来像:

ROI是每个主题内的一个感兴趣区域,因此所有主题都有所有 4 个 ROI。

我想计算一个 2-way ANCOVA 4(ROIs [a/b/c/d] - within) x 2 (Age [young/old] - between) + Gender [covariate] 来确定age和, 控制.ROIvalueGender

为此,我计算了:

工作正常并输出:

然后,我想运行多重比较以生成 p 值,我可以将其绘制在箱线图上以进行分析的可视化。

我正在使用emmeans_test

但收到错误:

我不知道为什么,因为当我删除协变量时,成对比较工作正常。它是否与用作协变量的分类变量有关?我被困住了,想确保我在图表中报告了适当的 p 值。