我想对数据集中的一堆自变量进行回归分析。预测变量很多,所以我不想全部写出来。是否有跨多个列的注释,所以我不必每次都键入?

我的尝试是这样做的(我的预测指标在第20至43列中):

modelAllHexSubscales = lm(HHdata$garisktot~HHdata[,20:43])


显然,这是行不通的,因为HHdata[,20:43]是数据矩阵,而我真的需要它来将数据视为HHdata[,20]+HHdata[,21]等。

最佳答案

首先粘贴列名称来生成公式。

f <- as.formula(paste('garisktot ~', paste(colnames(HHdata)[20:43], collapse='+')))
modelAllHexSubscales <- lm(f, HHdata)

08-24 17:00