Recurrence models apply historical seismicity information to seismic hazard analysis. These models that play an important role in the obtained hazard curve are determined by their parameters. Recurrence parameters estimation has some features that lie in missing-data problems category. Thus, the observed data cannot be used directly to estimate model parameters. Furthermore discussion about results reliability and probable conservatism is impossible. The present study aims at offering an approach for Gutenberg-Richter parameters (�� and ��-values) estimation and determine their variation. Applying the proposed method to analyses of the heterogeneous data sets of seismic catalog, one would calculate valid estimates for recurrence parameters. This method has the capability to reflect all known sources of variability. The results of the case study clearly demonstrate applicability and efficiency of the proposed method, which can easily be implemented not only in advanced but also in practical seismic hazard analyses.