R中混合模型ANOVA的结果矛盾

Question

我是混合模型方差分析的新手。我试图将混合模型ANOVA拟合到以下结构的某些数据上：我将许多主题划分为三个类别（因子之间）之一。对于每个对象，我在12个偏心率（因子内）下进行了多次测量。为了得到我的结果，我运行了以下命令：

aov(measurement ~ ref_group*eccentricity + Error(subject_id/eccentricity),
    data = data_clean_long)

aov_car(measurement~ ref_group*eccentricity + Error(subject_id/eccentricity),
        data = data_clean_long)

lmer(measurement~ ref_group * eccentricity + (1 | subject_id), 
     data = data_clean_long) %>% 
  anova()

这些函数返回不同的结果。此外，我有两个如上所述构成的数据集。对于一个数据集，两个功能是一致的，一个是不同的；对于另一个数据集，两个不同的功能是一致的。我在这里做蠢事吗？为什么我没有得到相同的结果？