Pulsar timing array experiments have recently found evidence for a stochasticgravitational wave (GW) background, which induces correlations among pulsar timing residualsdescribed by the Hellings and Downs (HD) curve. Standard calculations of the HD correlation andits variance assume an isotropic background. However, for a background of astrophysical origin, weexpect a higher GW spectral density in directions with higher galaxy number densities. In acompanion paper, we have developed a theoretical formalism to account for the anisotropies arisingfrom large-scale galaxy clustering, leading to a new contribution to the variance of the HDcorrelation. In this subsequent work, we provide numerical results for this novel effect. Weconsider a GW background resulting from mergers of supermassive black hole binaries, and relatethe merger number density to the overdensity of galaxies. We find that anisotropies due tolarge-scale galaxy clustering lead to a standard deviation of the HD correlation at most atpercent level, remaining well below the standard contributions to the HD variance. Hence, thiskind of anisotropies in the GW source distribution does not represent a substantial contaminationto the correlations of timing residuals in present and future PTA surveys. Suitable statisticalmethods to extract the galaxy clustering signal from PTA data will be investigated in the future.