The doubt in the calculation algorithm of the standardized precipitation index (SPI), which is widely preferred in the evaluation and monitoring of drought, still remains up-to-date because its calculation process is performed in the form of standardization or normalization with a default probability distribution. Therefore, the success of this index is directly affected by the choice of the probability distribution model. This study is based on the effect of three different parameter estimation methods on the calculation process, as well as the comparison of the SPI results calculated based on the default Gamma distribution and the distribution with the best ability to represent the 3-and 12-month consecutive summed rainfall data among the 15 candidate distributions namely Gamma (GAM), Generalized Extreme Value (GEV), Pearson Type III (P III), Log Pearson Type III (LP III), two-parameter Lognormal (LN2), three-parameter Lognormal (LN3), Generalized Logistic (GLOG), Extreme Value Type I (EVI), Generalized Pareto (GPAR), Weilbul (W), Normal (N), Exponential (EXP), Logistic (LOG), four-parameter Wakeby (WK4), and five-parameter Wakeby (WK5) distributions. Approximately 68.4% and 18.4% of the 3-month data considered had the best fit to the Weibull and Pearson III distribution, while approximately 24% and 18% of the 12-month data had the best fit to the Weibull and Logistic distribution. On the other hand, it was found that the default Gamma distribution calculated the extreme drought categories significantly more than the best-fit distribution model. In terms of parameter estimation methods, L-moments for 3-month series and maximum likelihood approaches for 12-month series were most dominant.