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

Question

我想按编号而不是名称指定列，而不是 lm(bp~height+age, data=mydata) 之类的东西。

我尝试了lm(mydata[[1]]~mydata[[2]]+mydata[[3]])，但问题在于，在拟合模型中，系数被命名为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循环）。