Fit a (frequentist) maximum likelihood model in Bambi

Question

Let's say I have a simple binomial regression model:

import pandas as pd
import bambi as bmb

df = pd.DataFrame(data={'group': ['A', 'B'], 'y': [3, 4], 'n': 10})

model = bmb.Model("p(y, n) ~ group", df, family="binomial")
idata = model.fit()

Is there a possibility to fit it using maximum likelihood estimate (frequentist without priors)? (I have big amounts of data and sampling is getting too slow.)