We present the results of a Bayesian analysis of solar neutrino data in terms of ν e → ν µ,τ and ν e → ν s oscillations, where ν s is a sterile neutrino. We perform a Rates Analysis of the rates of solar neutrino experiments, including the first SNO CC result, and spectral data of the CHOOZ experiment, and a Global Analysis that takes into account also the Super-Kamiokande day and night electron energy spectra. We show that the Bayesian analysis of solar neutrino data does not suffer any problem from the inclusion of the numerous bins of the CHOOZ and Super-Kamiokande electron energy spectra and allows to reach the same conclusions on the favored type of neutrino transitions and on the determination of the most favored values of the oscillation parameters in both the Rates and Global Analysis. Our Bayesian analysis shows that ν e → ν s transitions are strongly disfavored with respect to ν e → ν µ,τ transitions. In the case of ν e → ν µ,τ oscillations, the Large Mixing Angle region is favored by the data (86% probability), the LOW region has some small chance (13% probability), the Vacuum Oscillation region is almost excluded (1% probability) and the Small Mixing Angle region is practically excluded (0.01% probability). We calculate also the marginal posterior probability distributions for tan 2 ϑ and ∆m 2 in the case of ν e → ν µ,τ oscillations and we show that the data imply large mixing almost with certainty and large values of ∆m 2 are favored (2 × 10 −6 eV 2 < ∆m 2 < 10 −3 eV 2 with 86% probability). We present also the results of a standard least-squares analysis of solar neutrino data and we show that the standard goodness of fit test is not able to reject pure ν e → ν s transitions. The likelihood ratio test, which is insensitive to the number of bins of the CHOOZ and Super-Kamiokande energy spectra, allows to reject pure ν e → ν s transitions in favor of ν e → ν µ,τ transitions only in the Global Analysis.