Non-parametric reconstruction methods, such as Gaussian process (GP) regression, provide a model-independent way of estimating an underlying function and its uncertainty from noisy data. We demonstrate how GP-reconstruction can be used as a consistency test between a given data set and a specific model by looking for structures in the residuals of the data with respect to the model's best-fit. Applying this formalism to the Planck temperature and polarisation power spectrum measurements, we test their global consistency with the predictions of the base $\Lambda$CDM model. Our results do not show any serious inconsistencies, lending further support to the interpretation of the base $\Lambda$CDM model as cosmology's gold standard.

We modify the algorithm we proposed before in Aghamousa & Shafieloo (2015) on time delay estimation of the strong lens systems incorporating weighted cross-correlation and weighted summation of correlation coefficients. We show the high performance of this algorithm by applying it on Time Delay Challenge (TDC1) simulated data. We apply then our proposed method on the light curves of the lensed quasar SDSS J1001+5027 since this system has been well studied by other groups to compare our results with their findings. In this work we propose a new estimator namely "mirror" estimator along with a list of criteria for reliability test of estimation. Our mirror estimator results to $-117.1^{+7.1}_{-3.7}$ and $-117.1^{+7.2}_{-8.8}$ using simple Monte Carlo simulations and simulated light curves provided by Rathna Kumar et al. (2013) respectively. Although TDC1 simulations do not reflect the properties of SDSS J1001+5027 light curves, using these simulations results to smaller uncertainty which shows the higher quality observations can lead to substantially more precise time delay estimation. Our time delay estimation is in agreement with findings of the other groups for this strong lens system and the difference in the size of the error bars reflects the importance of appropriate light curve simulations.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

The observable time delays between the multiple images of strong lensing systems with time variable sources can provide us with some valuable information to probe the expansion history of the Universe. Estimation of these time delays can be very challenging due to complexities of the observed data where there are seasonal gaps, various noises and systematics such as unknown microlensing effects. In this paper we introduce a novel approach to estimate the time delays for strong lensing systems implementing various statistical methods of data analysis including the method of smoothing and cross-correlation. The method we introduce in this paper has been recently used in TDC0 and TDC1 Strong Lens Time Delay Challenges and has shown its power in reliable and precise estimation of time delays dealing with data with different complexities.

We present a novel approach to test the consistency of the cosmological models with multiband CMB data using a nonparametric approach. In our analysis we calibrate the REACT (Risk Estimation and Adaptation after Coordinate Transformation) confidence levels associated with distances in function space (\textit{confidence distances}) based on the Monte Carlo simulations in order to test the consistency of an assumed cosmological model with observation. To show the applicability of our algorithm, we confront Planck 2013 temperature data with concordance model of cosmology considering two different Planck spectra combination. In order to have an accurate quantitative statistical measure to compare between the data and the theoretical expectations, we calibrate REACT confidence distances and perform a bias control using many realizations of the data. Our results in this work using Planck 2013 temperature data put the best fit $\Lambda$CDM model at $95\% (\sim 2\sigma)$ confidence distance from the center of the nonparametric confidence set while repeating the analysis excluding the Planck $217 \times 217$ GHz spectrum data, the best fit $\Lambda$CDM model shifts to $70\% (\sim 1\sigma)$ confidence distance. The most prominent features in the data deviating from the best fit $\Lambda$CDM model seems to be at low multipoles $ 18 < \ell < 26$ at greater than $2\sigma$, $\ell \sim 750$ at $\sim1$ to $2 \sigma$ and $\ell \sim 1800$ at greater than $2\sigma$ level. Excluding the $217\times217$ GHz spectrum the feature at $\ell \sim 1800$ becomes substantially less significance at $\sim1$ to $2 \sigma$ confidence level. Results of our analysis based on the new approach we propose in this work are in agreement with other analysis done using alternative methods.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

Estimation of the angular power spectrum is one of the important steps in Cosmic Microwave Background (CMB) data analysis. Here, we present a nonparametric estimate of the temperature angular power spectrum for the Planck 2013 CMB data. The method implemented in this work is model-independent, and allows the data, rather than the model, to dictate the fit. Since one of the main targets of our analysis is to test the consistency of the $\Lambda$CDM model with Planck 2013 data, we use the nuisance parameters associated with the best-fit $\Lambda$CDM angular power spectrum to remove foreground contributions from the data at multipoles $\ell \geq50$. We thus obtain a combined angular power spectrum data set together with the full covariance matrix, appropriately weighted over frequency channels. Our subsequent nonparametric analysis resolves six peaks (and five dips) up to $\ell \sim1850$ in the temperature angular power spectrum. We present uncertainties in the peak/dip locations and heights at the $95\%$ confidence level. We further show how these reflect the harmonicity of acoustic peaks, and can be used for acoustic scale estimation. Based on this nonparametric formalism, we found the best-fit $\Lambda$CDM model to be at $36\%$ confidence distance from the center of the nonparametric confidence set -- this is considerably larger than the confidence distance ($9\%$) derived earlier from a similar analysis of the WMAP 7-year data. Another interesting result of our analysis is that at low multipoles, the Planck data do not suggest any upturn, contrary to the expectation based on the integrated Sachs-Wolfe contribution in the best-fit $\Lambda$CDM cosmology.

The Planck mission, designed for making measurements of the cosmic microwave background (CMB) radiation with unprecedented accuracy and angular resolution, is expected to release its first data set in the near future. For the first time in the CMB history, extensive measurements of the CMB polarization will be made available for the entire sky. Such precise and rich data are expected to contain a great wealth of information about the Universe. The information in the CMB data is conveniently represented in terms of angular power spectra for temperature and polarization. A proper estimation of these CMB power spectra from data is the first step in making inferences about the Universe and, in particular, cosmological parameters that govern the Universe. In this paper, we provide forecasts for the $TT$, $EE$, and $TE$ angular power spectra for the Planck mission. Our forecasts are made using synthetic data based on the best-fit $\Lambda$CDM model while conforming to the characteristics and parameters of the Planck mission. We have analysed such synthetic data sets using a nonparametric function estimation formalism that is otherwise asymptotically model-independent. Indeed, the power spectra estimated from such synthetic data turn out to be sufficiently close to the corresponding $\Lambda$CDM spectra. Our analysis further indicates that the three angular power spectra are determined sufficiently well for $2 \le l \lesssim l_{max}$, where $l_{max} = 2500 {($TT$),} 1377 {($EE$), and} 1727 {($TE$)}$ respectively. A signature of reionization is evident in the $EE$ and $TE$ power spectra in the form of a bump at low $l$s ...

What do the data, as distinguished from cosmological models, tell us about cosmological parameters that determined the model of the universe? In this paper, we address this question in the context of the WMAP angular power spectra for the cosmic microwave background radiation. Nonparametric methods are ideally suited for this purpose because they are model-independent by construction, and therefore allow inferences that are as data-driven as possible. Our analysis is based on a nonparametric fit to the WMAP 7-year power spectrum data, with uncertainties characterized in the form of a high-dimensional confidence set centered at this fit. For the purpose of making inferences about cosmological parameters, we have devised a sampling method to explore the projection of this confidence set around the nonparametric fit, into the space of seven cosmological parameters Omega_b, Omega_c, Omega_Lambda, Omega_k, H_0, n_s, tau). Our sampling method is justified by its computational simplicity, and validated by the fact that well-known degeneracies in this cosmological parameter space are correctly reproduced in our results. Our results show that cosmological parameters are not as tightly constrained by these data alone. However, incorporating additional prior information in the analysis (e.g., constraining the values of H_0 or Omega_k) leads to tighter confidence intervals on parameters that are consistent with uncertainties estimated by the mainstream parametric methods.

Using a nonparametric function estimation methodology, we present a comparative analysis of the WMAP 1-, 3-, 5-, and 7-year data releases for the CMB angular power spectrum with respect to the following key questions: (a) How well is the power spectrum determined by the data alone? (b) How well is the $\Lambda$CDM model supported by a model-independent, data-driven analysis? (c) What are the realistic uncertainties on peak/dip locations and heights? Our results show that the height of the power spectrum is well determined by data alone for multipole l approximately less than 546 (1-year), 667 (3-year), 804 (5-year), and 842 (7-year data). We show that parametric fits based on the $\Lambda$CDM model are remarkably close to our nonparametric fits in $l$-regions where data are sufficiently precise. In contrast, the power spectrum for an H$\Lambda$CDM model gets progressively pushed away from our nonparametric fit as data quality improves with successive data realizations, suggesting incompatibility of this particular cosmological model with respect to the WMAP data sets. We present uncertainties on peak/dip locations and heights at the 95% ($2 \sigma$) level of confidence, and show how these uncertainties translate into hyperbolic "bands" on the acoustic scale ($l_A$) and peak shift ($\phi_m$) parameters. Based on the confidence set for the 7-year data, we argue that the low-l up-turn in the CMB power spectrum cannot be ruled out at any confidence level in excess of about 10% ($\approx 0.12 \sigma$). Additional outcomes of this work are a numerical formulation for minimization of a noise-weighted risk function subject to monotonicity constraints, a prescription for obtaining nonparametric fits that are closer to cosmological expectations on smoothness, and a method for sampling cosmologically meaningful power spectrum variations from the confidence set of a nonparametric fit.