Mam ramkę danych z powiedzmy N + 2 kolumny. Pierwszy to tylko daty (używane głównie do kreślenia później), drugi to zmienna, której odpowiedź na pozostałe kolumny N chciałbym obliczyć. Myślę, że musi być coś takiego

df = data.frame(y = 1:10, x1 = runif(10), x2 = rnorm(10))
fit = lm(y~df[,2:3],data=df)

To nie działa. Próbowałem też i nie

fit = lm(y~sapply(colnames(df)[2:3],as.name),data=df)

jakieś pomysły?