The current accelerated expansion of the Universe remains one of the most intriguing topics in modern cosmology, driving the search for innovative statistical techniques. Recent advancements in machine learning have significantly enhanced its application across various scientific fields, including physics, and particularly cosmology, where data analysis plays a crucial role in problem-solving. In this work, a non-parametric regression method with Gaussian processes is presented along with several applications to reconstruct some cosmological observables, such as the deceleration parameter and the dark energy equation of state, in order to contribute some information that helps to clarify the behavior of the Universe. It was found that the results are consistent with λCDM and the predicted value of the Hubble parameter at redshift zero is H0=68.798±6.340(1σ)kms−1Mpc−1.