R语言自定义geeglm回归函数运行报unexpected ','错误排查

报错来自回归公式构造逻辑错误：geeglm的聚类id参数是独立于模型公式的函数入参，原代码错误把id参数拼接进了模型公式字符串，导致生成的公式末尾多出非法逗号，R解析公式时遇到不符合语法的逗号直接抛出错误。
原代码拼接出的错误公式内容为：

outcome ~ var1 + conf1+conf2 , id = id

R公式语法中~后仅允许出现变量项、交互项、函数项，不允许出现逗号分隔的参数定义，因此解析到逗号位置直接触发报错。

library(geepack)
# 示例测试数据
id<-rep(c(1:100),times=3)
df1<-data.frame(id)
df1$var1 <- sample(500:1000, length(df1$id))
df1$conf1 <- sample(500:1000, length(df1$id))
df1$conf2 <- sample(500:1000, length(df1$id))
df1$outcome <- sample(0:1, length(df1$id), replace=TRUE)

getReg<-function(varname, outcomename, confounderlist,  pat_id, data){
  data.sub <- na.omit(subset(data, select = c(outcomename, varname, confounderlist, pat_id)))
  # 仅构造符合语法的回归公式，不混入其他函数参数
  formula1<-as.formula(
    paste(outcomename,"~",varname,  "+", paste(confounderlist, collapse="+"))
  )
  # 聚类id作为geeglm的独立参数传入，不写入公式
  fit1<-geeglm(
    formula = formula1,
    data = data.sub,
    family = binomial(link="logit"),
    id = data.sub[[pat_id]]
  )
  # 显式返回拟合的模型对象
  return(fit1)
}

# 测试调用
model_res <- getReg(varname="var1", outcomename="outcome",  confounderlist=c("conf1", "conf2"), pat_id="id", data=df1)
# 查看模型结果
summary(model_res)

• 拆分公式构造与函数传参逻辑：模型公式仅保留「因变量~核心自变量+协变量」的合法结构，不再拼接id相关内容
• 用*as.formula()*将拼接的字符串显式转换为R公式对象，避免隐式转换带来的潜在兼容问题
• 聚类id通过data.sub[[pat_id]]从清洗后的数据集提取，作为geeglm的独立参数传入，完全和公式逻辑分离
• 补充显式返回语句return(fit1)，原函数未设置返回值，运行后无法直接获取拟合好的模型对象
• 补全geeglm的family参数名，避免位置传参可能带来的参数匹配错误

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