在 lm() 中使用列号而不是名称

Question

我想按数字而不是名称指定列，而不是像 lm(bp~height+age, data=mydata) 。

我尝试了 lm(mydata[[1]]~mydata[[2]]+mydata[[3]]) 但问题是，在拟合模型中，系数被命名为 < code>mydata[[2]]、mydata[[3]] 等，而我希望它们具有真实的列名称。

也许这是鱼与熊掌兼得的情况，但如果专家能建议这是否可能，我将不胜感激

Answer 1

lm(
    as.formula(paste(colnames(mydata)[1], "~",
        paste(colnames(mydata)[c(2, 3)], collapse = "+"),
        sep = ""
    )),
    data=mydata
)

您可以使用所需的索引数量来代替 c(2, 3) （不需要 for 循环）。

Answer 2

lm(mydata[,1] ~ ., mydata[-1])

我在 R 课程中发现的技巧是删除响应列，否则您会收到警告“本质上完美适合：摘要可能不可靠”。我不知道它为什么有效，它不符合文档。通常，我们保留响应栏。

托马斯之前的回答的简化版本：

lm(
    as.formula(paste(colnames(mydata)[1], "~ .")),
    data=mydata
)