No Arabic abstract
Bayesian graphical models are an efficient tool for modelling complex data and derive self-consistent expressions of the posterior distribution of model parameters. We apply Bayesian graphs to perform statistical analyses of Type Ia supernova (SN Ia) luminosity distance measurements from the joint light-curve analysis (JLA) data set. In contrast to the $chi^2$ approach used in previous studies, the Bayesian inference allows us to fully account for the standard-candle parameter dependence of the data covariance matrix. Comparing with $chi^2$ analysis results, we find a systematic offset of the marginal model parameter bounds. We demonstrate that the bias is statistically significant in the case of the SN Ia standardization parameters with a maximal 6 $sigma$ shift of the SN light-curve colour correction. In addition, we find that the evidence for a host galaxy correction is now only 2.4 $sigma$. Systematic offsets on the cosmological parameters remain small, but may increase by combining constraints from complementary cosmological probes. The bias of the $chi^2$ analysis is due to neglecting the parameter-dependent log-determinant of the data covariance, which gives more statistical weight to larger values of the standardization parameters. We find a similar effect on compressed distance modulus data. To this end, we implement a fully consistent compression method of the JLA data set that uses a Gaussian approximation of the posterior distribution for fast generation of compressed data. Overall, the results of our analysis emphasize the need for a fully consistent Bayesian statistical approach in the analysis of future large SN Ia data sets.
We analyze the three-year SDSS-II Superernova (SN) Survey data and identify a sample of 1070 photometric SN Ia candidates based on their multi-band light curve data. This sample consists of SN candidates with no spectroscopic confirmation, with a subset of 210 candidates having spectroscopic redshifts of their host galaxies measured, while the remaining 860 candidates are purely photometric in their identification. We describe a method for estimating the efficiency and purity of photometric SN Ia classification when spectroscopic confirmation of only a limited sample is available, and demonstrate that SN Ia candidates from SDSS-II can be identified photometrically with ~91% efficiency and with a contamination of ~6%. Although this is the largest uniform sample of SN candidates to date for studying photometric identification, we find that a larger spectroscopic sample of contaminating sources is required to obtain a better characterization of the background events. A Hubble diagram using SN candidates with no spectroscopic confirmation, but with host galaxy spectroscopic redshifts, yields a distance modulus dispersion that is only ~20 - 40% larger than that of the spectroscopically-confirmed SN Ia sample alone with no significant bias. A Hubble diagram with purely photometric classification and redshift-distance measurements, however, exhibit biases that require further investigation for precision cosmology.
The standard cosmology strongly relies upon the Cosmological Principle, which consists on the hypotheses of large scale isotropy and homogeneity of the Universe. Testing these assumptions is, therefore, crucial to determining if there are deviations from the standard cosmological paradigm. In this paper, we use the latest type Ia supernova compilations, namely JLA and Union2.1 to test the cosmological isotropy at low redshift ranges ($z<0.1$). This is performed through a Bayesian selection analysis, in which we compare the standard, isotropic model, with another one including a dipole correction due to peculiar velocities. We find that the Union2.1 sample favors the dipole-corrected model, but the opposite happens for the JLA. Nonetheless, the velocity dipole results are in good agreement with previous analyses carried out with both datasets. We conclude that there are no significant indications for large anisotropic signals from nearby supernova compilations, albeit this test should be greatly improved with the upcoming cosmological surveys.
We perform a model independent reconstruction of the cosmic expansion rate based on type Ia supernova data. Using the Union 2.1 data set, we show that the Hubble parameter behaviour allowed by the data without making any hypothesis about cosmological model or underlying gravity theory is consistent with a flat LCDM universe having H_0 = 70.43 +- 0.33 and Omega_m=0.297 +- 0.020, weakly dependent on the choice of initial scatter matrix. This is in closer agreement with the recently released Planck results (H_0 = 67.3 +- 1.2, Omega_m = 0.314 +- 0.020) than other standard analyses based on type Ia supernova data. We argue this might be an indication that, in order to tackle subtle deviations from the standard cosmological model present in type Ia supernova data, it is mandatory to go beyond parametrized approaches.
The large spectral bandwidth and wide field of view of the Australian SKA Pathfinder radio telescope will open up a completely new parameter space for large extragalactic HI surveys. Here we focus on identifying and parametrising HI absorption lines which occur in the line of sight towards strong radio continuum sources. We have developed a method for simultaneously finding and fitting HI absorption lines in radio data by using multi-nested sampling, a Bayesian Monte Carlo algorithm. The method is tested on a simulated ASKAP data cube, and is shown to be reliable at detecting absorption lines in low signal-to-noise data without the need to smooth or alter the data. Estimation of the local Bayesian evidence statistic provides a quantitative criterion for assigning significance to a detection and selecting between competing analytical line-profile models.
Much of the cosmological utility thus far extracted from Type Ia supernovae (SNe Ia) relies on the assumption that SN~Ia peak luminosities do not evolve significantly with the age (local or global) of their stellar environments. Two recent studies have provided conflicting results in evaluating the validity of this assumption, with one finding no correlation between Hubble residuals (HR) and stellar environment age, while the other claims a significant correlation. In this Letter we perform an independent reanalysis that rectifies issues with the statistical methods employed by both of the aforementioned studies. Our analysis follows a principled approach that properly accounts for regression dilution and critically (and unlike both prior studies) utilises the Bayesian-model-produced SN environment age estimates (posterior samples) instead of point estimates. Moreover, the posterior is used as an informative prior in the regression. We find the Pearson correlation between the HR and local (global) age to be in excess of $4sigma$ ($3sigma$). Assuming there exists a linear relationship between HR and local (global) age, we find a corresponding slope of $-0.035 pm 0.007,mathrm{mag,Gyr}^{-1}$ ($-0.036 pm 0.007,mathrm{mag,Gyr}^{-1}$). We encourage further usage of our approach to examine possible cosmological implications of the HR and age correlation.