R语言单因素ANOVA无P值问题：Estimated effects may be unbalanced报错排查

1. 提示信息解读 #

• "8 observations deleted due to missingness"：你的数据中存在8条含缺失值的记录（hatch或conc列有NA），aov()函数自动删除了这些行，这是常规清理操作，但需要确认缺失值的分布是否合理。
• "Estimated effects may be unbalanced"：将conc转为因子后，不同浓度组的样本量差异过大（从你提供的数据表格来看，部分浓度组仅1个样本），导致方差分析的设计不平衡，违反了ANOVA对组内样本量相对均衡的隐含要求。

2. 无P值的核心原因 #

当某个浓度组的样本量仅为1时，该组的组内变异无法计算（自由度为0），会导致整个模型的残差自由度不足，最终无法生成F值和P值。

3. 具体解决步骤 #

步骤1：检查数据缺失与分组样本量 #

先运行以下代码确认数据细节：

# 查看各列缺失值数量
colSums(is.na(CMTT_nano))
# 查看每个浓度组的样本数
table(CMTT_nano$conc)

步骤2：处理不平衡分组 #

• 若存在样本量为1的浓度组：结合实验逻辑，要么将其与邻近浓度组合并（比如把0.001和0.01组合并），要么直接剔除这些小样本组。
• 若缺失值随机分布且数量少，可保留当前删除操作；若缺失值集中在某组，考虑用组均值填充补全，或直接剔除该组。

步骤3：重新分析数据 #

• 调整分组后，重新运行ANOVA即可正常生成P值：

# 假设已处理好分组和缺失值
model <- aov(hatch ~ conc, data = CMTT_nano)
summary(model)

• 若分组不平衡问题无法解决，改用非参数检验（无需满足ANOVA的正态性和方差齐性假设）：

kruskal.test(hatch ~ conc, data = CMTT_nano)

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