Как автоматически включить все двухсторонние взаимодействия в модели glm в R

Предположим, вы используете R и имеете данные, хранящиеся в кадре данных, M. Тогда я знаю, что

g <- glm(Y ~ ., data=M) 

будет автоматически соответствовать модели, где Y - зависимые переменные, а все остальные столбцы M являются предикторами. Есть ли аналогично простой способ дополнительно включить каждое двустороннее взаимодействие?