在R中无法获取重复测量ANOVA模型残差的问题求助

我有一个R数据框，头部数据如下：

Sujet T E O P meanTR
1     1 0 0 0 0.97
1     1 0 0 0 1.44
2     0 1 0 1 0.94

各变量的含义：

• Sujet：取值范围是1到12
• T、E、O、P：都是二分类变量，取值0或1
• meanTR：数值型变量

我想要做方差分析，于是运行了下面的代码：

model_all <- aov(meanTR ~ E*O*P+ Error(Sujet/E*O*P), data = df)

接下来我想提取模型的残差来绘制诊断图，执行了：

res <- residuals(model_all)

但返回的结果是NULL，这让我很困惑。

我在网上找了方法，尝试用proj()函数处理模型：

model_all.pr <- proj(model_all)
res <- model_all.pr[[3]][, "Residuals"]

结果报了"subscript out of bound"（下标越界）的错误；我又试了：

res <- model_all.pr[[3]]["Residuals"]

这次返回的是NA，还是不行。

我搞不清楚问题出在哪，希望有人能帮我解决，最终目标是能顺利运行下面的代码来做残差诊断：

plot(res)
qqnorm(res)

内容的提问来源于stack exchange，提问作者HelpASisterOut