<p>We carry out a Bayesian model comparison analysis to determine the dark matter model that best describes the Milky Way rotation curve, using four recent compilations of rotation-curve velocities with the Gaia data. We compare these data using three different baryonic models along with seven different dark matter models. Within the adopted modeling framework and Gaia-based rotation curve datasets, we find that the Einasto profile provides the preferred phenomenological fit compared to the NFW profile across most combinations of datasets and baryonic models. We also find that the cored dark profiles are preferred over the (cuspy) NFW profile. We also test MOND using three different interpolating functions and find that, within the implementations considered, the models provide poorer fits than the preferred dark matter profiles. Finally, among the different baryonic models considered, none is decisively favored over the others.</p>