在R中，当存在多个因素时，如何将TukeyHSD结果导出到行/列

Question

我正在尝试将TukeyHSD结果导出到文件中。但是，我只能找到一个简单的比较方法，而不是我也想考虑另一个因素时。

例如，这是我尝试执行的方差分析的类型：

dat1 <- lm(Count~Treatment*Keio_gyrA, data=colony_counts2)
Anova(dat1)

为了查看这些类型的结果（按治疗方式分类）：

我唯一能找到的解决方案是：

mod1 <- aov(Count~Treatment*Keio_gyrA, data=colony_counts2)
tk1 <- TukeyHSD(mod1, 'Treatment', ordered = TRUE)
write.csv(as.data.frame(tk1$'Treatment:Keio_gyrA'), "colony counts ANOVA2.csv")

但是无论我做什么，它只会给我Count~Treatment，Count~Keio_gyrA或一个错误。

上面的错误给了我

  

TukeyHSD.aov（mod1，“ Treatment”，“ Keio_gyrA”，ordered = TRUE）中的错误：“未指定任何因素”

但是重新整理该语句并尝试以任何方式一起包含Treatment和Keio_gyrA似乎不起作用。

我还尝试了cat函数（下面的代码），虽然我可以获取.txt文件，但它不在列/行中，因此很难从中获取信息。

dat3 <- lm(Max_OD~Treatment:Keio_gyrA, data=max_OD2)
Anova(dat3)
out3 <- capture.output(TukeyHSD(aov(dat3)))
cat(out3, file="ANOVA Max OD.txt", sep=" ", append=TRUE)

任何帮助将不胜感激。谢谢。