Likelihood-free Model Selection in Cosmic Reionization with Three-dimensional Tomographic 21 cm Lightcone Images

We explore likelihood-free (aka simulation-based) Bayesian model selection to quantify model comparison analyses of reionisation scenarios. We iteratively train the 3D Convolutional Neural Network (CNN) on four toy EoR models based on 21cmFAST simulations with contrasting morphology to obtain summaries of the 21 cm lightcone. Within the pyDelfi framework, we replaced the Emcee sampler with MultiNest to integrate learnt posteriors and produce the Bayesian Evidence. We comfortably distinguish the model used to produce the mock data set in all cases. However, we struggle to produce accurate posterior distributions for outside-in reionisation models. After a variety of cross-checks and alternate analyses we discuss the flexibility of summarising models that differ from precisely the intended network training conditions as this should be more widely scrutinised before CNN can reliably analyse observed data.