我想对数据集中的一堆自变量进行回归分析。预测变量很多,所以我不想全部写出来。是否有跨多个列的注释,所以我不必每次都键入?
我的尝试是这样做的(我的预测指标在第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)