尝试获取重复测量方差分析的 AOV 结果时出错答案

【问题标题】：An error while trying to get AOV results of repeated measures ANOVA尝试获取重复测量方差分析的 AOV 结果时出错
【发布时间】：2020-10-14 13:58:12
【问题描述】：

我的数据是长格式的。我正在寻找时间对分数的影响。尝试了以下代码：

res.aov <- aov(score~time+Error(ID/time), data=dat)

显示如下错误：

Warning message:
In aov(score ~ time + Error(ID/time), data = dat) :
  Error() model is singular

我得到以下结果：

> summary(res.aov)

Error: ID
          Df Sum Sq Mean Sq F value Pr(>F)
time       1   36.9   36.91   2.664  0.112
Residuals 32  443.4   13.85               

Error: ID:time
          Df Sum Sq Mean Sq F value Pr(>F)
time       1   0.02  0.0168    0.01  0.919
Residuals 32  51.51  1.6096               

Error: Within
           Df Sum Sq Mean Sq F value Pr(>F)
Residuals 110  79.43  0.7221               
>

为什么我看不到内表的 p 值？以及如何处理这个错误？

【问题讨论】：

标签： r statistics data-analysis anova

【解决方案1】：

您遇到该错误是因为您的变量之一每个级别只有一个条目，这意味着无法估计其影响，例如：

dat =  data.frame(
ID = factor(letters[1:20]),
time=rep(1:5,4),
score=rnorm(20))

table(dat$ID)

a b c d e f g h i j k l m n o p q r s t 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

res.aov <- aov(score~time+Error(ID/time), data=dat)
Warning message:
In aov(score ~ time + Error(ID/time), data = dat) :
  Error() model is singular

您需要确定设计并运行正确的模型。如果你有20个不相关的数据点，即20个人，不同时间不同分数，你需要的是：

aov(score~time,data=dat)

如果您拥有的是重复测量，即在不同时间点测试相同分数的个人，您的数据将如下所示：

dat =  data.frame(
    ID = factor(rep(letters[1:4],each=5)),
    time=factor(rep(1:5,4)),
    score=rnorm(20))

table(dat$ID,dat$time)
   
    1 2 3 4 5
  a 1 1 1 1 1
  b 1 1 1 1 1
  c 1 1 1 1 1
  d 1 1 1 1 1

summary(aov(score~time+Error(ID/time),data=dat))

Error: ID
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals  3  2.182  0.7272               

Error: ID:time
          Df Sum Sq Mean Sq F value Pr(>F)
time       4  3.161  0.7903   0.654  0.635
Residuals 12 14.492  1.2077

还要注意，在上面的重复测量示例中，时间必须是一个因素，而不是连续的

【讨论】：