Fitting glm without explicit declaration of each covariate
This has already been answered in stackoverflow here and here.

The main idea is to create the formula from a string with as.formula:

xnam &lt;- paste("x", 1:25, sep="")
fmla &lt;- as.formula(paste("y ~ ", paste(xnam, collapse= "+")))