It is well known that sufficiently smooth, hyperbolic dynamical systems admit strong statistical descriptions e.g. limit laws such as a central limit theorem or large deviation principle. Given the existence of these laws one is led to the question of their stability: do nearby systems have similar statistical descriptions, and to what extent can one numerically approximate the statistics of any particular system? In this thesis such questions are investigated by building on the so-called functional analytic approach, and in particular the spectral perturbation theory of Keller and Liverani. For deterministic systems it is shown that the Keller-Liverani perturbation theory is compatible with the naive Nagaev-Guivarc’h method – the method used to obtain the aforementioned statistical limit laws – yielding a general framework for deducing the statistical stability of deterministic dynamical systems under a variety of perturbations. This theory is then applied to piecewise expanding maps in one and many dimensions, in addition to Anosov maps on tori. Of particular note is the development of new, efficient and rigorous numerical methods for the approximation of the statistical properties of multidimensional piecewise expanding maps and Anosov maps. In the second part of this thesis this program is begun again for random systems. Here there is no analogue of the Keller-Liverani perturbation theory, and so an appropriate random version of the theory is developed. This theory is then applied to smooth random expanding maps on the circle, and the stability of some basic statistical properties is deduced with respect to fiber-wise deterministic perturbations and a Fourier-analytic numerical method.