afex in R: Error: Empty cells in within-subjects design (i.e., bad data structure). But the data structure is fine
afex in R: Error: Empty cells in within-subjects design (i.e., bad data structure). But the data structure is fine
我正在尝试 运行 一个具有 1 个主体间变量和 2 个主体内变量的 ANCOVA，但我 运行 遇到了一个对我来说毫无意义的错误。我的数据如下所示：
Scan
ID
Region
ALFF
Age
Resp
1
20
AID
0.826
Adol
77.25
2
20
AID
1.116
Adol
73.18
1
22
AID
0.362
Adult
78.70
2
22
AID
0.849
Adult
72.58
1
20
MDM
0.826
Adol
79.25
2
20
MDM
1.116
Adol
71.18
1
22
MDM
0.778
Adult
79.70
2
22
MDM
0.291
Adult
73.58
我的 ANCOVA 代码是：
Full_Anova_ALFF<- AlFF_Resp %>% group_by(地区) %>% do(fit=aov_car(ALFF ~ Age扫描+错误（ID/ScanResp），数据=.））
当我 运行 时出现此错误：
转换为因数：年龄
错误：受试者内设计中的空单元格（即错误的数据结构）。
table(数据[c("扫描", "响应")])
Resp
Scan
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
Scan
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 是数字，我不知道为什么会出现此错误。没有空单元格！而作为错误消息的一部分输出的这个奇怪的 table 也很奇怪。它似乎将呼吸视为一个因素？但我明确地告诉过它呼吸是数字的。当我从模型中取出呼吸时，它 运行 完全没问题。然而，不幸的是，我确实需要包括呼吸。
有人知道出了什么问题吗？或者甚至只是我可以用来完成此操作的解决方法？
在此先感谢您的帮助！！
对于那些感兴趣的人，我想通了！我最终使用了线性混合模型而不是上面的模型，因为我的协变量处于扫描级别而不是主体间变量。
我的命令最终是：
Full_Anova_ALFF<- AlFF_Resp %>% group_by(地区) %>% do(fit=mixed(ALFF~Resp+Age*Scan+(1|ID)+ (1|扫描),数据=.))
希望这对以后的人有所帮助！
我正在尝试 运行 一个具有 1 个主体间变量和 2 个主体内变量的 ANCOVA，但我 运行 遇到了一个对我来说毫无意义的错误。我的数据如下所示：
| Scan | ID | Region | ALFF | Age | Resp |
|---|---|---|---|---|---|
| 1 | 20 | AID | 0.826 | Adol | 77.25 |
| 2 | 20 | AID | 1.116 | Adol | 73.18 |
| 1 | 22 | AID | 0.362 | Adult | 78.70 |
| 2 | 22 | AID | 0.849 | Adult | 72.58 |
| 1 | 20 | MDM | 0.826 | Adol | 79.25 |
| 2 | 20 | MDM | 1.116 | Adol | 71.18 |
| 1 | 22 | MDM | 0.778 | Adult | 79.70 |
| 2 | 22 | MDM | 0.291 | Adult | 73.58 |
我的 ANCOVA 代码是：
Full_Anova_ALFF<- AlFF_Resp %>% group_by(地区) %>% do(fit=aov_car(ALFF ~ Age扫描+错误（ID/ScanResp），数据=.））
当我 运行 时出现此错误：
转换为因数：年龄
错误：受试者内设计中的空单元格（即错误的数据结构）。 table(数据[c("扫描", "响应")])
Resp
| Scan | 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
| Scan | 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 是数字，我不知道为什么会出现此错误。没有空单元格！而作为错误消息的一部分输出的这个奇怪的 table 也很奇怪。它似乎将呼吸视为一个因素？但我明确地告诉过它呼吸是数字的。当我从模型中取出呼吸时，它 运行 完全没问题。然而，不幸的是，我确实需要包括呼吸。
有人知道出了什么问题吗？或者甚至只是我可以用来完成此操作的解决方法？
在此先感谢您的帮助！！
对于那些感兴趣的人，我想通了！我最终使用了线性混合模型而不是上面的模型，因为我的协变量处于扫描级别而不是主体间变量。
我的命令最终是：
Full_Anova_ALFF<- AlFF_Resp %>% group_by(地区) %>% do(fit=mixed(ALFF~Resp+Age*Scan+(1|ID)+ (1|扫描),数据=.))
希望这对以后的人有所帮助！