问题标签 [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.
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 进行成对比较......
r - 如何获得 r 中多个组的平均 SD 和 Pvalue?
你好我所有的数据框看起来像
我想分别根据类和类别计算(test1和test2之间)的平均SD和Pvalue
我使用 dplyr 计算平均值和 SD 并努力计算 Pvalue,因为我的数据集包含 1000 行和 4 个不同的类别和 8 个类
这是使用 dplyr 表示 mean 和 sd 后得到的结果
而SD只是手动输入以显示它的错误值我想要的输出是
提前致谢
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 下构建的
你以前处理过这个问题吗?你是怎么解决的?谢谢!
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 的数据集的代码:
r - 我可以将 anova_test (rstatix) 类的对象转换为 aov、anova 或 aovlist 类的对象吗?
重复测量方差分析的 DATANOVIA 指南中的一些样本数据:
我想使用report
来自出版商 EasyStats ( https://github.com/easystats/report ) 的软件包。它接受 、 或 类的aov
ANOVAanova
对象aovlist
。ANOVA 中的对象rstatix
属于类anova_test
,因此与report
.
有什么方法可以转换对象类,还是应该使用基 R 进行重复测量方差分析?
回复 1(已被删除):
class(res.aov) <- "aov"
或class(res.aov) <- c("aov", class(res.aov))
。您可以添加或更改对象的类属性
所以我尝试了这个:
然后你的建议:
类更改似乎已经奏效,但将对象传递给report
. 我该如何进行?
附带说明一下,如果有人可以澄清这些是否等效,我的问题可以通过将我的 ANOVA 传递给基础 R 来解决:
是
anova_test(response ~ A * B + Error(subject/(A * B)))
在rstatix
相当于
aov(response ~ A * B + Error(subject/(A * B)))
在stats
?
r - 使用 rstatix 在 ggplots 上的 p 值位置不正确
我无法使用 rstatix 将 p 值放置在 ggplot 的 y 轴上的正确位置。我可以让包作者博客上提供的示例正常工作,但是当我更改值时,位置不正确。工作版本:
但是当我更改值时,p 值的位置太高了。
我想导致问题的第二个数据帧存在差异,但我无法弄清楚。谢谢!
r - 如何映射数据框,这是一个 tidyeval 错误吗?
想要映射数据框中的列并针对固定列对每列执行 t 检验。所需的输出将是一个数据框,每行都是 t 检验结果 - 可以使用map_dfr
一次映射过程
深入研究 tidy eval,不确定这是否是 tidy eval 错误 - 非常感谢任何帮助!
(mtcars
作为玩具数据集)
遵循 tidyeval 指导:https ://tidyeval.tidyverse.org/dplyr.html 尝试查看是否引用/取消引用是问题,但没有骰子
r - 循环变量以生成余额表
我正在使用rstatix
. 我可以为每个变量生成我想要的结果,但无法遍历多个变量以一次性生成一个链接表。
这失败了:
r - 多个组和变量的卡方
我有一个包含多个组(2 到 6 个级别的因子)和二分变量(0、1)的庞大数据集。
示例数据
我想为所有变量的每个组做一个卡方。
我设法通过使用两个嵌套的 for 循环使其工作,但我确信有更好的解决方案
我想我需要某种应用功能,但我想不通
r - R中具有两路方差分析的分类协变量
这是我的数据集的随机样本:
看起来像:
ROI
是每个主题内的一个感兴趣区域,因此所有主题都有所有 4 个 ROI。
我想计算一个 2-way ANCOVA 4(ROIs [a/b/c/d] - within) x 2 (Age [young/old] - between) + Gender [covariate] 来确定age
和, 控制.ROI
value
Gender
为此,我计算了:
工作正常并输出:
然后,我想运行多重比较以生成 p 值,我可以将其绘制在箱线图上以进行分析的可视化。
我正在使用emmeans_test
:
但收到错误:
我不知道为什么,因为当我删除协变量时,成对比较工作正常。它是否与用作协变量的分类变量有关?我被困住了,想确保我在图表中报告了适当的 p 值。