Resumo:
O objetivo deste trabalho foi selecionar, sob a perspectiva bayesiana, genótipos de feijão-caupi (Vigna unguiculata) que reúnam alta adaptabilidade e estabilidade fenotípicas, no Estado do Mato Grosso do Sul. Foram utilizados dados de quatro experimentos, conduzidos em delineamento de blocos ao acaso, em que a produtividade de grãos de 20 genótipos de feijão-caupi semiprostrado foi avaliada. Para representar as distribuições a priori pouco informativas, utilizaram-se distribuições de probabilidade com grande variância; e, para representar distribuições a priori informativas, adotou-se o conceito de metanálise, com uso de informações de trabalhos anteriores. A comparação entre as distribuições a priori foi realizada por meio do fator de Bayes. A abordagem bayesiana proporciona maior acurácia na seleção de genótipos de feijão-caupi semiprostrado, com elevadas adaptabilidade e estabilidade fenotípicas avaliadas por meio da metodologia de Eberhart & Russell. Com base nas prioris informativas, os genótipos MNC99-507G-4, TE97-309G-24, MNC99-542F-7 e BR 17-Gurguéia são classificados como de alta adaptabilidade a ambientes favoráveis. Já os genótipos TE96-290-12G, MNC99-510F-16, MNC99-508G-1, MNC99-541F-21, MNC99-542F-5 e MNC99-547F-2 apresentam alta adaptabilidade a ambientes desfavoráveis.
Termos para indexação:
Vigna unguiculata; fator de Bayes; interação genótipo x ambiente; metanálise; priori informativa.
Abstract:
The objective of this work was to select, under the Bayesian perspective, cowpea (Vigna unguiculata) genotypes that meet high phenotypic adaptability and stability, in the state of Mato Grosso do Sul, Brazil. Data from four experiments, conducted in a randomized complete block design, were used, in which grain yield of 20 semiprostrate cowpea genotypes was evaluated. To represent non-informative prior distributions, probability distributions with high variance were used; and, to represent informative prior distributions, a metanalysis concept was adopted using information from previous studies. The comparison between the prior distributions was done using the Bayes factor. The Bayesian approach provides greater accuracy in the selection of semiprostrate cowpea genotypes, with high phenotypic adaptability and stability assessed by the Eberhart & Russell methodology. Based on the informative priors, the MNC99-507G-4, TE97-309G-24, MNC99-542F-7, and BR 17-Gurguéia genotypes are classified as with high adaptability to favorable environments. The TE96-290-12G, MNC99-510F-16, MNC99-508G-1, MNC99-541F-21, MNC99-542F-5, and MNC99-547F-2 genotypes have high adaptability to unfavorable environments.
Index terms:
Vigna unguiculata; Bayes fator; genotype x environment interaction; metanalysis; informative prior.
Introdução
O feijão-caupi [Vigna unguiculata (L.) Walp.] é uma das fontes alimentares mais importantes para regiões tropicais e subtropicais do planeta. Atualmente, o Brasil é o terceiro maior produtor mundial da cultura, que é cultivada, sobretudo, nas regiões Norte (55,8 mil hectares) e Nordeste (1,2 milhão de hectare) do País, onde constitui o prato básico das classes de menor poder aquisitivo (Santos et al., 2014SANTOS, J.A. da S.; SOARES, C.M.G.; CORRÊA, A.M.; TEODORO, P.E.; RIBEIRO, L.P.; ABREU, H.K.A. de. Agronomic performance and genetic dissimilarity among cowpea (Vigna unguiculata (L.) Walp.) genotypes. Global Advanced Research Journal of Agricultural Science, v.3, p.271-277, 2014.).
O feijão-caupi tem sido cultivado em três regiões distintas do Brasil: Centro-oeste, Norte e Nordeste (Torres et al., 2015TORRES, F.E.; SAGRILO, E.; TEODORO, P.E. ; RIBEIRO, L.P.; CARGNELUTTI FILHO, A. Número de repetições para avaliação de caracteres em genótipos de feijão-caupi. Bragantia, v.74, p.161-168, 2015. DOI: 10.1590/1678-4499.0393.
https://doi.org/10.1590/1678-4499.0393....
). Isso acarreta dificuldades na recomendação de cultivares, uma vez que as adaptadas a uma determinada condição de cultivo podem não apresentar bom desempenho em outra condição ambiental. Essa resposta diferencial dos genótipos aos diferentes ambientes de cultivo é comumente avaliada por meio da interação genótipo x ambiente (GxE).
Nesse sentido, alguns trabalhos têm sido conduzidos com o objetivo de selecionar genótipos de feijão-caupi superiores quanto à adaptabilidade e à estabilidade produtiva. As metodologias mais comumente utilizadas para esse fim são as de modelos mistos (REML/BLUP); a de Eberhart & Russel (1966)EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
https://doi.org/10.2135/cropsci1966.0011...
, baseada na regressão linear simples; e as de Yates & Cochran (1938)YATES, F.; COCHRAN, W.G. The analysis of group of experiments. The Journal of Agricultural Science, v.28, p.556-580, 1938. DOI: 10.1017/S0021859600050978.
https://doi.org/10.1017/S002185960005097...
(tradicional) e de Wricke (1965)WRICKE, G. Zur berechnung der okovalenz bei sommerweizen und hafer. Zeitschrift für Pflanzenzüchtung, v.52, p.127-138, 1965., ambas baseadas na análise de variância (Freire Filho et al., 2002FREIRE FILHO, F.R.; RIBEIRO, V.Q.; ROCHA, M. de M.; LOPES, Â.C. de A. Adaptabilidade e estabilidade da produtividade de grãos de linhagens de caupi de porte enramador. Revista Ceres, v.49, p.383-393, 2002.; Rocha et al., 2007ROCHA, M. de M.; FREIRE FILHO, F.R.; RIBEIRO, V.Q.; CARVALHO, H.W.L. de; BELARMINO FILHO, J.; RAPOSO, J.A.A.; ALCÂNTARA, J. dos P.; RAMOS, S.R.R.; MACHADO, C. de F. Adaptabilidade e estabilidade produtiva de genótipos de feijão-caupi de porte semi-ereto na Região Nordeste do Brasil. Pesquisa Agropecuária Brasileira, v.42, p.1283-1289, 2007. DOI: 10.1590/S0100-204X2007000900010.
https://doi.org/10.1590/S0100-204X200700...
; Barros et al., 2013BARROS, M.A.; ROCHA, M. de M.; GOMES, R.L.F.; SILVA, K.J.D. e; NEVES, A.C. das. Adaptabilidade e estabilidade produtiva de feijão-caupi de porte semiprostrado. Pesquisa Agropecuária Brasileira, v.48, p.403-410, 2013. DOI: 10.1590/S0100-204X2013000400008.
https://doi.org/10.1590/S0100-204X201300...
; Torres et al., 2015TORRES, F.E.; SAGRILO, E.; TEODORO, P.E. ; RIBEIRO, L.P.; CARGNELUTTI FILHO, A. Número de repetições para avaliação de caracteres em genótipos de feijão-caupi. Bragantia, v.74, p.161-168, 2015. DOI: 10.1590/1678-4499.0393.
https://doi.org/10.1590/1678-4499.0393....
).
Apesar da grande quantidade de metodologias disponíveis para estudo da interação GxE, poucas fazem uso de informações a priori na estimação dos parâmetros. Nascimento et al. (2011)NASCIMENTO, M.; SILVA, F.F.; SÁFADI, T.; NASCIMENTO, A.C.M.C.; FERREIRA, R. de P.; CRUZ, C.D. Abordagem bayesiana para avaliação da adaptabilidade e estabilidade de genótipos de alfafa. Pesquisa Agropecuária Brasileira, v.46, p.26-32, 2011. DOI: 10.1590/S0100-204X2011000100004.
https://doi.org/10.1590/S0100-204X201100...
formularam uma abordagem bayesiana para o método de Eberhart & Russel (1966)EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
https://doi.org/10.2135/cropsci1966.0011...
, que permitiu maior acurácia na seleção de genótipos de alface, em comparação ao método tradicional. Esta abordagem utiliza três conceitos básicos: uma informação inicial (probabilidade a priori), geralmente assumida como uma lei de probabilidade conjunta sobre os parâmetros, antes de se obter a informação particular das amostras y1,...,yn da variável aleatória; o modelo probabilístico da variável aleatória resposta y, com o qual se obtém a verossimilhança da amostra; e o teorema de Bayes, que combina a informação a priori e a função de verossimilhança, e gera uma distribuição a posteriori (Cotes et al., 2006COTES, J.M.; CROSSA, J.; SANCHES, A.; CORNELIUS, P.A. A Bayesian approach for assessing the stability of genotypes. Crop Science, v.46, p.2654-2665, 2006. DOI: 10.2135/cropsci2006.04.0227.
https://doi.org/10.2135/cropsci2006.04.0...
). Assim, no contexto bayesiano, os parâmetros passam a ser interpretados como variáveis aleatórias, com uma lei de probabilidade (distribuição a priori) que reflete uma informação inicial (ou a falta de informação) sobre eles, independentemente do que os dados possam mostrar (Ramos Molina et al., 2012RAMOS MOLINA, L.M.; SANCHE, A.; COTES TORRES, J.M. Inferência bayesiana na análise de testes regionais de arroz em dois sistemas de cultivo. Revista Facultad Nacional de Agronomía, v.64, p.5883-5891, 2012.). Outra peculiaridade dessa análise é que a homogeneidade da variância e a adequação dos resíduos à distribuição normal, pressupostos para muitas metodologias, não são requeridos na abordagem bayesiana (Cotes et al., 2006COTES, J.M.; CROSSA, J.; SANCHES, A.; CORNELIUS, P.A. A Bayesian approach for assessing the stability of genotypes. Crop Science, v.46, p.2654-2665, 2006. DOI: 10.2135/cropsci2006.04.0227.
https://doi.org/10.2135/cropsci2006.04.0...
). Embora a inferência bayesiana seja um procedimento estatístico robusto e tenha muitas aplicações possíveis, sua utilização no melhoramento genético de plantas ainda é pequena.
O objetivo deste trabalho foi selecionar, sob a perspectiva bayesiana, genótipos de feijão-caupi (Vigna unguiculata) que reúnam alta adaptabilidade e estabilidade fenotípicas, no Estado do Mato Grosso do Sul.
Material e Métodos
Foram conduzidos quatro ensaios de valor de cultivo e uso com 18 linhagens e duas cultivares de feijão-caupi, em 2005 e 2006, nos municípios de Aquidauana, Chapadão do Sul e Dourados (Tabela 1). Utilizou-se o delineamento experimental de blocos ao acaso, com 20 tratamentos e quatro repetições. A unidade experimental foi constituída por quatro fileiras de 5,0 m de comprimento, espaçadas em 0,5 m entre fileiras e 0,25 m entre plantas dentro da fileira. Em cada unidade experimental, a produtividade de grãos foi avaliada nas duas fileiras centrais, tendo-se corrigido a umidade dos grãos para 13%, e extrapolada para kg ha-1.
As 18 linhagens de feijão-caupi avaliadas foram: MNC99-505G-11, MNC99-507G-4, MNC99-507G-8, MNC99-510G-8, MNC99-510F-16, TE96-290-12G, TE97-309G-18, TE97-304G-4, TE97-304G-12, TE97-309G-24, MNC99-508G-1, MNC99-541F-15, MNC99-541F-18, MNC99-541F-21, MNC99-542F-5, MNC99-542F-7, MNC99-547F-2 e CNC x 409-11F-P2. Além das linhagens, foram avaliadas as cultivares BRS Paraguaçu e BR 17-Gurguéia. Tanto as linhagens quanto as cultivares são oriundas do programa de melhoramento genético de feijão-caupi da Embrapa Meio-Norte.
Para verificar o efeito da interação GxE, o fator de Bayes foi calculado para os dados de produtividade de grãos ajustados a diferentes modelos. Essa abordagem foi utilizada com o intuito de manter toda a apresentação no contexto bayesiano, visto que o fator de Bayes é consistente com a análise de variância clássica (Rouder et al., 2012ROUDER, J.N.; MOREY, R.D.; SPECKMAN, P.L.; PROVINCE, J.M. Default Bayes factors for ANOVA designs. Journal of Mathematical Psychology, v.56, p.356-374, 2012. DOI: 10.1016/j.jmp.2012.08.001.
https://doi.org/10.1016/j.jmp.2012.08.00...
). Após passarem pela abordagem bayesiana, os dados foram submetidos às análises de adaptabilidade e estabilidade, por meio da metodologia de Eberhart & Russell (1966)EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
https://doi.org/10.2135/cropsci1966.0011...
.
O modelo de regressão linear de Eberhart & Russell (1966)EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
https://doi.org/10.2135/cropsci1966.0011...
adotado foi Yij = β0i + βliIj + Ψij, em que Yij é a média observada do genótipo i no ambiente j; β0i é o coeficiente linear referente ao i-ésimo genótipo; βli é o coeficiente de regressão do genótipo i; Ij é o índice ambiental j; e Ψij são os erros aleatórios, compostos pelo desvio da regressão do genótipo i no ambiente j e pelo erro médio associado à média. O índice ambiental foi estimado de acordo com a equação com
em que é a média geral; é a média do ambiente j; e n é o número de ambientes.
De acordo com o método de Eberhart & Russell (1966)EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
https://doi.org/10.2135/cropsci1966.0011...
, a adaptabilidade do genótipo foi medida pelo parâmetro β1i, enquanto a estabilidade de comportamento foi avaliada por meio da variância dos desvios da regressão (s2di) e pelo coeficiente de determinação (R2), que, segundo Cruz & Regazzi (2007)CRUZ, C.D.; REGAZZI, A.J. Modelos biométricos aplicados ao melhoramento genético. Viçosa: Imprensa Universitária, 2007. 480p., é uma medida auxiliar de avaliação da estabilidade. Quando a s2di é significativa e o R2 é superior a 80%, considera-se que o R2indica previsibilidade aceitável. Nessa abordagem frequentista, as hipóteses de interesse, ou seja, são avaliadas pelo testes t e F, respectivamente.
Para a análise bayesiana, foram considerados apenas os genótipos avaliados em ao menos um dos seguintes trabalhos: Freire Filho et al. (2002)FREIRE FILHO, F.R.; RIBEIRO, V.Q.; ROCHA, M. de M.; LOPES, Â.C. de A. Adaptabilidade e estabilidade da produtividade de grãos de linhagens de caupi de porte enramador. Revista Ceres, v.49, p.383-393, 2002., Carvalho et al. (2006)CARVALHO, H.W.L. de; BRITO NETO, J.; FREIRE FILHO, F.R.; ROCHA, M. de M.; RIBEIRO, V.Q.; RODRIGUES, A.R. dos S.; RIBEIRO, S.S.; NOGUEIRA, L.C. Adaptabilidade e estabilidade de linhagens de feijão-caupi de porte prostrado em Sergipe e Alagoas. In: CONGRESSO NACIONAL DE FEIJÃO-CAUPI, 2006, Teresina. Anais eletrônicos... Teresina: Embrapa Meio Norte, 2006. Disponível em: <Disponível em: http://www.cpamn.embrapa.br/anaisconac2006/listaresumo.htm >. Acesso em: 3 fev. 2015.
http://www.cpamn.embrapa.br/anaisconac20...
e Valadares et al. (2010)VALADARES, R. de N.; MOURA, M. da C.C.L.; SILVA, A.F.A. da; SILVA, L.S. da; VASCONCELOS, M. da C.C. de A.; SILVA, R.G. Adaptabilidade e estabilidade fenotípica em genótipos de feijão-caupi (Vigna unguiculata (L.) Walp.) de porte ereto/semi-ereto nas mesorregiões leste e sul maranhense. Agropecuária Científica no Semi-árido, v.6, p.21-27, 2010., utilizados como referência para a especificação das distribuições a priori sob um contexto de metanálise (Tabela 2).
Estimativas da média (β0i) e dos parâmetros de adaptabilidade (β1i) e estabilidade (σ2di) obtidas com a metodologia de Eberhart & Russel (1966)EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
https://doi.org/10.2135/cropsci1966.0011... , nos trabalhos de Freire Filho et al. (2002)FREIRE FILHO, F.R.; RIBEIRO, V.Q.; ROCHA, M. de M.; LOPES, Â.C. de A. Adaptabilidade e estabilidade da produtividade de grãos de linhagens de caupi de porte enramador. Revista Ceres, v.49, p.383-393, 2002., Carvalho et al. (2006)CARVALHO, H.W.L. de; BRITO NETO, J.; FREIRE FILHO, F.R.; ROCHA, M. de M.; RIBEIRO, V.Q.; RODRIGUES, A.R. dos S.; RIBEIRO, S.S.; NOGUEIRA, L.C. Adaptabilidade e estabilidade de linhagens de feijão-caupi de porte prostrado em Sergipe e Alagoas. In: CONGRESSO NACIONAL DE FEIJÃO-CAUPI, 2006, Teresina. Anais eletrônicos... Teresina: Embrapa Meio Norte, 2006. Disponível em: <Disponível em: http://www.cpamn.embrapa.br/anaisconac2006/listaresumo.htm >. Acesso em: 3 fev. 2015.
http://www.cpamn.embrapa.br/anaisconac20... e Valadares et al. (2010)VALADARES, R. de N.; MOURA, M. da C.C.L.; SILVA, A.F.A. da; SILVA, L.S. da; VASCONCELOS, M. da C.C. de A.; SILVA, R.G. Adaptabilidade e estabilidade fenotípica em genótipos de feijão-caupi (Vigna unguiculata (L.) Walp.) de porte ereto/semi-ereto nas mesorregiões leste e sul maranhense. Agropecuária Científica no Semi-árido, v.6, p.21-27, 2010..
Na abordagem bayesiana, ao se considerar o modelo estatístico Yij= β0i + βliIj + Ψij e se assumir que cada observação Yij tem distribuição a função de verossimilhança para cada genótipo i é dada por:
Para a estimação dos parâmetros de adaptabilidade e estabilidade, é necessário atribuir distribuições a priori para os parâmetros. Para foram consideradas as seguintes distribuições: esta última, uma gama inversa com média e variância iguais a respectivamente.
Ao se assumir a independência entre os parâmetros dessas distribuições, a distribuição a priori conjunta para cada genótipo é dada por:
Para fazer inferências sobre os parâmetros de interesse, é necessário obter suas distribuições marginais a posteriori. Ao denotar o vetor de parâmetros para cada genótipo i por em que p=1, 2, 3, a distribuição marginal a posteriori para o parâmetro θpi é obtida pela seguinte integral: isto é, a integral em relação a todos os parâmetros do vetor, exceto o p-ésimo componente.
Na maioria dos casos, essas integrais são complexas se não apresentam soluções exatas. Para contornar esse problema, utiliza-se outra metodologia: uma amostra da distribuição conjunta a posteriori é obtida ao se utilizar as cadeias de Markov e o método de Monte Carlo para a determinação dos momentos de interesse das distribuições marginais (Cassella & George, 1992CASSELLA, G.; GEORGE, E.I. Explaining the Gibbs sampler. The American Statistician, v.46, p.167-174, 1992.). No presente trabalho, a metodologia foi implementada no programa R (The R Foundation, 2010THE R FOUNDATION. R: the R project for statistical computing. Vienna, 2010. Available at: <Available at: https://www.r-project.org >. Accessed on: 15 Apr. 2015.
https://www.r-project.org...
), e a amostra da distribuição conjunta foi obtida por meio da função MCMCregress do pacote MCMCpack. Para avaliar a influência de informações a priori na estimação dos parâmetros de adaptabilidade e estabilidade, optou-se pela utilização de dois modelos diferentes (Nascimento et al., 2011NASCIMENTO, M.; SILVA, F.F.; SÁFADI, T.; NASCIMENTO, A.C.M.C.; FERREIRA, R. de P.; CRUZ, C.D. Abordagem bayesiana para avaliação da adaptabilidade e estabilidade de genótipos de alfafa. Pesquisa Agropecuária Brasileira, v.46, p.26-32, 2011. DOI: 10.1590/S0100-204X2011000100004.
https://doi.org/10.1590/S0100-204X201100...
).
No primeiro modelo (M1), foram consideradas distribuições a priori informativas, cujas informações eram provenientes da aplicação da técnica de metanálise, que se caracteriza pela utilização de informações dos trabalhos anteriores de Freire Filho et al. (2002)FREIRE FILHO, F.R.; RIBEIRO, V.Q.; ROCHA, M. de M.; LOPES, Â.C. de A. Adaptabilidade e estabilidade da produtividade de grãos de linhagens de caupi de porte enramador. Revista Ceres, v.49, p.383-393, 2002., Carvalho et al. (2006)CARVALHO, H.W.L. de; BRITO NETO, J.; FREIRE FILHO, F.R.; ROCHA, M. de M.; RIBEIRO, V.Q.; RODRIGUES, A.R. dos S.; RIBEIRO, S.S.; NOGUEIRA, L.C. Adaptabilidade e estabilidade de linhagens de feijão-caupi de porte prostrado em Sergipe e Alagoas. In: CONGRESSO NACIONAL DE FEIJÃO-CAUPI, 2006, Teresina. Anais eletrônicos... Teresina: Embrapa Meio Norte, 2006. Disponível em: <Disponível em: http://www.cpamn.embrapa.br/anaisconac2006/listaresumo.htm >. Acesso em: 3 fev. 2015.
http://www.cpamn.embrapa.br/anaisconac20...
e Valadares et al. (2010)VALADARES, R. de N.; MOURA, M. da C.C.L.; SILVA, A.F.A. da; SILVA, L.S. da; VASCONCELOS, M. da C.C. de A.; SILVA, R.G. Adaptabilidade e estabilidade fenotípica em genótipos de feijão-caupi (Vigna unguiculata (L.) Walp.) de porte ereto/semi-ereto nas mesorregiões leste e sul maranhense. Agropecuária Científica no Semi-árido, v.6, p.21-27, 2010.. Dessa forma, todos os genótipos avaliados nos experimentos apresentados foram considerados para a análise bayesiana (Tabela 2). As informações foram inseridas na análise por meio dos valores assumidos para os parâmetros das distribuições a priori, denominados de hiperparâmetros.
Esses valores foram dados pela média e pela variância da amostra composta pelas estimativas dos parâmetros obtidas das referências citadas (Tabela 2), o que resultou nas distribuições: em que são as médias das estimativas dos β0i; são as médias das estimativas dos é a variância dos valores médios de é variância dos valores médios de e αi e βi são os valores obtidos pela resolução do sistema:
No segundo modelo (M2), foram utilizadas distribuições a priori pouco informativas, que representam distribuições de probabilidade com grande variância. Foram adotadas as seguintes distribuições: ~GamaInv(αi=0,00001; βi=5.000). A Figura 1 apresenta as distribuições pouco informativas obtidas por meio dos hiperparâmetros fornecidos, de acordo com Nascimento et al. (2011)NASCIMENTO, M.; SILVA, F.F.; SÁFADI, T.; NASCIMENTO, A.C.M.C.; FERREIRA, R. de P.; CRUZ, C.D. Abordagem bayesiana para avaliação da adaptabilidade e estabilidade de genótipos de alfafa. Pesquisa Agropecuária Brasileira, v.46, p.26-32, 2011. DOI: 10.1590/S0100-204X2011000100004.
https://doi.org/10.1590/S0100-204X201100...
.
Distribuições pouco informativas utilizadas como informação a priori: A, normal (μ=0 e =100.000); e B, gama inversa (a = 0,00001 e b = 5.000).
A comparação entre M1 e M2, ou seja, entre as distribuições a priori informativas e não informativas, foi feita pelo fator de Bayes (Kass & Raftery, 1995KASS, R.E.; RAFTERY, A.E. Bayes factors. Journal of the American Statistical Association, v.90, p.773-795, 1995. DOI: 10.1080/01621459.1995.10476572.
https://doi.org/10.1080/01621459.1995.10...
). O cálculo do fator de Bayes foi realizado por meio da função Bayes Factor do pacote MCMCpack. Conforme Jeffreys (1961)JEFFREYS, H. Theory of probability. 3rd ed. Oxford: Claredon, 1961. 447p., o fator de Bayes pode ser interpretado da seguinte forma: FBij<1 mostra evidência a favor do modelo j; 1≤FBij<3 mostra evidência moderada a favor do modelo i; 3≤FBij<10 mostra evidência substancial a favor do modelo i; 10≤FBij<30 mostra evidência forte a favor do modelo i; 30≤FBij<100 mostra evidência muito forte a favor do modelo i; e FBij≥100 mostra evidência decisiva a favor do modelo i.
Com relação ao parâmetro de estabilidade (), as amostras de sua distribuição marginal foram obtidas indiretamente, uma vez que esse parâmetro representa uma função de . Ao se obter valores para indiretamente, em cada iteração, obtêm-se valores para , por meio da seguinte expressão: em que QMR é o quadrado médio do resíduo fornecido pela análise da variância; e r é número de repetições do experimento.
As hipóteses de interesse foram testadas pela construção de intervalos de credibilidade para os parâmetros, obtidos diretamente da distribuição marginal a posteriori dos parâmetros. Assim, seja o vetor de parâmetros, o intervalo de credibilidade para θi, com probabilidade de cobertura δ, foi dado por:
em que θ* e θ* representam, respectivamente, os limites inferiores e superiores do intervalo de credibilidade.
Como o Gibbs sampler é um algoritmo iterativo, foi necessário verificar sua convergência, que, no presente trabalho, foi constatada pela aplicação dos critérios de Heidelberger & Welch (1983)HEIDELBERGER, P.; WELCH, P.D. Simulation run length control in the presence of an initial transient. Operations Research, v.31, p.1109-1144, 1983. DOI: 10.1287/opre.31.6.1109.
https://doi.org/10.1287/opre.31.6.1109....
, Geweke (1992)GEWEKE, J. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments. In: BERNARDO, J.M.; BERGER, J.O.; DAVID, A.P.; SMITH, A.F.M. (Ed.). Bayesian statistics. New York: Oxford University, 1992. p.625-631. e Raftery & Lewis (1992)RAFTERY, A.E.; LEWIS, S.M. Comment: one long run with diagnostics: implementation strategies for Markov chain Monte Carlo. Statistical Science, v.7, p.493-497, 1992. DOI: 10.1214/ss/1177011143.
https://doi.org/10.1214/ss/1177011143....
, implementados no pacote Bayesian Output Analysis (BOA) do programa R (The R Foundation, 2010THE R FOUNDATION. R: the R project for statistical computing. Vienna, 2010. Available at: <Available at: https://www.r-project.org >. Accessed on: 15 Apr. 2015.
https://www.r-project.org...
).
Na análise bayesiana de adaptabilidade e estabilidade, foram consideradas 110.000 iterações no algoritmo Gibbs sampler para cada parâmetro do modelo de regressão adotado, com um período de aquecimento ("burn-in") de 10.000 iterações. Para obter-se uma amostra não correlacionada, considerou-se um espaçamento entre os pontos amostrados de duas iterações ("thinning"), o que resultou em uma amostra final de tamanho 50.000. As amostras representam amostras das distribuições marginais a posteriori de cada parâmetro, sob as quais foi feita a inferência de cada parâmetro.
Resultados e Discussão
O fator de Bayes obtido da comparação entre o modelo completo [M1 (numerador): bloco/ambiente + genótipos + ambientes + GxA] e o modelo restrito [M2 (denominador): bloco/ambiente + genótipos + ambientes] foi de 8.630. Esse resultado é indicativo de ocorrência de resposta diferenciada dos genótipos de acordo com o ambiente, e está de acordo com as diferenças observadas entre os ambientes em termos de altitude, latitude, longitude, precipitação e temperatura (Tabela 1). Efeitos significativos de genótipos, ambientes e da interação GxE foram relatados em diversos trabalhos com genótipos de feijão-caupi de porte semiprostrado avaliados em ensaios multiambientes, na região Meio-Norte do Brasil (Freire Filho et al., 2002FREIRE FILHO, F.R.; RIBEIRO, V.Q.; ROCHA, M. de M.; LOPES, Â.C. de A. Adaptabilidade e estabilidade da produtividade de grãos de linhagens de caupi de porte enramador. Revista Ceres, v.49, p.383-393, 2002.; Rocha et al., 2007ROCHA, M. de M.; FREIRE FILHO, F.R.; RIBEIRO, V.Q.; CARVALHO, H.W.L. de; BELARMINO FILHO, J.; RAPOSO, J.A.A.; ALCÂNTARA, J. dos P.; RAMOS, S.R.R.; MACHADO, C. de F. Adaptabilidade e estabilidade produtiva de genótipos de feijão-caupi de porte semi-ereto na Região Nordeste do Brasil. Pesquisa Agropecuária Brasileira, v.42, p.1283-1289, 2007. DOI: 10.1590/S0100-204X2007000900010.
https://doi.org/10.1590/S0100-204X200700...
; Barros et al., 2013BARROS, M.A.; ROCHA, M. de M.; GOMES, R.L.F.; SILVA, K.J.D. e; NEVES, A.C. das. Adaptabilidade e estabilidade produtiva de feijão-caupi de porte semiprostrado. Pesquisa Agropecuária Brasileira, v.48, p.403-410, 2013. DOI: 10.1590/S0100-204X2013000400008.
https://doi.org/10.1590/S0100-204X201300...
; Torres et al., 2015TORRES, F.E.; SAGRILO, E.; TEODORO, P.E. ; RIBEIRO, L.P.; CARGNELUTTI FILHO, A. Número de repetições para avaliação de caracteres em genótipos de feijão-caupi. Bragantia, v.74, p.161-168, 2015. DOI: 10.1590/1678-4499.0393.
https://doi.org/10.1590/1678-4499.0393....
).
Em relação à convergência entre as cadeias simuladas, o fator de dependência de Raftery & Lewis (1992)RAFTERY, A.E.; LEWIS, S.M. Comment: one long run with diagnostics: implementation strategies for Markov chain Monte Carlo. Statistical Science, v.7, p.493-497, 1992. DOI: 10.1214/ss/1177011143.
https://doi.org/10.1214/ss/1177011143....
forneceu valores menores que cinco, e o valor p do critério de Geweke (1992)GEWEKE, J. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments. In: BERNARDO, J.M.; BERGER, J.O.; DAVID, A.P.; SMITH, A.F.M. (Ed.). Bayesian statistics. New York: Oxford University, 1992. p.625-631.) foi sempre maior que o nível de significância pré-fixado (5%). Portanto, ambos os critérios indicaram a convergência das cadeias geradas pelo algoritmo Gibbs sampler, confirmada pelo critério de Heidelberger & Welch (1983)HEIDELBERGER, P.; WELCH, P.D. Simulation run length control in the presence of an initial transient. Operations Research, v.31, p.1109-1144, 1983. DOI: 10.1287/opre.31.6.1109.
https://doi.org/10.1287/opre.31.6.1109....
, que verifica se os valores das cadeias são provenientes de distribuição estacionária.
As estimativas dos parâmetros de adaptabilidade e estabilidade foram obtidas pelo cálculo da média a posteriori (Tabela 3). No M1 (prioris informativas), dos 20 genótipos avaliados com a metodologia bayesiana, apenas os genótipos MNC99-507G-4, TE97-309G-24, MNC99-542F-7 e BR 17-Gurguéia foram classificados como de adaptabilidade específica a ambientes favoráveis (LI>1). Os genótipos TE96-290-12G, MNC99-510F-16, MNC99-508G-1, MNC99-541F-21, MNC99-542F-5 e MNC99-547F-2 apresentaram adaptabilidade específica a ambientes desfavoráveis (LS<1). Os demais genótipos foram classificados como de adaptabilidade e estabilidade geral, com valores dos parâmetros de adaptabilidade e estabilidade dentro dos índices de confiabilidade de 95%. Entretanto, com a análise dos genótipos por meio do M2 (prioris não informativas), observou-se que nenhum genótipo apresentou adaptabilidade específica a ambientes favoráveis (LI>1). Desse modo, é possível inferir que a utilização do modelo frequentista, em que informações a priori não são levadas em conta, tende a atribuir aos genótipos adaptabilidade e estabilidade a ambientes diversos, o que compromete a recomendação acurada para o cultivo.
Na comparação entre o M1 e o M2, foram observadas diferenças de magnitudes semelhantes entre os parâmetros. Portanto, é necessário que se informe qual dos dois modelos apresenta maior qualidade de ajuste, a partir do cálculo do fator de Bayes (Nascimento et al., 2011NASCIMENTO, M.; SILVA, F.F.; SÁFADI, T.; NASCIMENTO, A.C.M.C.; FERREIRA, R. de P.; CRUZ, C.D. Abordagem bayesiana para avaliação da adaptabilidade e estabilidade de genótipos de alfafa. Pesquisa Agropecuária Brasileira, v.46, p.26-32, 2011. DOI: 10.1590/S0100-204X2011000100004.
https://doi.org/10.1590/S0100-204X201100...
). Com exceção do genótipo MNC99-541F-21, todos os demais apresentaram FB>10 (Tabela 4), o que indica que o uso de prioris informativas fornece resultados mais acurados. Apesar de a informação a priori para cada genótipo ter sido baseada em apenas um trabalho, resultados similares a esses foram relatados em estudos semelhantes com avaliação da produtividade de matéria seca de genótipos de alfafa (Medicago sativa L.) (Nascimento et al., 2011) e da produtividade de grãos de milho-pipoca (Zea mays L.) (Couto et al., 2015COUTO, M.F.; NASCIMENTO, M.; AMARAL, A.T. do; SILVA, F.F. e; VIANA, A.P.; VIVAS, M. Eberhart and Russel's Bayesian method in the selection of popcorn cultivars. Crop Science, v.55, p.571-577, 2015. DOI: 10.2135/cropsci2014.07.0498.
https://doi.org/10.2135/cropsci2014.07.0...
). Assim, espera-se que os resultados sejam ainda mais precisos, em trabalhos com maior disponibilidade de informação a priori, obtida pela metanálise.
Valores do fator de Bayes obtidos na comparação entre os modelos com prioris informativas (i) e não informativas (j), para os genótipos de feijão-caupi (Vigna unguiculata) semiprostrado.
Conclusões
-
A abordagem bayesiana proporciona maior acurácia na seleção de genótipos de feijão-caupi (Vigna unguiculata) semiprostrado com alta adaptabilidade e estabilidade fenotípicas.
-
Com base nas prioris informativas, os genótipos MNC99-507G-4, TE97-309G-24, MNC99-542F-7 e BR 17-Gurguéia são classificados como de adaptabilidade específica a ambientes favoráveis.
-
Os genótipos TE96-290-12G, MNC99-510F-16, MNC99-508G-1, MNC99-541F-21, MNC99-542F-5 e MNC99-547F-2 apresentam adaptabilidade específica a ambientes desfavoráveis.
Agradecimentos
À Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes) e à Fundação Arthur Bernardes (Funarbe), pela concessão de bolsas.
Referências
- BARROS, M.A.; ROCHA, M. de M.; GOMES, R.L.F.; SILVA, K.J.D. e; NEVES, A.C. das. Adaptabilidade e estabilidade produtiva de feijão-caupi de porte semiprostrado. Pesquisa Agropecuária Brasileira, v.48, p.403-410, 2013. DOI: 10.1590/S0100-204X2013000400008.
» https://doi.org/10.1590/S0100-204X2013000400008. - CARVALHO, H.W.L. de; BRITO NETO, J.; FREIRE FILHO, F.R.; ROCHA, M. de M.; RIBEIRO, V.Q.; RODRIGUES, A.R. dos S.; RIBEIRO, S.S.; NOGUEIRA, L.C. Adaptabilidade e estabilidade de linhagens de feijão-caupi de porte prostrado em Sergipe e Alagoas. In: CONGRESSO NACIONAL DE FEIJÃO-CAUPI, 2006, Teresina. Anais eletrônicos... Teresina: Embrapa Meio Norte, 2006. Disponível em: <Disponível em: http://www.cpamn.embrapa.br/anaisconac2006/listaresumo.htm >. Acesso em: 3 fev. 2015.
» http://www.cpamn.embrapa.br/anaisconac2006/listaresumo.htm - CASSELLA, G.; GEORGE, E.I. Explaining the Gibbs sampler. The American Statistician, v.46, p.167-174, 1992.
- COTES, J.M.; CROSSA, J.; SANCHES, A.; CORNELIUS, P.A. A Bayesian approach for assessing the stability of genotypes. Crop Science, v.46, p.2654-2665, 2006. DOI: 10.2135/cropsci2006.04.0227.
» https://doi.org/10.2135/cropsci2006.04.0227. - COUTO, M.F.; NASCIMENTO, M.; AMARAL, A.T. do; SILVA, F.F. e; VIANA, A.P.; VIVAS, M. Eberhart and Russel's Bayesian method in the selection of popcorn cultivars. Crop Science, v.55, p.571-577, 2015. DOI: 10.2135/cropsci2014.07.0498.
» https://doi.org/10.2135/cropsci2014.07.0498. - CRUZ, C.D.; REGAZZI, A.J. Modelos biométricos aplicados ao melhoramento genético. Viçosa: Imprensa Universitária, 2007. 480p.
- EBERHART, S.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, v.6, p.36-40, 1966. DOI: 10.2135/cropsci1966.0011183X000600010011x.
» https://doi.org/10.2135/cropsci1966.0011183X000600010011x. - FREIRE FILHO, F.R.; RIBEIRO, V.Q.; ROCHA, M. de M.; LOPES, Â.C. de A. Adaptabilidade e estabilidade da produtividade de grãos de linhagens de caupi de porte enramador. Revista Ceres, v.49, p.383-393, 2002.
- GEWEKE, J. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments. In: BERNARDO, J.M.; BERGER, J.O.; DAVID, A.P.; SMITH, A.F.M. (Ed.). Bayesian statistics. New York: Oxford University, 1992. p.625-631.
- HEIDELBERGER, P.; WELCH, P.D. Simulation run length control in the presence of an initial transient. Operations Research, v.31, p.1109-1144, 1983. DOI: 10.1287/opre.31.6.1109.
» https://doi.org/10.1287/opre.31.6.1109. - JEFFREYS, H. Theory of probability. 3rd ed. Oxford: Claredon, 1961. 447p.
- KASS, R.E.; RAFTERY, A.E. Bayes factors. Journal of the American Statistical Association, v.90, p.773-795, 1995. DOI: 10.1080/01621459.1995.10476572.
» https://doi.org/10.1080/01621459.1995.10476572. - NASCIMENTO, M.; SILVA, F.F.; SÁFADI, T.; NASCIMENTO, A.C.M.C.; FERREIRA, R. de P.; CRUZ, C.D. Abordagem bayesiana para avaliação da adaptabilidade e estabilidade de genótipos de alfafa. Pesquisa Agropecuária Brasileira, v.46, p.26-32, 2011. DOI: 10.1590/S0100-204X2011000100004.
» https://doi.org/10.1590/S0100-204X2011000100004. - RAFTERY, A.E.; LEWIS, S.M. Comment: one long run with diagnostics: implementation strategies for Markov chain Monte Carlo. Statistical Science, v.7, p.493-497, 1992. DOI: 10.1214/ss/1177011143.
» https://doi.org/10.1214/ss/1177011143. - RAMOS MOLINA, L.M.; SANCHE, A.; COTES TORRES, J.M. Inferência bayesiana na análise de testes regionais de arroz em dois sistemas de cultivo. Revista Facultad Nacional de Agronomía, v.64, p.5883-5891, 2012.
- ROCHA, M. de M.; FREIRE FILHO, F.R.; RIBEIRO, V.Q.; CARVALHO, H.W.L. de; BELARMINO FILHO, J.; RAPOSO, J.A.A.; ALCÂNTARA, J. dos P.; RAMOS, S.R.R.; MACHADO, C. de F. Adaptabilidade e estabilidade produtiva de genótipos de feijão-caupi de porte semi-ereto na Região Nordeste do Brasil. Pesquisa Agropecuária Brasileira, v.42, p.1283-1289, 2007. DOI: 10.1590/S0100-204X2007000900010.
» https://doi.org/10.1590/S0100-204X2007000900010. - ROUDER, J.N.; MOREY, R.D.; SPECKMAN, P.L.; PROVINCE, J.M. Default Bayes factors for ANOVA designs. Journal of Mathematical Psychology, v.56, p.356-374, 2012. DOI: 10.1016/j.jmp.2012.08.001.
» https://doi.org/10.1016/j.jmp.2012.08.001. - SANTOS, J.A. da S.; SOARES, C.M.G.; CORRÊA, A.M.; TEODORO, P.E.; RIBEIRO, L.P.; ABREU, H.K.A. de. Agronomic performance and genetic dissimilarity among cowpea (Vigna unguiculata (L.) Walp.) genotypes. Global Advanced Research Journal of Agricultural Science, v.3, p.271-277, 2014.
- THE R FOUNDATION. R: the R project for statistical computing. Vienna, 2010. Available at: <Available at: https://www.r-project.org >. Accessed on: 15 Apr. 2015.
» https://www.r-project.org - TORRES, F.E.; SAGRILO, E.; TEODORO, P.E. ; RIBEIRO, L.P.; CARGNELUTTI FILHO, A. Número de repetições para avaliação de caracteres em genótipos de feijão-caupi. Bragantia, v.74, p.161-168, 2015. DOI: 10.1590/1678-4499.0393.
» https://doi.org/10.1590/1678-4499.0393. - VALADARES, R. de N.; MOURA, M. da C.C.L.; SILVA, A.F.A. da; SILVA, L.S. da; VASCONCELOS, M. da C.C. de A.; SILVA, R.G. Adaptabilidade e estabilidade fenotípica em genótipos de feijão-caupi (Vigna unguiculata (L.) Walp.) de porte ereto/semi-ereto nas mesorregiões leste e sul maranhense. Agropecuária Científica no Semi-árido, v.6, p.21-27, 2010.
- WRICKE, G. Zur berechnung der okovalenz bei sommerweizen und hafer. Zeitschrift für Pflanzenzüchtung, v.52, p.127-138, 1965.
- YATES, F.; COCHRAN, W.G. The analysis of group of experiments. The Journal of Agricultural Science, v.28, p.556-580, 1938. DOI: 10.1017/S0021859600050978.
» https://doi.org/10.1017/S0021859600050978.
Datas de Publicação
-
Publicação nesta coleção
Out 2015
Histórico
-
Recebido
13 Abr 2015 -
Aceito
11 Ago 2015