The number of interaction terms increases dramatically, when the number of variables increases. This makes it almost impossible to produce model results in a reasonable time when ncol(X) is large. Is it possible to change this default value ncol(X) >= 20 to something smaller?