Seleção precoce para produção de grãos em café arábica pela avaliação de caracteres morfológicos

Early selection for grain yield in Arabica coffee using morphological traits evaluation

GENÉTICA

Seleção precoce para produção de grãos em café arábica pela avaliação de caracteres morfológicos

Early selection for grain yield in Arabica coffee using morphological traits evaluation

Alexsandro Lara Teixeira^I; Flávia Maria Avelar Gonçalves^II; Juliana Costa de Rezende^III; Samuel Pereira de Carvalho^II; Antônio Alves Pereira^III; Bráulio Fabiano Xavier de Moraes^II; Luíz Gustavo Vieira Teixeira^II

^IEmbrapa Rondônia, BR 364, Km 5,5, Zona Rural, Caixa Postal 127, CEP 76815‑800 Porto Velho, RO. E‑mail: [e-mail]

^IIUniversidade Federal de Lavras, Caixa Postal 3037, CEP 37200‑000 Lavras, MG. E‑mail: [e-mail], [e-mail], [e-mail], [e-mail]

IIIEpamig, Campus da Ufla, Rodovia Lavras‑Ijaci, Km 02, Caixa Postal 176, CEP 37200‑000 Lavras, MG. E‑mail: [e-mail], [e-mail]

RESUMO

O objetivo deste trabalho foi avaliar a viabilidade da seleção precoce em café arábica, por meio do emprego de correlações entre caracteres morfológicos e de produção de grãos. Foram avaliados 269 acessos de Coffea arabica, entre cultivares, híbridos e genótipos selvagens, oriundos do banco de germoplasma de café da Epamig, em Patrocínio, MG. Utilizou-se o delineamento em blocos ao acaso, com duas repetições, no espaçamento 3,5x0,8 m, com parcelas de dez plantas. Os caracteres morfológicos foram avaliados em 2006, 12 meses após a implantação da cultura, com as plantas ainda na fase juvenil. Os dados de produção de grãos foram referentes às safras de 2008/2009 e 2009/2010. O comprimento do primeiro ramo plagiotrópico, aos 12 meses de idade, apresentou alta correlação e alto efeito direto, em sentido favorável à produção de grãos, tendo sido o principal caráter morfológico responsável pela variação na produção de grãos. O vigor vegetativo, aos 12 meses de idade, foi a variável secundária de maior importância para explicar variações na produção de grãos. Esses caracteres podem ser utilizados efetivamente na seleção precoce para produção de grãos em café arábica.

Termos para indexação:Coffea arabica, análise de trilha, comprimento de ramo, melhoramento do cafeeiro, vigor vegetativo.

ABSTRACT

The objective of this work was to assess the viability of early selection in Arabica coffee using the correlation between morphological traits and grain yield. Two hundred and sixty‑nine accessions of Coffea arabica, among cultivars, hybrids, and wild genotypes, from Epamig's coffee germplasm bank, Patrocínio, MG, Brazil, were evaluated. A randomized complete block design was used, with two replicates, at a spacing of 3.5x0.8 m, with ten plants per plot. Morphological traits were evaluated in 2006, 12 months after planting, when plants were still at the juvenile phase. Grain yield data are from the 2008/2009 and 2009/2010 crop seasons. The length of the first plagiotropic stem, at 12 months of age, was highly correlated and showed a high direct effect in a favorable direction toward grain yield, being the main morphological trait responsible for variations in grain yield. The vegetative vigor at 12 months of age was the most important secondary variable in explaining the variations in grain yield. These traits can be used effectively in the early selection for grain yield in Arabica coffee.

Index terms:Coffea arabica, path analysis, stem length, coffee breeding, vegetative vigor.

Introdução

A busca por cultivares mais produtivas tem sido o principal foco dos programas de melhoramento do cafeeiro no Brasil (Carvalho et al., 1988, 2010; Martinez et al., 2007). O emprego da seleção indireta, com uso de caracteres correlacionados à produção de grãos, é uma estratégia que vem sendo utilizada em programas de melhoramento, pois permite maximizar os ganhos com a seleção. Essa ferramenta proporciona, aos melhoristas, uma melhor orientação na escolha das variáveis a serem utilizadas no momento da seleção (Cruz et al., 2004).

Estimativas de correlações positivas entre produtividade e caracteres morfológicos, em cafeeiro, têm sido encontradas em diversos trabalhos (Carvalho et al., 2010). Apesar de auxiliar no entendimento de um caráter complexo, a análise da correlação não permite determinar a importância dos efeitos diretos e indiretos dos caracteres que a compõem. Assim, não é possível constatar se sua estimativa foi estabelecida por verdadeiras relações de causa e efeito, e a interpretação direta de suas magnitudes pode resultar em equívocos na estratégia de seleção, pois a alta correlação entre dois caracteres pode ser resultado do efeito de um terceiro caráter, ou de um grupo de caracteres, sobre eles (Wright, 1921).

A análise de trilha consiste no estudo dos efeitos diretos e indiretos de caracteres sobre uma variável principal, cujas estimativas são obtidas por meio de equações de regressão, em que as variáveis são previamente padronizadas. Embora a correlação seja uma característica intrínseca a dois caracteres em dada condição experimental, sua decomposição depende do conjunto de caracteres estudados, que normalmente são avaliados pelo conhecimento prévio de suas importâncias e de possíveis inter‑relações expressas em diagramas de trilha (Cruz et al., 2004). Esta técnica tem sido amplamente utilizada no melhoramento de diversas culturas; entretanto, são poucos os trabalhos desta natureza na cultura do cafeeiro.

Ao analisar os efeitos diretos e indiretos de alguns caracteres morfológicos sobre a produtividade de grãos em café, Severino et al. (2002) constataram que o vigor vegetativo e a seca de ponteiros são bons critérios de avaliação do potencial produtivo. Já em estudos com capim‑elefante, Daher et al. (2004) verificaram que caracteres morfológicos, como número de perfilhos por metro linear e diâmetro de perfilhos a 10 cm do solo, foram capazes de melhor explicar o potencial de produção de matéria seca, tendo atuado de forma direta e inversamente proporcional sobre a variável principal, de acordo com as condições ambientais durante o crescimento. Segundo estes autores, a análise de trilha pode auxiliar o melhorista na escolha dos caracteres a serem utilizados na seleção para produtividade de grãos, com vistas à obtenção de maiores ganhos com a seleção.

A identificação de caracteres com maior efeito direto no sentido favorável à seleção é primordial para que a resposta correlacionada, por meio da seleção indireta, seja eficiente. Além disso, em se tratando de plantas perenes, a existência de correlação entre caracteres na fase juvenil e na fase produtiva pode tornar viável a realização da seleção precoce (Freitas et al., 2007).

O objetivo deste trabalho foi avaliar a viabilidade de utilização da seleção precoce em café arábica (Coffea arabica L.), por meio do emprego de correlações entre caracteres morfológicos e de produção de grãos.

Material e Métodos

O experimento foi realizado em Patrocínio, MG (18º56'38"S e 46º59'34"W, a 972 m de altitude média). O clima da região, conforme a classificação de Köppen, é do tipo Cwa, tropical de altitude, com temperaturas médias anuais de 20,2ºC e precipitação pluvial média de 1.620 mm, com 65 a 70% desse total concentrados no período de dezembro a março.

Foram avaliados 269 acessos de C. arabica, representados por cultivares, híbridos e alguns genótipos selvagens introduzidos de outros países, oriundos do banco de germoplasma de café, instalado no ano de 2005, na fazenda experimental da Empresa de Pesquisa Agropecuária de Minas Gerais, em Patrocínio, MG. Utilizou-se o delineamento em blocos ao acaso, com duas repetições, no espaçamento de 3,5x0,8 m, com parcelas de dez plantas. Os caracteres morfológicos foram avaliados em 2006, 12 meses após a implantação da cultura, com as plantas ainda na fase juvenil. Os dados de produção de grãos foram referentes às safras de 2008/2009 e 2009/2010. Os experimentos foram conduzidos de acordo com as recomendações de adubação para a cultura do cafeeiro. Foram adotadas as práticas de manejo usualmente empregadas na cultura.

Os caracteres avaliados foram: vigor vegetativo, tendo-se atribuído notas de 1 a 10 para sua avaliação, em que 1 foi utilizado para plantas depauperadas, e 10 para plantas muito vigorosas; altura da planta (cm), medida do colo da planta até a gema apical do caule; número de pares de ramos plagiotrópicos (NPRP), avaliado por meio da contagem de todos os pares de ramos laterais primários que apresentaram comprimento superior a 5 cm; comprimento do primeiro ramo plagiotrópico (CRP), avaliado por meio da medição do primeiro ramo plagiotrópico acima do colo da planta, em cm; número de nós do primeiro ramo plagiotrópico (NRP), obtido pela contagem de todos os nós do ramo; comprimento do quarto par de folhas (CPF), em cm; largura do quarto par de folhas (LPF), em cm; e produção de grãos (café beneficiado), tendo-se determinado a produtividade de grãos por parcela, com posterior conversão para sacas por hectare.

Com a finalidade de atender às pressuposições da análise de variância, foram realizados os testes de homogeneidade e normalidade de variância do erro e, após verificação de que os dados atendiam às pressuposições da análise de variância, foram obtidas as estimativas de correlações genotípicas e do coeficiente de determinação genotípico.

Foi avaliado o grau de multicolinearidade da matriz de correlação das médias X'X (Montgomery & Peck, 1981), e, para detectar as variáveis que contribuíram para o surgimento da multicolinearidade, efetuou-se a análise dos autovalores e autovetores. De acordo com Montgomery & Peck (1981), valores de NC<100 não constituem problema sério (multicolinearidade fraca). Se 100<NC<1.000, a multicolinearidade é de moderada a forte, e, quando NC>1.000, há fortes indícios de multicolinearidade severa. Para contornar os efeitos da multicolinearidade sem que houvesse necessidade de descarte de variáveis, foi utilizado o método de regressão em crista (Hoerl & Kennard, 1970), em que se adicionou uma variável (k=0,15) à diagonal da matriz X'X. Com esse procedimento, foi possível obter um número de condição (NC) menor que 100, que caracteriza multicolinearidade fraca e não compromete a análise de trilha. Posteriormente, foram realizados os desdobramentos das correlações genotípicas em efeitos diretos e indiretos por meio da análise de trilha (Wright, 1921).

Para análise dos dados, adotou-se um diagrama causal em cadeia (Figura 1). O critério adotado na elaboração do diagrama foi baseado na biologia reprodutiva do cafeeiro, tendo-se estipulado, como variáveis primárias, as diretamente relacionadas com a produção de grãos, e, como variáveis secundárias, o restante. Dessa forma, a produção de grãos foi considerada a variável principal, com NPRP, CRP e NRP, como variáveis primárias, e vigor, altura, CPF e LPF, como variáveis secundárias.

Para a avaliação dos efeitos diretos e indiretos dos caracteres secundários sobre cada caráter primário, o modelo apropriado e os sistemas de equações derivados para cada variável explicativa foram utilizados. O estudo dos efeitos dos componentes secundários sobre a variável principal produção foi realizado por meio dos desdobramentos da correlação entre a variável principal e as variáveis secundárias (Cruz et al., 2004).

Os dados obtidos referentes aos caracteres morfológicos e à produtividade de grãos foram submetidos à análise de variância, com a significância dos efeitos verificada pelo teste F, a 5% de probabilidade. A acurácia seletiva (ˆr_ˆgg), determinada por meio da expressão: ˆr_ˆgg = (1‑1/F)^0,5, em que F é o valor do teste F de Snedecor para o efeito de genótipo (Resende & Duarte, 2007), foi estimada para aferir a precisão experimental. As análises de variância foram realizadas com o pacote estatístico SAS (SAS Institute, 1999). As estimativas de herdabilidade no sentido amplo foram obtidas pelas variâncias genéticas (s²_g) e fenotípicas (s²_f), por meio da equação h²_a = s²_g/s²_f. As correlações genotípicas e a análise de trilha foram realizadas com uso do programa Genes (Cruz, 2006).

Resultados e Discussão

Todas as características estudadas apresentaram diferenças significativas entre genótipos (Tabela 1). As estimativas de acurácia foram de alta magnitude (65< ˆr_ˆgg<90), o que indica boa precisão experimental. O uso da acurácia, como medida de precisão experimental, conforme Resende & Duarte (2007), tem a vantagem de não depender da magnitude da média, o que proporciona maior segurança na utilização da expressão fenotípica como indicador da variação genotípica. Valores de acurácia acima de 70% indicam alta precisão experimental.

A média da análise conjunta para produção de grãos foi de 21,70 sacas por hectare. A herdabilidade observada para o caráter foi baixa (h² = 45%). A baixa média registrada para produção já era esperada, em razão da alta variabilidade genética presente entre os acessos, que apresentaram produtividades de 2,70 até 63,80 sacas por hectare. Entretanto, os valores obtidos para herdabilidade são aceitáveis. Valores acima de 70% para este parâmetro são encontrados em experimentos com cafeeiros em gerações mais avançadas ou mesmo em ensaios de valor de cultivo e uso (Carvalho et al., 1988).

Com relação aos caracteres morfológicos, a média observada para vigor vegetativo (7,86) pode ser considerada alta, o que evidencia, no geral, que as plantas sob avaliação apresentavam boas condições vegetativas. Quanto à herdabilidade, o caráter altura de planta foi o que apresentou maior valor (75%), o que mostra que, para essa característica, é possível obter maiores ganhos genéticos com a seleção (Tabela 1).

As maiores correlações foram observadas entre a produção de grãos e as variáveis comprimento do primeiro ramo plagiotrópico, vigor vegetativo, largura do quarto par de folhas e número de nós do primeiro ramo plagiotrópico (Tabela 2). Portanto, plantas de maior vigor aos 12 meses de idade e que apresentem o primeiro ramo plagiotrópico longo e com grande número de nós tendem a ser mais produtivas. O vigor vegetativo também apresentou alta correlação com número de pares de ramos plagiotrópicos, comprimento do primeiro ramo plagiotrópico e número de nós do primeiro ramo plagiotrópico. Estes resultados estão de acordo com Freitas et al. (2007), que constataram que a seleção precoce, quanto ao comprimento de ramos plagiotrópicos aos 12 meses de idade, pode ser utilizada para identificação de genótipos superiores. Estes autores, no entanto, não realizaram a avaliação direta da produção de grãos e basearam suas conclusões apenas na capacidade fotossintética das plantas.

Estimativas de correlação semelhantes foram obtidas por Severino et al. (2002), para vigor vegetativo e produção de grãos (0,90), em dados acumulados durante os três primeiros anos de produção. Esses resultados são indicativos de que a seleção precoce pode ser utilizada na identificação de genótipos com alto potencial produtivo e, ainda, auxiliar o melhorista nas tomadas de decisão em etapas iniciais dos programas de melhoramento.

Ao se analisar as correlações por meio do diagrama causal em cadeia (Figura 1), todas as três variáveis primárias apresentaram coeficientes de correlação acima de 0,65, com a variável primária produção (Tabela 3). Já entre os componentes primários e secundários, a variável secundária vigor apresentou alta correlação genotípica com as variáveis primárias, com estimativas acima de 0,83 (Tabela 4). Para a variável altura, foi observada alta correlação apenas com o caráter CRP. Esse resultado assemelha-se ao encontrado por Carvalho et al. (2010), que observaram alta correlação entre altura da planta e comprimento dos ramos plagiotrópicos (0,92), em plantas já em fase de produção. É interessante ressaltar que, para o caráter altura de planta, os programas de melhoramento têm preconizado plantas de porte médio, em razão das práticas culturais de plantios adensados e colheita mecanizada (Androcioli Filho, 2002), atualmente comuns.

Algumas variáveis, apesar de terem apresentado alta associação com a variável principal, podem não ser a causa determinante das variações sobre o caráter de interesse. Nessas situações, a concentração de esforços na seleção desta variável poderá não resultar em ganhos satisfatórios na variável principal (Cruz et al., 2004). Assim, para um melhor entendimento da inter‑relação entre os caracteres, realizou-se o desdobramento dos coeficientes de correlação (análise de trilha), com o objetivo de fornecer uma medida da influência de cada causa e seu efeito. Os desdobramentos foram realizados em conformidade com o diagrama causal em cadeia (Figura 1).

As magnitudes e os sinais dos efeitos diretos e indiretos das variáveis primárias sobre a variável principal produção encontram-se na Tabela 3. Ao verificar os efeitos diretos e indiretos por meio das vias de associação (Figura 1), é possível observar que apenas a variável CRP teve efeito direto de elevada magnitude. As variáveis NPRP e NRP mostraram pequeno efeito direto sobre a produção, tendo atuado de forma indireta. Já a variável CRP, apresentou efeito direto expressivo sobre a produção. Isso mostra que, apesar de todas essas variáveis primárias apresentarem alta correlação com a produção de grãos (Tabela 2), apenas a variável CRP tem efeito direto considerável sobre esta. Assim, no processo de seleção precoce, o comprimento do primeiro ramo plagiotrópico, aos 12 meses de idade, pode ser utilizado pelo melhorista na identificação de genótipos promissores.

Situações contrárias à exposta podem ocorrer, como observado por Carvalho et al. (2010), que verificaram baixas correlações da produtividade de grãos com as variáveis número e comprimento de ramos plagiotrópicos. Nesse caso, seria interessante o emprego da análise de trilha para verificar as relações de causa e efeito desses caracteres, já que, biologicamente, esses componentes devem estar estreitamente relacionados com a produção de grãos (Cruz et al., 2004).

O coeficiente de determinação do modelo da análise de trilha, em questão, foi elevado (R² = 0,88), o que mostra que as variações na variável principal são, em grande parte, explicadas por esse arranjo causal (Tabela 3). Cabe destacar que, em determinadas situações, valores maiores que o valor absoluto 1,0 podem ocorrer em consequência do modelo utilizado na estimação das variâncias e das covariâncias, que determinam a correlação (Cruz et al., 2004). Por meio da matriz de correlação entre as variáveis primárias e a principal, registrou-se um diagnóstico de multicolinearidade fraca (NC = 23,12), o que não denota problemas nas análises. Na presença de multicolinearidade, as variâncias associadas aos estimadores dos coeficientes de trilha podem atingir valores excessivamente altos, o que os torna pouco confiáveis. Além disso, as estimativas dos parâmetros podem assumir valores absurdos ou sem nenhuma coerência com o fenômeno biológico estudado (Cruz et al., 2004; Rodrigues et al., 2010).

No desdobramento das correlações entre as estimativas dos efeitos diretos e indiretos das quatro variáveis secundárias sobre as três variáveis primárias, conforme o diagrama causal em cadeia (Figura 1), verificou-se que o caráter vigor teve alta magnitude de efeitos diretos sobre as três variáveis primárias, com valores acima de 0,51, tendo chegado a 0,90 sobre NPRP (Tabela 4). O caráter vigor apresentou altas correlações e altos efeitos diretos sobre os três caracteres primários, o que o indica que pode ser utilizado na seleção destes caracteres e, por consequência, na seleção para produção de grãos. A variável secundária altura mostrou efeito direto elevado sobre a variável CRP. Já os caracteres CPF e LPF, apresentaram baixa correlação com NPRP e NRP; porém, seus efeitos indiretos via vigor superaram a magnitude dos seus respectivos efeitos residuais (Tabela 4). A variável CPF foi a que apresentou o maior efeito direto negativo, tendo atuado em sentido desfavorável sobre a variável NRP (‑0,33). Os coeficientes de determinação (R²) da análise de trilha, neste caso, foram maiores que 0,70, o que torna possível inferir que o modelo causal explicou bem as variações encontradas para as três variáveis primárias. O número de condição (NC) obtido foi de 42,01, o que indica multicolinearidade fraca entre as variáveis.

Nas análises do efeito total, pode-se concluir que, entre os caracteres secundários, o vigor foi o que apresentou maior relação com produção (Tabela 5), com efeito direto total de elevada magnitude (0,57). Verifica-se, também, que a influência desse caráter sobre os componentes primários foi diferente para NPRP e CRP. A variável LPF apresentou efeito indireto considerável via vigor (0,37), o que ressalta a importância da variável vigor sobre a variável principal. Severino et al. (2002) também observaram que, além da alta correlação, há grande efeito direto positivo entre o vigor vegetativo e a produção de grãos. Esse fato explica porque o vigor vegetativo está entre as características mais utilizadas para estimação da capacidade produtiva de cafeeiros (Carvalho et al., 2010).

Por meio das variáveis secundárias, pode-se inferir que serão efetivos os ganhos indiretos para produção de grãos, via resposta correlacionada por seleção quanto a vigor. Esse caráter tem como vantagem adicional ter efeito direto positivo sobre CRP, que, entre os componentes primários, é o de maior efeito direto sobre a produção de grãos.

Ao se compararem os resultados obtidos entre as estimativas de correlação genotípica e as do coeficiente de trilha, percebe-se que é possível identificar quais são as variáveis que exercem maior influência sobre o caráter em questão, após os desdobramentos da correlação em efeitos diretos e indiretos. Além disso, a utilização da seleção precoce torna o processo mais dinâmico e eficiente, e concentra esforços nas variáveis que proporcionarão maiores ganhos com a seleção indireta.

Conclusão

O comprimento do primeiro ramo plagiotrópico e o vigor vegetativo, aos 12 meses de idade, podem ser utilizados efetivamente na seleção precoce quanto à produção de grãos, em café arábica.

Recebido em 18 de fevereiro de 2011 e aprovado em 23 de julho de 2012

• ANDROCIOLI FILHO, A. Café adensado: espaçamentos e cuidados no manejo da lavoura. Londrina: Iapar, 2002. 32p. (Iapar. Circular, 121).
• CARVALHO, A.; FAZUOLI, L.C.; MAZZAFERA, P. Produtividade de progênies derivadas de hibridação dos cultivares Laurina e Mundo Novo. Bragantia, v.47, p.213‑222, 1988.
• CARVALHO, A.M. de; MENDES, A.N.G.; CARVALHO, G.R.; BOTELHO, C.E.; GONÇALVES, F.M.A.; FERREIRA, A.D. Correlação entre crescimento e produtividade de cultivares de café em diferentes regiões de Minas Gerais, Brasil. Pesquisa Agropecuária Brasileira, v.45, p.269‑275, 2010.
• CRUZ, C.D. Programa Genes: estatística experimental e matrizes. Viçosa: UFV, 2006. 285p.
• CRUZ, C.D.; REGAZZI, A.J.; CARNEIRO, P.C.S. Análise de trilha. In: CRUZ, C.D.; REGAZZI, A.J.; CARNEIRO, P.C.S. (Ed.). Modelos biométricos aplicados ao melhoramento Viçosa: UFV, 2004. p.180‑193.
• DAHER, R.F.; PEREIRA, A.V.; PEREIRA, M.G.; LÉDO, F.J. da S.; AMARAL JUNIOR, A.T. do; ROCABADO, J.M.A.; FERREIRA, C.F.; TARDIN, F.D. Análise de trilha de caracteres forrageiros do capim‑elefante (Pennisetum purpureum Schum.). Ciência Rural, v.34, p.1531‑1535, 2004.
• FREITAS, Z.M.T.S. de; OLIVEIRA, F.J. de; CARVALHO, S.P. de; SANTOS, V.F. dos; SANTOS, J.P. de O. Avaliação de caracteres quantitativos relacionados com o crescimento vegetativo entre cultivares de café arábica de porte baixo. Bragantia, v.66, p.267‑275, 2007.
• HOERL, A.E.; KENNARD, R.W. Ridge regression: applications to nonorthogonal problems. Technometrics, v.12, p.69‑82, 1970.
• MARTINEZ, H.E.P.; AUGUSTO, H.S.; CRUZ, C.D.; PEDROSA, A.W.; SAMPAIO, N.F. Crescimento vegetativo de cultivares de café (Coffea arabica L.) e sua correlação com a produção em espaçamentos adensados. Acta Scientiarum. Agronomy, v.29, p.481‑489, 2007.
• MONTGOMERY, D.C.; PECK, E.A. Introduction to linear regression analysis New York: John Wiley, 1981. 504p.
• RESENDE, M.D.V. de; DUARTE, J.B. Precisão e controle de qualidade em experimentos de avaliação de cultivares. Pesquisa Agropecuária Tropical, v.37, p.182‑194, 2007.
• RODRIGUES, G.B.; MARIM, B.G.; SILVA, D.J.H. da; MATTEDI, A.P.; ALMEIDA, V. de S. Análise de trilha de componentes de produção primários e secundários em tomateiro do grupo Salada. Pesquisa Agropecuária Brasileira, v.45, p.155‑162, 2010.
• SAS INSTITUTE. SAS/STAT user's guide Version 8. Cary: SAS Institute, 1999.
• SEVERINO, L.S.; SAKIYAMA, N.S.; PEREIRA, A.A.; MIRANDA, G.V.; ZAMBOLIM, L.; BARROS, U.V. Associações da produtividade com outras características agronômicas de café (Coffea arabica L. "Catimor"). Acta Scientiarum. Agronomy, v.24, p.1467‑1471, 2002.
• WRIGHT, S. Correlation and causation. Journal of Agricultural Research, v.20, p.557‑585, 1921.

Datas de Publicação

• Publicação nesta coleção
  17 Out 2012
• Data do Fascículo
  Ago 2012

Histórico

• Recebido
  18 Fev 2011
• Aceito
  23 Jul 2012