问题标签 [emmeans]

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 回答
595 浏览

r - 有没有办法将 emmeans() 中的数据输出为干净的表格格式?

我想从 R 中的 emmeans() 函数的均值和对比中制作一个表格。有没有简单的方法可以做到这一点?

0 投票
0 回答
351 浏览

r - 使用 clm 的 emmeans 对组进行多重比较:UseMethod 中的错误(“vcoc”):没有适用于“vcov”的方法应用于“clm”类的对象

我正在尝试使用 emmean/lsmeans 进行事后测试,以进行组的多重比较。事实上,我正在尝试重新创建一个在这里找到的示例:https ://rcompanion.org/handbook/G_11.html 。

因此,我有一个看起来像这样的模型:

然后我正在尝试探索交互效果的组分离:

但我收到以下错误:

我也对 emmeans 进行了同样的尝试,但出现了同样的错误。

有人可以指导我出了什么问题吗?:)

0 投票
1 回答
652 浏览

r - emmeans - 控制与治疗超过一个因素

我正在使用对对照组进行自定义比较。trt.vs.ctrl如果我只对比较一个因素感兴趣,那么该方法对我来说非常有效,但是当我将比较设置为更复杂时(即控制组由 2+ 变量的特定组合描述)失败(或失败) )。

下面的示例代码。假设使用pigs数据,我想将所有饮食与低百分比的鱼类饮食进行比较。请注意在nd数据框中,“鱼”只有 9% 与之关联。但是,当我运行 时emmeans,该函数不会对嵌套进行处理,虽然控制是正确的,但处理组还包括鱼和百分比的各种值。这意味着调整是错误的。

所以我能想到的两种方法:

  1. 在这种情况下,我该如何emmeans拾取嵌套,或者
  2. 如何dunnettx手动进行调整(=我可以使用调整“无”,然后取出我真正想要的测试,并自己调整 p 值?)。

感谢你的帮助。

使用建议include非常有效。在这里发布我的代码,以防其他人将来遇到同样的问题。

0 投票
1 回答
644 浏览

r - 在 emmip 图中更改线条粗细

有人可以帮我理解为什么我似乎无法用随附的代码更改附图的线条粗细吗?我想让线条更粗,特别是对于虚线类型的组。添加 scale_size_manual 线不会改变任何东西(尤其是线的粗细)。

我不是 ggplot2 方面的专家,只是把它搞砸了,所以毫无疑问我做错了什么。谢谢

0 投票
2 回答
37 浏览

r - 向 data.frame 添加一个新列,其因子取决于另一个 data.frame 的条件

我有 2 个不同的 data.frames。我想在组 data.frame(LEVEL 和 SPECIES)给定的条件下将 grouping$.group 列添加到物候 data.frame。我尝试了使用 by= 的 merge() 函数,但它一直给我“fix.by(by.y, y) 中的错误:'by' 必须指定一个唯一有效的列”。抱歉,这似乎是一件非常容易的事情。我是初学者。。

0 投票
1 回答
121 浏览

r - R - 计划比较中完全关闭的自由度(使用 afex - aov_ez - emmeans - 对)?

在我们的任务中,参与者对与厌恶刺激配对的可能性不同的刺激(轨迹)提供主观恐惧评级(因此 DV 是恐惧评级)。有两组、三个块和三个轨迹,导致 2(组:实验/轭控制)x 3(块:GEN1-3)x 3(轨迹:G1/G2/G3)RM ANOVA。为此,我使用以下代码:

并得到这些结果:

我有一个先验假设要测试,我需要将评级与组内的不同轨迹进行比较,因此我使用代码运行计划的比较

并得到以下结果:

我在这里只显示两个组的第一个块 (GEN1),但其他块中的 DF 完全相同。

有人可以向我解释这些自由度是如何计算的,它们是否正确?我似乎找不到直接适用于我的问题的答案。另外,大约一年前我进行了类似的分析,但没有得到这么大的 DF?

先感谢您!

0 投票
1 回答
391 浏览

r - emmeans 如何计算用于比较均值的置信区间

我正在寻找更多关于如何emmeans计算置信区间的背景和文档,这些置信区间在以下小插图中概述的均值的图形比较中使用: https ://cran.r-project.org/web/packages/emmeans/vignettes/comparisons.html #图形

在“图形比较”部分有一个我所指的例子。特别是,我对用于比较均值的带有箭头的红线感兴趣。

上面写着:

如果一个平均值的箭头与另一组的箭头重叠,则差异并不“显着”。

但是用于比较均值的红线区间是如何计算的呢?这是在某处记录的吗?

0 投票
2 回答
812 浏览

r - R:使用 emmeans 包一次运行多个事后测试

我正在研究一个以几种不同类型的蛋白质作为列的数据集。它看起来像这样这是简化的,原始数据集包含 100 多种蛋白质。我想看看在考虑随机效应(=id)时,蛋白质的浓度是否因治疗而异。我设法一次运行多个重复的方差分析。但我也想根据治疗对所有蛋白质进行成对比较。我首先想到的是使用 emmeans 包,但我在编码时遇到了麻烦。

所以我尝试了不同的方法

我不明白为什么即使我在方差分析模型中使用了“cbind(Protein1, Protein2, Protein3)”。R 仍然只给我一个结果,而不是像下面这样

我该如何编码或者我应该尝试不同的包/功能?

我一次运行一种蛋白质没有问题。但是,由于我要运行 100 多种蛋白质,因此将它们逐一编码非常耗时。

任何建议表示赞赏。谢谢!

0 投票
1 回答
755 浏览

r - R中的包emmeans不返回效果大小

我正在关注本教程以及?eff_size从包emmeans到计算eff_size()下面的回归模型。

但我得到了错误:need an object with call component来自eff_size()电话。我错过了什么吗?

0 投票
1 回答
444 浏览

r - r - 多级重复测量 ANCOVA 的 emmeans 成对分析

我正在通过 AOV 包使用多级框架构建重复测量 ANCOVA。我有一个连续响应变量、两个因子预测变量和 3 个连续协变量。我的模型脚本如下:

我的数据集片段显示了它的格式:

我的条件嵌入在参与者 ID 的错误术语中,因为这是一个完全重复测量模型。

我正在尝试对这些值进行成对分析。我的输出提供综合 F 检验:

通过以下方式计算 emmeans 时:

我得到一个合理的输出。

但是,将其用于协变量时:

我得到以下输出:

我在这里做错了什么,成对比较不适用于协变量?