Crear automáticamente fórmulas para todos los modelos lineales posibles

Digo que tengo un conjunto de entrenamiento en un marco de datostrain con columnasColA, ColB, ColC, etc. Una de estas columnas designa una clase binaria, digamos columnaClass, con valores "sí" o "no".

Estoy probando algunos clasificadores binarios, por ejemplo

library(klaR)
mynb <- NaiveBayes(Class ~ ColA + ColB + ColC, train)

Me gustaría ejecutar el código anterior en un bucle, generando automáticamente todas las combinaciones posibles de columnas en la fórmula, es decir

mynb <- append(mynb, NaiveBayes(Class ~ ColA, train)
mynb <- append(mynb, NaiveBayes(Class ~ ColA + ColB, train)
mynb <- append(mynb, NaiveBayes(Class ~ ColA + ColB + ColC, train)
...
mynb <- append(mynb, NaiveBayes(Class ~ ColB + ColC + ColD, train)
...

Cómo puedo generar automáticamente fórmulas para cada posible modelo lineal que involucre columnas de un marco de datos?

Respuestas a la pregunta(3)

Su respuesta a la pregunta