用于显示带有均值、SE 和显着性字母（大写字母和小写字母）的双向方差分析的 R 包

Question

我正在寻找一个可以帮助我显示双向方差分析结果的 R 包。具体来说，我需要在输出中包含以下元素：

平均值（非估计值，原始值） 标准误差 (SE)（未估计，原始） 用字母表示含义，使用大写和小写字母（例如“A”、“a”、“B”、“b”）

适合<- aov(log(response) ~ treatments(4 levels) * species(4 types), data = sample_data) tukeyHSD(fit)

桌子应该看起来像 在此输入图片描述

大写字母按列比较物种，小写字母比较处理。 是否有提供此功能的 R 包，或者可以帮助我实现此功能的功能/包的任何组合？交互作用不显着，因此应该删除交互作用吗？

Emmeans 只给出估计平均值，MultcompView 只给出小字母（可能是我错了），最终的表格字母太多。

Answer 1

您可以通过关注基本步骤，以更简洁的方法实现您想要的输出。这是该过程的简化版本：

执行双向方差分析：

fit <- aov(log(response) ~ treatments * species, data = sample_data)

使用“emmeans”获取估计边际均值并执行成对比较：

library(emmeans)
emm <- emmeans(fit, ~ treatments + species)

使用 multcompView 生成具有意义的字母：

library(multcompView)
tukey_letters <- cld(emm, Letters = letters)

将结果合并到最终表中：

final_table <- data.frame(
  Treatment = tukey_letters$treatments,
  Species = tukey_letters$species,
  Mean = summary(emm)$emmean,
  SE = summary(emm)$SE,
  Letters = tukey_letters$.group
)

print(final_table)