RESUMO
Os objetivos deste trabalho foram determinar o tamanho de amostra necessário para a estimação da média de caracteres e avaliar as relações lineares existentes entre caracteres de plantas de Crotalaria spectabilis. Foi conduzido um experimento com C. spectabilis na safra 2014/2015, sendo os seguintes caracteres avaliados em 100 plantas: altura da planta (AP), número de ramos produtivos (NRP), massa seca da parte aérea (MSPA), massa seca de vagens (MSVG), número de vagens (NVG), peso médio de cada vagem (PMVG), número de sementes (NS), número médio de sementes por vagem (NSVG), massa de cem grãos (MCG) e peso total de sementes (PTS). Foram calculadas medidas de tendência central, de dispersão e de distribuição, foi verificada a normalidade dos dados e determinado o tamanho de amostra necessário para a estimação da média de cada caractere, assumindo-se erros de estimação iguais a 1, 2, …, 20% da média. Depois, foram calculados os coeficientes de correlação linear de Pearson e estimados os efeitos diretos e indiretos das variáveis explicativas sobre PTS por meio de análise de trilha. Em C. spectabilis, são necessárias 121 plantas para estimação da média dos caracteres NRP, MSPA, MSVG, NVG, NS e PTS com erro máximo de 10% da média. Já para a estimação da média de AP, PMVG, NSVG e MCG nesse mesmo nível de precisão, são necessárias 6 plantas. A variável NS apresentou elevada correlação e efeito direto sobre PTS, indicando relação de causa e efeito e possibilidade de seleção indireta. As demais variáveis explicativas apresentaram efeitos indiretos sobre PTS via NS.
Palavras-chave
precisão experimental; planejamento experimental; análise de trilha; seleção indireta
ABSTRACT
The aims of this research were to determine the sample size needed to estimate the mean of traits and evaluate the linear relationships between traits of Crotalaria spectabilis plants. An experiment with C. spectabilis was conducted in the 2014/2015 season, and the following characters were evaluated in 100 plants: plant height (AP), number of productive branches (NRP), dry matter of the aerial part (MSPA), dry matter of the pods (MSVG), number of pods (NVG), weight of each pod (PMVG), number of seeds (NS), average number of seeds per pod (NSVG), mass of hundred grains (MCG), and total weight of seeds (PTS). Measures of central tendency, dispersion and distribution were calculated, the normality of the data was checked, and the sample size needed to estimate the mean of each trait was determined, assuming estimation errors equal to 1, 2, …, 20% of the mean. Then, the Pearson’s linear correlation coefficients were calculated, and the direct and indirect effects of the explanatory variables on PTS were estimated by path analysis. In C. spectabilis, 121 plants are required to the mean estimation of NRP, MSPA, MSVG, NVG, NS, and PTS traits, with a maximum estimation error of 10% of the mean. In this same level of accuracy, 6 plants are needed for the estimation of the mean of AP, PMVG, NSVG, and MCG. The variable NS showed high correlation and direct effect on PTS, indicating a relation of cause and effect and the possibility of indirect selection. The other explanatory variables had indirect effects on PTS via NS.
Key words
experimental precision; experimental planning; path analysis; indirect selection
INTRODUÇÃO
A Crotalaria spectabilis Roth pertence à família Fabaceae e juntamente com as demais espécies do gênero Crotalaria são popularmente conhecidas como crotalária, guizo ou chocalho de cascavel, sendo utilizadas como cultura de cobertura ou de adubação verde (Albuquerque et al. 2000Albuquerque, M. C. F., Rodrigues, T. J. D. e Mendonça, E. A. F. (2000). Absorção de água por sementes de Crotalaria spectabilis Roth determinada em diferentes temperaturas e disponibilidade hídrica. Revista Brasileira de Sementes, 22, 206-215.; Leal et al. 2012Leal, M. A. A., Guerra, J. G. M., Peixoto, R. T. G. e Almeida, D. L. (2012). Desempenho de crotalária cultivada em diferentes épocas de semeadura e de corte. Revista Ceres, 59, 386-391. http://dx.doi.org/10.1590/S0034-737X2012000300014.
http://dx.doi.org/10.1590/S0034-737X2012...
). As espécies de crotalária apresentam elevada capacidade de produção de biomassa e fixação de nitrogênio (Teodoro et al. 2011Teodoro, R. B., Oliveira, F. L., Silva, D. M. N., Fávero, C. e Quaresma, M. A. L. (2011). Aspectos agronômicos de leguminosas para adubação verde no Cerrado do Alto Vale do Jequitinhonha. Revista Brasileira de Ciência do Solo, 35, 635-640. http://dx.doi.org/10.1590/S0100-06832011000200032.
http://dx.doi.org/10.1590/S0100-06832011...
, 2015Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
). Ainda, de acordo com Araújo et al. (2015)Araújo, L. S., Cunha, P. C. R., Silveira, P. M., Sousa Netto, M. e Oliveira, F. C. (2015). Potencial de cobertura do solo e supressão de tiririca (Cyperus rotundus) por resíduos culturais de plantas de cobertura. Revista Ceres, 62, 483-488. http://dx.doi.org/10.1590/0034-737X201562050009.
http://dx.doi.org/10.1590/0034-737X20156...
, a cobertura morta com resíduos culturais de C. spectabilis apresenta interferência na dinâmica populacional da tiririca (Cyperus rotundus), com redução do número e da fitomassa seca dessa planta daninha. Contudo, atualmente, existem poucas informações na literatura sobre a importância da C. spectabilis, o que intensifica a necessidade de pesquisas sobre a cultura.
Para se obterem resultados confiáveis em pesquisas envolvendo a C. spectabilis e demais culturas agrícolas, é indispensável o correto dimensionamento do tamanho da amostra (número de plantas) que deverá ser utilizado. Nesse sentido, Bussab e Morettin (2011)Bussab, W. O., Morettin, P. A. (2011). Estatística básica. 7. ed. São Paulo: Saraiva. destacam que o tamanho de amostra é inversamente proporcional ao erro permitido a priori pelo pesquisador, sendo maior de acordo com o aumento da variabilidade dos dados e da confiabilidade desejada. Estudos de dimensionamento amostral já foram realizados para a avaliação de caracteres em plantas de cobertura como Crambe abyssinica (Cargnelutti Filho et al. 2010Cargnelutti Filho, A., Toebe, M., Silveira, T. R., Casarotto, G., Haesbaert, F. M. e Lopes, S. J. (2010). Tamanho de amostra e relações lineares de caracteres morfológicos e produtivos de crambe. Ciência Rural, 40, 2262-2267. http://dx.doi.org/10.1590/S0103-84782010001100003.
http://dx.doi.org/10.1590/S0103-84782010...
), nabo forrageiro (Cargnelutti Filho et al. 2014Cargnelutti Filho, A., Facco, G., Lúcio, A. D., Toebe, M., Burin, M., Fick, A. L. e Neu, I. M. M. (2014). Tamanho de amostra para a estimação da média de caracteres morfológicos e produtivos de nabo forrageiro. Ciência Rural, 44, 223-227. http://dx.doi.org/10.1590/S0103-84782014000200005.
http://dx.doi.org/10.1590/S0103-84782014...
), tremoço branco (Burin et al. 2014Burin, C., Cargnelutti Filho, A., Toebe, M., Alves, B. M. e Fick, A. L. (2014). Dimensionamento amostral para a estimação da média de caracteres de tremoço branco. Comunicata Scientiae, 5, 205-212.) e aveia preta (Cargnelutti Filho et al. 2015Cargnelutti Filho, A., Toebe, M., Alves, B. M., Burin, C., Santos, G. O., Facco, G. e Neu, I. M. M. (2015). Dimensionamento amostral para avaliar caracteres morfológicos e produtivos de aveia preta em épocas de avaliação. Ciência Rural, 45, 9-13. http://dx.doi.org/10.1590/0103-8478cr20140504.
http://dx.doi.org/10.1590/0103-8478cr201...
). Em Crotalaria juncea e Crotalaria spectabilis, foram determinados tamanhos de amostra para a estimação da média de massas fresca e seca e da produtividade (Teodoro et al. 2015Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
).
Por sua vez, o estudo das associações entre caracteres é importante em programas de melhoramento pela possibilidade da seleção indireta de plantas com maior produção de massa, de grãos ou nos casos em que a variável de interesse é de difícil mensuração (Cruz e Regazzi 1997Cruz, C. D. e Regazzi, A. J. (1997). Modelos biométricos aplicados ao melhoramento genético. 2. ed. Viçosa: UFV.). Para o estudo da associação entre pares de caracteres, pode-se utilizar o coeficiente de correlação linear de Pearson (r), que assume valores no intervalo de −1 ≤ r ≤1 (Ferreira 2009Ferreira, D. F. (2009). Estatística básica. 2. ed. Lavras: UFLA.). Segundo o autor, o coeficiente de correlação mede a força e a direção da associação, podendo esta ser positiva, negativa, ou, ainda, existir ausência de relacionamento linear entre pares de variáveis.
Quando um conjunto de variáveis (caracteres) é avaliado simultaneamente, a quantificação e a interpretação da magnitude de correlações específicas podem resultar em inferências equivocadas, devido a efeitos de outras variáveis não consideradas na correlação (Cruz e Regazzi 1997Cruz, C. D. e Regazzi, A. J. (1997). Modelos biométricos aplicados ao melhoramento genético. 2. ed. Viçosa: UFV.). Assim, a análise de trilha é uma técnica multivariada importante, que permite o desdobramento dos coeficientes de correlação em efeitos diretos e indiretos sobre a variável principal (Wright 1921Wright, S. (1921). Correlation and causation. Journal of Agricultural Research, 20, 557-585., 1923Wright, S. (1923). The theory of path coefficients: a reply to Niles’s criticism. Genetics, 8, 239-255., 1934Wright, S. (1934). The method of path coefficients. The Annals of Mathematical Statistics, 5, 161-215.). Essa técnica pode contribuir para a verificação das relações de causa e efeito de caracteres explicativos sobre a variável principal, indicando quais caracteres podem ser utilizados na seleção indireta (Cruz e Regazzi 1997Cruz, C. D. e Regazzi, A. J. (1997). Modelos biométricos aplicados ao melhoramento genético. 2. ed. Viçosa: UFV.; Cruz e Carneiro 2006Cruz, C. D. e Carneiro, P. C. S. (2006). Modelos biométricos aplicados ao melhoramento genético. v. 2. 2. ed. Viçosa: UFV.). Estudos de relações lineares por meio de correlações e análise de trilha já foram desenvolvidos em culturas de cobertura e de produção de grãos, como em crambe (Cargnelutti Filho et al. 2010Cargnelutti Filho, A., Toebe, M., Silveira, T. R., Casarotto, G., Haesbaert, F. M. e Lopes, S. J. (2010). Tamanho de amostra e relações lineares de caracteres morfológicos e produtivos de crambe. Ciência Rural, 40, 2262-2267. http://dx.doi.org/10.1590/S0103-84782010001100003.
http://dx.doi.org/10.1590/S0103-84782010...
).
Diante do exposto, os objetivos deste trabalho foram determinar o tamanho de amostra necessário para a estimação da média e avaliar as relações lineares existentes entre caracteres de plantas de Crotalaria spectabilis.
MATERIAL E MÉTODOS
Foi conduzido um experimento (ensaio de uniformidade) com Crotalaria spectabilis na safra 2014/2015, em área experimental da Universidade Federal do Pampa, campus Itaqui, situada no município de Itaqui, Estado do Rio Grande do Sul, nas coordenadas geográficas lat 29°09′S e long 56°33′W e altitude de 74 m. Pela classificação de Köppen, o clima da região é do tipo Cfa, subtropical úmido, com verões quentes e sem estação seca definida, sendo o solo classificado como Plintossolo Háplico (Embrapa 2013Empresa Brasileira de Pesquisa Agropecuária (2013). Sistema brasileiro de classificação de solos. Brasília: Embrapa.). A semeadura foi realizada em linha, no dia 16/10/2014, utilizando-se 33 sementes por metro linear e espaçamento entre linhas de 0,45 m, num ensaio de 8,1 m de comprimento × 8,1 m de largura (65,61 m2). A adubação de base foi de 25 kg∙ha−1 de N, 100 kg∙ha−1 de P2O5 e 100 kg∙ha−1 de K2O. Todos os tratos culturais foram realizados de modo uniforme dentro da área amostral, também considerada como sendo experimento em branco sem tratamentos, conforme descrito por Storck et al. (2011)Storck, L., Garcia, D. C., Lopes, S. J. e Estefanel, V. (2011). Experimentação vegetal. 3. ed. Santa Maria: UFSM..
Em 100 plantas aleatoriamente identificadas dentro do ensaio, foram realizadas 2 avaliações (colheitas), sendo a primeira realizada no dia 18/05/2015 (214 dias após a semeadura — DAS) e a segunda realizada no dia 02/06/2015 (229 DAS), em função da deiscência das vagens. Em cada planta avaliada, foram mensurados os seguintes caracteres: altura da planta (AP, em cm, avaliada aos 229 DAS), número de ramos produtivos (NRP, em unidades, avaliado aos 229 DAS), massa seca da parte aérea (MSPA, em g, corresponde à massa seca da parte aérea sem vagens após a colheita aos 229 DAS adicionada da massa seca de vagens obtida pela soma das avaliações após as colheitas aos 214 e 229 DAS), massa seca de vagens (MSVG, em g, soma das avaliações após as colheitas aos 214 e 229 DAS), número de vagens (NVG, em unidades — soma das avaliações aos 214 e 229 DAS), peso médio de cada vagem (PMVG, em g, obtido pela divisão MSVG/NVG), número de sementes (NS, em unidades, soma das avaliações aos 214 e 229 DAS), número médio de sementes por vagem (NSVG, em unidades, obtido pela divisão NS/NV), peso total de sementes (PTS, em g, soma das avaliações após as colheitas aos 214 e 229 DAS) e massa de 100 grãos (MCG, em g, obtida pela relação 100 × PTS/NS).
Para cada caractere mensurado, foram calculados os valores mínimo e máximo, média, mediana, variância, desvio padrão, erro padrão, coeficiente de variação, curtose, assimetria e o p-valor do teste de normalidade de Kolmogorov-Smirnov (Massey 1951Massey, F. J. (1951). The Kolmogorov-Smirnov test for goodness of fit. Journal of the American Statistical Association, 46, 68-78.). Após, foi determinado o tamanho de amostra necessário para a estimação da média de cada caractere, assumindo-se erros de estimação (em % da média) iguais a 1, 2, …, 20% da estimativa da média, com grau de confiança (1 − α) de 95%. Para isso, utilizou-se a equação (Bussab e Morettin 2011Bussab, W. O., Morettin, P. A. (2011). Estatística básica. 7. ed. São Paulo: Saraiva.):
η = (t2α/2 s2)/(erro de estimação)2
onde: tα/2 é o valor crítico da distribuição t de Student, cuja área à direita é igual a α/2, com (η − 1) graus de liberdade; α = 5% de probabilidade de erro; s2 representa a estimativa da variância.
A seguir, fixou-se η em 100 plantas (número de plantas avaliadas), e foi calculado o erro de estimação em percentagem da estimativa da média (X) para cada caractere considerando-se o total de plantas avaliadas, por meio da equação:
erro de estimação =
onde: s é a estimativa do desvio padrão.
Após a determinação do tamanho de amostra, foram construídos histogramas de frequência e gráficos de dispersão entre os 10 caracteres avaliados. Foi calculada a matriz de coeficientes de correlação linear de Pearson (r) entre os caracteres, sendo a significância dos coeficientes verificada por meio do teste t de Student, a 5% de probabilidade de erro, conforme descrito por Ferreira (2009)Ferreira, D. F. (2009). Estatística básica. 2. ed. Lavras: UFLA.. A seguir, foi realizado o diagnóstico de multicolinearidade na matriz X’X, que representa a matriz de correlação entre os 9 caracteres explicativos (AP, NRP, MSPA, MSVG, NVG, PMVG, NS, NSVG e MCG), conforme critérios estabelecidos por Montgomery e Peck (1982)Montgomery, D. C. e Peck, E. A. (1982). Introduction to linear regression analysis. New York: John Wiley & Sons.. Assim, se o número de condição (NC) obtido pela razão entre o maior e o menor autovalor da matriz de correlação X’X foi NC ≤ 100, considerou-se multicolinearidade fraca entre os caracteres explicativos; 100 < NC < 1.000, multicolinearidade moderada a severa; e NC ≥ 1.000, multicolinearidade severa. No caso de multicolinearidade severa, procedeu-se à eliminação de caracteres explicativos altamente correlacionados, conforme recomendações de Cruz e Carneiro (2006)Cruz, C. D. e Carneiro, P. C. S. (2006). Modelos biométricos aplicados ao melhoramento genético. v. 2. 2. ed. Viçosa: UFV..
Em seguida, foi realizada a análise de trilha de PTS em função dos caracteres explicativos, em único diagrama causal, conforme metodologia proposta por Wright (1921Wright, S. (1921). Correlation and causation. Journal of Agricultural Research, 20, 557-585., 1923Wright, S. (1923). The theory of path coefficients: a reply to Niles’s criticism. Genetics, 8, 239-255., 1934)Wright, S. (1934). The method of path coefficients. The Annals of Mathematical Statistics, 5, 161-215. e descrita por Cruz e Carneiro (2006)Cruz, C. D. e Carneiro, P. C. S. (2006). Modelos biométricos aplicados ao melhoramento genético. v. 2. 2. ed. Viçosa: UFV. e Cruz e Regazzi (1997)Cruz, C. D. e Regazzi, A. J. (1997). Modelos biométricos aplicados ao melhoramento genético. 2. ed. Viçosa: UFV., utilizando-se, para isso, o sistema de equações normais:
onde: X’X representa a matriz de correlação entre caracteres explicativos; X’Y representa o vetor de correlações de cada caractere explicativo com PTS; é o estimador dos efeitos diretos dos caracteres explicativos sobre PTS.
As análises estatísticas foram realizadas com auxílio do programa GENES (Cruz 2013Cruz, C. D. (2013). GENES — a software package for analysis in experimental statistics and quantitative genetics. Acta Scientiarum. Agronomy, 35, 271-276. http://dx.doi.org/10.4025/actasciagron.v35i3.21251.
http://dx.doi.org/10.4025/actasciagron.v...
) e do Microsoft® Office Excel.
RESULTADOS E DISCUSSÃO
Os valores de média e mediana foram similares entre si para os caracteres AP, PMVG, NSVG e MCG (Tabela 1). Já os caracteres NRP, MSPA, MSVG, NVG, NS e PTS apresentaram valores maiores de médias em relação às medianas, indicando padrão assimétrico positivo. Isso se deve ao fato de que algumas das 100 plantas avaliadas apresentaram elevados escores desses caracteres (Figura 1) e, por consequência, contribuíram para o deslocamento assimétrico à direita (Bussab e Morettin 2011Bussab, W. O., Morettin, P. A. (2011). Estatística básica. 7. ed. São Paulo: Saraiva.). Além disso, em 8 dos 10 caracteres avaliados, houve comportamento leptocúrtico, indicando maiores concentrações de plantas na região central da distribuição dos dados e menor número de plantas nas extremidades.
Valores mínimo e máximo, média, mediana, variância, desvio padrão, erro padrão, coeficiente de variação, curtose, assimetria e valor de p do teste de normalidade de Kolmogorov-Smirnov para os caracteres: altura da planta, número de ramos produtivos, massa seca da parte aérea, massa seca de vagens, número de vagens, peso médio de cada vagem, número de sementes, número médio de sementes por vagem, massa de cem grãos e peso total de sementes, avaliados em 100 plantas de Crotalaria spectabilis.
Representação gráfica da distribuição de frequência (na diagonal) e das relações entre os caracteres altura da planta (AP), número de ramos produtivos (NRP), massa seca da parte aérea (MSPA), massa seca de vagens (MSVG), número de vagens (NVG), peso médio de cada vagem (PMVG), número de sementes (NS), número médio de sementes por vagem (NSVG), massa de cem grãos (MCG) e peso total de sementes (PTS) avaliados em 100 plantas de Crotalaria spectabilis.
A média de produção de massa seca de parte aérea (27,317 g∙planta−1 — Tabela 1) foi inferior à média obtida por Teodoro et al. (2015)Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
, de 61,24 g∙ planta−1. O peso total de sementes (produtividade) no presente estudo foi de 8,591 g∙planta−1, sendo que Teodoro et al. (2015)Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
obtiveram valores médios de 4,45 g∙planta−1. Destaca-se que o número de plantas avaliadas no estudo desenvolvido por Teodoro et al. (2015)Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
foi inferior (n = 45 plantas), a densidade de semeadura foi cerca de 50% inferior, sem utilização de fertilizantes, e as condições ambiental e de solo foram distintas, podendo contribuir com as diferenças verificadas em relação às produções de massa seca e de sementes.
Os desvios de assimetria e de curtose contribuíram para o afastamento da distribuição normal (p ≤ 0,036) dos caracteres NRP, MSVG e NVG (Tabela 1). Os demais caracteres apresentaram bom ajuste à distribuição normal (p ≥ 0,072), mesmo apresentando afastamentos de assimetria e curtose. Em relação à variabilidade dos dados, verificou-se que os caracteres AP, PMVG, NSVG e MCG, que apresentaram valores próximos de médias e medianas e melhor ajuste à distribuição normal (p ≥ 0,368), apresentaram menor coeficiente de variação (7,819% ≤ CV ≤ 12,271%). Os demais caracteres (NRP, MSPA, MSVG, NVG, NS e PTS) apresentaram comportamento assimétrico positivo, o que pode ter contribuído para o aumento da variabilidade (45,946% ≤ CV ≤ 55,452%). Com base nos valores de coeficientes de variação e de distribuição dos dados, espera-se a necessidade de maior tamanho de amostra (número de plantas) para a avaliação de NRP, MSPA, MSVG, NVG, NS e PTS em relação ao requerido para a avaliação de AP, PMVG, NSVG e MCG.
O tamanho de amostra para estimação da média de cada caractere, com semiamplitude do intervalo de confiança igual a 1% da estimativa da média e grau de confiança de 95%, oscilou entre 241 e 12.107 plantas, respectivamente, para a estimação da média de MCG e do NS (Tabela 2). Ainda nesse mesmo nível de precisão, menor tamanho de amostra foi necessário para a mensuração de AP, PMVG, NSVG e MCG (número de plantas ≤ 593), e maior tamanho de amostra foi necessário para a mensuração de NRP, MSPA, MSVG, NVG, NS e PTS (número de plantas ≥ 8.311). Esses resultados estão de acordo com a maior variabilidade constatada para NRP, MSPA, MSVG, NVG, NS e PTS (Tabela 1). Segundo Teodoro et al. (2015)Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
, para a avaliação das massas fresca e seca e produção de grãos de C. Juncea e C. Spectabilis nesse nível de precisão, são necessárias de 1.108 a 8.510 plantas. Na cultura de crambe, Cargnelutti Filho et al. (2014)Cargnelutti Filho, A., Facco, G., Lúcio, A. D., Toebe, M., Burin, M., Fick, A. L. e Neu, I. M. M. (2014). Tamanho de amostra para a estimação da média de caracteres morfológicos e produtivos de nabo forrageiro. Ciência Rural, 44, 223-227. http://dx.doi.org/10.1590/S0103-84782014000200005.
http://dx.doi.org/10.1590/S0103-84782014...
verificaram a necessidade de avaliação de 1.516 a 15.405 plantas para a estimação da média de caracteres com erro máximo de 1% da média. Em nabo forrageiro, Cargnelutti Filho et al. (2014)Cargnelutti Filho, A., Facco, G., Lúcio, A. D., Toebe, M., Burin, M., Fick, A. L. e Neu, I. M. M. (2014). Tamanho de amostra para a estimação da média de caracteres morfológicos e produtivos de nabo forrageiro. Ciência Rural, 44, 223-227. http://dx.doi.org/10.1590/S0103-84782014000200005.
http://dx.doi.org/10.1590/S0103-84782014...
verificaram necessidade de avaliação de 217 a 5.751 plantas, dependendo do caractere, para estimação da média com erro de estimação de 2% da média e 95% de confiança.
Considerando-se a dificuldade de mensuração de 12.107 plantas para a estimação dos caracteres de C. Spectabilis com erro de, no máximo, 1% da média, foram determinados tamanhos de amostra para condições de maiores erros de estimação (2, 3, …, 20% da média). Nesse sentido, se o pesquisador desejar estimar a média de AP, PMVG, NSVG e MCG com erros de estimação de, no máximo, 5 e 10% da média, deverá mensurar, respectivamente, 24 e 6 plantas (Tabela 2). Para a mensuração de NRP, MSPA, MSVG, NVG, NS e PTS com erros de estimação de, no máximo, 5 e 10% da média, deverá mensurar, respectivamente, 484 e 121 plantas. Com a utilização de 100 plantas nas avaliações, o erro cometido oscilou entre 1,55 e 11,00%, respectivamente, para MCG e NS. Para a avaliação de massas fresca e seca e produção de sementes de C. spectabilis com erro de estimação de 5 e 10% da média, Teodoro et al. (2015)Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101.
http://dx.doi.org/10.3923/ja.2015.98.101...
recomendaram a avaliação de, respectivamente, 197 e 49 plantas. Já Cargnelutti Filho et al. (2010)Cargnelutti Filho, A., Toebe, M., Silveira, T. R., Casarotto, G., Haesbaert, F. M. e Lopes, S. J. (2010). Tamanho de amostra e relações lineares de caracteres morfológicos e produtivos de crambe. Ciência Rural, 40, 2262-2267. http://dx.doi.org/10.1590/S0103-84782010001100003.
http://dx.doi.org/10.1590/S0103-84782010...
e Cargnelutti Filho et al. (2014)Cargnelutti Filho, A., Facco, G., Lúcio, A. D., Toebe, M., Burin, M., Fick, A. L. e Neu, I. M. M. (2014). Tamanho de amostra para a estimação da média de caracteres morfológicos e produtivos de nabo forrageiro. Ciência Rural, 44, 223-227. http://dx.doi.org/10.1590/S0103-84782014000200005.
http://dx.doi.org/10.1590/S0103-84782014...
indicaram a avaliação de 155 e 231 plantas, respectivamente, para a estimação da média de caracteres de crambe e de nabo forrageiro, com erro de estimação máximo de 10% da média. Em tremoço branco, Burin et al. (2014)Burin, C., Cargnelutti Filho, A., Toebe, M., Alves, B. M. e Fick, A. L. (2014). Dimensionamento amostral para a estimação da média de caracteres de tremoço branco. Comunicata Scientiae, 5, 205-212. verificaram a necessidade de avaliação de 192 plantas para a estimação da média de todos os caracteres com amplitude do intervalo de confiança de 95% de, no máximo, 20% da média.
Tamanho de amostra (número de plantas) para a estimação da média dos caracteres: altura da planta, número de ramos produtivos, massa seca da parte aérea, massa seca de vagens, número de vagens, peso médio de cada vagem, número de sementes, número médio de sementes por vagem, massa de cem grãos e peso total de sementes para os erros de estimação iguais a 1, 2, ..., 20% da estimativa da média e a semiamplitude do intervalo de confiança com base na avaliação de 100 plantas de Crotalaria spectabilis
O coeficiente de correlação linear de Pearson (r) entre os pares de caracteres oscilou de r = −0,192ns a r = 0,991*, respectivamente, para as correlações entre NRP e NSVG e entre NS e PTS (Tabela 3). Entre os caracteres explicativos (AP, NRP, MSPA, MSVG, NVG, PMVG, NS, NSVG e MCG), correlações de maior magnitude foram verificadas entre MSPA, MSVG, NVG e NS (0,865* ≤ r ≤ 0,984*). Já o caractere principal PTS apresentou maiores correlações com NS (r = 0,991*), MSVG (r = 0,967*), NVG (r = 0,956*) e MSPA (r = 0,871*) e correlações positivas e significativas de menor intensidade com AP (r = 0,380*) e NRP (r = 0,214*). Assim, é provável que plantas com maior número de sementes e de vagens e maiores massas seca de vagens e de parte aérea apresentem maior peso total de sementes.
O diagnóstico de multicolinearidade na matriz de correlação entre caracteres explicativos (AP, NRP, MSPA, MSVG, NVG, PMVG, NS, NSVG e MCG) indicou número de condição de 1.113,08, considerado severo de acordo com Montgomery e Peck (1982)Montgomery, D. C. e Peck, E. A. (1982). Introduction to linear regression analysis. New York: John Wiley & Sons.. Após a eliminação de MSVG, altamente correlacionada com MSPA, NVG e NS (Tabela 3), e de NSVG, houve redução do número de condição para 107,32. Dessa forma, a análise de trilha do PTS foi realizada em função dos caracteres explicativos AP, NRP, MSPA, NVG, PMVG, NS e MCG (Tabela 4) e, portanto, o modelo causal foi: PTS = Resíduo.
Coeficientes de correlação linear de Pearson (r) entre os caracteres: altura da planta, número de ramos produtivos, massa seca da parte aérea, massa seca de vagens, número de vagens, peso médio de cada vagem, número de sementes, número médio de sementes por vagem, massa de cem grãos e peso total de sementes avaliados em 100 plantas de Crotalaria spectabilis
Os efeitos diretos de AP (−0,010), NRP (−0,012), MSPA (0,004), NVG (0,024), PMVG (0,006) e MCG (0,116) sobre o PTS foram de pequena magnitude e sentidos distintos (Tabela 4) em relação aos observados nas correlações (Tabela 3). Além disso, as elevadas correlações lineares positivas do PTS com MSPA (r = 0,871*) e NVG (r = 0,956*) e as correlações de menor magnitude com AP (r = 0,380*) e NRP (r = 0,214*) se devem a efeitos indiretos dessas variáveis sobre o NS. Dessa forma, esse caractere foi o único que apresentou elevado efeito direto positivo sobre o PTS (0,981), com magnitude e direção similares às observadas para a correlação entre NS e PTS (r = 0,991*), indicando relação de causa e efeito. A elevada predição da análise de trilha (R2 = 0,995) e o reduzido efeito residual (0,073) indicam elevada confiabilidade nos resultados obtidos. Em crambe, Cargnelutti Filho et al. (2010)Cargnelutti Filho, A., Toebe, M., Silveira, T. R., Casarotto, G., Haesbaert, F. M. e Lopes, S. J. (2010). Tamanho de amostra e relações lineares de caracteres morfológicos e produtivos de crambe. Ciência Rural, 40, 2262-2267. http://dx.doi.org/10.1590/S0103-84782010001100003.
http://dx.doi.org/10.1590/S0103-84782010...
verificaram que, para a seleção indireta de plantas com maior produção de grãos, deverá ser avaliado o número de frutos por planta.
Efeitos diretos e indiretos dos caracteres: altura da planta, número de ramos produtivos, massa seca da parte aérea, número de vagens, peso médio de cada vagem, número de sementes e massa de cem grãos sobre o peso total de sementes com base na avaliação de 100 plantas de Crotalaria spectabilis.
Pode-se inferir com base no presente estudo que, para a seleção indireta de plantas com maior PTS, deve-se mensurar o número de sementes (Tabelas 3,4). Caracteres como AP, NRP, MSPA e NVG atuam indiretamente sobre o PTS, via NS. Outros caracteres como PMVG e MCG não apresentam associação linear e efeitos diretos sobre o PTS e, portanto, não devem ser considerados na seleção indireta visando à obtenção de plantas com maior produção de sementes.
CONCLUSÃO
Em Crotalaria spectabilis, são necessárias 121 plantas para a estimação da média dos caracteres NRP, MSPA, MSVG, NVG, NS e PTS com erro de estimação máximo de 10% da média. Já para a estimação de AP, PMVG, NSVG e MCG nesse mesmo nível de precisão, são necessárias 6 plantas.
A variável NS apresentou elevada correlação e efeito direto sobre PTS, indicando relação de causa e efeito e possibilidade de seleção indireta. As demais variáveis explicativas apresentaram efeitos indiretos sobre PTS via NS.
AGRADECIMENTOS
Ao Programa de Educação Tutorial do Ministério da Educação, ao Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), à Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS) e à Fundação Universidade Federal do Pampa (UNIPAMPA) pela concessão de bolsas. À empresa Piraí Sementes pela concessão das sementes de C. spectabilis para fins de pesquisa.
REFERÊNCIAS
- Albuquerque, M. C. F., Rodrigues, T. J. D. e Mendonça, E. A. F. (2000). Absorção de água por sementes de Crotalaria spectabilis Roth determinada em diferentes temperaturas e disponibilidade hídrica. Revista Brasileira de Sementes, 22, 206-215.
- Araújo, L. S., Cunha, P. C. R., Silveira, P. M., Sousa Netto, M. e Oliveira, F. C. (2015). Potencial de cobertura do solo e supressão de tiririca (Cyperus rotundus) por resíduos culturais de plantas de cobertura. Revista Ceres, 62, 483-488. http://dx.doi.org/10.1590/0034-737X201562050009
» http://dx.doi.org/10.1590/0034-737X201562050009 - Burin, C., Cargnelutti Filho, A., Toebe, M., Alves, B. M. e Fick, A. L. (2014). Dimensionamento amostral para a estimação da média de caracteres de tremoço branco. Comunicata Scientiae, 5, 205-212.
- Bussab, W. O., Morettin, P. A. (2011). Estatística básica. 7. ed. São Paulo: Saraiva.
- Cargnelutti Filho, A., Facco, G., Lúcio, A. D., Toebe, M., Burin, M., Fick, A. L. e Neu, I. M. M. (2014). Tamanho de amostra para a estimação da média de caracteres morfológicos e produtivos de nabo forrageiro. Ciência Rural, 44, 223-227. http://dx.doi.org/10.1590/S0103-84782014000200005
» http://dx.doi.org/10.1590/S0103-84782014000200005 - Cargnelutti Filho, A., Toebe, M., Alves, B. M., Burin, C., Santos, G. O., Facco, G. e Neu, I. M. M. (2015). Dimensionamento amostral para avaliar caracteres morfológicos e produtivos de aveia preta em épocas de avaliação. Ciência Rural, 45, 9-13. http://dx.doi.org/10.1590/0103-8478cr20140504
» http://dx.doi.org/10.1590/0103-8478cr20140504 - Cargnelutti Filho, A., Toebe, M., Silveira, T. R., Casarotto, G., Haesbaert, F. M. e Lopes, S. J. (2010). Tamanho de amostra e relações lineares de caracteres morfológicos e produtivos de crambe. Ciência Rural, 40, 2262-2267. http://dx.doi.org/10.1590/S0103-84782010001100003
» http://dx.doi.org/10.1590/S0103-84782010001100003 - Cruz, C. D. (2013). GENES — a software package for analysis in experimental statistics and quantitative genetics. Acta Scientiarum. Agronomy, 35, 271-276. http://dx.doi.org/10.4025/actasciagron.v35i3.21251
» http://dx.doi.org/10.4025/actasciagron.v35i3.21251 - Cruz, C. D. e Carneiro, P. C. S. (2006). Modelos biométricos aplicados ao melhoramento genético. v. 2. 2. ed. Viçosa: UFV.
- Cruz, C. D. e Regazzi, A. J. (1997). Modelos biométricos aplicados ao melhoramento genético. 2. ed. Viçosa: UFV.
- Empresa Brasileira de Pesquisa Agropecuária (2013). Sistema brasileiro de classificação de solos. Brasília: Embrapa.
- Ferreira, D. F. (2009). Estatística básica. 2. ed. Lavras: UFLA.
- Leal, M. A. A., Guerra, J. G. M., Peixoto, R. T. G. e Almeida, D. L. (2012). Desempenho de crotalária cultivada em diferentes épocas de semeadura e de corte. Revista Ceres, 59, 386-391. http://dx.doi.org/10.1590/S0034-737X2012000300014
» http://dx.doi.org/10.1590/S0034-737X2012000300014 - Massey, F. J. (1951). The Kolmogorov-Smirnov test for goodness of fit. Journal of the American Statistical Association, 46, 68-78.
- Montgomery, D. C. e Peck, E. A. (1982). Introduction to linear regression analysis. New York: John Wiley & Sons.
- Storck, L., Garcia, D. C., Lopes, S. J. e Estefanel, V. (2011). Experimentação vegetal. 3. ed. Santa Maria: UFSM.
- Teodoro, P. E., Silva Junior, C. A., Ribeiro, L. P., Silva, F. A., Corrêa, C. C. G., Zanuncio, A. S. e Torres, F. E. (2015). Sample dimension for estimation of biomass and yield of sunn (Crotalaria juncea L.) and showy rattlebox (C. spectabilis Roth.). Journal of Agronomy, 14, 98-101. http://dx.doi.org/10.3923/ja.2015.98.101
» http://dx.doi.org/10.3923/ja.2015.98.101 - Teodoro, R. B., Oliveira, F. L., Silva, D. M. N., Fávero, C. e Quaresma, M. A. L. (2011). Aspectos agronômicos de leguminosas para adubação verde no Cerrado do Alto Vale do Jequitinhonha. Revista Brasileira de Ciência do Solo, 35, 635-640. http://dx.doi.org/10.1590/S0100-06832011000200032
» http://dx.doi.org/10.1590/S0100-06832011000200032 - Wright, S. (1921). Correlation and causation. Journal of Agricultural Research, 20, 557-585.
- Wright, S. (1934). The method of path coefficients. The Annals of Mathematical Statistics, 5, 161-215.
- Wright, S. (1923). The theory of path coefficients: a reply to Niles’s criticism. Genetics, 8, 239-255.
Datas de Publicação
-
Publicação nesta coleção
12 Jan 2017 -
Data do Fascículo
Jan-Mar 2017
Histórico
-
Recebido
20 Dez 2015 -
Aceito
07 Jun 2016