Maximum likelihood estimation and uncertainty quantification for Gaussian process approximation of deterministic functions