Bayesian global analysis of neutrino oscillation data

We perform a Bayesian analysis of current neutrino oscillation data. When estimating the oscillation parameters we find that the results generally agree with those of the ¬χ² method, with some differences involving s²₂₃ and CP-violating effects. We discuss the additional subtleties caused by the circular nature of the CP-violating phase, and how it is possible to obtain correlation coefficients with s²₂₃. When performing model comparison, we find that there is no significant evidence for any mass ordering, any octant of s²₂₃ or a deviation from maximal mixing, nor the presence of CP-violation.