R - 创建具有多列和多行的交叉表

Question

我习惯了 SPSS，我真的很喜欢使用自定义表格来报告调查数据。如果我能在 R 中做类似的事情，我会很高兴。

我想做的是一个包含多行和多列的表格，其中包含列百分比和计数（N - 百分比的基数）

这是一个调查数据的示例代码：

set.seed(321)
ID <- seq(1:200)
Age <- sample(c("18-34", "35-59"), 200, replace = TRUE)
Sex <- sample(c("Male", "Female"), 200, replace = TRUE)
TOTAL <- rep(c("TOTAL"), 200)
Edu <- sample(c("Lower", "Middle", "Higher"), 200, replace = TRUE)
PurchaseInt <- sample(c("Definitely yes", "Somewhat yes", "Somewhat not", "Definitely not"),200, replace=TRUE)
Relevance <- sample(c("Definitely fits my needs", "Somewhat fits my needs", "Somewhat does not fit", "Definitely does not fit"),200, replace=TRUE)

DF <- data.frame(ID,TOTAL,Sex,Age,Edu,PurchaseInt,Relevance)
head(DF)
  ID TOTAL    Sex   Age    Edu    PurchaseInt               Relevance
1  1 TOTAL   Male 35-59  Lower Definitely yes  Somewhat fits my needs
2  2 TOTAL   Male 35-59 Higher   Somewhat not Definitely does not fit
3  3 TOTAL   Male 18-34 Higher Definitely yes   Somewhat does not fit
4  4 TOTAL Female 18-34  Lower   Somewhat not Definitely does not fit
5  5 TOTAL Female 18-34 Higher Definitely yes   Somewhat does not fit
6  6 TOTAL Female 18-34 Higher Definitely not Definitely does not fit

# Simple table, 1 variable by 1 variable, no N (BASE) BAD TABLE :(
prop.table(table(DF$PurchaseInt, DF$Sex),2)

                 Female Male
  Definitely not   0.28 0.30
  Definitely yes   0.25 0.28
  Somewhat not     0.29 0.24
  Somewhat yes     0.17 0.18

我真的很想得到这样的东西（从 SPSS 完成）：

我意识到将计数与 col 百分比结合起来可能会更加棘手。对我来说至关重要的是找到在一个表中报告多行和多列（尤其是多列）的可能性，因为这有助于数据分析很多。

Answer 1

这有两个部分：首先，您要创建表格，然后要报告它。您按不同的边距拆分，然后将其放在一个表中，这有点奇怪；我也不确定你是如何得到这些数字的，它们是指列百分比吗？如果是这样，我会用你的随机种子得到不同的。

无论如何，这是第 1 部分，它会为您获取数据。

# a useful function
table_by <- function(row_var, col_var = NULL) {
  # the repeated t() below ensures you have a 4 x 1 matrix
  tbl <- if (is.null(col_var)) t(t(table(DF[[row_var]]))) else table(DF[[row_var]], DF[[col_var]])
  tbl <- prop.table(tbl, 2)
  tbl <- round(tbl, 2) * 100
  tbl
}

col12 <- rbind(table_by("PurchaseInt", "Sex"), table_by("Relevance", "Sex"))
col34 <- rbind(table_by("PurchaseInt", "Age"), table_by("Relevance", "Age"))
col56 <- rbind(table_by("PurchaseInt", "Edu"), table_by("Relevance", "Edu"))
percent_rows <- cbind(col12, col34, col56)
whole_table <- cbind(
  rbind(table_by("PurchaseInt"), table_by("Relevance")),
  percent_rows
)

# should be the data you want
whole_table

对于第二部分，你可以使用我的huxtable 包——还有其他的：

library(huxtable)
wt_hux <- as_hux(whole_table, add_colnames = TRUE, add_rownames = TRUE)
number_format(wt_hux)[-2,] <- "%.0f%%"
number_format(wt_hux)[2,]  <- "%.0f"
wt_hux[1, 1:2] <- c("", "Total")
wt_hux[2, 1]   <- "Total"
wt_hux <- insert_row(wt_hux, c("", "Total", "Sex", "", "Age", "", "Edu", "", ""))
colspan(wt_hux)[1, c(3, 5, 7)] <- c(2, 2, 3)
align(wt_hux)[1, c(3, 5, 7)] <- "center"
wt_hux <- insert_column(wt_hux, c("", "", "Total", "PurchaseInt", "", "", "", "Relevance", "", "", ""))
rowspan(wt_hux)[c(4, 8), 1] <- 4

bottom_border(wt_hux)[c(1, 6, 10), ] <- 1 # for example

# should look roughly the way you want. You can print it to PDF or HTML:
wt_hux

Answer 2

非常感谢您的回答！很有帮助。我还没有遇到huxtable。

通过一些修改，我可以让它按照我想要的方式工作。代码如下：

# a useful function
table_by <- function(row_var, col_var = NULL) {
  # the repeated t() below ensures you have a 4 x 1 matrix
  tbl <- if (is.null(col_var)) t(t(table(DF[[row_var]]))) else table(DF[[row_var]], DF[[col_var]])
  tbl <- prop.table(tbl, 2)
  tbl <- round(tbl, 2) * 100
  tbl
}

# HERE I also added a table showing counts for demographics
col12 <- rbind(table(DF$TOTAL), table_by("PurchaseInt", "TOTAL"), table_by("Relevance", "TOTAL"))
col34 <- rbind(table(DF$Sex), table_by("PurchaseInt", "Sex"), table_by("Relevance", "Sex"))
col56 <- rbind(table(DF$Age), table_by("PurchaseInt", "Age"), table_by("Relevance", "Age"))
col78 <- rbind(table(DF$Edu), table_by("PurchaseInt", "Edu"), table_by("Relevance", "Edu"))

# should be the data you want
whole_table <- cbind(col12, col34, col56,col78)
whole_table

library(huxtable)
wt_hux <- as_hux(whole_table, add_colnames = TRUE, add_rownames = TRUE)
number_format(wt_hux)[-2,] <- "%.0f%%"
number_format(wt_hux)[2,]  <- "%.0f"
wt_hux[1, 1:2] <- c("", "Total")
wt_hux[2, 1]   <- "Total"
wt_hux <- insert_row(wt_hux, c("", "Total", "Sex", "", "Age", "", "Edu", "", ""))
colspan(wt_hux)[1, c(3, 5, 7)] <- c(2, 2, 3)
align(wt_hux)[1, c(3, 5, 7)] <- "center"
wt_hux <- insert_column(wt_hux, c("", "", "Total", "PurchaseInt", "", "", "", "Relevance", "", "", ""))
rowspan(wt_hux)[c(4, 8), 1] <- 4

bottom_border(wt_hux)[c(2,3, 7), ] <- 1 # for example

# should look roughly the way you want. You can print it to PDF or HTML:
wt_hux

想知道这是否可以包装成一个函数。我还不是很擅长编写 R 函数，但我很喜欢一个衬垫（因为我可能需要许多具有不同行和列的此类表，因为这只是一个示例）。

干杯， 格热谢克。