We develop a conceptual framework and methodology to study the drivers of the quenching of galaxies, including the drivers of galactic conformity. The framework is centred on the statistic $\Delta$, which is defined as the difference between the observed star-formation state of a galaxy, and a prediction of its state based on an empirical model of quenching. In particular, this work uses the average quenching effects of stellar mass and local density to construct an empirical model of quenching. $\Delta$ is therefore a residual which reflects the effects of drivers of quenching not captured by stellar mass and local density, or so-called 'hidden variables'. Through a toy model, we explore how the statistical properties of $\Delta$ can be used to learn about the internal and external hidden variables which control the quenching of a sample of galaxies. We then apply this analysis to a sample of local galaxies and find that, after accounting for the average quenching effects of stellar mass and local density, $\Delta$ remains correlated out to separations of 3 Mpc. Furthermore, we find that external hidden variables remain important for driving the residual quenching of low-mass galaxies, while the residual quenching of high-mass galaxies is driven mostly by internal properties. These results, along with a similar analysis of a semi-analytical mock catalogue, suggest that it is necessary to consider halo-related properties as candidates for hidden variables. A preliminary halo-based analysis indicates that much of the correlation of $\Delta$ can be attributed to the physics associated with individual haloes.

We re-examine the observational evidence for large-scale (4 Mpc) galactic conformity in the local Universe, as presented in Kauffmann et al. (2013). We show that a number of methodological features of their analysis act to produce a misleadingly high amplitude of the conformity signal. These include a weighting in favour of central galaxies in very high-density regions, the likely misclassification of satellite galaxies as centrals in the same high-density regions, and the use of medians to characterize bimodal distributions. We show that the large-scale conformity signal in Kauffmann et al. clearly originates from a very small number of central galaxies in the vicinity of just a few very massive clusters, whose effect is strongly amplified by the methodological issues that we have identified. Some of these 'centrals' are likely misclassified satellites, but some may be genuine centrals showing a real conformity effect. Regardless, this analysis suggests that conformity on 4 Mpc scales is best viewed as a relatively short-range effect (at the virial radius) associated with these very large neighbouring haloes, rather than a very long-range effect (at tens of virial radii) associated with the relatively low-mass haloes that host the nominal central galaxies in the analysis. A mock catalogue constructed from a recent semi-analytic model shows very similar conformity effects to the data when analysed in the same way, suggesting that there is no need to introduce new physical processes to explain galactic conformity on 4 Mpc scales.