在一项调查中，我有40个答案（Q1_Biodiv，Q1_Poll，...），我想分析每个答案的社会描述符。 我想写出数据的描述性结果，并将其与方差分析的汇总结果结合起来以得到几个答案。 然后，我想将所有内容写在同一张Excel工作表中。
我的数据是这样的：
Q1_Biodiv Q1_Poll gender2 edu3 job3 country2
2 2 1 1 1 1
3 3 1 1 1 1
2 2 1 1 2 1
2 5 1 1 1 1
2 5 1 1 1 2
5 5 2 1 2 2
3 3 1 1 2 2
5 5 2 1 1 1
5 5 2 1 1 1
3 3 1 1 2 1
3 3 1 1 1 2
2 2 1 1 1 2
我尝试了以下代码： 我使用图书馆（WriteXLS）和图书馆（扫帚）
db1$job <- factor(db1$job)
job <- tabular( Q1_Biodiv+Q1_Poll ~ job* (mean+sd), data=db1
db1$education <- factor(db1$education)
edu<- tabular( Q1_Biodiv+Q1_Poll ~ education* (mean+sd), data=db1 )
db1$gender <- factor(db1$gender)
ge <- tabular( Q1_Biodiv+Q1_Poll ~ gender* (mean+sd), data=db1)
db1$country <- factor(db1$country)
co <- tabular( Q1_Biodiv+Q1_Poll~country* (mean+sd), data=db1)
####### ANOVA
av11 <- anova(lm(Q1_Biodiv~job+edu+gender+country,data=db1))
av11<-tidy(av11)
# save everything to the same excel sheet in the right order
# Dosen't work with this
WriteXLS(av1, "descriptives.xlsx")
有人能帮助我完成整个任务吗？？？ 提前非常感谢！