问题标签 [ancova]

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 投票
0 回答
17 浏览

permutation - 为什么使用 set.seed() 不能得到恒定的结果?

我使用 aovp() 在 Rstudio 中进行置换 ANCOVA,并使用 set.seed() 来获得恒定的结果。

但是每次我重启Rstudio,结果还是会变!

这是我的代码:</p>

以及ftreatment(factor)的p值结果:</p>

  1. 第一次: p < 2e-16 ***
  2. 第二次:p = 0.00220 **

我应该怎么做才能得到恒定的结果?

0 投票
1 回答
106 浏览

r - R 中的 afex:错误:受试者内部设计中的空单元格(即错误的数据结构)。但是数据结构很好

我正在尝试运行具有 1 个主体间变量和 2 个主体内变量的 ANCOVA,但我遇到了一个对我来说毫无意义的错误。我的数据如下所示:

扫描 ID 地区 阿尔法夫 年龄 响应
1 20 援助 0.826 阿多尔 77.25
2 20 援助 1.116 阿多尔 73.18
1 22 援助 0.362 成人 78.70
2 22 援助 0.849 成人 72.58
1 20 MDM 0.826 阿多尔 79.25
2 20 MDM 1.116 阿多尔 71.18
1 22 MDM 0.778 成人 79.70
2 22 MDM 0.291 成人 73.58

我的 ANCOVA 代码是:

Full_Anova_ALFF<- AlFF_Resp %>% group_by(Region) %>% do(fit=aov_car(ALFF ~ Age Scan+Error(ID/Scan Resp), data = .))

当我运行它时出现此错误:

转换为因子:年龄

错误:受试者内部设计中的空单元格(即错误的数据结构)。表(数据[c(“扫描”,“响应”)])

扫描 X77.25 X73.1777777766667 X63.1944444433333 X70.3333333333333 X78.7 X72.5833333333333
X1 1 0 0 0 1 0
X2 0 1 0 0 0 1
X3 0 0 1 0 0 1
X4 0 0 0 1 0 0
扫描 X72.4833333333333 X78.25 X65.1833333333333 X71.9166666666667 X57.333333335 X65.55
X1 0 0 1 0 0 1
X2 0 0 0 1 0 0
X3 1 0 0 0 0 0
X4 0 1 0

Scan 是因子变量,而 resp 是数字,我不知道为什么会发生此错误。没有空单元格!而这个作为错误信息的一部分输出的奇怪表格也很奇怪。它似乎将呼吸视为一个因素?但我已经明确告诉它呼吸是数字的。当我从模型中取出呼吸时,它运行得非常好。然而,不幸的是,我确实需要包括呼吸。

有人知道出了什么问题吗?或者甚至只是我可以用来完成这项工作的解决方法?

在此先感谢您的帮助!!

0 投票
0 回答
158 浏览

r - WRS2“`.rowNamesDF<-`(x, value = value) 中的错误:'row.names' 长度无效”

我正在使用 WRS2 包在我的数据集上使用分类自变量(条件)、连续协变量(年龄)和连续因变量进行稳健的 ancovas。我在这里附上了数据集。

但是,当使用 WRS2 中的 ancova 函数时,我不断收到以下错误Error in `.rowNamesDF<-`(x, value = value) : invalid 'row.names' length

通过搜索,这似乎与使用 tibble 而不是数据框或重复的变量名称有关。这些都不是问题,因为我在导入后将数据转换为数据框,并且所有变量名称都是唯一的(我只将感兴趣的变量复制到单独的电子表格中并重命名它们)

但是,我的所有变量都不会发生错误,只有其中一些变量会发生错误,请参见下面的 MWE

起初我认为问题可能是因为缺少数据,所以我创建了一个过滤器来排除所有丢失的情况,但这并没有帮助解决问题。然后我尝试手动删除所有缺少数据的行,但它仍然抛出相同的错误。

我能找到的任何答案似乎都无法解释为什么它适用于某些变量而不适用于其他变量。

0 投票
0 回答
42 浏览

r - 混合MANCOVA \ MANCOVA与R中的重复测量

有谁能够帮我?

我试图运行一个混合 Mancova,其中包含 2 个主体内因素和 1 个主体间因素以及 1 个协变量。有6个因变量。哪个包或功能是正确的?我试过 stats::manova() 函数

但是当我将它与 SPSS 结果进行比较时,R 输出不正确。

非常感谢!

0 投票
0 回答
9 浏览

simulation - 多重插补的 Ancova 模拟

我正在尝试模拟 ANCOVA 分析的数据。我不知道从哪里开始。我需要3 个组的前测后测变量。我假设数据需要来自多元正态分布。

我预计每次比较(Grp1 与 grp3、grp2 与 grp 3 和 grp2 与 grp1)的差异为 5%(SD=10%),因为我遵循的是顺序设计。

我对平均差异(sd)和组均值感兴趣。

我将需要这些数据进行多重插补。

请帮忙!!

0 投票
0 回答
20 浏览

model - 重复测量 ANCOVA

我想运行重复测量 ANCOVA,但不知道如何正确编程...

我的重复测量模型如下: model <- aov(value ~ (Condition) + Error(pp_no/Condition), data = data_SCR_resh)

如何添加我的连续协变量(lum_diff)?...

一切都会有帮助的,娜塔莉

0 投票
1 回答
75 浏览

r - R WRS ancova 函数:if (nuhat < 2) stop("自由度必须大于或等于 2") 中的错误

我想从 WRS2 包 ( https://rdrr.io/cran/WRS2/man/ancova.html ) 运行一个强大的 ANCOVA。“比率”是我的 DV,“诱饵”是预测变量,“会话”是协变量。

我使用以下功能:

这是我的数据集:

在此处输入图像描述

但是,我收到以下错误消息:

if (nuhat < 2) stop("自由度必须大于或等于 2") 中的错误:需要 TRUE/FALSE 的缺失值另外:警告消息:1:In min(sub[vecn >= 12 ]) : min 没有非缺失参数;返回 Inf 2: In max(sub[vecn >= 12]) : max 没有非缺失参数;返回-Inf

任何想法我做错了什么?

0 投票
0 回答
41 浏览

r - 在 R 中提取 ANCOVA 的系数

我正在尝试可视化 ANCOVA 的预测输出

型号为:

m1 <- lm(PD~SR*Type, data = df(带交互项)

m2 <- lm(PD~SR+Type, data = df(无交互项)

我的示例代码

(例子)

在此处输入图像描述

我只想为所有行设置灰色,但上面的代码不起作用,因为我使用了geom_line. 这就是为什么我想得到interceptand slope,然后用它geom_abline来制作只用一种颜色填充的多条回归线。有人对此有任何建议吗?

或者如果我的理解有误,请指导我。让我们假设 ANCOVA 的交互项很重要,所以我期望在这种情况下每种类型 ( a, b, c) 的截距/斜率是不同的。如果我在这里为每种类型提取模型的截距和斜率并构建它们的直方图以便从其他方法(即平均方法)评估我的结果是否有意义?

0 投票
0 回答
44 浏览

r - sm.ancova 的事后检验(非参数 ancova)

我正在尝试在我的数据集上使用非参数 ancova。

目前,我正在sm为该功能使用包sm.ancova

但是,我不确定我应该对数据使用哪种事后测试。

会有推荐的统计测试或软件包吗?

0 投票
0 回答
16 浏览

spss - 运行 MANOVA,其中 Pillai's/Wilks 不显着,但其中一个 DV 在我的受试者间效应输出表中非常重要

我是一个统计新手,我的教授告诉我要为我正在检查的东西运行 MANOVA。基本上,我想看看种族和某个象限分组之间是否存在交互作用,这些结果变量是总体测量 (ders_tot) 的子量表。

ANCOVA(一个 DV)已经发现种族和 ders_tot 的象限分组之间存在交互作用。

我的 MANOVA 输出显示,使用 Pillai's/Wilks 没有显着性(两者的 p = .098),但在 SPSS 中还有一个自动生成的受试者间效应表,表明一个特定结果变量的强交互显着性( p = .003)。其他 DV 远没有显着性(有些高达 p = .27 或 p = .66)。

我的 MANOVA 显着性(或缺乏显着性)是否被高度不显着的变量严重扭曲?我是否仍然“允许”对 MANOVA 中包含的一个表明重要意义的特定变量进行分析?我也有数据可视化/图表输出,为分析该特定变量提供了强有力的理由。

(编辑:以下问题已修复)

[另外,我注意到我的一个协变量总是在 SPSS 中以 1 df 运行,而它应该是 2。我已经三重检查了变量类型并添加了标签等等,但无法让它运行适当地。当我在 R 中运行相同的分析时,df = 2。这不会影响我的信号。发现很多,但这让我发疯!]