It may be a basic question, but I could not seem to find a solution anywhere. If we have a data frame with 100 factors (call them `a1` to `a100`), how could a linear model be entered in R? I understand you could write 

```
lm(y~ a1*...*a100)

```

but if the names are long, it would take a long time to write them all out. Is there a faster way? For example, by referencing columns or something similar? Somewhat related, if I get a data table with a column name that involves parentheses (e.g. `y-max()`), how could I enter that? It reads as a function in R, but it is not.

I apologize if this has already been asked, but I could not seem to find an answer.

Thank you all in advance

---Edit---

Thank you for the answers. However, if I did want higher-order interaction terms, how would I accomplish that? Would I need to write a script or is there a smarter way?