Genotype × environment (G×E) interaction is an important source of variation in soybean yield, which can significantly influence selection in breeding programs. This study aimed to select superior soybean genotypes for performance and yield stability, from data from multi-environment trials (METs), through GGE biplot analysis that combines the main effects of the genotype (G) plus the genotype-by-environment (G×E) interaction. As well as, through path analysis, determine the direct and indirect influences of yield components on soybean grain yield, as a genotype selection strategy. Eight soybean genotypes from the breeding program of Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA) were evaluated in field trials using a randomized block experimental design, in an 8 x 8 factorial scheme with four replications in eight different environments of the Cerrado of Northeastern Brazil during two crop seasons. Phenotypic performance data were measured for the number of days to flowering (NDF), height of first pod insertion (HPI), final plant height (FPH), number of days to maturity (NDM), mass of 100 grains (M100) and grain yield (GY). The results revealed that the variance due to genotype, environment, and G×E interaction was highly significant (P < 0.001) for all traits. The ST820RR, BRS 333RR, BRS SambaíbaRR, M9144RR and M9056RR genotypes exhibited the greatest GY stability in the environments studied. However, only the BRS 333RR genotype, followed by the M9144RR, was able to combine good productive performance with high yield stability. The study also revealed that the HPI and the NDM are traits that should be prioritized in the selection of soybean genotypes due to the direct and indirect effects on the GY.