Data-driven reduced order modeling for parametric PDE eigenvalue problems using Gaussian process regression