fit_bt_model function

Fit a Bradley–Terry model with sirt and fallback to BradleyTerry2