Cómo incluir automáticamente todas las interacciones bidireccionales en un modelo glm en R

Supongamos que está utilizando R y tiene datos almacenados en un dataframe, M Entonces sé que

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

se ajustará automáticamente a un modelo donde Y es las variables dependientes y todas las otras columnas de M son los predictores. ¿Existe una forma análogamente simple de incluir, además, cada interacción bidireccional?

Puede hacer interacciones bidireccionales simplemente usando .*. y arbitrarias interacciones de n -way escribiendo .^n . formula(g) le dirá la versión expandida de la fórmula en cada uno de estos casos.