Como você deve ter adivinhado, esta questão é de grande interesse para a equipe do LIGO. Simultaneamente com a publicação do artigo que você mencionou anunciando a descoberta, a equipe do LIGO enviou vários documentos complementares com mais detalhes sobre a descoberta e previsões. Um deles aborda sua pergunta:
A taxa de fusões binárias de buracos negros inferidas a partir de observações avançadas do LIGO em torno de GW150914
O método de estimativa da taxa de eventos considera o GW150914 e outro evento significativamente mais fraco (e menos estatisticamente significativo). Eles consideram vários modelos de como a taxa de eventos pode depender das propriedades do sistema e perguntam o que as observações do GW150914 e do outro evento candidato implicam na taxa geral. Os resultados variam de modelo para modelo, mas eles escolheram modelos que achavam que poderiam apoiar comportamento astrofisicamente plausível. Como resumido em seu resumo:
2 - 53G p c- 3y r- 16 - 400G p c- 3y r- 12 - 400G p c- 3y r- 1
Observe que o artigo é enviado, não publicado, ou seja, ainda está em revisão por pares. Falando como alguém com experiência em tais cálculos, alguns aspectos do método parecem duvidosos para mim, então acho que vale a pena conferir o artigo em poucas semanas para revisões. Não é preciso uma metodologia sofisticada para ver que a ordem de magnitude aqui (de alguns a ~ 100 por gigaparsec cúbico por ano) está no estádio certo. Mas o artigo apresenta uma metodologia que poderia fazer estimativas e previsões mais detalhadas e mais precisas à medida que os dados se acumulam, por isso é importante garantir que a metodologia seja sólida.