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