自动化数据框的交叉制表和重要性测试？

Question

这似乎很容易，但我找不到它。我有一个id，sex，...，age，Q1，...，Q2等数据框。

# EXAMPLE CODE
big.data= data.frame(ID= c(1,2,3,4,5,6,7), sex= c(1,1,1,2,2,2,NA), 
age= c(25,25,35,NA,35,45,45), Q1= c(1,2,NA,2,2,2,2), Q2= c(1,NA,2,2,3,3,3) )

我想选择性别...年龄作为类别，Q1 ... Q2作为感兴趣的变量。我希望R根据感兴趣的变量（本例中为四个表）为每个类别生成CSV表，这些表与信息（不一定是格式或布局）相似（通过回答问题的列数和百分比）：

按性别分列的Q2表显示了对于每个对Q2的响应的所有性别响应之间的整个表的卡方检验和显着性检验的结果，其中p

如何在R中轻松完成？