Though debated, the existence of claimed large-scale anomalies in the CMB is not totally dismissed. In parallel to the debate over their statistical significance, recent work focussed on masks and secondary anisotropies as potential sources of these anomalies. In this work we investigate simultaneously the impact of the method used to account for masked regions and the impact of the ISW effect, which is the large-scale secondary anisotropy most likely to affect the CMB anomalies. In this sense, our work is an update of both Francis & Peacock 2010 and Kim et al. 2012. Our aim is to identify trends in CMB data with different mask treatments. We reconstruct the ISW signal due to 2MASS and NVSS galaxies. We account for missing data using the sparse inpainting technique of Abrial et al. 2008 and sparse inpainting of the CMB, LSS and ISW and find that it constitutes a bias-free reconstruction method suitable to study large-scale statistical isotropy and the ISW effect. We focus on three large-scale CMB anomalies: the low quadrupole, the quad/oct alignment, and the octopole planarity. After sparse inpainting, the low quadrupole becomes more anomalous, whilst the quad/oct alignment becomes less anomalous. The significance of the low quadrupole is unchanged after subtraction of the ISW effect, while the trend is that the quad/oct alignment has reduced significance, yet other hypotheses remain possible as well (e.g. exotic physics). Our results also suggest that both of these anomalies may be due to the quadrupole alone. The octopole planarity significance is also reduced after inpainting and after ISW subtraction, however, we do not find that it was very anomalous to start with. In the spirit of reproducible research, we make all codes and resulting products which constitute main results of this paper public here: http://www.cosmostat.org/anomaliesCMB.html.

One of the main challenges of modern cosmology is to understand the nature of dark energy. The Integrated Sachs-Wolfe (ISW) effect is sensitive to dark energy and presents an independent signature of dark energy in the absence of modified gravity and curvature. The ISW effect occurs on large scales, where cosmic variance is high and where there are large amounts of missing data in the CMB and large scale structure maps due to Galactic confusion. Moreover, existing methods in the literature often make strong assumptions about the statistics of the underlying fields or estimators. Together these effects can severely limit signal extraction. We review literature on the ISW effect, comparing statistical subtleties between existing methods, and identifying several limitations. We propose a novel method to detect and measure the ISW signal. This method assumes only that the primordial CMB field is Gaussian. It is based on a sparse inpainting method to reconstruct missing data and uses a bootstrap technique to avoid assumptions about the statistics of the estimator. It is a complete method, using three complementary statistical methods. We apply our method to Euclid-like simulations and show we can expect a \sim 7\sigma model-independent detection of the ISW signal with WMAP7-like data, even with missing data. Other tests return \sim 4.7\sigma detection levels for a Euclid-like survey, with levels independent from whether the galaxy field is normally or lognormally distributed. We apply our method to the 2 Micron All Sky Survey (2MASS) and WMAP7 CMB data and find detections in the 1.1 - 2.0\sigma range, as expected from simulations. As a by-product, we reconstruct the full-sky temperature ISW field due to 2MASS data. We have presented a novel technique, based on sparse inpainting and bootstrapping, which accurately detects and reconstructs the ISW effect.