My question is a more detailed version of the one found here, which elicited some good information but the question was never really answered.
From table 4 in a 2010 paper we see the estimated rate of BH-BH mergers ranged from $10^{-4}$ to $0.3 \rm \,Mpc^{-3} Myr^{-1}$. This corresponds to $10^{-10}$ to $3\times10^{-7} \rm\,Mpc^{-3} yr^{-1}$. They write:
For BH–BH inspirals, horizon distances of ...2187 Mpc are assumed. These distances correspond to a choice of ... 10M for BH mass.
In the left panel of figure 4 found in this 2016 paper by LIGO there is a plot of BH mass vs horizon distance. For 2015-2016 sensitivity levels and 10M mass, the horizon distance is 300 Mpc. So we need to multiply the above numbers by $300^3=2.7\times10^{7}$. This gives a range of $2.7\times10^{-3}$ to $8.1$ detectable mergers per year, given the prior assumptions that a typical signal would come from black holes with ~10 solar masses. Also, in the discovery paper we find:
We present the analysis of 16 days of coincident observations between the two LIGO detectors from September 12 to October 20, 2015.
Therefore the duration of the experiment was $16/365$ years. This gives an expected rate ranging from $1.18\times10^{-4}$ to $3.55\times10^{-1}$ mergers during the experiment. As far as I can tell, this can be taken as an upper bound, since not all of these events will really be detected (due to non-optimal location, coincident noise, etc).
Are there other factors to take into account?