of the EdS model.

The results on the clustering of QSOs derived with this particular approach are consistent with the ones obtained with the scrambling of the redshifts. In the following all the computations will be carried out with the latter method.

First, we calculate the TPCF integrated over a sphere, $\bar{\xi}(r)$, as a function of the sphere radius $r$, for the three sub-samples separately (DSS, USNO, GSC) and the total sample (AERQS). Fig. \[totflat\] reports the results for the EdS universe, while Fig. \[totlambda\] refers to a $\Lambda$ universe. The error bars in Fig. \[totflat\] and \[totlambda\] represent the 1$\sigma$ interval for $\bar{\xi}(r)$ and are obtained by assuming a Poisson distribution [@gehrels].

In order to investigate the possible presence of a spurious clustering signal at large scales, we have computed the angular TPCF binned in intervals of 3 degrees (corresponding to $\sim 
14.6 h^{-1}$ Mpc comoving). Fig. \[angclust\] shows the absence of any significant bias on the large scales sampled by AERQS, up to 150 $h^{-1}$ Mpc.

The signal shown in Fig. \[totflat\] and \[totlambda\] for separations smaller than $15 h^{-1}$ Mpc is due to 25 and 28 QSO pairs for the EdS and $\Lambda$ models, respectively. For a 