Is there, either for the lm() function or for some other function for linear regression, an **argument** such that the reference group can be set to **always** be the biggest group rather than the alphabetical/numerical default in lm()?

As this's often done in stats, I'm thinking I somehow keep missing it when I search the documentation or that I'm looking in the wrong places. Any help would be appreciated! 

Below, even when in a UDF, is what I'd like **NOT** **to have to keep doing**. 

```
mtcars # load dataset 
mtcars <- mtcars[1:31, ]  # remove a now so that there is a single biggest group
lm(mpg ~ gear+carb+disp, data = mtcars ) # carb's group 1 is the reference by default 
mtcars$carb <- as.factor(mtcars$carb) 
mtcars <- within(mtcars, carb <- relevel(carb, ref = "4")) # set carb's group 4 as the reference
lm(mpg ~ gear+carb+disp, data = mtcars ) 

```