Genotype-by-environment (G × E) interaction is a critical concept in plant breeding and agronomic research, as it reflects the differential responses of genotypes to varying environmental conditions. The presence of significant G × E interaction poses challenges for the identification and selection of superior genotypes, given that their performance may fluctuate across environments. To assess and interpret G × E interactions effectively, stability statistics play a pivotal role in quantifying genotype adaptability and performance consistency. Over the past few decades, numerous stability statistics have been developed for the analysis G × E interactions. These statistics are generally classified into two broad categories: univariate and multivariate approaches. For greater conceptual clarity, this chapter presents all statistics within the frameworks of independent and model-based approaches. Furthermore, a practical example is provided to guide researchers in implementing these methods using R software.