The code is freely available (https://github.com/singmann/gumbel-reanalysis), so feel free to include it in your package. I think the main question is how to hand-over a multinomial response to brms. It is trivial if data is disaggreagted, but less so with aggregated data that can speed-up fitting considerably.