In this study, the decay of earthquake aftershock sequences of some major earthquakes in different tectonic regimes in the Iranian plateau is discussed. The studied earthquakes are Rigan [2010], Ahar-Varzaghan [2012], Goharan [2013], Sefidsang [2017] and Ezgeleh [2017]. The spatial and temporal windows are considered based on the method proposed by Gardner and Knopoff [1974] to compute decay parameters for each sequence. The decay rates of sequences were compared to well-known models to find the best fit for each sequence. The results showed that the modified Omori is the best fit for Ahar-Varzaghan and Ezgeleh sequences, for Rigan and Sefidsang sequences the modified Omori and the Kisslinger ones found as the best fits. The values of the p parameter of the Reasenberg and Shcherbakov models were larger compared to the Omori model, but the parameter of the Kisslinger model was slightly smaller compared to the Omori one. The c parameter showed an inverse relation to the threshold magnitude. The correlation between the p and c parameters and also the and the Gutenberg and Richter (G-R) parameters were investigated. In addition, we made use of a graphical method to analyze the seismic sequence of the Ezgeleh earthquake during 13 months after the main event. The graphical method was successful to estimate the occurrence of an event with an approximate magnitude of M=6.4 in the sequence.