Species–area curve is widely used to describe biological diversity across spatial scales. However, there are major ecological and statistical limitations in its application and inference, such as sampling methods, species-area scales, and statistical functions that must be considered in diversity assessment. This study aimed to explain relations among sampling methods, species-area curve scales and statistical functions in the Perk forested catchment, Zagros forests, western Iran. We studied plant diversity during the peak of the growing season, from early April to late May, 2019. The stratified random sampling was applied to collect species richness data on the same sized plots (4m2). The mean of data from noncontiguous plots was used to construct Scheiner's type IIIB curve. A comprehensive set of species-area models including five convex and three sigmoid models were tested. Akaike information criterion (AIC) and Bayesian information criterion (BIC) were used to choose the best-fitted function model. Furthermore, the Log Error of Extrapolation (LEE) was used to assess the reliability of extrapolation. Both AIC and BIC results suggested that the Lomolino and Cumulative Weibull functions as the best-fit curves for extrapolating species richness beyond sampling area. However, the LEE showed the superiority of the Rational function. This discrepancy between LEE and both AIC and BIC introduced models originates from the LEE limitation of application to tenfold of the sampling area. This study showed that the sampling scheme and the study scale are the two important factors in the shape of statistical function (convex or sigmoid) which is using in constructing species-area curves. In contrast to the prevailing use of power function, the results showed that in a broader scale like the landscape scale, the Lomolino and Cumulative Weibull functions are superior in estimating species richness in Zagros forests, western Iran.