The Advanced LIGO and Advanced Virgo gravitational wave detectors have detected a population of binary black hole mergers in their first two observing runs. For each of these events we have been able to associate a potential sky location region represented as a probability distribution on the sky. Thus, at this point we may ask the question of whether this distribution agrees with the isotropic model of the universe or if there is any evidence of anisotropy. We perform Bayesian model selection between an isotropic and a simple anisotropic model, taking into account the anisotropic selection function caused by the underlying antenna patterns and sensitivity of the interferometers over the sidereal day. Using the ten binary mergers, we find a Bayes factor of $\sim~16:1$ in favour of the isotropic model over a simple pixel-based distribution on the sky. This suggests that the available data from the first two observing runs is better described with an isotropic model. However, the first detections were mostly poorly localised in the sky (before the Advanced Virgo joined the network), spanning large portions of the sky and hampering detection of potential anisotropy. Thus, it would be appropriate to repeat this analysis with events from the ongoing third LIGO observational run and with a more sophisticated cosmological model.