How to automatically include all 2-way interactions in a glm model in R
===

Suppose you are using `R` and have data stored in a data frame, `M`. Then I know that 

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

will automatically fit a model where `Y` is the dependent variables and all other columns of `M` are the predictors. Is there an analogously simple way to additionally include every two way interaction? 