公式传参引号困境：R语言自定义lmer函数无法接收非字符串公式

你的第二个函数有两个关键错误导致它无法正常工作，咱们一步步拆解：

错误1：错误处理了rhs参数的表达式捕获 #

你用quote(rhs)时，只是捕获了符号rhs本身，而不是用户传入的实际表达式（Species + (1 + Petal.Length | Petal.Width)）。要捕获传入的真实表达式，得用substitute(rhs)，再把它转换成字符串，才能和因变量拼接成完整公式。

错误2：lapply调用lmer的语法错误 #

第二个函数里的这行代码：

models <- lapply(models, function(x) lmer, data=df, REML=FALSE)

并没有实际调用lmer函数，只是把lmer函数本身当成返回值放进了列表。正确写法是直接把lmer作为lapply的第三个参数（lapply会自动把每个公式传递给它），或者在匿名函数里显式调用lmer(x, data=df, REML=FALSE)。

修正后的函数 #

下面是修复后的get_pvals_rhs函数，解决了上述问题：

library(lme4)
library(lazyeval)

get_pvals_rhs <- function(df, cols, rhs) {
  cols <- substitute(cols)
  col_pos <- setNames(as.list(seq_along(df)), names(df))
  pos <- eval(cols, col_pos)
  
  # 捕获传入的rhs表达式并转为字符串
  rhs_str <- deparse(substitute(rhs))
  formulas <- lapply(pos, function(x) formula(paste(colnames(df[x]), "~", rhs_str)))
  
  # 正确调用lmer，传递每个公式
  models <- lapply(formulas, lmer, data=df, REML=FALSE)
  
  tvals <- lapply(models, function(x) data.frame(coef(summary(x)))[c(2), ]$t.value)
  pvals <- lapply(tvals, function(x) { 2 * (1 - pnorm(abs(x))) })
  
  return(unlist(pvals))
}

测试验证 #

现在运行修正后的函数，就能得到和第一个函数完全一致的结果：

fails_fixed <- get_pvals_rhs(iris, Sepal.Length:Sepal.Width, Species + (1 + Petal.Length | Petal.Width))
fails_fixed

# 和第一个函数的结果对比验证
works <- get_pvals_qrhs(iris, Sepal.Length:Sepal.Width, "Species + (1 + Petal.Length | Petal.Width)")
all.equal(fails_fixed, works) # 会返回TRUE，说明结果一致

额外小技巧：如果想避免字符串拼接的繁琐，可以用reformulate函数直接组合公式，可读性更好，把生成公式的行替换成下面这行即可：

formulas <- lapply(pos, function(x) reformulate(rhs_str, response = colnames(df[x])))

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