Abstract: Appropriate selection of Ground Motion Prediction Equation (GMPE) is one of the most important parts of Probabilistic Seismic Hazard Analysis (PSHA). This relationship is a mathematical equation which calculates the ground motion parameters as a function of earthquake magnitude, distance, site conditions and possibly other parameters using various regression methods. Each of these GMPEs uses the preliminary data of the selected earthquakes. In this study, considering the fundamental role of these equations in PSHA, an alternative method was proposed to select suitable equation at every intensity (earthquake magnitude) and distance (site distance to fault) according to the preliminary data aggregation at every intensity and distance. The results show that the use of this method makes a significant difference to the PSHA results over selecting one equation or using logic tree. Also, a practical example of this new method was described in Iran as one of the world’s earthquake-prone areas.