如何修正RStudio中ANOVA分析的NA/NaN/Inf错误？

问题诊断 #

报错NA/NaN/Inf in 'y'和警告NAs introduced by coercion说明：

• 因变量Freq_RTC在转为数值型时生成了NA值，导致线性模型拟合失败
• 你的数据转换逻辑未处理全原始数据中的非数值格式问题

分步解决 #

1. 排查原始数据格式
   先确认变量类型和原始内容，定位问题根源：

# 查看数据集整体结构
str(CCR)
# 查看Freq_RTC的前几行原始值
head(CCR$Freq_RTC)
# 统计Freq_RTC的唯一值，排查特殊字符
table(CCR$Freq_RTC)

常见问题是欧洲格式数字（用.做千分位、,做小数点，如1.234,56），仅删除.会留下逗号，转数值时会生成NA。

2. 修正数据转换逻辑
   针对千分位/小数点格式问题，先删除千分位的.，再将小数点的,替换为.：

# 处理Freq_RTC：清理格式后转数值型
CCR$Freq_RTC <- as.numeric(gsub(",", ".", gsub("\\.", "", CCR$Freq_RTC)))
# 同理处理Spesa_RTC
CCR$Spesa_RTC <- as.numeric(gsub(",", ".", gsub("\\.", "", CCR$Spesa_RTC)))

3. 清理转换产生的NA值
   检查NA数量并选择处理方式：

# 统计两个变量的NA数量
sum(is.na(CCR$Freq_RTC))
sum(is.na(CCR$Spesa_RTC))

# 方式1：删除含NA的行（适合NA占比低的场景）
CCR_clean <- na.omit(CCR)
# 方式2：用均值填充NA（适合NA占比不高且数据分布均匀的场景）
CCR$Freq_RTC[is.na(CCR$Freq_RTC)] <- mean(CCR$Freq_RTC, na.rm = TRUE)
CCR$Spesa_RTC[is.na(CCR$Spesa_RTC)] <- mean(CCR$Spesa_RTC, na.rm = TRUE)

4. 确认ANOVA变量类型
   ANOVA要求自变量为定性变量（因子类型），若Spesa_RTC是数值型，需转为因子：

CCR$Spesa_RTC <- as.factor(CCR$Spesa_RTC)

5. 重新执行ANOVA分析
   用清理后的数据集运行代码：

anova3 <- aov(Freq_RTC ~ Spesa_RTC, data = CCR_clean) # 或用处理后的CCR
model.tables(anova3, type = "mean")
summary(anova3)
etaSquared(anova3)

内容的提问来源于stack exchange，提问作者Lavinia Coletta