Bayesian Parameter Shift Rule in Variational Quantum Eigensolvers