Quantifying changes in dipole polarity reversal frequency using a combination of geodynamo simulations and stochastic models