AMMI-Bayesian perspective in the selection of pre-cultivars of carioca beans in Agreste-Sertão of Pernambuco, Brazil.

The productivity of beans is greatly influenced by the different edaphoclimatic conditions in the Agreste-Sertão region, requiring the identification of adapted and stable genotypes to minimize the effects of the interaction between genotypes per environments (GxE). The objective of this work was to analyze the adaptability and stability of carioca bean pre-cultivars in three municipalities in the Agreste-Sertão of Pernambuco using the AMMI model in its Bayesian version BAMMI and compare the results with the frequentist approach. According to the results, the BAMMI analysis showed better predictive capacity, as well as better performance in the study of adaptability and stability. The cultivar BRS Notável stood out in terms of main effect and stability. Adaptability of genotypes to specific locations was also observed, enabling the use of the positive effect of the GxE interaction, which was more evident with the BAMMI model. From this work, the flexibility of BAMMI model to deal with data resulting from multi-environmental experiments can be seen, overcoming limitations of the standard analysis of the AMMI model.

Bayesian factor analytic model: An approach in multiple environment trials.

One of the main challenges in plant breeding programs is the efficient quantification of the genotype-by-environment interaction (GEI). The presence of significant GEI may create difficulties for breeders in the selection and recommendation of superior genotypes for a wide environmental network. Among the diverse statistical procedures developed for this purpose, we highlight those based on mixed models and factor analysis that are called factor analytic (FA) models. However, some inferential issues are related to the factor analytic model, such as Heywood cases that make the model non-identifiable. Moreover, the representation of the loads and factors in the conventional biplot does not involve any measurement of uncertainty. In this work, we propose dealing with the FA model using the Bayesian framework with direct sampling of factor loadings via spectral decomposition; this guarantees identifiability in the estimation process and eliminates the need for the rotationality of factor loadings or imposition of any ad hoc constraints. We used simulated and real data to illustrate the method's application in multi-environment trials (MET) and to compare it with traditional FA mixed models on controlled unbalancing. In general, the Bayesian FA model was robust under different simulated unbalanced levels, presenting the superior predictive ability of missing data when compared to competing models, such as those based on FA mixed models. In addition, for some scenarios, the classical FA mixed model failed in estimating the full FA model, illustrating the parametric problems of convergence in these models. Our results suggest that Bayesian factorial models might be successfully used in plant breeding for MET analysis.

A Bayesian Shrinkage Approach for AMMI Models.

Linear-bilinear models, especially the additive main effects and multiplicative interaction (AMMI) model, are widely applicable to genotype-by-environment interaction (GEI) studies in plant breeding programs. These models allow a parsimonious modeling of GE interactions, retaining a small number of principal components in the analysis. However, one aspect of the AMMI model that is still debated is the selection criteria for determining the number of multiplicative terms required to describe the GE interaction pattern. Shrinkage estimators have been proposed as selection criteria for the GE interaction components. In this study, a Bayesian approach was combined with the AMMI model with shrinkage estimators for the principal components. A total of 55 maize genotypes were evaluated in nine different environments using a complete blocks design with three replicates. The results show that the traditional Bayesian AMMI model produces low shrinkage of singular values but avoids the usual pitfalls in determining the credible intervals in the biplot. On the other hand, Bayesian shrinkage AMMI models have difficulty with the credible interval for model parameters, but produce stronger shrinkage of the principal components, converging to GE matrices that have more shrinkage than those obtained using mixed models. This characteristic allowed more parsimonious models to be chosen, and resulted in models being selected that were similar to those obtained by the Cornelius F-test (α = 0.05) in traditional AMMI models and cross validation based on leave-one-out. This characteristic allowed more parsimonious models to be chosen and more GEI pattern retained on the first two components. The resulting model chosen by posterior distribution of singular value was also similar to those produced by the cross-validation approach in traditional AMMI models. Our method enables the estimation of credible interval for AMMI biplot plus the choice of AMMI model based on direct posterior distribution retaining more GEI pattern in the first components and discarding noise without Gaussian assumption as requested in F-based tests or deal with parametric problems as observed in traditional AMMI shrinkage method.