在R中创建双变量表格报错：非数值参数传入二元运算符的解决方法

错误原因 #

你遇到的non numeric argument to binary operator错误，核心问题有两个：

1. 公式写法错误：subgroup1+subgroup2~... 里的+是二元运算符，要求两侧为数值型，但你的亚组变量是分类变量，触发了类型不匹配的错误。
2. 代码语法不完整：write.csv、summary、univariateTable的调用都未闭合括号，直接导致语法解析失败。

修正方案（基于tableone包，这是做基线表最常用的工具包） #

先确认加载tableone包，再根据你的亚组类型选择对应方法：

情况1：两个亚组是交叉分层维度（比如性别+年龄组）

先把两个亚组变量合并成一个交叉分层变量，再进行分析：

library(tableone)

# 合并两个亚组变量为一个分层变量（假设你的数据集叫data）
data$combined_subgroup <- interaction(data$subgroup1, data$subgroup2)

# 创建分层表格：vars参数填你要分析的所有变量
tab <- CreateTableOne(
  vars = c("var1", "var2", "var3"), # 替换为你的实际变量名
  strata = "combined_subgroup",     # 用合并后的分层变量
  data = data
)

# 导出为CSV
write.csv(summary(tab), "stratified_basetable.csv", row.names = FALSE)

情况2：两个亚组是同一分组的不同水平（比如本来就是一个分组变量的两组）

直接用原始分组变量作为分层参数即可：

library(tableone)

tab <- CreateTableOne(
  vars = c("var1", "var2"),
  strata = "subgroup", # 假设原始分组变量是subgroup，包含你说的两个亚组水平
  data = data
)

write.csv(summary(tab), "basetable.csv", row.names = FALSE)

额外注意事项 #

• 确保亚组变量是因子型：如果是字符型，先转成因子避免分析错误：data$subgroup <- as.factor(data$subgroup)
• 检查变量名拼写：确保代码里的变量名和数据集中的完全一致
• 如果你用的是其他包的univariateTable函数，核心逻辑不变：不要在公式左边用+连接分类变量，而是通过strata参数指定分层变量，或者先合并分层维度

内容的提问来源于stack exchange，提问作者idkwtfisgoingon