r - 在 R 中使用变量获取 lmList()

Question

我有作为 data.table 读入的数据。它看起来像这样:

d <- data.table(TYPE = c('A', 'A', 'A', 'B', 'B', 'C', 'C'),
                CLASS = c(1, 2, 3, 1, 4, 2, 1),
                NUM = c(48, 23, 32, 56, 75, 45, 76))

如果我运行 fits <- lmList(NUM ~ TYPE | CLASS, data = d)但是，如果我尝试将列名存储到变量中(以减少硬编码)并像这样运行它，它就会运行:

Y <- 'NUM'
X <- 'TYPE'
grp <- 'CLASS'
fits <- lmList(Y ~ X | grp, data = d)

我收到这个错误:

Warning messages:
1: In storage.mode(v) <- "double" : NAs introduced by coercion
2: In lmList(Y ~ X | grp, data = d) :
  Fitting failed for 1 group(s), probably because a factor only had one level:
  contrasts can be applied only to factors with 2 or more levels

我也试过fits <- lmList(as.name(Y) ~ as.name(X) | as.name(grp), data = d)我得到这个错误:

Error in model.frame.default(formula = as.name(Y) ~ as.name(X) + as.name(grp),  : 
  object is not a matrix

我尝试的其他方法是:

fits <- lmList(eval(paste(Y ~ X | grp, data = d)))

结果是:

Error in UseMethod("lmList") : 
  no applicable method for 'lmList' applied to an object of class "character"

我在这里做错了什么？欢迎任何帮助或建议。谢谢!

更新
这是一个link至 lmList供引用。

Answer 1

我知道这晚了，但我刚遇到这个问题。我正在为 future 的搜索者回答这个问题。您必须将公式视为 paste() 中的文本字符串。尝试:

fits <- lmList(as.formula(paste(Y, " ~ ", X, " | grp")), data = d)