Linear discriminant analysis scores for the classification of G. selincuoensis in northern Tibet using otolith shape variables including shape indices and Fourier descriptors. The first discriminant axis explained 32.6% variation and the second axis explained 18.6% variation (Wilks’ λ = 0.263).