问题标签 [ancova]

我是 R 的初学者，我想从对 200 多个结果变量执行的 ANCOVA 测试中提取边际均值。当我只stack()对一个结果变量使用时效果很好，但是当我同时使用stack()and时出现错误lapply()。

这里我使用内置数据集“iris”来显示问题。数据集“iris”在 Species 中具有三个级别，我将其用Petal.Width作协变量、Species预测变量，并将前三列变量用作结果变量。

我的目的是同时提取对应结果变量的多个边际均值，而不是逐个提取。

在我提出前一个问题（如何使用 effect() 函数提取多个变量的边际均值？ ）之后，上面的代码在 @StupidWolf 的帮助下运行良好。然后在执行以下代码时出现错误：

错误是 Error in rep.int(factor(names(x), unique(names(x))), lengths(x)) : invalid 'times' value。

但是，当我只提取一个结果变量Sepal.Length的边际均值时（例如），我可以使用以下代码提取边际均值：

我不知道如何正确使用stack()和lapply()来提取边际手段。

非常感谢！

艾拉

我根据三年观察的空气质量数据创建了以下图表，并想知道这些斜率在两个时间段（2018 年 3 月至 2019 年 3 月至 6 月的平均值与 2020 年 3 月至 6 月的平均值）是否不同：

此处显示了我的数据框的快照：

该图是使用以下代码制作的：

我知道我需要首先消除序列相关的影响（因为自变量是时间序列），但我不确定如何做到这一点。我应该使用该date列吗？或者我应该使用虚拟列Date_2来做到这一点？此列只是Month.Date用于创建一系列数字和连续的 x 值的串联。

我使用该gls()函数来执行此操作，并且相信我已将该date列指定为我的序列相关。

我的尝试显示在这里：

输出：

当谈到线性回归输出时，我有点生疏，不知道如何解释这个。

此外，我想检查我的模型结构是否正确，以实现我想要的输出。

对此的任何帮助将不胜感激。

-TL;博士-

1. 我想在两条线上运行 ANCOVA，以确定Period变量的斜率是否不同。
2. 我想消除序列相关的影响，因为自变量是时间序列。

实现这一目标的最有效方法是什么？

如有必要，可以提供更多信息。

我正在用一个连续变量（DENSITY）和一个因子（SEASON）在 R 中做一个 ANCOVA。当我检查模型假设时，我得到一个名为：“恒定杠杆：残差与因子水平”的图，而不是“残差与杠杆”图。

图片：我得到了什么

图片：我想要什么

如何获得“残差 vs 杠杆”图？为什么我教科书中完全相同的代码会给出另一个 autoplot() 输出？

在此先感谢您的帮助！

我有一个基本的前后试验设计。两个随机组和两个测试，每个组中的每个参与者在干预之前（这里是 V1）和一个后（V2）。我对此完全陌生，并且对此进行了很多阅读，并且基于一些来源，有人建议将预测试作为协变量的 ANCOVA 测试是最合适的。

所以，我建模如下：

并检查残差的正态性并使用 Levene 检验来检验 V2 和 Group 之间的相关性。

对于某个感兴趣的变量，我得到了以下结果-

我随后进行了事后处理Turkey Test，发现它在那里也很重要。

我有一些问题：

1. 这是正确的方法吗？
2. 为什么我的 V1（预测试）协变量具有如此高的显着性水平，这意味着什么？（我假设随机化本质上意味着基线组之间没有差异）。
3. 基于此，我可以得出结论，在这个特定方面，两组之间确实存在差异吗？

我正在尝试创建这个

但在这样做时，R 会说以下内容：

在 plot.xy(xy, type, ...) 中：强制引入的 NAs`

我将在下面附上我的整个代码：

这是我的数据：

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

看起来像：

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 值。

这是我正在尝试做的事情的简要概述。我的数据集是从一项实验中收集的，该实验比较了 3 个条件（对照组、第 1 组和第 2 组）的性能。我的第一个任务是比较 gr1 和 control 之间是否存在差异。我对此进行了独立的 t 检验，结果并不显着（所以组平均值之间没有 sig 差异？）。我的下一个任务是包含一个协变量。该协变量有两个级别：高（分数范围从 4 到 7）和低（分数范围从 0 到 3）。我的代码如下所示：

我为 gr1 vs control 做了这个。输出显示条件显着，协变量显着，但交互作用不显着。我很困惑为什么条件现在很重要。我也只用性能〜条件进行了分析，但这并不重要。任何人都可以为我降低它。

我试图找到一些数据的估计边际均值，但不断收到错误消息。

返回

但是没有重复的列。我的所有列都没有命名为 std.error

检查错误表明这一点

我有一个基于分类变量 a[x] 预测 a0 的模型。现在，我需要添加另一个连续变量，从而将模型转换为 ANCOVA。很高兴听到任何建议。

在调整了一些线性模型后，我首先要测试回归斜率的同质性。第二步，这是我的疑问，我想使用事后测试来两个两个比较斜率。

这是一个从https://www.datanovia.com/en/lessons/ancova-in-r/修改的示例

获取数据

绘制回归线

检查回归斜率的同质性

在这里，我们可以看到统计上不显着的 p 值 4.15e-01

事后测试 emmeans_test 执行成对比较以确定哪些组不同。不过，我想采用多重比较程序来确定哪些B（斜率）与其他哪些不同。

有这个功能吗？提前致谢。