feval using constant-only model?

Hi there - it seems that feval won't evaluate a glm produced by GeneralizedLinearModel.fit if it only includes a constant (y ~ 1). Granted, this is an easy model to solve manually, but I'm trying to make code that's universal and will accept whatever model with whatever distribution. Is there any way to have feval solve a constant-only model?

Thanks, Alden