A soja apresenta grande importância socioeconômica no cenário nacional. Para a Região Meio-Norte brasileira, apesar das suas condições edafoclimáticas favoráveis para a produção de grãos, ainda são escassos, na literatura, trabalhos objetivando a seleção de genótipos de soja adaptados e com comportamento produtivo estável. O objetivo deste trabalho foi selecionar genótipos superiores quanto a adaptabilidade, estabilidade e produtividade em ambientes dos estados do Maranhão e Piauí, a partir de duas diferentes metodologias. Foram avaliados 16 genótipos de soja, em oito ambientes, formados pela combinação local x safra agrícola. Em todos os ensaios, empregou-se o delineamento em blocos casualizados, com três repetições. Avaliou-se os seguintes caracteres: número de dias para maturação, altura da planta, peso de cem sementes e produtividade de grãos. Foram realizadas análises de variância individuais, conjunta e análise de adaptabilidade e estabilidade a partir da produtividade de grãos, utilizando duas metodologias: GGE Biplot e REML/BLUP. Pela análise conjunta, a interação genótipos x ambientes foi altamente significativa para todos os caracteres avaliados. Os genótipos com maior adaptabilidade e estabilidade produtiva pela metodologia REML/BLUP são: BRS8980 IPRO, BRASBT13-0528, M8372 IPRO e BRASBT13-0621. Observou-se concordância com os resultados obtidos pela metodologia GGE Biplot, a qual considerou os genótipos M8372 IPRO, BRASBT13-0528 e BRASBT13-0621 como os mais próximos do ideal. Pela facilidade de identificação e seleção de genótipos superiores quanto a adaptabilidade, estabilidade e produtividade, a metodologia REML/BLUP mostrou-se mais vantajosa que a GGE Biplot.