We use gwpopulation [1]} and dynesty [2]} within the Bilby framework [3]} to sample the likelihoods in Eq. (REF ) and Eq. (REF ).
Parameter estimation for each event is also performed using Bilby [4]}, following the settings used to produce GWTC-1 results [5]}.
GW data from LIGO and Virgo are obtained from the Gravitational Wave Open Science Center [6]}.
