This thesis presents the simulation of approximated Gaussian process autoregressive models. Gaussian process models are a Bayesian nonparametric regression method, the main advantage of which is the quantification of uncertainty in closed-form. However, the closedform solution for the marginal likelihood results in a cubic computational complexity with respect to the …