We use the Iocco et al. (2015) compilation of 2,780 circular velocity measurements to analyze the Milky Way rotation curve. We find that the error bars for individual measurements are non-gaussian, and hence instead derive median statistics binned central circular velocity values and error bars from these data. We use these median statistics central values and error bars to fit the data to simple, few parameter, rotation curve functions. These simple functions are unable to adequately capture the significant small scale spatial structure in these data and so provide poor fits. We introduce and use the Gaussian Processes (GP) method to capture this small scale structure and use it to derive Milky Way rotation curves from the binned median statistics circular velocity data. The GP method rotation curves have significant small-scale spatial structure superimposed on a broad rise to galactocentric radius $R\approx7$ kpc and a decline at larger $R$. We use the GP method median statistics rotation curve to measure the Oort $A$ and $B$ constants and other characteristic rotation curve quantities. We study correlations in the residual circular velocities (relative to the GP method rotation curve). Along with other evidence for azimuthal asymmetry of the Milky Way circular rotation velocity field, we find that larger residual circular velocities seem to favor parts of spiral arms.

Zavarygin et al. (2018) compiled a list of 15 deuterium abundance measurements, discarded two because the remaining 13 measurements are then consistent with gaussianity, and found that the weighted mean baryon density (Omega_b h^2) determined from the 13 measurements is mildly discrepant (1.6sigma) with that determined from the Planck 2015 cosmic microwave background anisotropy data in a flat cosmogony. We find that a median statistic central estimate of Omega_b h^2 from all 15 deuterium abundance measurements is a more accurate estimate, is very consistent with Omega_b h^2 estimated from Planck 2015 data in a flat cosmogony, but is about 2sigma lower than that found in a closed cosmogonical model from the Planck 2015 data.