We present a Bayesian hierarchical framework for an end-to-end analysis of peculiar velocity surveys, which makes explicit the inference problem of constraining cosmological parameters from redshift-independent distance indicators. We demonstrate our method for a Fundamental Plane-based survey. The essence of our approach is to work closely with observables (angular size, surface brightness, redshift, etc), through which we bypass the use of summary statistics by working with the probability distributions. The hierarchical approach improves upon the usual analysis in several ways. In particular, it allows a consistent analysis without having to make prior assumptions about cosmology during the calibration phase. Moreover, calibration uncertainties are automatically accounted for in parameter estimation. Results are presented for a new, fully analytic posterior marginalised over all latent variables, which we expect to allow for more principled analyses in upcoming surveys. A maximum a posteriori estimator is also given for peculiar velocities derived from Fundamental Plane data.

The standard model of cosmology is based on two unknown dark components that are uncoupled from each other. In this paper we investigate whether there is evidence for an interaction between these components of cold dark matter (CDM) and dark energy (DE). In particular, we focus on a minimal extension and reconstruct the interaction history at low-redshifts non-parametrically using a variation of the commonly used principal component analysis. Although we focus on the interaction in the dark sector, any significant deviation from the standard model that changes the expansion history of the Universe, should leave imprints detectable by our analysis. Thus, detecting signatures of interaction could also be indicative of other non-standard phenomena even if they are not the results of the interaction. It is thus interesting to note that the results presented in this paper do not provide support for the interaction in the dark sector, although the uncertainty is still quite large. In so far as interaction is present but undetectable using current data, we show from a Fisher forecast that forthcoming LSST and DESI surveys will be able to constrain a DM-DE coupling at $20\%$ precision --- enough to falsify the non-interacting scenario, assuming the presence of a modest amount of interaction.

Parameters that quantify the acceleration of cosmic expansion are conventionally determined within the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) model, which fixes spatial curvature to be homogeneous. Generic averages of Einstein's equations in inhomogeneous cosmology lead to models with non-rigidly evolving average spatial curvature, and different parametrizations of apparent cosmic acceleration. The timescape cosmology is a viable example of such a model without dark energy. Using the largest available supernova data set, the JLA catalogue, we find that the timescape model fits the luminosity distance-redshift data with a likelihood that is statistically indistinguishable from the standard spatially flat $\Lambda$ cold dark matter cosmology by Bayesian comparison. In the timescape case cosmic acceleration is non-zero but has a marginal amplitude, with best-fitting apparent deceleration parameter, $q_0=-0.043^{+0.004}_{-0.000}$. Systematic issues regarding standardization of supernova light curves are analysed. Cuts of data at the statistical homogeneity scale affect light curve parameter fits independent of cosmology. A cosmological model dependence of empirical changes to the mean colour parameter is also found. Irrespective of which model ultimately fits better, we argue that as a competitive model with a non-FLRW expansion history, the timescape model may prove a useful diagnostic tool for disentangling selection effects and astrophysical systematics from the underlying expansion history.