Best-fit parameters obtained after performing the likelihood analysis for the SBG catalogue on our datasets. The value of f_aniso and θ allowing to maximize the likelihood for each dataset are plotted against each other for the astrophysical model considered in Fig. 8 and two GMF models (see legend). The 1σ ellipse reported in Caccianiga et al. (2019) from Auger data is shown. The individual distributions of the different quantities plotted are shown on top of the coordinate axis.