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.

We present an original phenomenological model to describe the evolution of galaxy number counts, morphologies, and spectral energy distributions across a wide range of redshifts (0.2<z<15) and stellar masses [Log10 M/Msun >6]. Our model follows observed mass and luminosity functions of both star-forming and quiescent galaxies, and reproduces the redshift evolution of colors, sizes, star-formation and chemical properties of the observed galaxy population. Unlike other existing approaches, our model includes a self-consistent treatment of stellar and photoionized gas emission and dust attenuation based on the BEAGLE tool. The mock galaxy catalogs generated with our new model can be used to simulate and optimize extragalactic surveys with future facilities such as the James Webb Space Telescope (JWST), and to enable critical assessments of analysis procedures, interpretation tools, and measurement systematics for both photometric and spectroscopic data. As a first application of this work, we make predictions for the upcoming JWST Advanced Deep Extragalactic Survey (JADES), a joint program of the JWST/NIRCam and NIRSpec Guaranteed Time Observations teams. We show that JADES will detect, with NIRCam imaging, thousands of galaxies at z>6, and tens at z>10 at m_AB<30 (5-sigma) within the 236 arcmin^2 of the survey. The JADES data will enable accurate constraints on the evolution of the UV luminosity function at z>8, and resolve the current debate about the rate of evolution of galaxies at z>8. Ready to use mock catalogs and software to generate new realizations are publicly available as the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) package.

We compile an updated list of 38 measurements of the Hubble parameter $H(z)$ between redshifts $0.07 \leq z \leq 2.36$ and use them to place constraints on model parameters of constant and time-varying dark energy cosmological models, both spatially flat and curved. We use five models to measure the redshift of the cosmological deceleration-acceleration transition, $z_{\rm da}$, from these $H(z)$ data. Within the error bars, the measured $z_{\rm da}$ are insensitive to the model used, depending only on the value assumed for the Hubble constant $H_0$. The weighted mean of our measurements is $z_{\rm da} = 0.72 \pm 0.05\ (0.84 \pm 0.03)$ for $H_0 = 68 \pm 2.8\ (73.24 \pm 1.74)$ km s$^{-1}$ Mpc$^{-1}$ and should provide a reasonably model-independent estimate of this cosmological parameter. The $H(z)$ data are consistent with the standard spatially-flat $\Lambda$CDM cosmological model but do not rule out non-flat models or dynamical dark energy models.

We construct error distributions for a compilation of 232 Large Magellanic Cloud (LMC) distance moduli values from de Grijs (2014) that give an LMC distance modulus of (m-M)_{0}=18.49 plus/minus 0.13 mag (median and 1 sigma symmetrized error). Central estimates found from weighted mean and median statistics are used to construct the error distributions. The weighted mean error distribution is non-Gaussian --- flatter and broader than Gaussian --- with more (less) probability in the tails (center) than is predicted by a Gaussian distribution; this could be the consequence of unaccounted-for systematic uncertainties. The median statistics error distribution, which does not make use of the individual measurement errors, is also non-Gaussian --- more peaked than Gaussian --- with less (more) probability in the tails (center) than is predicted by a Gaussian distribution; this could be the consequence of publication bias and/or the non-independence of the measurements. We also construct the error distributions of 247 SMC distance moduli values from de Grijs (2015). We find a central estimate of (m-M)_{0}=18.94 plus/minus 0.14 mag (median and 1 sigma symmetrized error), and similar probabilities for the error distributions.

We construct the error distribution of $\rm{^{7}Li}$ abundance measurements for 66 observations (with error bars) used by Spite12 that give $\rm{A(Li)}=2.21 \pm 0.065$ (median and 1$\sigma$ symmetrized error). This error distribution is somewhat non-Gaussian, with larger probability in the tails than is predicted by a Gaussian distribution. The 95.4% confidence limits are 3.0$\sigma$ in terms of the quoted errors. We fit the data to four commonly used distributions: Gaussian, Cauchy, Student's t, and double exponential with the center of the distribution found with both weighted mean and median statistics. It is reasonably well described by a widened $n=8$ Student's $t$ distribution. Assuming Gaussianity, the observed A(Li) is 6.5$\sigma$ away from that expected from standard Big Bang nucleosynthesis given the $Planck$ observations Coc2014. Accounting for the non-Gaussianity of the observed A(Li) error distribution reduces the discrepancy to 4.9$\sigma$, which is still significant.

We present median statistics central values and ranges for 12 cosmological parameters, using 582 measurements (published during 1990-2010) collected by Croft & Dailey (2011). On comparing to the recent Planck collaboration Ade et al. 2013 estimates of 11 of these parameters, we find good consistency in nine cases.

Weighted mean and median statistics techniques are used to combine 23 independent lower redshift, $z<1.04$, Hubble parameter, $H(z)$, measurements and determine binned forms of $H(z)$. When these are combined with 5 higher redshift, $1.3\leqslant z \leqslant 2.3$, $H(z)$ measurements the resulting constraints on cosmological parameters, of three cosmological models, that follow from the weighted-mean binned data are almost identical to those derived from analyses using the 28 independent $H(z)$ measurements. This is consistent with what is expected if the lower redshift measurements errors are Gaussian. Plots of the binned weighted-mean $H(z)/(1+z)$ versus $z$ data are consistent with the presence of a cosmological deceleration-acceleration transition at redshift $z_{\rm da}=0.74 \pm 0.05$ \citep{farooq3}, which is expected in cosmological models with present-epoch energy budget dominated by dark energy as in the standard spatially-flat $\Lambda$CDM cosmological model.