We compute the two-dimensional correlation functions of the binary black hole coalescence detections in LIGO-Virgo first and second observation runs .
The observed correlation function is compared to two reference functions obtained from artificial maps to test for any possible excess correlation of the binary black hole coalescence events at different angular scales .
No excess correlation at any angular scale is found .
The power-law slope of the correlation function is estimated to be \gamma = 2.24 \pm 0.33 at 3- \sigma confidence level , a value consistent with a random distribution of sources .