广义估计方程自己工作但不在函数内R

问题描述

我正在尝试编写一个函数来使用 geepack 包运行 GEE。它“单独”工作正常，但不在函数内，请参见下面的示例：

library(geepack)
library(pstools)

df <- data.frame(study_id = c(1:20),leptin = runif(20),insulin = runif(20),age = runif(20,min = 20,max = 45),sex = sample(c(0,1),size = 20,replace = TRUE))

#Works
geepack::geeglm(leptin ~ insulin + age + sex,id = study_id,data = df)

#Doesn't work
model_function_covariates_gee <- function(x,y) {
  M1 <- paste0(x,"~",y,"+ age + sex")
  M1_fit <- geepack::geeglm(M1,data = df)
  s <- summary(M1_fit)
  return(s)
}

model_function_covariates_gee("leptin","insulin")

错误信息：

Error in mcall$formula[3] <- switch(match(length(sformula),c(0,2,3)),: 
  incompatible types (from language to character) in subassignment type fix

有谁知道这是为什么？我已经摆弄过它，但无法改变它。提前致谢。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）