1

我正在尝试将我的 TukeyHSD 结果导出到文件中。但是，我能找到的唯一解决方案是当它是一个简单的比较时，而不是当我还想考虑另一个因素时。

例如，这是我想要做的 ANOVA 类型：

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）：'which' 未指定任何因素“

但是重新排列语句并试图以任何方式将Treatmentand包含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)

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

4

0 回答 0