我正在寻找一种简单的方法将作为变量名的函数参数传递给R函数中的公式。
测试数据集:
set.seed(4892)
df.pass <- data.frame("alfa"=sample(1:9, 100, replace=T), "beta"=sample(1:9, 100, replace=T),
"theta"=sample(1:9, 100, replace=T), "out"=runif(100, 0, 1))
示例分析(测试交互模型是否不同)将被用于函数:
lrtest(glm(out~alfa*beta, family = binomial("logit"), df.pass),
glm(out~alfa + beta, family = binomial("logit"), df.pass))
如果目标是使用任意变量名称和数据集创建解决上述问题的泛型函数invinteract
,那么将变量名称从function()
参数传递给公式术语的最简单方法是什么?对应于out
,alfa
和beta
的位置?
将原始变量名称插入到公式中不起作用,因为R会尝试将名称计算为对象而不会发现任何内容。
将字符串变量名直接插入公式中也不起作用。
是否有必要使用paste()
重建公式,还是有更直接的方式?
答案 0 :(得分:1)
glm
也接受字符串而不是公式。因此,您可以这样做:
mytest <- function(DF, y, x1, x2) {
lrtest(glm(sprintf("%s ~ %s * %s", y, x1, x2), family = binomial("logit"), DF),
glm(sprintf("%s ~ %s + %s", y, x1, x2), family = binomial("logit"), DF))
}
mytest(df.pass, "out", "alfa", "beta")