The next generation 21 cm surveys open a new window onto the early stages of cosmic structure formation and provide new insights about the Epoch of Reionization (EoR). However, the non-Gaussian nature of the 21 cm signal along with the huge amount of data generated from these surveys will require more advanced techniques capable to efficiently extract the necessary information to constrain the Reionization History of the Universe. In this paper we present the use of Bayesian Neural Networks (BNNs) to predict the posterior distribution for four astrophysical and cosmological parameters. Besides achieving state-of-the-art prediction performances, the proposed methods provide accurate estimation of parameters uncertainties and infer correlations among them. Additionally, we demonstrate the advantages of Normalizing Flows (NF) combined with BNNs, being able to model more complex output distributions and thus capture key information as non-Gaussianities in the parameter conditional density distribution for astrophysical and cosmological dataset. Finally, we propose novel calibration methods employing Normalizing Flows after training, to produce reliable predictions, and we demonstrate the advantages of this approach both in terms of computational cost and prediction performances.

Upcoming experiments such as Hydrogen Epoch of Reionization Array (HERA) and Square Kilometre Array (SKA) are intended to measure the 21cm signal over a wide range of redshifts, representing an incredible opportunity in advancing our understanding about the nature of cosmic Reionization. At the same time these kind of experiments will present new challenges in processing the extensive amount of data generated, calling for the development of automated methods capable of precisely estimating physical parameters and their uncertainties. In this paper we employ Variational Inference, and in particular Bayesian Neural Networks, as an alternative to MCMC in 21 cm observations to report credible estimations for cosmological and astrophysical parameters and assess the correlations among them.

These are the proceedings of the First Astrostatistics School: Bayesian Methods in Cosmology, held in Bogot\'a D.C., Colombia, June 9-13, 2014. The first astrostatistics school has been the first event in Colombia where statisticians and cosmologists from some universities in Bogot\'a met to discuss the statistic methods applied to cosmology, especially the use of Bayesian statistics in the study of Cosmic Microwave Background (CMB), Baryonic Acoustic Oscillations (BAO), Large Scale Structure (LSS) and weak lensing.

There is strong observational evidence for the presence of large-scale magnetic fields MF in galaxies and clusters, with strength $\sim \mu$G and coherence lenght on the order of Kpc. However its origin remains as an outstanding problem. One of the possible explanations is that they have been generated in the early universe. Recently, it has been proposed that helical primordial magnetic fields PMFs, could be generated during the EW or QCD phase transitions, parity-violating processes and predicted by GUT or string theory. Here we concentrate on the study of two mechanisms to generate PMFs, the first one is the $\nu$MSM which triggers instability in the Maxwell's equations and leads to the generation of helical PMFs. The second one is the usual electroweak baryogenesis scenario. Finally, we calculate the exact power spectra of these helical PMFs and we show its role in the production of gravitational waves finding a scale-invariant on large scales and an oscillatory motion (damping) for $k\eta \gg 1$