我使用2 x 2受试者之间的方差分析来计算性别（男性，女性）和自我效能（低，高）的主要影响，以及性别和自我效能之间的相互作用。我想计算平均值的标准误差。设计是不平衡的。我计算了调整后的平均值，但我无法计算标准误差。
> anova <- aov(RT ~Gender*efficacy, data=data)
> summary(anova, type="III")
Df Sum Sq Mean Sq F value Pr(>F)
Gender 1 14.3 14.289 4.285 0.0399 *
efficacy 1 2.1 2.053 0.616 0.4337
Gender:efficacy 1 0.2 0.167 0.050 0.8231
Residuals 173 576.9 3.334
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> model.tables(anova, type="means", se=T)
Design is unbalanced - use se.contrast() for se's