R：并排几个模型的 Anova 输出

Question

我正在尝试以与我通常为回归表所做的类似方式获得几个方差分析的良好发布准备输出。代码看起来有点像：

head(iris)
library(car)

# run model 1
lm1 <- lm(Sepal.Length ~ Species, data = iris)
summary(lm1)
a1 <- Anova(lm1)
# run model 2
lm2 <- lm(Petal.Width ~ Species, data = iris)
summary(lm2)        # long format (regression type output)
a2 <- Anova(lm2)    # short format (anova type output, which I need)

# what I usually do for regression outputs:
library(stargazer)
stargazer(lm1, lm2, type = "html", out ="iris.html")
# which yields a nice table with models side by side, including dependent variable names for each model; 
# this one I'd export as .html to Word and process there

# trying a similar joint table for anova type output:
stargazer(a1, a2, type = "html", out ="iris2.html")
# ... but that yields 2 separated tables and they can't be distinguished by dependent variables etc

# same problem here:
table <- rbind(a1, a2)
write.csv(as.data.frame(table), file = "iris2.csv")
# when I do this, I only get the same two tables underneath each other with their columnwise headers, 
# but without the distinguishing dependent variables 

由于我必须一遍又一遍地处理更多模型，因此我希望 Word 中的后期处理尽可能少。 我知道使用 LaTeX 有非常好的解决方案，但不幸的是，由于合著者，这些完全不可能。我使用了“xtable”、“pixiedust”和“export”包，但没有得到我想要的结果。

希望有人能帮忙， 提前致谢！

Answer 1

如果您有兴趣定期生成 HTML 和 Latex 表格，我建议您查看library(kableExtra)。这是关于HTML 和Latex 用法的优秀教程。

以下是如何创建 html 表格的示例。你会得到以下结果：

# Convert ANOVA results into dataframes allows for easier name manipulation
a1_new <- data.frame(a1)
a2_new <- data.frame(a2) 

# Putting all into one dataframe/table
anova_results <- data.frame(cbind(c("Species", "Residuals", "Species", "Residuals"), 
                                  rbind(a1_new, a2_new))) 
colnames(anova_results) <- c("", "Sum Sq", "Df", "F value", "Pr(>F)")
row.names(anova_results) <- NULL

# create HTML table using kableExtra
library(kableExtra)
anova_results %>% kable("html", digits=2) %>% 
  kable_styling(bootstrap_options = "striped", full_width = F) %>% 
  pack_rows(., "Sepal Length", 1, 2) %>% # groups rows with label
  pack_rows(., "Petal Width", 3, 4) # groups rows with label

如果你想要并排的值，这里有一种完成它的 hacky 方法...

anova_results2 <- data.frame(rbind(c("Sum Sq", "Df", "F value", "Pr(>F)", 
                                     "Sum Sq", "Df", "F value", "Pr(>F)"), 
                                   cbind(round(a1_new, 2), round(a2_new,2)))) 
colnames(anova_results2) <- c("", "", "", "","", "", "", "")
row.names(anova_results2)[1] <- ""

anova_results2 %>% kable("html") %>% 
  kable_styling(bootstrap_options = "striped", full_width = F) %>%
  add_header_above(c("", "Sepal Length" = 4, "Petal Width" = 4))

注意：为了获得正常工作的 Latex 表，您需要进行不同的修改。我确实认为，从长远来看，如果您想创建出版质量表，在 RMarkdown 中使用 Latex 可能是最好的方法。