Equações de predição de valores energéticos de alimentos obtidas utilizando meta-análise e componentes principais

Prediction equations of energetic values of feedstuffs obtained using meta-analysis and principal components

Equações de predição de valores energéticos de alimentos obtidas utilizando meta-análise e componentes principais

Prediction equations of energetic values of feedstuffs obtained using meta-analysis and principal components

Flávia Cristina Martins Queiroz Mariano^I,1; Renato Ribeiro de Lima^I; Paulo Borges Rodrigues^II; Renata Ribeiro Alvarenga^II; Germano Augusto Jerônimo do Nascimento^III

^IDepartamento de Ciências Exatas, Universidade Federal de Lavras (UFLA), CP 3037, 37200-000, Lavras, MG, Brasil. E-mail: [e-mail]

^IIDepartamento de Zootecnia, UFLA, Lavras, MG, Brasil

^IIIDepartamento de Zootecnia, Universidade Federal do Ceará (UFC), Fortaleza, CE, Brasil

RESUMO

Neste estudo, foi proposta a utilização da análise de componentes principais, na formação de grupos homogêneos de artigos científicos, a serem considerados em uma meta-análise. Nessa meta-análise, foram utilizados resultados referentes à composição química e energia metabolizável aparente corrigida pelo balanço de nitrogênio (EMAn) de alimentos para aves, na obtenção de equações de predição da EMAn. Foram considerados 293 experimentos com resultados já publicados. Dados provenientes de ensaios metabólicos foram utilizados na validação das equações de predição obtidas, sendo que estas apresentaram resultados semelhantes às disponíveis na literatura. No procedimento de meta-análise, a formação de grupos homogêneos de resultados experimentais, que é uma das maiores dificuldades, foi facilitada com a utilização de componentes principais, uma vez que não houve a necessidade de determinar variáveis ou fatores a serem considerados nessa classificação. Assim, tem-se uma forma rápida e eficiente de definir tais grupos.

Palavras-chave: composição química, energia metabolizável, validação de equações.

ABSTRACT

The purpose of this study was to develop a meta-analysis study by using the principal components analysis to obtain homogeneous groups of experimental results. In the process of the meta-analysis, it was considered data from 293 experiments carried out in Brazil. Prediction equations were obtained to estimate the nitrogen-corrected apparent metabolizable energy (AMEn) of poultry feedstuffs. Data from metabolic trials were used to validate the prediction equations obtained, which were similar than other available equations in the literature. One of the problems in meta-analysis is the determination of the homogeneous groups of experiments and, this problem was eliminated by using principal components, since there was no need to establish variables or factors to be considered in this classification.

Key words: chemical composition, metabolizable energy, validation of equations.

INTRODUÇÃO

A diversidade de alimentos e subprodutos utilizados na alimentação de frangos de corte exige o conhecimento preciso da composição química e valores de energia metabolizável (EM) desses alimentos, possibilitando a formulação de dietas equilibradas nutricional e economicamente. Uma forma rápida de determinar os valores de EM consiste na utilização de equações de predição, estabelecidas em função da composição química deles, normalmente de fácil e rápida obtenção.

Existem, à disposição na literatura, vários trabalhos de pesquisa, em que se estabeleceram equações de predição para os valores energéticos de uma série de alimentos (ZHAO et al., 2008; WAN et al., 2009). Porém, os resultados obtidos utilizando-se essas equações não têm sido satisfatórios. Com o intuito de melhorá-los, pode-se utilizar o princípio da meta-análise. Um dos objetivos da meta-análise é a obtenção de novos resultados, a partir da síntese de n observações distintas, porém relacionadas entre si (LOVATTO et al., 2007). Assim, avalia-se um número maior de trabalhos, o que garante a obtenção de resultados mais consistentes.

Segundo LOVATTO et al. (2007), algumas etapas necessárias no desenvolvimento de uma meta-análise são: a) definição do objetivo do trabalho, do tema da pesquisa; b) sistematização das informações; c) codificação dos dados; d) filtragem dos dados; e) análise dos dados incluídos na base; f) formação de grupos; g) ponderações; e h) escolha de um modelo estatístico. Embora os estudos considerados envolvam um mesmo assunto, na meta-análise, é importante que grupos homogêneos destes sejam formados e a variabilidade existente entre eles seja modelada e incorporada na análise. No entanto, a formação de grupos homogêneos é uma das maiores dificuldades na meta-análise, podendo ser feita por combinação de níveis de fatores, análises de agrupamento, etc.

NASCIMENTO et al. (2009, 2011a, 2011b) estabeleceram equações de predição de EMAn utilizando o princípio da meta-análise, definindo 24 grupos homogêneos pela combinação de três fatores, que influenciavam na variabilidade da EMAn: sexo (macho, fêmea e mistos); idade (1^a e 2^a semanas de vida; 3^a e 4^a semanas; 5^a e 6^a semanas; acima ou indefinido); e metodologia empregada (coleta total - CT; alimentação forçada + CT).

Neste trabalho, foi proposta a utilização da análise de componentes principais como uma forma mais simples e direta de agrupamento dos resultados experimentais mais homogêneos, a serem considerados na meta-análise. Essa meta-análise foi desenvolvida para a obtenção de equações de predição para energia metabolizável aparente corrigida pelo balanço de nitrogênio (EMAn) em função da composição química de alimentos para aves. Tais equações obtidas foram validadas, utilizando-se dados provenientes de ensaios metabólicos, e, posteriormente, foram comparadas com outras já publicadas na literatura.

MATERIAL E MÉTODOS

Os dados utilizados na obtenção das equações de predição, via meta-análise, foram catalogados por NASCIMENTO et al. (2009, 2011a, 2011b). Esses dados referem-se a valores de energia metabolizável aparente corrigida pelo balanço de nitrogênio (EMAn) e composição química (proteína bruta-PB, extrato etéreo-EE, matéria mineral-MM, fibra bruta-FB, fibra em detergente neutro-FDN e fibra em detergente ácido-FDA) de alimentos concentrados energéticos e proteicos de origem vegetal, usualmente utilizados na formulação das rações avícolas. Neste estudo, foram utilizados diferentes resultados experimentais, realizados no Brasil, classificados em três categorias de alimentos: 197 trabalhos envolvendo equações de predição da EMAn de alimentos energéticos, 67 de soja+subprodutos (proteicos) e 293 de concentrados em geral (energéticos + proteicos).

Na formação de grupos homogêneos de resultados experimentais, foi proposta a utilização da técnica multivariada de componentes principais como método de agrupamento. Com o uso dessa técnica, o número de grupos a serem formados não é explicitamente definido como no caso em que os grupos são formados por combinação de níveis de fatores. Assim, para identificar qual é o número mais adequado de grupos homogêneos foram considerados 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 25, 30, 40 e 50 grupos. Além disso, foi considerada a situação na qual não se utiliza a meta-análise, ou seja, sem a formação de grupos. Portanto, foram realizadas 45 meta-análises, isto é, 15 critérios para a formação de grupos para cada uma das três categorias de alimentos.

Os passos desde a formação de grupos até a obtenção das equações de predição, em cada meta-análise, foram: (i) efetuar a análise de componentes principais, considerando as variáveis EMAn, PB, EE, MM, FB, FDN e FDA, obtendo-se o valor do primeiro componente principal (PC1) para cada um dos resultados experimentais; (ii) ordenar os valores referentes aos PC1; (iii) definir o número de grupos; (iv) estabelecer os quantis que delimitam os grupos homogêneos. Por exemplo: definindo-se 5 grupos, tem-se os quantis 0,20; 0,40; 0,60 e 0,80. Abaixo do quantil 0,20, estão os 20% dos resultados experimentais com os menores escores de PC1, que compõem o primeiro grupo. Entre o quantil 0,20 e 0,40, estão os próximos 20% dos resultados experimentais com menores valores de PC1, compondo o segundo grupo e assim por diante; (v) identificar os resultados experimentais dentro de cada grupo estabelecido; (vi) calcular a variância da EMAn observada nos resultados experimentais de cada grupo; (vii) obter o fator de ponderação que, neste trabalho, foi definido como o inverso da variância dentro de cada grupo; (viii) obter a equação de regressão estimada, utilizando o método dos mínimos quadrados ponderados, sendo os fatores de ponderação definidos em (vii).

Para o ajuste das equações, foi proposto o modelo de regressão linear múltipla , em que: Y_i refere-se ao valor da EMAn do alimento, no i-ésimo estudo; X_i1, ..., X_i6 representam as variáveis de composição química do alimento, no i-ésimo estudo, sendo respectivamente PB, MM, EE, FB, FDN e FDA e ε_i é o erro associado na i-ésima observação, assumido normal e independentemente distribuído, com média 0 e variância δ_i². O método stepwise foi adotado, sendo o critério de Akaike (AIC) utilizado na inclusão e exclusão de variáveis no modelo. Após a aplicação do stepwise, se alguma variável foi não significativa a 10%, pelo teste t, esta foi retirada do modelo.

Para efeito de seleção das equações obtidas neste trabalho que apresentaram melhores ajustes, utilizou-se, simultaneamente, o valor de AIC com os valores dos critérios de validação de cada meta-análise efetuada.

As equações que apresentaram melhores ajustes neste estudo, bem como equações de predição da EMAn de alimentos para aves recomendadas na literatura (RODRIGUES et al., 2001; 2002, NAGATA et al., 2004; ZONTA et al., 2004; NASCIMENTO et al., 2009, 2011a, 2011b) foram submetidas à validação. Nesta validação, foram utilizados dados provenientes de dois ensaios metabólicos, conduzidos por ALVARENGA et al. (2011). Nesses ensaios, foram utilizados nove ingredientes proteicos e nove energéticos, fornecidos a pintos machos da linhagem Cobb 500^®, utilizando o método tradicional de coleta total de excretas.

Para efeito de comparação, foram confrontados os resultados das equações selecionadas neste estudo e na literatura (Tabelas 1 e 2), considerando diferentes critérios de validação, descritos por TEDESCHI (2006). Um critério consistiu no ajuste do modelo de regressão linear simples Y_i = β₀ + β₁X_i + ε_i, em que: Y_i é a EMAn do alimento i obtida no ensaio metabólico; X_i é o valor predito da EMAn, utilizando-se uma equação selecionada; e ε_i é o erro tal que ε_i~ N(0,δ²). As hipóteses H₀:β₀ = 0 e H₀:β₁= 1 foram verificadas utilizando-se o teste t de Student, o que corresponde à hipótese de igualdade entre valores preditos e observados. Outro critério foi o coeficiente de correlação de concordância (CCC) entre os valores observados (Y_i) e os valores preditos de EMAn (X_i), proposto por LIN (1989). O CCC verifica, simultaneamente, se o modelo apresenta resultado mais acurado e preciso.

Toda rotina computacional utilizada neste trabalho foi desenvolvida no software livre R (R Development Core Team, 2011).

RESULTADOS E DISCUSSÃO

Equações de predição para EMAn selecionadas, obtidas neste estudo

Para a categoria dos alimentos energéticos, o primeiro componente principal (PC1) de cada trabalho explicou 69% da variação total. Em ambas as categorias dos alimentos proteicos (soja+subprodutos) e na dos concentrados em geral, 52% da variação total foi explicada pelo uso dos PC1. Quanto maior a porcentagem de explicação da variabilidade pelo PC1, melhor será o processo de formação de grupos de resultados experimentais homogêneos.

Para os alimentos energéticos, os valores dos critérios utilizados na seleção e validação das equações de predição estão representados na figura 1. À medida que foi utilizado um maior número de grupos, os valores de AIC aumentaram rapidamente, até próximo de cinco grupos (Figura 1a). A partir daí, observa-se um aumento menos expressivo com uma tendência de se estabilizar. Optou-se por utilizar o AIC, pois, apesar do critério do coeficiente de determinação (R²) (Figura 1b) ser largamente utilizado para a escolha do melhor modelo, ele é muito influenciado pelo número de parâmetros do modelo estatístico considerado. Verifica-se que, na figura 1c, com o aumento do número de grupos, principalmente a partir de 20 grupos, os valores de CCC diminuem. Isso não é desejado, pois, se a partir de uma equação obtêm-se valores preditos mais próximos dos valores observados de EMAn, o CCC tende a um. Na figura 1d, têm-se os p-valores do teste t, correspondentes aos resultados da validação, utilizando o ajuste da equação linear simples (H₀:β₀ = 0 e H₀: β₁= 1). Nota-se que, nas situações em que foram considerados 1, 3, 6, 8 e 9 grupos, essas hipóteses foram rejeitadas. Assim, analisando-se concomitantemente o AIC, o CCC e os p-valores, conclui-se que as equações de predição da EMAn que tiveram os melhores ajustes para os concentrados energéticos foram aquelas obtidas com a definição de 10 e 20 grupos, pois, com esses números de grupos, para o conjunto de dados utilizados, têm-se melhores resultados em relação aos critérios utilizados.

Nas outras duas situações, nas quais consideram os 293 resultados experimentais de alimentos concentrados e os 67 de soja+subprodutos, foram realizados procedimentos semelhantes aos dos alimentos energéticos. Coincidentemente, as equações de predição com 10 e 20 grupos no processo de estimação, para os dados utilizados, foram consideradas as melhores. Assim, as equações que apresentaram melhores ajustes, segundo o critério AIC e os resultados da validação, estão apresentadas nas quatro primeiras linhas da tabela 1, para alimentos energéticos, e, nas quatro primeiras linhas da tabela 2, para soja + subprodutos.

Comparação de equações de predição via validação com ensaio metabólico

Alimentos energéticos e concentrados

As equações de predição para EMAn de alimentos energéticos, utilizadas na comparação, foram aquelas descritas na tabela 1. Observa-se que essas equações são obtidas considerando-se diferentes bases de dados (apenas alimentos energéticos ou concentrados em geral), mas que foram utilizadas na validação da EMAn de alimentos energéticos.

Na tabela 1, observa-se que, com exceção das equações 5 e 6, as equações consideradas para alimentos energéticos não apresentam diferenças significativas entre os valores ajustados e os observados no ensaio metabólico. Ao analisar o coeficiente de correlação de concordância (CCC), nota-se que, até a equação 9, utilizada nessa comparação, todas possuem valores semelhantes de CCC, sendo todos próximos de um. Também, verifica-se que, apesar de as equações 1, 2, 3, 8 e 9 não apresentarem, em sua maioria, coeficientes de determinação (R²) mais altos, estas possuem os maiores p-valores para os parâmetros da regressão linear simples. Então, com a análise simultânea desses critérios, conclui-se que as equações 1, 2, 3, 4, 7, 8 e 9 descreveram com acurácia e precisão os valores da EMAn obtidos no ensaio metabólico para alimentos energéticos.

Portanto, as equações obtidas neste trabalho apresentaram resultados semelhantes às disponíveis na literatura (RODRIGUES et al., 2001; NASCIMENTO et al., 2009, 2011a). Em NASCIMENTO et al. (2009, 2011a), assim como neste estudo, teve-se um baixo custo na obtenção de tais equações, pois foi utilizada meta-análise ao invés de executarem experimentos específicos como em RODRIGUES et al. (2001). Porém, o desenvolvimento da meta-análise foi facilitado neste estudo, uma vez que não houve a necessidade de definir fatores mais adequados para a formação de grupos homogêneos de resultados experimentais, como em NASCIMENTO et al. (2009, 2011a).

Soja+subprodutos e concentrados

As equações de predição para EMAn de soja+subprodutos, utilizadas na comparação, estão apresentadas na tabela 2. A comparação foi desenvolvida entre as diferentes equações de predição obtidas e aquelas adquiridas em trabalhos de outros autores, análogo ao feito para os alimentos energéticos e concentrados.

Pelos valores dos critérios apresentados na tabela 2, observa-se que as únicas equações de predição para EMAn para soja+subprodutos, em que as hipóteses estabelecidas foram rejeitadas durante a validação, foram as equações 2 e 9. As demais não apresentaram diferenças significativas entre os valores ajustados e os observados no ensaio metabólico. Verifica-se, também, que, por possuírem valores próximos de CCC tendendo a um, as equações 2, 3, 4, 6 e 7 apresentam os melhores ajustes dentre as equações comparadas. No entanto, as equações 3, 4, 6 e 7 não apresentaram os maiores coeficientes de determinação (R²), porém possuem maiores p-valores. Contudo, se fosse considerado só o coeficiente de determinação como critério único de seleção, o melhor ajuste considerado seria o da equação 8, a qual possui valor de CCC igual a 0,21. No decorrer de todo este trabalho, verifica-se que o R² não é indicado a ser usado como critério único de seleção, pois o seu resultado não refletiria o que foi concluído na análise concomitante dos demais critérios. Portanto, diante desses resultados, recomenda-se, sempre que possível, a utilização conjunta dos critérios de seleção e validação de modelos. Assim, fazendo a análise simultânea desses critérios, conclui-se que as equações 3, 4, 6 e 7 descreveram com maior acurácia e precisão os valores da EMAn obtidos no ensaio metabólico para soja+subprodutos.

Portanto, as equações obtidas neste trabalho apresentaram resultados semelhantes às disponíveis em NASCIMENTO et al. (2009, 2011b). Análogo ao discutido na comparação para alimentos energéticos, a formação de grupos homogêneos foi facilitada através da técnica multivariada de componentes principais, uma vez que não é muito simples determinar quais os fatores que causam variabilidade entre os experimentos e que devam ser considerados na identificação de resultados experimentais homogêneos.

As equações ajustadas e selecionadas neste trabalho, cujos valores de R² foram altos, estão sendo testadas por meio de ensaios metabólicos realizados com frangos de corte em fase de crescimento, conforme trabalho realizado por ALVARENGA et al. (2011).

CONCLUSÃO

As equações de predição obtidas neste trabalho apresentaram resultados semelhantes a outras propostas na literatura. Porém, a formação de grupos homogêneos de trabalhos, que é uma das maiores dificuldades na meta-análise, foi facilitada com a utilização da técnica de componentes principais. Assim, a utilização dessa técnica é viável por ser uma forma mais simples e direta de agrupamento dos trabalhos mais homogêneos, sem a necessidade da determinação de fatores coerentes para esse agrupamento como proposto em outros trabalhos.

AGRADECIMENTO

Os autores agradecem o apoio financeiro do Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) e da Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

Recebido para publicação 06.10.11

Aprovado em 07.06.12

Devolvido pelo autor 04.07.12

CR-6152

• ALVARENGA, R.R. et al. Energetic values of feedstuffs for broilers determined with in vivo assays and prediction equations. Animal Feed Science Technology, Amsterdam, v.168, Issues 3-4, p.257-266, set. 2011. Disponível em: <http://www.sciencedirect.com/science/article/pii/S0377840111002276>. Acesso em: 22 set. 2011. doi:10.1016/j.anifeedsci.2011.04.092.
• LIN, L. A concordance correlation coefficient to evaluate reprodutibility. Biometrics, Washington, v.45, n.1, p.255-268, mar. 1989. Disponível em: <http://www.jstor.org/stable/2532051>. Acesso em: 13 set. 2009.
• LOVATTO, P.A. et al. Meta-análise em pesquisas científicas: enfoque em metodologias. Revista Brasileira de Zootecnia, Viçosa, v.36, supl., p.285-294, jul. 2007. Disponível em: <http://www.scielo.br/scielo.php?pid=S1516-35982007001000026&script=sci_arttext>. Acesso em: 11 out. 2008. doi: 10.1590/S1516-35982007001000025.
• NAGATA, A.K. et al. Energia metabolizável de alguns alimentos energéticos para frangos de corte, determinada por ensaios metabólicos e por equações de predição. Ciência e Agrotecnologia, Lavras, v.28, n.3, p.668-677, maio/jun. 2004. Disponível em: <http://www.scielo.br/scielo.php?pid=S1413-70542004000300025&script=sci_arttext>. Acesso em: 20 jun. 2009. doi: 10.1590/S1413-70542004000300024.
• NASCIMENTO, G.A.J. et al. Equações de predição para estimar os valores energéticos de alimentos concentrados de origem vegetal para aves utilizando a metanálise. Revista Brasileira de Zootecnia, Viçosa, v.38, n.7, p.1265-1271, dez. 2009. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1516-35982009000700015>. Acesso em: 10 jan. 2010. doi: 10.1590/S1516-35982009000700014.
• NASCIMENTO, G.A.J. et al. Equações de predição para estimar valores da energia metabolizável de alimentos concentrados energéticos para aves utilizando meta-análise. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, Belo Horizonte, v.63, n.1, p.222-230, fev. 2011a. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102-09352011000100032>. Acesso em: 10 set. 2011. doi: 10.1590/S0102-09352011000100032.
• NASCIMENTO, G.A.J. et al. Equações de predição para estimar os valores da EMAn de alimentos protéicos para aves utilizando a meta-análise. Revista Brasileira de Zootecnia, Viçosa, v.40, n.10, p.2172-2177, out. 2011b. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1516-35982011001000016>. Acesso em: 15 dez. 2011. doi: 10.1590/S1516-35982011001000016.
• R DEVELOPMENT CORE TEAM. R: a language and enviroment for statistical computing. Viena: R Foundation for Statistical Computing, 2011. Software.
• RODRIGUES, P.B. et al. Valores energéticos do milheto, do milho e subprodutos do milho, determinados com frangos de corte e galos adultos. Revista Brasileira de Zootecnia, Viçosa, v.30, n.6, p.1767-1778, nov./dez. 2001. Disponível em: <http://www.scielo.br/scielo.php?pid=S1516-35982001000700015&script=sci_arttext>. Acesso em: 23 jun. 2009. doi: 10.1590/S1516-35982001000700015.
• RODRIGUES, P.B. et al. Valores energéticos da soja e subprodutos da soja, determinados com frangos de corte e galos adultos. Revista Brasileira de Zootecnia, Viçosa, v.31, n.4, p.1771-1782, jul./ago. 2002. Disponível em: <http://www.scielo.br/scielo.php?pid=S1516-35982002000700020&script=sci_arttext>. Acesso em: 23 jun. 2009. doi: 10.1590/S1516-35982002000700020.
• TEDESCHI. L.O. Assessment of the adequacy of mathematical models. Agricultural systems, Amsterdam, v.89, n.2/3, p.225-247, set. 2006. Disponível em: <http://www.sciencedirect.com/science/article/pii/S0308521X05002568>. Acesso em: 11 out. 2008. doi:10.1016/j.agsy.2005.11.004.
• WAN, H.F. et al. Prediction of true metabolizable energy from chemical composition of wheat milling by-products for ducks. Poultry Science, Champaign, v.88, p.92-97, jan. 2009. Disponível em: <http://ps.fass.org/content/88/1/92.full.pdf+html>. Acesso em: 20 dez. 2010. doi: 10.3382/ps.2008-00160.
• ZHAO, F. et al. Predicting metabolizable energy of normal corn from its chemical composition in adult pekin ducks. Poultry Science, Champaign, v.87, p.1603-1608, ago. 2008. Disponível em: <http://ps.fass.org/content/87/8/1603.full.pdf+html>. Acesso em: 20 dez. 2010. doi: 10.3382/ps.2007-00494.
• ZONTA, M.C.M. et al. Energia metabolizável de ingredientes protéicos determinada pelo método de coleta total e por equações de predição. Ciência e Agrotecnologia, Lavras, v.28, n.6, p.1400-1407, nov./dez. 2004. Disponível em: <http://www.scielo.br/scielo.php?pid=S1413-70542004000600024&script=sci_arttext>. Acesso em: 20 jun. 2009. doi: 10.1590/S1413-70542004000600023.

1

Autor para correspondência.

Datas de Publicação

• Publicação nesta coleção
  31 Jul 2012
• Data do Fascículo
  Set 2012

Histórico

• Recebido
  06 Out 2011
• Aceito
  07 Jun 2012