Hi Shabnam,
This sounds a bit like rstanarm issue #551, so I would try to write out the formula explicitly (will be very long in your case) or at least avoid reformulate() and try to use dot notation (e.g., y ~ .), but I’m not sure if this will actually help (given that rstanarm issue). Alternatively, it could be something like Error message from projpred with binomial regression model fitted in rstanarm, so do you have NAs in your dataset?
In any case, I currently think this is related to rstanarm, so using brms or something else should avoid this.