Descrição do crescimento de ovinos Santa Inês utilizando modelos não-lineares selecionados por análise multivariada

Description of Santa Ines sheep growth using non-linear models selected by multivariate analysis

Milton Rezende TEIXEIRA NETO Jurandir Ferreira da CRUZ Helder Henrique Neves FARIA Emanuely Santos SOUZA Paulo Luiz Souza CARNEIRO Carlos Henrique Mendes MALHADO Sobre os autores

RESUMO

O objetivo deste estudo foi descrever o crescimento de ovinos Santa Inês utilizando modelos não-lineares selecionados por análise multivariada. Dezessete modelos não-lineares foram utilizados para descrever o crescimento de cordeiros do nascimento aos 210 dias de idade. Os modelos foram classificados por meio de análise multivariada utilizando o coeficiente de determinação ajustado, critério de informação de Akaike, critério de informação Bayesiano, erro quadrático médio de predição, coeficiente de determinação de predição, desvio médio absoluto dos resíduos e percentual de convergência como avaliadores da qualidade do ajuste. Os modelos Richards, Reynolds e Weibull não convergiram e os modelos Schnute, Meloun II, IV e VI, Gamito, Michaelis Menten e Michaelis Menten Modificado não apresentaram padrão de crescimento tradicional. Os modelos Mitscherlich, Gompertz, Logístico, Meloun I, III e V, Brody e Von Bertalanffy convergiram e foram utilizados na análise de agrupamento. Dentre os quatroclusters formados, o primeiro, composto pelos modelos Mitscherlich, Meloun I e Brody, apresentou os melhores valores para os critério de informação de Akaike e bayesiano. Os parâmetros da curva foram influenciados pelo sexo e o tipo de parto (P<0,05) e o ponto inflexão ocorreu aos 94 dias de idade. Os modelos Mitscherlich, Meloun I e Brody são adequados e equivalentes para descrição do crescimento de ovinos Santa Inês.

análise de agrupamento; cordeiro; curva de crescimento; ponto de inflexão; taxa de crescimento

SUMMARY

The objective of this study was to describe the growth of Santa Ines sheep using nonlinear models selected by multivariate analysis. Seventeen nonlinear models were used to describe the growth of sheep from birth to 210 days of age. The models were classified by means of multivariate analysis using the adjusted coefficient of determination, the Akaike information criterion, the Bayesian information criterion, the mean squared error of prediction, the coefficient of determination, mean absolute deviation and percentage of convergence of prediction as evaluators of adjustment quality. The Richards, Reynolds and Weibull models did not converge and the Schnute, Meloun II, IV and VI, Gamito, Michaelis Menten and Modified Michaelis Menten models did not show any traditional growth pattern. Models Mitscherlich, Gompertz, Logistic, Meloun I, III and V, Brody and Von Bertalanffy converged and were used in grouping analysis. Of the four formed clusters formed, the first, composed by models Mitscherlich, Meloun I and Brody, showed the best values for all the evaluators considered. The curve parameters were influenced by gender and type of birth (P<0.05) and the inflection point occurred at 94 days of age. The Mitscherlich, Meloun I and Brody models are adequate and equivalent to describe the growth of Santa Ines sheep.

grouping analysis; growth curve; growth rate; lamb; point of inflection

INTRODUÇÃO

Diferentes modelos não-lineares têm sido utilizados para descrever o crescimento de ovinos, sendo os modelos de Richards, Gompertz, Brody, Logístico e Von Bertalanffy, os mais comuns (LÔBO et al., 2006LÔBO, R.N.B.; VILLELA, L.C.V.; LÔBO, A.M.B.O.; PASSOS, J.R.S.; OLIVEIRA, A.A. Parâmetros genéticos de características estimadas da curva de crescimento de ovinos da raça Santa Inês. Revista Brasileira de Zootecnia , v.35, p.1012-1019, 2006. Supl.; MALHADO et al., 2009MALHADO, C.H.M.; CARNEIRO, P.L.S.; AFFONSO, P.R.A.M.; SOUZA JÚNIOR, A.A.O.; SARMENTO, J.L.R. Growth curves in Dorper sheep crossed with the local Brazilian breeds, Morada Nova, Rabo Largo, and Santa Inês. Small Ruminant Research , v.84, p.16-21, 2009.; CHALH & EL GAZZAH, 2014CHALH, A.; EL GAZZAH, M. Variogram investigation of covariance shape within longitudinal data with possible use of a krigeage technique as an interpolation tool: Sheep growth data as an example. Irish Journal of Agricultural and Food Research , v.53, p.51–64, 2014.). Alternativamente, os modelos Schnute, Mitscherlich, Gamito e Meloun I, II, III e IV têm sido propostos (SILVEIRA et al., 2011)SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M; MUNIZ, J.A. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural , v.41, n.4, p.692-698, 2011..

Os modelos Schnute, Mitscherlich e Meloun I e II, utilizados anteriormente para descrever o crescimento da estrutura auricular de fetos (MELOUN & MILITKY, 1996MELOUN, M.; MILITKY J. Sbírka úloh, Statistické zpracování experimentálnich dat . Pardubice: Universita Pardubice. 1996.), devido à grande flexibilidade podem ser adequados para descrição de curvas de crescimento de animais.

Na maioria dos estudos sobre curvas de crescimento animal o coeficiente de determinação ajustado tem sido utilizado como critério de seleção do modelo mais adequado. No entanto o critério de informação de Akaike, critério de informação Bayesiano, erro quadrático médio de predição e o coeficiente de determinação de predição são avaliadores importantes da qualidade do ajuste do modelo que podem ser considerados (GÓMEZ et al., 2008GÓMEZ, D.A.A.; MUÑOZ, M.F.C.; BETANCUR, L.F.R. Modelación de las funciones de crecimiento aplicadas a la producción animal. Revista Colombiana de Ciencias Pecuárias , v.21, p.39-58, 2008.).

Quanto maior o número de critérios considerados, mais segura é a indicação dos melhores modelos. Por outro lado, quando o número de avaliadores é grande, a escolha dos modelos pode se constituir em um processo complexo (SILVEIRA et al., 2011SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M; MUNIZ, J.A. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural , v.41, n.4, p.692-698, 2011.). A utilização de técnicas de agrupamento com base nas médias de cada avaliador permite a identificação dos melhores modelos de forma mais fácil (SILVEIRA et al., 2012SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M. Classificação multivariada de modelos de crescimento para grupos genéticos de ovinos de corte. Revista Brasileira de Saúde e Produção Animal [online], v.13, n1, p.62-73, 2012.).

Este estudo teve como objetivo identificar, por meio de análise multivariada e avaliadores de qualidade, o modelo não-linear mais adequado para descrever o crescimento de ovinos da raça Santa Inês.

MATERIAL E MÉTODOS

Os dados foram coletados no período de setembro de 2011 a maio de 2012, na Fazenda União, localizada no município de Ribeirão do Largo-BA (15º27'33'' S e 40º44'28''). Foram utilizados 89 ovinos puros da raça Santa Inês, sendo 43 machos e 46 fêmeas. Os animais foram criados em sistema semi-intensivo e mantidos com as mães até a idade de 90 dias, momento em que ocorreu o desmame.

A alimentação das mães foi constituída de pastagem de grama estrela e as crias receberam suplementação com alimento concentrado (16% PB) em quantidade equivalente a 1,0% peso vivo/dia.

A pesagem dos animais foi realizada sistematicamente em intervalos de 30 dias, desde o nascimento até os 210 dias de idade. O controle parasitário foi realizado com base em exames parasitológicos de fezes (OPG) conforme programa profilático da propriedade.

Para determinação dos parâmetros da curva de crescimento foram utilizados os modelos não-lineares de Richards, Brody, Gompertz, Logístico, Von Bertalanffy, Reynolds, Weibull, Mitscherlich, Schnute, Gamito, Michaelis Menten, Michaelis Menten Modificado e Meloun I, II, III, IV e V (Tabela 1). Para todos os modelos utilizados, o parâmetro β1representou o peso adulto, ou peso assintótico, do animal; o parâmetro β3, a taxa de maturidade, ou velocidade de crescimento; o parâmetro β4, o ponto de inflexão; e o parâmetro β2, uma constante de integração.

Tabela 1
Modelos de regressão não-linear utilizados para descrição da curva de crescimento dos cordeiros

Todos os modelos foram ajustados pelo método da Máxima Verossimilhança via algoritmo de otimização de Gauss-Newton utilizando o PROC MODEL ® versão 9.1 (SAS, 2003SAS. Statistical Analysis System user’s guide . Version 9.1 ed. Cary, USA: SAS Institute, 2003.). Os valores dos parâmetros foram definidos após identificação daqueles que apresentavam menor número de interações para convergência.

Para avaliar a qualidade do ajuste dos modelos foram utilizados como avaliadores o coeficiente de determinação ajustado (R2aj), critério de informação de Akaike (AIC), critério de informação bayesiano (BIC), erro quadrático médio de predição (MEP), coeficiente de determinação de predição (R2p), percentual de convergência (C%) e desvio médio absoluto dos resíduos (DMA), conforme sugerido por Silveira et al. (2011)SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M; MUNIZ, J.A. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural , v.41, n.4, p.692-698, 2011. e Sarmento et al. (2006)SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006..

Após constatação da convergência, os valores dos avaliadores de cada modelo foram calculados a partir da média dos valores dos avaliadores obtidos para cada animal. Posteriormente, os dados foram submetidos à análise de agrupamento (PROC CLUSTER, SAS, versão 9.1).

Depois de selecionado o modelo, foi avaliada a influência do sexo e tipo de parto sobre os parâmetros da curva (GLM, SAS, versão 9.1), bem como foram estimadas as correlações de Pearson entre os parâmetros e a taxa de crescimento absoluto (TCA). Para os modelos que não apresentam ponto de inflexão utilizou-se o método da máxima curvatura geométrica para determinar a idade em que os animais desaceleram o ganho de peso, conforme sugerido por Cecon et al. (2008)CECON, P.R.; SILVA, F.F.; FERREIRA, A.; FERRÃO, R.G.; CARNEIRO, A.P.S.; DETMANN, E.; FARIA, P.N.; MORAIS, T.S.S. Análise de medidas repetidas na avaliação de clones de café ‘Conilon’. Pesquisa Agropecuária Brasileira , v.43, p.1171-1176, 2008..

RESULTADOS E DISCUSSÃO

Os modelos Richards, Reynolds e Weibull não convergiram e os modelos Schnute, Meloun II, IV e VI, Gamito, Michaelis Menten e Michaelis Menten Modificado não apresentaram padrão de crescimento tradicional e, por essa razão seus parâmetros não foram considerados. Os modelos Mitscherlich, Gompertz, Logístico, Meloun I, III e V, Brody e Von Bertalanffy mostraram-se adequados, sendo o Meloun V o único modelo com quatro parâmetros (Tabela 2).

Tabela 2
Valores de peso assintótico (β1), constante de integração (β2), taxa de maturidade (β3) e ponto de inflexão (β4), estimados por diferentes modelos de crescimento

Com exceção do modelo Meloun V, todos os demais modelos superestimaram o peso ao nascimento. O peso assintótico estimado (β1) variou de 37,56 a 66,54kg, sendo que os modelos Mitscherlich, Meloun I e III e Brody apresentaram estimativas semelhantes.

O peso assintótico representa à estimativa do peso à maturidade e não o peso máximo que o animal pode alcançar (BROW et al., 1976BROW, J.E.; FITZHUGH JÚNIOR, H.A.; CARTWRIGHT, T.C.A. A comparison of nonlinear models for describing weigt-age relationships in cattle.Journal of Animal Science , v.42, p.810-818, 1976.) e a variação verificada na sua estimativa nos diferentes modelos pode ser atribuída à dificuldade dos modelos não-lineares em ajustar o peso inicial, comprometendo o processo de estimação desse parâmetro, fato esse também verificado por Cruz et al. (2009)CRUZ, G.R.B.; COSTA, R.G.; RIBEIRO, M.N. Curva de crescimento de caprinos mestiços no estado da Paraíba. Revista Brasileira de Ciências Agrárias , v.4, n.2, p.204-210, 2009.. A estimativa do β1 pode muitas vezes ser sujeita a erros porque o conjunto de dados ajustado por modelos não-lineares de crescimento não inclui o crescimento do animal até a maturidade (TORAL, 2008TORAL, F.L.B. Número e intervalo de pesagens para estimação de parâmetros de curvas de crescimento em bovinos. Revista Brasileira de Zootecnia , v.37, n.12, p.2120–2128, 2008.).

Por outro lado, a estimativa do peso assintótico para uma determinada raça pode variar em função de outros fatores. Guedes et al. (2004)GUEDES, M.H.P.; MUNIZ, J.A.; PEREZ, J.R.O.; SILVA, F.F.; AQUINO, L.H.; SANTOS, C.L. Estudo das curvas de crescimento de cordeiros das raças Santa Inês e Bergamácia considerando heterogeneidade de variâncias. Ciência Agrotécnica , v.28, p.381-388, 2004. utilizando os modelos Brody, Von Bertalanffy, Gompertz e Logístico para descrever o crescimento de animais Santa Inês criados em sistema intensivo encontraram valores superiores aos verificados no presente estudo (65,28; 51,06; 43,18 e 37,24kg, respectivamente). Nesse caso, a diferença verificada entre os dois estudos pode ser devida ao manejo dos animais, uma vez que estimativa do peso assintótico pode sofrer influencia da raça e do sistema de criação (MALHADO et al., 2008MALHADO, C.H.M.; CARNEIRO, P.L.S.; SANTOS, P.F.; AZEVEDO, D.M.M.; SOUZA, J.C.; AFFONSO, P.R.M. Curva de crescimento em ovinos mestiços Santa Inês x Texel criados no Sudoeste do Estado da Bahia. Revista Brasileira de Saúde e Produção Animal [online], v.9, n.2, p.210-218, 2008.).

A taxa de maturidade (β3) apresentou grande amplitude, cujos valores variaram de 0,0063 a 70,3911. Considerando que a taxa de maturidade indica a velocidade de crescimento para atingir o peso assintótico, os valores elevados de β3 podem ser interpretados como indicadores de precocidade.

Todos os modelos apresentaram bons ajustes em relação aos coeficientes de determinação ajustado (R2aj) e de predição (R2p), com médias superiores a 0,97. O critério de informação de Akaike (AIC) apresentou menores médias para Brody e Meloun I, enquanto que, para critério de informação Bayesiano (BIC) o menor valor foi verificado para o modelo Brody. O erro quadrático médio de predição (MEP) apresentou menor e maior média para os modelos Meloun I e Meloun V, respectivamente. Com base no percentual de convergência (C%), os modelos Gompertz, Logístico e Von Bertalanffy foram superiores aos demais, convergindo a curva de todos os animais. No entanto, esses modelos apresentaram os maiores valores de desvio médio absoluto (Tabela 3).

Tabela 3
Valores do critério de informação de Akaike (AIC), critério de informação Bayesiano (BIC), coeficiente de determinação ajustado (R2aj), erro quadrático médio de predição (MEP), coeficiente de determinação de predição (R2p), percentual de convergência (C%) e desvio médio absoluto dos resíduos (DMA) de acordo com os modelos estudados

Os modelos com melhor ajuste apresentam maiores valores de coeficiente de determinação (R2aj e R2p) e percentual de convergência e menores AIC, BIC, MEP DMA (SILVA et al., 2011SILVA, F.L.; ALENCAR, M.M.; FREITAS, A.R.; PACKER, I.U.; MOURÃO, G.B. Curvas de crescimento em vacas de corte de diferentes tipos biológicos.Pesquisa Agropecuária Brasileira , v.46, n.3, p.262-271, 2011.; SARMENTO et al., 2006SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006.). Em estudo com ovinos F1 Santa Inês x Dorper Silveira et al., (2011)SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M; MUNIZ, J.A. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural , v.41, n.4, p.692-698, 2011. encontraram médias semelhantes as do presente estudo para R2aj (0,9914), R2p(0,9992), AIC (25,9349) e BIC (27,9102) e inferior para MEP (0,6708).

Todos os avaliadores de qualidade, com exceção do MEP e DMA, apresentaram baixo coeficiente de variação (CV). O alto CV obtido para o MEP deveu-se ao fato de este avaliador ser calculado após a remoção sucessiva de observações do conjunto de dados originais, avaliando a capacidade do modelo de se ajustar na ausência de informações (NOBRE et al., 1987NOBRE, P.R.C.; ROSA, A.N; SILVA, A.N.; EVANGELISTA, S.R.M. Curvas de crescimento de gado Nelore ajustadas para diferentes frequências de pesagens.Pesquisa Agropecuária Brasileira , v.22, p. 1027-1037, 1987.). Quanto ao DMA, a dificuldade dos modelos para estimar o peso ao nascimento proporcionou diferentes estimativas do peso nas subsequentes idades, contribuindo para aumentar o CV.

A amplitude observada para o C%, parâmetro o qual expressa a quantidade de animais que o modelo conseguiu convergir a curva de crescimento, foi resultante do menor número de curvas traçadas pelo modelo com quatro parâmetros (Meloun V), evidenciando a dificuldade de se trabalhar com modelos mais parametrizados.

A importância de cada avaliador da qualidade de ajuste depende do CV, sendo que baixos valores indicam equivalência dos modelos em relação ao avaliador. Considerando que as diferenças entre os modelos para os avaliadores R2aj e R2p foram irrelevantes, o AIC e BIC contribuíram estrategicamente para definição do modelo com melhor ajuste.

Dentre os quatro clusters formados, o primeiro, composto pelos modelos Mitscherlich, Meloun I e Brody, apresentou as melhores médias para os avaliadores AIC, BIC, R2aj, MEP, R2p e DMA. O segundo cluster, composto pelos modelos Meloun III, Von Bertalanffy e Gompertz, apresentou o maior valor para o percentual de convergência (Tabela 4).

Tabela 4
Clusters com respectivos modelos e médias dos avaliadores de qualidade do ajuste

Apesar de alguns modelos não-lineares serem apontados como adequados para descrever o crescimento de ovinos Santa Inês, existe divergências quanto à definição do melhor modelo. Enquanto Lobo et al. (2006)LÔBO, R.N.B.; VILLELA, L.C.V.; LÔBO, A.M.B.O.; PASSOS, J.R.S.; OLIVEIRA, A.A. Parâmetros genéticos de características estimadas da curva de crescimento de ovinos da raça Santa Inês. Revista Brasileira de Zootecnia , v.35, p.1012-1019, 2006. Supl.recomendam o modelo Richards, Guedes et al. (2004)GUEDES, M.H.P.; MUNIZ, J.A.; PEREZ, J.R.O.; SILVA, F.F.; AQUINO, L.H.; SANTOS, C.L. Estudo das curvas de crescimento de cordeiros das raças Santa Inês e Bergamácia considerando heterogeneidade de variâncias. Ciência Agrotécnica , v.28, p.381-388, 2004. e Sarmento et al. (2006)SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006., indicam o modelo Gompertz, e Silva et al. (2012)SILVA, L.S.A.; FRAGA, A.B.; SILVA, F.L. Growth curve in Santa Inês sheep. Small Ruminant Research , v.105, p.182-185, 2012. o modelo logístico. Essa dificuldade de precisar o melhor modelo pode ser em função da desuniformidade dos estudos quanto a quantidade de dados, oscilação nos valores dos pesos, número de pesagens por animal e idade da última pesagem (TORAL, 2008TORAL, F.L.B. Número e intervalo de pesagens para estimação de parâmetros de curvas de crescimento em bovinos. Revista Brasileira de Zootecnia , v.37, n.12, p.2120–2128, 2008.).

No entanto, vale salientar que a dificuldade na definição do melhor modelo deve-se também às particularidades da população avaliada, uma vez que os ovinos Santa Inês, segundo Teixeira Neto et al. (2015)TEIXEIRA NETO, M. R.; CRUZ, J. F.; CARNEIRO, P. L. S.; MALHADO, C. H. M.; BARBOSA, J. A.; SOUZA, L. E. B. Diversidade fenotípica de linhagens de ovinos Santa Inês por meio de análise multivariada. Revista Brasileira de Saúde e Produção Animal [online], v.16, n.4, p. 784-795, 2015., apresentam elevada variabilidade fenotípica.

Na escolha do modelo que melhor descreve o crescimento devem ser atendidos três pontos de igual importância: a interpretação biológica dos parâmetros, a qualidade de ajuste e o percentual de convergência (FITZHUGH JÚNIOR, 1976FITZHUGH JÚNIOR, H.A. Analysis of growth curves and strategies for altering their shape. Journal of Animal Science , v.42, n.4, p.1036-1051, 1976.). Nas situações em que os modelos são similares em relação aos parâmetros e que o modelo com melhor ajuste não é o que apresenta maior percentual de convergência, a escolha do modelo mais adequado fica a critério do pesquisador, priorizando o que melhor descreve o crescimento ou o que retrata o crescimento do maior número de animais.

De acordo com Cavalcante et al. (2013)CAVALCANTE, D.H.; CAMPELO, J.E.G.; SOUSA JÚNIOR, S.C.; SOUSA, G.G.T.; ARAÚJO, J.I.M.; ARAÚJO, A.C.; FONSECA, W.J.L.; BARROS JÚNIOR, C.P.; ARAÚJO, A.M. Modelos não paramétricos para ajustes de curva de crescimento em caprinos Sem Raça Definida (SRD). Revista Acadêmica: Ciências Agrárias e Ambientais, v.11, n.3, p.283-289, 2013., o melhor modelo a ser utilizado é aquele que melhor se adapta e apresenta resultados mais adequados ao caso de estudo.

No presente estudo, ao priorizar a qualidade de ajuste, os modelos Mitscherlich, Meloun I e Brody mostraram-se equivalentes. No entanto, o modelo Brody apresentou valores ligeiramente melhores para AIC, BIC e R2aj, sendo, por essa razão, escolhido para o cálculo das correlações e avaliação da influência do sexo e tipo de parto sobre os parâmetros da curva, bem como para descrever a taxa de crescimento absoluto (TCA) dos animais.

O efeito da interação sexo x tipo de parto não foi significativo (P>0,05). Animais oriundos de partos simples e duplos apresentaram parâmetros de crescimento diferentes (P<0,05), sendo que aqueles oriundos de partos simples apresentaram maior velocidade de crescimento. Por outro lado, os animais oriundos de parto duplo apresentaram maior peso à maturidade (Tabela 5).

Tabela 5
Peso assintótico (β1) constante de integração (β2) e taxa de maturidade (β3) estimados pelo modelo não-linear de Brody, de acordo com o sexo e tipo de parto para ovinos da raça Santa Inês

Esse resultados corroboram com os achados por Sarmento et al. (2006)SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006., quando verificaram que animais com maiores taxas de crescimento têm menor probabilidade de atingir maiores pesos à maturidade que aqueles que crescem mais lentamente no início da vida.

Machos e fêmeas apresentaram taxa de crescimento semelhante (P>0,05). Entretanto, os machos apresentaram maior peso à maturidade (P<0,05). De forma diferente ao verificado no presente estudo, Souza et al. (2011)SOUZA, L.A.; CARNEIRO, P.L.S.; MALHADO, C.H.M.; PAIVA, S.R.; CAÍRES, D.N.; BARRETO, D.L.F. Curvas de crescimento em ovinos da raça morada nova criados no estado da Bahia. Revista Brasileira de Zootecnia , v.40, n.8, p.1700-1705, 2011. avaliando o crescimento de ovinos Morada Nova pelo modelo Von Bertalanffy verificaram que o sexo influenciou os parâmetros β1 e β3.

No presente estudo, verificou-se correlações negativas do parâmetro β3 com β1 e β2 (-0,76 e -0,50; P<0,01). Por outro lado, os parâmetros β1 e β2 apresentaram correlação positiva elevada entre si (0,75). Do Ó et al. (2012)DO Ó, A.O.; RÊGO NETO, A.A.; SANTOS, G.V.; SARMENTO, J.L.R.; BIAGIOTTI, D.; SOUSA, J.E.R. Curva de crescimento de ovinos Santa Inês no Vale do Gurgueia. Revista Brasileira de Saúde e Produção Animal [online], v.13, n.4, p.912-922, 2012. e Lôbo et al. (2006)LÔBO, R.N.B.; VILLELA, L.C.V.; LÔBO, A.M.B.O.; PASSOS, J.R.S.; OLIVEIRA, A.A. Parâmetros genéticos de características estimadas da curva de crescimento de ovinos da raça Santa Inês. Revista Brasileira de Zootecnia , v.35, p.1012-1019, 2006. Supl. encontraram correlação negativa entre β1 e β3 de -0,26 e -0,87, em ovinos Santa Inês, respectivamente.

De acordo com Sarmento et al. (2006)SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006., uma possível conseqüência biológica da correlação negativa entre β2 e β3 seria que cordeiros desmamados mais leves (valores altos de β2) tendem apresentar menores taxas de crescimento relativo (quociente entre ganho de peso e peso) e, conseqüentemente, menores taxas de maturidade (menores valores de β3), do que aqueles desmamados mais pesados. Por outro lado, a correlação positiva entre β1 e β2poderia indicar que os cordeiros mais pesados ao desmame atingiriam maiores pesos à maturidade.

A relação biológica mais importante para uma curva está entre os parâmetros β1 e β3 (McMANUS et al., 2003McMANUS, C.; EVANGELISTA, C.; FERNANDES, L.A.C.; MIRANDA, R.M.; MORENO-BERNAL, F.E.; SANTOS, N.R. Curvas de crescimento de ovinos Bergamácia criados no Distrito Federal. Revista Brasileira de Zootecnia , v.32, p.1207- 1212, 2003.). O interessante seria que animais de produção de carne apresentassem maior velocidade de crescimento e maturidade o mais precocemente possível, reduzindo o tempo que o animal necessitaria para atingir a idade reprodutiva, melhorando, assim, a taxa de retorno da atividade.

O ponto de inflexão da curva ocorreu aos 94 dias de idade, momento em que os animais apresentaram TCA de 0,158kg/dia (Figura 1). A identificação do ponto de inflexão imediatamente após o desmame sugere a necessidade de ajuste no manejo nutricional para reduzir os efeitos da mudança de dieta sobre a TCA. Segundo Souza et al. (2013)SOUZA, L.D.A.; CARNEIRO, P.L.S.; MALHADO, C.H.M.; SILVA, F.F.; SILVEIRA, F.G.D. Traditional and alternative nonlinear models for estimating the growth of Morada Nova sheep. Revista Brasileira de Zootecnia , v.42, n.9, p.651-655, 2013., esta informação é de importância estratégica para os produtores, auxiliando na determinação da idade mais econômica para abate, oferecendo a melhor relação custo / benefício.

Figura 1
Taxa de crescimento absoluto, (A), ponto de máxima curvatura (B) e curva de crescimento estimada pelo modelo Brody (C) de ovinos Santa Inês

As informações sobre o desempenho animal expressas pelo TCA torna possível a determinação da idade ideal de abate. Segundo Silva et al. (2011)SILVA, F.L.; ALENCAR, M.M.; FREITAS, A.R.; PACKER, I.U.; MOURÃO, G.B. Curvas de crescimento em vacas de corte de diferentes tipos biológicos.Pesquisa Agropecuária Brasileira , v.46, n.3, p.262-271, 2011., após o ponto inflexão o consumo animal passa a ser destinado à mantença, ocorre redução da conversão alimentar e o custo para manutenção do animal é elevado. Desta forma o TCA possibilita a identificação do momento, em torno do qual o animal, com ganho de peso satisfatório, deve ser abatido ou descartado.

O ponto de inflexão e o TCA, evidentemente, podem ser influenciados por fatores ambientais e/ou genéticos. Estudos anteriores mostram ponto de inflexão e TCA de 200 dias de idade e 0,150kg/dia (SILVA et al., 2012SILVA, L.S.A.; FRAGA, A.B.; SILVA, F.L. Growth curve in Santa Inês sheep. Small Ruminant Research , v.105, p.182-185, 2012.) e de 52 dias de idade e 0,170kg/dia (SARMENTO et al., 2006SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006.), respectivamente. Além das diferenças ambientais, a variabilidade genética existente dentro da raça Santa Inês certamente contribui para ocorrência dessas variações.

Assim, a análise multivariada dos avaliadores de qualidade de ajuste se constitui em um facilitador na identificação de modelos não lineares para descrição da curva de crescimento em ovinos Santa Inês. Os modelos Mitscherlich, Meloun I e Brody são adequados e equivalentes para descrição do crescimento, sendo que, o modelo Brody revela uma notável redução do ganho de peso após o desmame.

REFERÊNCIAS

  • BROW, J.E.; FITZHUGH JÚNIOR, H.A.; CARTWRIGHT, T.C.A. A comparison of nonlinear models for describing weigt-age relationships in cattle.Journal of Animal Science , v.42, p.810-818, 1976.
  • CAVALCANTE, D.H.; CAMPELO, J.E.G.; SOUSA JÚNIOR, S.C.; SOUSA, G.G.T.; ARAÚJO, J.I.M.; ARAÚJO, A.C.; FONSECA, W.J.L.; BARROS JÚNIOR, C.P.; ARAÚJO, A.M. Modelos não paramétricos para ajustes de curva de crescimento em caprinos Sem Raça Definida (SRD). Revista Acadêmica: Ciências Agrárias e Ambientais, v.11, n.3, p.283-289, 2013.
  • CECON, P.R.; SILVA, F.F.; FERREIRA, A.; FERRÃO, R.G.; CARNEIRO, A.P.S.; DETMANN, E.; FARIA, P.N.; MORAIS, T.S.S. Análise de medidas repetidas na avaliação de clones de café ‘Conilon’. Pesquisa Agropecuária Brasileira , v.43, p.1171-1176, 2008.
  • CHALH, A.; EL GAZZAH, M. Variogram investigation of covariance shape within longitudinal data with possible use of a krigeage technique as an interpolation tool: Sheep growth data as an example. Irish Journal of Agricultural and Food Research , v.53, p.51–64, 2014.
  • CRUZ, G.R.B.; COSTA, R.G.; RIBEIRO, M.N. Curva de crescimento de caprinos mestiços no estado da Paraíba. Revista Brasileira de Ciências Agrárias , v.4, n.2, p.204-210, 2009.
  • DO Ó, A.O.; RÊGO NETO, A.A.; SANTOS, G.V.; SARMENTO, J.L.R.; BIAGIOTTI, D.; SOUSA, J.E.R. Curva de crescimento de ovinos Santa Inês no Vale do Gurgueia. Revista Brasileira de Saúde e Produção Animal [online], v.13, n.4, p.912-922, 2012.
  • FITZHUGH JÚNIOR, H.A. Analysis of growth curves and strategies for altering their shape. Journal of Animal Science , v.42, n.4, p.1036-1051, 1976.
  • GÓMEZ, D.A.A.; MUÑOZ, M.F.C.; BETANCUR, L.F.R. Modelación de las funciones de crecimiento aplicadas a la producción animal. Revista Colombiana de Ciencias Pecuárias , v.21, p.39-58, 2008.
  • GUEDES, M.H.P.; MUNIZ, J.A.; PEREZ, J.R.O.; SILVA, F.F.; AQUINO, L.H.; SANTOS, C.L. Estudo das curvas de crescimento de cordeiros das raças Santa Inês e Bergamácia considerando heterogeneidade de variâncias. Ciência Agrotécnica , v.28, p.381-388, 2004.
  • LÔBO, R.N.B.; VILLELA, L.C.V.; LÔBO, A.M.B.O.; PASSOS, J.R.S.; OLIVEIRA, A.A. Parâmetros genéticos de características estimadas da curva de crescimento de ovinos da raça Santa Inês. Revista Brasileira de Zootecnia , v.35, p.1012-1019, 2006. Supl.
  • MALHADO, C.H.M.; CARNEIRO, P.L.S.; SANTOS, P.F.; AZEVEDO, D.M.M.; SOUZA, J.C.; AFFONSO, P.R.M. Curva de crescimento em ovinos mestiços Santa Inês x Texel criados no Sudoeste do Estado da Bahia. Revista Brasileira de Saúde e Produção Animal [online], v.9, n.2, p.210-218, 2008.
  • MALHADO, C.H.M.; CARNEIRO, P.L.S.; AFFONSO, P.R.A.M.; SOUZA JÚNIOR, A.A.O.; SARMENTO, J.L.R. Growth curves in Dorper sheep crossed with the local Brazilian breeds, Morada Nova, Rabo Largo, and Santa Inês. Small Ruminant Research , v.84, p.16-21, 2009.
  • McMANUS, C.; EVANGELISTA, C.; FERNANDES, L.A.C.; MIRANDA, R.M.; MORENO-BERNAL, F.E.; SANTOS, N.R. Curvas de crescimento de ovinos Bergamácia criados no Distrito Federal. Revista Brasileira de Zootecnia , v.32, p.1207- 1212, 2003.
  • MELOUN, M.; MILITKY J. Sbírka úloh, Statistické zpracování experimentálnich dat . Pardubice: Universita Pardubice. 1996.
  • NOBRE, P.R.C.; ROSA, A.N; SILVA, A.N.; EVANGELISTA, S.R.M. Curvas de crescimento de gado Nelore ajustadas para diferentes frequências de pesagens.Pesquisa Agropecuária Brasileira , v.22, p. 1027-1037, 1987.
  • SARMENTOS, J.L.R.; REZAZZI, A.J.; SOUZA, W.H.; TORRES, R.A.; BREDA, F.C.; MENEZES, G.R.O. Estudo da curva de crescimento de ovinos Santa Inês.Revista Brasileira de Zootecnia , v.35, p.435-442, 2006.
  • SAS. Statistical Analysis System user’s guide . Version 9.1 ed. Cary, USA: SAS Institute, 2003.
  • SILVA, F.L.; ALENCAR, M.M.; FREITAS, A.R.; PACKER, I.U.; MOURÃO, G.B. Curvas de crescimento em vacas de corte de diferentes tipos biológicos.Pesquisa Agropecuária Brasileira , v.46, n.3, p.262-271, 2011.
  • SILVA, L.S.A.; FRAGA, A.B.; SILVA, F.L. Growth curve in Santa Inês sheep. Small Ruminant Research , v.105, p.182-185, 2012.
  • SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M; MUNIZ, J.A. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural , v.41, n.4, p.692-698, 2011.
  • SILVEIRA, F.G.; SILVA, F.F.; CARNEIRO, P.L.S.; MALHADO, C.H.M. Classificação multivariada de modelos de crescimento para grupos genéticos de ovinos de corte. Revista Brasileira de Saúde e Produção Animal [online], v.13, n1, p.62-73, 2012.
  • SOUZA, L.A.; CARNEIRO, P.L.S.; MALHADO, C.H.M.; PAIVA, S.R.; CAÍRES, D.N.; BARRETO, D.L.F. Curvas de crescimento em ovinos da raça morada nova criados no estado da Bahia. Revista Brasileira de Zootecnia , v.40, n.8, p.1700-1705, 2011.
  • SOUZA, L.D.A.; CARNEIRO, P.L.S.; MALHADO, C.H.M.; SILVA, F.F.; SILVEIRA, F.G.D. Traditional and alternative nonlinear models for estimating the growth of Morada Nova sheep. Revista Brasileira de Zootecnia , v.42, n.9, p.651-655, 2013.
  • TEIXEIRA NETO, M. R.; CRUZ, J. F.; CARNEIRO, P. L. S.; MALHADO, C. H. M.; BARBOSA, J. A.; SOUZA, L. E. B. Diversidade fenotípica de linhagens de ovinos Santa Inês por meio de análise multivariada. Revista Brasileira de Saúde e Produção Animal [online], v.16, n.4, p. 784-795, 2015.
  • TORAL, F.L.B. Número e intervalo de pesagens para estimação de parâmetros de curvas de crescimento em bovinos. Revista Brasileira de Zootecnia , v.37, n.12, p.2120–2128, 2008.

Datas de Publicação

  • Publicação nesta coleção
    Abr 2016

Histórico

  • Recebido
    05 Jul 2015
  • Aceito
    14 Dez 2015
UFBA - Universidade Federal da Bahia Avenida Adhemar de Barros nº 500 - Ondina , CEP 41170-110 Salvador-BA Brasil, Tel. 55 71 32836725, Fax. 55 71 32836718 - Salvador - BA - Brazil
E-mail: rbspa@ufba.br