r - lapply 函数将单个和 + 参数传递给 LM

Question

我无法尝试将“+”参数传递给 lm。

我下面的两行代码适用于单个参数，例如:

model_combinations=c('.', 'Long', 'Lat', 'Elev')


lm_models = lapply(model_combinations, function(x) {
                               lm(substitute(Y ~ i, list(i=as.name(x))), data=climatol_ann)})

但如果我在 model_combinations 列表的末尾添加“Lat+Elev”，则相同的代码会失败，如下所示:

model_combinations=c('.', 'Long', 'Lat', 'Elev', 'Lat+Elev') 

Error in eval(expr, envir, enclos) : object 'Lat+Elev' not found

我已经扫描了帖子，但无法找到解决方案。

Answer 1

我通常发现使用 reformulate 通过字符串操作构造公式比尝试使用 substitute() 修改表达式更可靠/更容易理解，例如

model_combinations <- c('.', 'Long', 'Lat', 'Elev', 'Lat+Elev')
model_formulas <- lapply(model_combinations,reformulate,
                         response="Y")
lm_models <- lapply(model_formulas,lm,data=climatol_ann)

因为 reformulate 在字符串级别工作，如果元素本身是非原子的(例如 Lat+Elev)，它不会有问题。这里唯一棘手的情况是，如果您的 data 参数或变量是在某些不容易找到的环境中构造的，但传递显式 data 参数通常可以避免问题。

(您也可以使用 as.formula(paste(...)) 或 as.formula(sprintf(...)); 重新公式化() 只是一个方便的包装器。)