Probability plotting position is generally used for the graphical analysis of the annual maximum quantile and the estimation of exceedance probability to display the fitness between sample and an appropriate probability distribution. In addition, it is used to apply a specific goodness of fit test. Plotting position formula to define the probability plotting position has been studied in many researches. Especially, the GEV distribution which is an important probability distribution to analyze the frequency of hydrologic data was popular. In this study, the theoretical reduced variates are derived using the mean value of order statistics to derived an appropriate plotting position formula for the GEV distribution. In addition, various forms of plotting position formula considering various sample sizes and coefficients of skewness related with shape parameters are applied. The parameters of plotting position formulas are estimated using the genetic algorithm. The accuracy of derived plotting position formula is estimated by the errors between the theoretical reduced variates and those by various plotting position formulas including the derived ones in this study. As a result, the errors by derived plotting position formula is the smallest at the range of shape parameter with -0.25～0.10.