Single Gaussian process method for arbitrary tokamak regimes with a statistical analysis