1

我正在尝试运行具有 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（“扫描”，“响应”）]）

 Resp
扫描 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 
 Resp
扫描 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 是数字，我不知道为什么会发生此错误。没有空单元格！而这个作为错误信息的一部分输出的奇怪表格也很奇怪。它似乎将呼吸视为一个因素？但我已经明确告诉它呼吸是数字的。当我从模型中取出呼吸时，它运行得非常好。然而，不幸的是，我确实需要包括呼吸。

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

在此先感谢您的帮助！！

4

1 回答 1

1

对于那些感兴趣的人，我想通了！我最终使用了线性混合模型而不是上面的模型，因为我的协变量处于扫描级别而不是主体间变量。

我的命令最终是：

Full_Anova_ALFF<- AlFF_Resp %>% group_by(Region) %>% do(fit=mixed(ALFF~Resp+Age*Scan+(1|ID)+(1|Scan),data=.))

希望这对将来的人有所帮助！

于 2021-05-12T15:17:25.630 回答