Можно использовать ярлык для формулы вlm()

m <- matrix(rnorm(100), ncol=5)
lm(m[,1] ~ m[,2:5]

здесь было бы так же, как

lm(m[,1] ~ m[,2] + m[,3] + m[,4] + m[,5]

но в случае, когда переменные не одного уровня (по крайней мере, это мое предположение на данный момент), это не работает, и я получаю ошибку:

Error in model.frame.default(formula = hm[, 1] ~ hm[, 2:4], drop.unused.levels = TRUE) : 
  invalid type (list) for variable 'hm[, 2:4]'

Данные (хм):

     N cor.distance switches  time
1   50   0.04707842        2 0.003
2  100  -0.10769441        2 0.004
3  200  -0.01278359        2 0.004
4  300   0.04229509        5 0.008
5  500  -0.04490092        6 0.010
6 1000   0.01939561        4 0.007

Есть ли еще какой-нибудь способ, чтобы избежать написания длинной формулы?