title: Using column numbers not names in lm()
body: <p>Instead of something like <code>lm(bp~height+age, data=mydata)</code> I would like to specify the columns by number, not name.</p>
  
  <p>I tried <code>lm(mydata[[1]]~mydata[[2]]+mydata[[3]])</code> but the problem with this is that, in the fitted model, the coefficients are named <code>mydata[[2]]</code>, <code>mydata[[3]]</code> etc, whereas I would like them to have the real column names. </p>
  
  <p>Perhaps this is a case of not having your cake and eating it, but if the experts could advise whether this is possible I would be grateful</p>