为何paste0会破坏ANOVA汇总表格式？如何保留格式追加SST？

使用ezANOVA()生成的ANOVA汇总表原本格式规整，但用paste0()拼接SST值后，输出会变成字符串形式，原表格格式完全被破坏。

原正常输出代码 #

anova_weather_score = ezANOVA(weather_summary, dv = .(Score), wid = .(SubjectID), within = .(Backlight, TargetColor), between = .(ColorVision), detailed = TRUE)
Mn = mean(weather_summary$Score)
SST = sum((weather_summary$Score - Mn)^2)
print(anova_weather_score)

破坏格式的错误代码 #

anova_weather_score = ezANOVA(weather_summary, dv = .(Score), wid = .(SubjectID), within = .(Backlight, TargetColor), between = .(ColorVision), detailed = TRUE)
Mn = mean(weather_summary$Score)
SST = sum((weather_summary$Score - Mn)^2)
anova_weather_score = paste0(anova_weather_score, "SST = ", SST)
print(anova_weather_score)

解决方案 #

方法1：分开打印，保留原表格格式

先输出规整的ANOVA表格，再单独打印SST信息，互不影响：

anova_weather_score = ezANOVA(weather_summary, dv = .(Score), wid = .(SubjectID), within = .(Backlight, TargetColor), between = .(ColorVision), detailed = TRUE)
Mn = mean(weather_summary$Score)
SST = sum((weather_summary$Score - Mn)^2)

# 打印原ANOVA表格
print(anova_weather_score)
# 用cat和sprintf控制SST的输出格式
cat(sprintf("SST = %.4f\n", SST))

方法2：将SST加入结果列表，统一输出

ezANOVA()返回的是列表结构，直接给列表新增SST元素，打印时原表格格式会保留：

anova_weather_score = ezANOVA(weather_summary, dv = .(Score), wid = .(SubjectID), within = .(Backlight, TargetColor), between = .(ColorVision), detailed = TRUE)
Mn = mean(weather_summary$Score)
SST = sum((weather_summary$Score - Mn)^2)

# 给结果列表新增SST元素
anova_weather_score$SST = SST

# 打印整个列表，原表格格式不变，同时显示SST
print(anova_weather_score)

方法3：将SST整合进ANOVA表格（更融合的格式）

把SST作为一行添加到ANOVA结果表格中，实现完全统一的格式：

anova_weather_score = ezANOVA(weather_summary, dv = .(Score), wid = .(SubjectID), within = .(Backlight, TargetColor), between = .(ColorVision), detailed = TRUE)
Mn = mean(weather_summary$Score)
SST = sum((weather_summary$Score - Mn)^2)

# 获取ANOVA结果表格（列表中的$ANOVA元素）
anova_table = anova_weather_score$ANOVA
# 构造SST行，列数与原表格匹配
sst_row = data.frame(
  Effect = "Total Sum of Squares",
  DFn = NA, DFd = NA, F = NA, p = NA, `p<.05` = "", SS = SST,
  stringsAsFactors = FALSE
)
# 合并表格
anova_weather_score$ANOVA = rbind(anova_table, sst_row)

# 打印包含SST的完整表格
print(anova_weather_score)

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