如何在 R 中以不相等的样本大小运行 2x2 混合方差分析？答案

【问题标题】：How to run 2x2 mixed anova with unequal sample size in R?如何在 R 中以不相等的样本大小运行 2x2 混合方差分析？
【发布时间】：2021-02-26 11:42:28
【问题描述】：

我正在尝试使用 R 运行具有不等样本大小的 2x2 混合 ANOVA。数据是从 30 个具有两种不同条件（即 2 个水平内因子）的个体中收集的，并通过 k 均值聚类分析（即，第 1 组的编号 = 11，第 2 组的编号 = 19）。这是我的代码示例和输出：

summary(aov(JH ~ Box+Group+Box:Group+Error(P/Box), data = d3))

Error: P
          Df Sum Sq Mean Sq F value Pr(>F)
Group      1  0.027 0.02715    1.56   0.22
Box:Group  1  0.001 0.00078    0.04   0.83
Residuals 27  0.470 0.01741               

Error: P:Box
          Df   Sum Sq  Mean Sq F value Pr(>F)
Box        1 0.000000 1.00e-07    0.00   0.97
Group      1 0.000022 2.17e-05    0.24   0.63
Box:Group  1 0.000032 3.24e-05    0.35   0.56
Residuals 27 0.002488 9.21e-05

与我使用另一个数据集运行具有相同样本大小的 2x2 混合 ANOVA 时的输出不同（下面附上的输出），我在错误：P 下有另一个交互效果。

summary(aov(HipLoadingK ~ Box*Sex+Error(P/Box), data = TW))

Error: P
          Df  Sum Sq  Mean Sq F value  Pr(>F)   
Sex        1 0.02578 0.025779   8.038 0.00841 **
Residuals 28 0.08980 0.003207                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Error: P:Box
          Df   Sum Sq  Mean Sq F value Pr(>F)  
Box        1 0.003454 0.003454   4.271 0.0481 *
Box:Sex    1 0.000066 0.000066   0.082 0.7765  
Residuals 28 0.022645 0.000809

只要我为每个组间/组内主效应和交互效应解释适当的结果，这个结果就可以了吗？或者我应该编辑我的代码以获得正确的输出？如果是这样，请告诉我应该在我的代码中添加或编辑什么。

【问题讨论】：

你的问题分两部分，解释部分是否可以，属于交叉验证。编码部分，你的输出看起来很奇怪我不明白你在这里没有可重复的例子
见stackoverflow.com/questions/5963269/…
@StupidWolf 感谢您考虑我的问题。我刚刚发现我的分配代码导致了错误。修复代码后，它工作了！

标签： r anova

【解决方案1】：

当我将数据集从短格式更改为长格式时，我不小心为每个参与者设置了不同的组名顺序。例如，P01 是 X 条件下的 A 组，但他或她的组是 Y 条件下的 B。在我修复代码后，它运行良好。

【讨论】：