Acessibilidade / Reportar erro

Estimação de parâmetros da cinética de trânsito de partículas em bovinos sob pastejo por diferentes seqüências amostrais

Estimation of solid transit kinetics parameters in cattle under grazing by different sample series

Resumos

Objetivou-se, neste estudo, avaliar a capacidade de ajustamento de um modelo não-linear no processo de estimação dos parâmetros da cinética de trânsito de partículas de bovinos em pastejo, empregando-se diferentes seqüências de amostragem fecal. Foram utilizados cinco novilhos F1 Limousin x Nelore, fistulados no esôfago e rúmen, sob pastejo de Brachiaria decumbens, com suplementação concentrada, durante o período das águas. O experimento constituiu de três períodos experimentais, conduzido em delineamento em blocos casualizados. Empregou-se como indicador o cromo mordante, produzido a partir de amostras de extrusa. Ajustou-se às curvas de excreção fecal do indicador o modelo duplo exponencial tempo-dependente, empregando-se as seguintes seqüências amostrais: SEQ 1 - 22 coletas (seqüência amostral completa); SEQ 2 e SEQ 3 - 17 coletas (redução de pontos de coleta nas fases ascendente e descendente da curva, respectivamente); e SEQ 4 e SEQ 5 - 13 e 10 coletas fecais (redução do número de coletas em todo o perfil da curva). As seqüências reduzidas foram produzidas a partir da omissão de pontos específicos da seqüência amostral completa (SEQ 1). A comparação entre as estimativas dos parâmetros cinéticos e da excreção fecal; a análise descritiva do número de iterações necessário à convergência do modelo e do coeficiente de determinação; e a avaliação de falta de ajustamento não apontaram diferenças entre as seqüências. A análise residual apontou, contudo, melhorias quanto ao comportamento gráfico e perfil de corridas de sinal dos resíduos com a redução do número de coletas para SEQ 4 e SEQ 5. Em função de pequena perda de eficiência, avaliada frente à variância residual, observada em SEQ 5, recomenda-se o emprego de 13 coletas fecais (SEQ 4) para avaliação da curva de excreção fecal do indicador em estudos desta natureza.

dose pulso; excreção fecal; indicadores externos; modelos matemáticos; taxa de passagem


The objective of this study was to evaluate the fitness capacity of a non-linear model in the process parameter estimation of the solid transit kinetics of bovine at pasture, using different sequences of fecal sampling. Five F1 Limousin X Nellore steers, fistulated in esophagus and rumen on Brachiaria decumbens grazing and concentrate supplemented during the rainy season period were used. The experiment was developed in three experimental periods; in a randomized completely blocks design. Chromium mordant, produced from extrusa samples, was used as marker. The double exponential time-dependent model was adjusted to the curves of fecal excretion of marker, and the following sample sequences were used: SEQ 1 - 22 samples (complete sample sequence); SEQ 2 and SEQ 3 - 17 samples (reduction of collection points in the ascendant and descendent phases of the curve, respectively); and SEQ 4 and SEQ 5 - 13 and 10 fecal samples (reduction of the number of samples in the whole profile of the curve). The reduced sequences were produced starting from the omission of specific points of the complete sample sequence (SEQ 1). The comparison among the estimates of the kinetic parameters and fecal excretion; the descriptive analysis of the necessary number of iterations to the convergence of the model and of the determination coefficient; and the evaluation of lack of fitness showed no differences among the sequences. However, the residual analysis pointed improvements in the graphic behavior and runs of signal profile of the residues with the reduction of the number of collections in SEQ 4 and SEQ 5. In function of small efficiency loss, evaluated from the residual variance, observed in SEQ 5, the employment of 13 fecal collections (SEQ 4) for evaluation of the fecal excretion curve of the marker in similar studies was recommended.

fecal excretion; external markers; mathematical models; pulse dose; rate of passage


Estimação de Parâmetros da Cinética de Trânsito de Partículas em Bovinos sob Pastejo por Diferentes Seqüências Amostrais

Edenio Detmann1 1 Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail: detmann@alunos.ufv.br , Paulo Roberto Cecon2 1 Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail: detmann@alunos.ufv.br , Mário Fonseca Paulino3 1 Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail: detmann@alunos.ufv.br , Joanis Tilemahos Zervoudakis1 1 Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail: detmann@alunos.ufv.br , Sebastião de Campos Valadares Filho3 1 Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail: detmann@alunos.ufv.br , Cláudio Vieira de Araújo1 1 Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail: detmann@alunos.ufv.br

RESUMO - Objetivou-se, neste estudo, avaliar a capacidade de ajustamento de um modelo não-linear no processo de estimação dos parâmetros da cinética de trânsito de partículas de bovinos em pastejo, empregando-se diferentes seqüências de amostragem fecal. Foram utilizados cinco novilhos F1 Limousin x Nelore, fistulados no esôfago e rúmen, sob pastejo de Brachiaria decumbens, com suplementação concentrada, durante o período das águas. O experimento constituiu de três períodos experimentais, conduzido em delineamento em blocos casualizados. Empregou-se como indicador o cromo mordante, produzido a partir de amostras de extrusa. Ajustou-se às curvas de excreção fecal do indicador o modelo duplo exponencial tempo-dependente, empregando-se as seguintes seqüências amostrais: SEQ 1 - 22 coletas (seqüência amostral completa); SEQ 2 e SEQ 3 - 17 coletas (redução de pontos de coleta nas fases ascendente e descendente da curva, respectivamente); e SEQ 4 e SEQ 5 - 13 e 10 coletas fecais (redução do número de coletas em todo o perfil da curva). As seqüências reduzidas foram produzidas a partir da omissão de pontos específicos da seqüência amostral completa (SEQ 1). A comparação entre as estimativas dos parâmetros cinéticos e da excreção fecal; a análise descritiva do número de iterações necessário à convergência do modelo e do coeficiente de determinação; e a avaliação de falta de ajustamento não apontaram diferenças entre as seqüências. A análise residual apontou, contudo, melhorias quanto ao comportamento gráfico e perfil de corridas de sinal dos resíduos com a redução do número de coletas para SEQ 4 e SEQ 5. Em função de pequena perda de eficiência, avaliada frente à variância residual, observada em SEQ 5, recomenda-se o emprego de 13 coletas fecais (SEQ 4) para avaliação da curva de excreção fecal do indicador em estudos desta natureza.

Palavras-chave: dose pulso, excreção fecal, indicadores externos, modelos matemáticos, taxa de passagem

Estimation of Solid Transit Kinetics Parameters in Cattle Under Grazing by Different Sample Series

ABSTRACT - The objective of this study was to evaluate the fitness capacity of a non-linear model in the process parameter estimation of the solid transit kinetics of bovine at pasture, using different sequences of fecal sampling. Five F1 Limousin X Nellore steers, fistulated in esophagus and rumen on Brachiaria decumbens grazing and concentrate supplemented during the rainy season period were used. The experiment was developed in three experimental periods; in a randomized completely blocks design. Chromium mordant, produced from extrusa samples, was used as marker. The double exponential time-dependent model was adjusted to the curves of fecal excretion of marker, and the following sample sequences were used: SEQ 1 - 22 samples (complete sample sequence); SEQ 2 and SEQ 3 - 17 samples (reduction of collection points in the ascendant and descendent phases of the curve, respectively); and SEQ 4 and SEQ 5 - 13 and 10 fecal samples (reduction of the number of samples in the whole profile of the curve). The reduced sequences were produced starting from the omission of specific points of the complete sample sequence (SEQ 1). The comparison among the estimates of the kinetic parameters and fecal excretion; the descriptive analysis of the necessary number of iterations to the convergence of the model and of the determination coefficient; and the evaluation of lack of fitness showed no differences among the sequences. However, the residual analysis pointed improvements in the graphic behavior and runs of signal profile of the residues with the reduction of the number of collections in SEQ 4 and SEQ 5. In function of small efficiency loss, evaluated from the residual variance, observed in SEQ 5, the employment of 13 fecal collections (SEQ 4) for evaluation of the fecal excretion curve of the marker in similar studies was recommended.

Key Words: fecal excretion, external markers, mathematical models, pulse dose, rate of passage

Introdução

As transformações digestivas são determinadas por atributos intrínsecos do alimento e por sua interações com os processos cinéticos. Neste enfoque, a expressão quantitativa dos processos cinéticos de digestão e passagem torna-se necessária para estimar mais precisamente a quantidade e composição dos nutrientes digeridos e sua subseqüente eficiência de utilização pelo animal (ELLIS et al., 1994).

A estimação da cinética de trânsito de partículas em ruminantes tem sido realizada, normalmente, por intermédio do emprego de indicadores complexados com a fração fibrosa do alimento, em procedimento de dose única, com amostragens fecais subseqüentes em intervalos de tempo pré-definidos, visando a caracterização da curva de excreção deste indicador, a qual, é submetida ao ajuste de modelos não-lineares por processos iterativos.

No entanto, ao se trabalhar com animais em pastejo, algumas considerações devem ser tecidas sobre o método de estimação descrito acima. Para a realização das coletas fecais, há a necessidade de contenção dos animais. Isto leva à ampliação do estresse, o qual pode alterar o comportamento natural de pastejo do animal (OWENS e HANSON, 1992; BURNS et al., 1994), podendo levar a alterações desconhecidas sobre as estimativas obtidas. Assim, é desejável que, ao trabalhar-se com animais em pastejo, estes se comportem o mais próximo possível do normal, sem grandes interferências, o que poderia ser alcançado com menor número de contenções (FRANCE et al., 1988).

Em estudos com animais confinados, tem-se utilizado, para estimar os parâmetros cinéticos de trânsito, número de coletas variando, normalmente, entre 20 (OLIVEIRA et al., 1999) e 26 (MOORE et al., 1992). Por outro lado, em estudos com animais a pasto, empregam-se, normalmente, de 17 (LIRA, 2000) a 25 (SOARES et al., 1999) coletas que podem ser consideradas elevadas, dado o curto intervalo de tempo no qual são realizadas (120 a 144 horas).

Objetivou-se, neste estudo, avaliar a capacidade de ajustamento de um modelo não-linear para estimar a cinética de trânsito de partículas de animais em pastejo a diferentes intensidades e seqüências de amostragem fecal.

Material e Métodos

O presente trabalho foi realizado na Fazenda Experimental de Felixlândia, em Felixlândia-MG, da Empresa de Pesquisa Agropecuária de Minas Gerais (EPAMIG), no período de dezembro de 1997 a março de 1998.

A área experimental, com 5,35 hectares e topografia plana, coberta uniformemente com a gramínea Brachiaria decumbens, foi redividida, por intermédio de cerca eletrificada, em cinco piquetes de 1,07 hectares, providos de bebedouro e comedouro coberto.

Foram utilizados cinco bovinos F1 Limousin x Nelore, castrados, com idade e peso médios iniciais de 22 meses e 396 kg, fistulados no esôfago e rúmen.

O experimento constou de três períodos experimentais de 31 dias, sendo os sete primeiros de adaptação, conduzidos durante a parte central do período chuvoso. Ao início de cada período experimental, os animais foram sorteados, com relação aos suplementos e piquetes.

Foram fornecidos os seguintes suplementos: sal mineral; suplementos constituídos de milho e farelo de soja, em nível de 1 e 2 kg/dia; e suplementos constituídos de farelo de trigo e farelo de soja, fornecido, em nível de 1 e 2 kg/dia. Todos foram balanceados para se atingir o nível de 20% de proteína bruta, com base na matéria natural, sendo fornecidos diariamente aos animais às 10 h. Os animais em todos os tratamentos tiveram acesso irrestrito à água e mistura mineral.

A partir de amostras de extrusa, coletadas no terceiro dia de cada período experimental, produziu-se fibra mordantada com dicromato de sódio, adaptando-se os procedimentos de ÚDEN et al. (1980). Alíquotas de 100 g foram acondicionadas em cartuchos de papel, as quais foram colocadas diretamente no rúmen de cada animal no 26o dia do período experimental, às 8 h. As amostras de fezes foram coletadas entre 0 e 144 horas após a dosagem, num total de 22 coletas, sendo imediatamente colocadas em sacos plásticos, identificadas e congeladas a -20ºC.

Posteriormente, as amostras fecais e o mordante foram analisados por espectrofotometria de absorção atômica, a fim de quantificar o teor de cromo, conforme metodologia descrita por WILLIANS et al. (1962).

Às curvas de excreção do indicador ajustou-se, por intermédio do procedimento estatístico iterativo Gauss-Newton, o modelo não-linear duplo exponencial tempo-dependente simplificado, proposto por DHANOA et al. (1985), o qual tem demonstrado melhor capacidade de ajustamento sob condições tropicais em confinamento (OLIVEIRA et al., 1999) e a pasto (LIRA, 2000), dado por:

em que é concentração fecal estimada do indicador no tempo t; A e B, parâmetros biologicamente indefinidos; k1 e k2, taxas de passagem pelo primeiro e segundo compartimentos, respectivamente; e t, tempo após o fornecimento do indicador.

O modelo foi considerado satisfatoriamente ajustado quando os seguintes critérios foram supridos: (1) as estimativas dos parâmetros cinéticos apresentaram valores biologicamente coerentes (DHANOA et al., 1985); e (2) convergência para uma solução única a partir do fornecimento de diferentes valores iniciais (DHANOA et al., 1985; SOUZA, 1998). De forma a melhorar o diagnóstico do ajuste, empregaram-se, de forma auxiliar, os testes t, para as estimativas dos parâmetros, e F, para a regressão. Estes testes forneceram diagnósticos apenas aproximados, uma vez que, sob condições não-lineares, os resíduos ordinários não apresentam distribuição aproximadamente normal (REGAZZI, 1999).

As seqüências amostrais empregadas para o ajustamento possuem as estruturas temporais descritas na Tabela 1. A SEQ 1 constitui-se de todas as amostragens realizadas (22 coletas), sendo as demais estruturadas a partir da omissão dos valores de concentração de tempos preestabelecidos durante o processo de ajustamento. Objetivou-se em SEQ 2 e SEQ 3 (17 coletas) observar a capacidade de ajustamento do modelo à redução amostral na fase ascendente e descendente da curva de excreção, respectivamente. A SEQ 4 (13 coletas) constitui a fusão de SEQ 2 e SEQ 3, com menor intensidade de coletas em todo o período amostral, sendo este o mesmo objetivo proposto em SEQ 5, porém, a redução foi mais severa (10 coletas). Para efeito de comparação, SEQ 1 foi considerada como grupo controle, servindo, portanto, como referência às estimativas obtidas pelas demais seqüências.

Após o ajustamento das curvas de excreção do indicador, procedeu-se à determinação da excreção fecal, adaptando-se a equação descrita por HOLLEMAN e WHITE (1989), definido-se com limites inferior e superior para integração da função t=0 e t=300 horas (DETMANN et al., 2000b), respectivamente:

em que EF é excreção fecal (kg/h); D, dose de cromo (mg); e 0ò300 Ct dt, integral da equação ajustada para concentração do indicador nas fezes, em função do tempo pós-dosagem (mg/kg × h). A título de comparação, os valores de excreção fecal foram avaliados em termos de valores absolutos (kg/dia) e relativos (g/100kgPV), com o objetivo de eliminar possível interferência de variações de peso entre os animais.

As seqüências propostas foram comparadas entre si, considerando-se como referência base a série amostral completa (SEQ 1), segundo os seguintes critérios: (1) comparação entre as médias das estimativas dos parâmetros k1, k2 e excreção fecal; (2) análise residual por intermédio de corridas de sinal e plotagem de resíduos "studentizados"; (3) eficiência relativa em função da variância residual; (4) ocorrência de falta de ajustamento; (5) número de iterações necessárias para convergência à solução do modelo; e (6) análise descritiva do coeficiente de determinação.

A análise de variância para estimativas dos parâmetros cinéticos e da excreção fecal seguiu o modelo:

em que Yijkl é valor observado relativo ao suplemento i, ao período j, à seqüência k e à observação l; m, média geral; Si,efeito do suplemento i (i = 1, 2, 3, 4 e 5); Pj, efeito do período ou bloco j ( j = 1, 2 e 3); eij, efeito residual das parcelas; Ak, efeito da seqüência amostral k (k = 1, 2, 3, 4 e 5); SAik, efeito da interação do suplemento i com a seqüência amostral k; e eijkl, erro aleatório, associado a cada observação, pressuposto normal e independentemente distribuído, com média zero e variância s². Para realização das comparações entre médias, empregou-se o teste Dunnet, adotando-se SEQ 1 como referência. Adotou-se a = 0,05.

Para comparação das variâncias residuais dos modelos ajustados propôs-se neste estudo o índice de eficiência média, descrito abaixo:

em que EMk é eficiência média da seqüência k, para k ¹ 1; s²ij1l, estimativa da variância para o suplemento i, no período j, para a seqüência 1 (SEQ 1), na repetição l, aninhada aos efeitos principais; s²ijkl, estimativa da variância para o suplemento i, no período j, para a seqüência k, sendo k ¹ 1, na repetição l, aninhada aos efeitos principais; e n = 15; sendo k fixo e ¹ 1 e i, j e l variáveis.

Para verificação da falta de ajustamento dos modelos ajustados às diferentes seqüências, adaptou-se o procedimento descrito por DRAPER e SMITH (1966). Análises prévias não revelaram efeitos significativos da suplementação sobre os parâmetros cinéticos (DETMANN et al., 2000a). Assim, aos valores obtidos para todos os animais, dentro de cada período, ajustou-se um modelo para cada seqüência, de forma a permitir a avaliação da existência de vício na estimativa do resíduo da regressão frente ao resíduo puro. Empregou-se o teste F, adotando-se a = 0,05.

Os resíduos ordinários foram definidos como a diferença entre os valores observados e estimados da concentração fecal do indicador no tempo "t". Para melhor interpretação e comparação entre seqüências, procedeu-se à "studentização" dos resíduos, segundo a equação (CORDEIRO e PAULA, 1989; SOUZA, 1998):

em que é resíduo "studentizado" para a observação i; , resíduo ordinário observação i; , desvio-padrão residual; e , elemento i da matriz de projeção , sendo a matriz jacobiana da função resposta . Os resíduos foram plotados em gráficos individuais, para cada seqüência, e avaliados quanto à dispersão.

Segundo DRAPER e SMITH (1966), quando a série temporal de um conjunto de resíduos é conhecida, pode-se avaliar o ajustamento do modelo por intermédio de grupos de resíduos consecutivos de mesmo sinal, os quais são denominados corridas de sinal. Uma vez que as seqüências propostas são compostas de números diferentes de observações, foi proposta, neste estudo, a avaliação das corridas de sinal por intermédio da razão entre o número de corridas e o número de observações, valor este que deve ser o maior possível. Outra avaliação residual, por intermédio de corridas de sinal, foi realizada segundo proposta de BARD (1974), por intermédio da equação: m = 2n1n2/(n1+n2 ) + 1, em que m é número esperado de corridas de sinal, sob a suposição de aleatoriedade dos erros, e n1 e n2 são número de resíduos negativos e positivos, respectivamente.

O coeficiente de determinação foi estimado como o quadrado do coeficiente de correlação entre os valores observados e preditos para a variável dependente (SOUZA, 1998; REGAZZI, 1999).

Todas as análises estatísticas foram realizadas por intermédio dos programas SAEG, versão 5.0, e SAS, versão 6.12.

Resultados e Discussão

A grande maioria dos métodos estatísticos empregados em análises residuais de modelos de regressão pressupõe que os erros são normal e independentemente distribuídos, com média zero e variância s². Se observadas falhas nestas pressuposições, pode-se afirmar que o modelo é incoerente ao fenômeno estudado (DRAPER e SMITH, 1966; MYERS, 1990). Embora sob condições não-lineares tenha-se comportamento distinto, com esperança de resíduos diferente de zero e presença de correlação destes com os valores preditos da variável dependente (CORDEIRO e PAULA, 1989; SOUZA, 1998), a análise residual assume importância semelhante em modelagem não-linear, funcionando como instrumento fundamental para o diagnóstico de atipicidades, como erros de especificação de resposta e heterocedasticidade (SOUZA, 1998). Devido à comprovada eficiência do modelo empregado neste trabalho em estimar os parâmetros relacionados à cinética de trânsito do trato gastrointestinal de ruminantes (DHANOA et al., 1985; OLIVEIRA et al., 1999; e LIRA, 2000), foram interpretadas as possíveis condições anormais observadas na análise residual como resultado da menor adequação do conjunto de dados em ser submetido ao ajuste do modelo proposto. Tanto para análise gráfica, como para o avaliação do perfil de corridas de sinal, omitiram-se as informações pertinentes ao tempo = 0 (VIEIRA et al., 1997).

Os gráficos demonstrando a plotagem dos resíduos "studentizados", em função da variável independente tempo para as diferentes seqüências amostrais, podem ser observados na Figura 1. Em termos ideais, espera-se que o conjunto de projeções horizontais dos resíduos plotados contra uma variável constituinte do modelo se aproxime de uma linha reta, ou seja, apresente-se constante, independente do nível desta variável em questão (DRAPER e SMITH, 1966). Observa-se, contudo, que para todas as seqüências avaliadas houve tendência de redução da dispersão dos resíduos, em relação à nulidade, com o aumento do tempo pós-aplicação do indicador. Este comportamento, denominado "efeito funil", detecta possível condição de heterocedasticidade (MYERS, 1990), o que reflete certo grau de correlação entre os resíduos, comumente observado em séries temporais (DRAPER e SMITH, 1966).


Primariamente, procurou-se buscar seqüências amostrais mais simples, que apresentassem características de ajustamento semelhantes à série completa (SEQ 1), considerada aqui como controle. Contudo, ao contrário do esperado, a redução do número de coletas fecais parece ter exercido efeito positivo sobre a dispersão residual, tornando a distribuição mais próxima da constância, como observado em SEQ 4 e SEQ 5.

Os resíduos podem desenvolver tendências e outros afastamentos da aleatoriedade que implicam em modificações no modelo (BARD, 1974), podendo ser sugestivo de erros de especificação de respostas (SOUZA, 1998). A maior falta de agregação entre os valores preditos e observados na fase ascendente da curva, semelhante ao encontrado em SEQ 1, SEQ 2 e SEQ 3, pode levar, em muitos casos, a estimativas viesadas da taxa constante de valor mais elevado (k2). DHANOA et al. (1985) observaram comportamento semelhante, ao trabalharem com o modelo tempo-independente de GROVUM e WILLIANS (1973), e sugeriram que, nestes casos, a transformação dos dados para uma escala logarítmica permitiria a estabilização da variância com pouco ou nenhum comprometimento sobre as estimativas dos parâmetros. SOUZA (1998) sugeriu ser a homocedasticidade condição desejada e afirmou que transformações não-lineares podem ser aplicadas à variável independente e à função resposta, no caso de independência dos resíduos, com este objetivo.

A descrição dos modelos quanto à freqüência de ocorrência de corridas de sinal é exposta na Tabela 2. Observa-se que a redução do número de coletas tornou mais elevado o número de corridas de sinal para os modelos ajustados, o que é demonstrado pela elevação na ocorrência de razões corridas de sinal/pontos amostrais de maior valor.

O sucesso de qualquer método de ajustamento de curvas pode ser julgado pela habilidade do modelo em descrever certa associação entre variáveis sem super ou subestimar alguma seção do perfil estudado (DHANOA et al., 1995; VIEIRA et al., 1997; e LÓPEZ et al., 1999). Um modelo ajustado que tem baixo número de corridas pode apresentar este problema (LÓPEZ et al., 1999), podendo-se ainda inferir sob estas circunstâncias que a aleatoriedade dos resíduos é suspeita (BARD, 1974). Assim, quanto maior o número de corridas de sinal em um perfil, maior o grau de aleatoriedade dos resíduos, tornando melhor a descrição do perfil (DHANOA et al., 1995). Neste enfoque, semelhantemente à análise de dispersão residual a redução no número de pontos amostrais implicou em melhorias no ajustamento do modelo aos dados, sendo SEQ 2 e SEQ 3 semelhantes entre si, porém inferiores à SEQ 4 e SEQ 5, que apresentaram as melhores freqüências.

O comportamento das seqüências com relação às corridas de sinal mostrou-se semelhante em relação ao número esperado de corridas, sob a suposição de aleatoriedade dos resíduos (Tabela 2), com os melhores resultados observados em SEQ 4 e SEQ 5. Estes resultados corroboram o perfil gráfico mais estabilizado dos resíduos observados nestas seqüências. Embora alguns testes estatísticos, como o de "Swed e Eisenhart", estejam disponíveis para possibilitar inferências sob o número de corridas de sinal, optou-se por sua não aplicação, devido à natureza aproximada, e não exata, do nível probabilístico observado em séries temporais (DRAPER e SMITH, 1966). Conclui-se, contudo, que, descritivamente, a redução no número de coletas contribuiu para que maior parte dos modelos ajustados apresentasse estrutura residual próxima à aleatoriedade.

Os resultados obtidos para o índice de eficiência média, em relação à seqüência controle, e para o número médio de iterações necessários à convergência do modelo são expressos na Tabela 3. A eficiência é uma medida que reflete o grau relativo de precisão de uma seqüência em relação à seqüência controle (SEQ 1). Observa-se que, apesar dos valores mais elevados obtidos para SEQ 3 e SEQ 5, o que indica alguma perda de eficiência, infere-se que os valores se encontram extremamente próximos à unidade, mostrando um comportamento semelhante entre seqüências no tocante à variação aleatória.

Grande número de iterações pode ser indicativo de modelos inapropriados ou mal-acondicionados (DHANOA et al., 1995). Alguns trabalhos conduzidos com modelos de cinéticas de trânsito e digestão têm sugerido como índice satisfatório valores menores que 20 (DHANOA et al., 1985, 1995; e LÓPEZ et al., 1999). Para todas as seqüências avaliadas, o número de iterações foi considerado extremamente satisfatório e semelhantes entre si, sempre com valores médios inferiores a 10.

A partição da variância residual da regressão, em falta de ajuste e erro puro, permite a detecção da ocorrência de vício no processo de estimação, que serve como indicativo de modelo inapropiado (SNEDECOR e COCHRAN, 1976; LÓPEZ et al., 1999). No entanto, em nenhum dos modelos conjuntos ajustados para cada seqüência, observou-se a ocorrência de falta de ajuste (P>0,05), corroborando a adequação do modelo empregado ao fenômeno estudado e mostrando que neste enfoque a redução do número de coletas não afetou o processo estimativo.

O coeficiente de determinação expressa a bondade de ajustamento ou adequação de ajuste do modelo à seqüência de dados (SOUZA, 1998). Como observado na Tabela 4, os valores obtidos para todas as seqüências foram extremamente elevados, demonstrando alto grau de associação entre os valores de concentração fecal do indicador preditos e observados, o que permite concluir que, independente do número de coletas, o modelo pôde explicar a maior parte da variação total observada sobre a variável em questão. No entanto, segundo LÓPEZ et al. (1999), se o R² está localizado muito próximo à unidade, o que indica alto nível de significância associado à razão de variâncias, este pode ser melhor utilizado como medida global de ajustamento do que como medida de comparação entre modelos.

Comportamento semelhante aos demais índices foi observado sobre as estimativas dos parâmetros da cinética de trânsito e da excreção fecal (Tabela 5). Nenhuma diferença entre seqüências foi detectada (P>0,05), o que leva novamente a inferir que a redução no número de coletas não afeta o processo de obtenção de estimativas da cinética de trânsito em ruminantes. Os valores obtidos nesta pesquisa para k1 e k2 agregam aos descritos por DHANOA et al. (1985), OLIVEIRA et al. (1999) e LIRA (2000). Segundo ELLIS et al. (1994), a razão k2/k1 deve exceder 1,5 para que a estimativa do parâmetro taxa para o segundo compartimento seja confiável e se estabeleça uma distinção real da existência de dois compartimentos, comportamento obtido satisfatoriamente para todas as seqüências avaliadas.

O processo de estimação de excreção fecal por intermédio de modelos matemáticos e cálculo diferencial e integral é recente (HOLLEMAN e WHITE, 1989) e tem sido pouco empregado em estudos de digestão. Avaliações recentes têm buscado aferir esta técnica (SUSMEL et al., 1996; DETMANN et al., 2001), a qual apresenta a grande vantagem de simplificar os procedimentos experimentais, unindo em um única metodologia a estimação de parâmetros cinéticos e da excreção fecal (DETMANN et al., 2001). Os resultados obtidos neste trabalho mostram que, semelhantemente aos parâmetros cinéticos, a redução do número de coletas fecais não compromete este procedimento.

Conclusões

A redução do número de coletas fecais na técnica dose pulso não altera as estimativas obtidas para os parâmetros da cinética de trânsito do trato gastrointestinal de bovinos em pastejo, a adequação do ajuste e o número de iterações necessárias para a convergência dos modelos. Por outro lado, a redução do número de pontos amostrais parece exercer efeitos positivos sobre o comportamento residual, melhorando o grau de aleatoriedade dos resíduos em série temporal.

Recomenda-se, em virtude de pequena perda de eficiência da seqüência com 10 coletas e na ausência de maiores estudos, o emprego de 13 coletas fecais, representadas aqui pela SEQ 4, na aplicação de técnica semelhante, o que possibilitará a redução do nível de interferência sobre o comportamento de animais manejados a pasto.

Recebido em: 15/04/00

Aceito em: 19/11/00

2 Professor, DPI-UFV, bolsista do CNPq.

3 Professor DZO-UFV, bolsista do CNPq.

  • BARD, Y. 1974. Nonlinear parameter estimation. New York: Academic Press. 341p.
  • BURNS, J.C., POND, K.R., FISHER, D.S. 1994. Measurement of forage intake. In: FAHEY JR., G.C. (Ed.) Forage quality, evaluation, and utilization Winsconsin: American Society of Agronomy. p.494-532.
  • CORDEIRO, G.M., PAULA, G.A. 1989. Modelos de regressăo para análise de dados univariados. Rio de Janeiro: IPMA/CNPq. 353p.
  • DETMANN, E., PAULINO, M.F., ZERVOUDAKIS, J.T. et al. Estimaçăo da excreçăo fecal por intermédio de modelos matemáticos em novilhos mestiços, suplementados, a pasto. In: REUNIĂO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 37, Viçosa, 2000. Anais.. Viçosa: SBZ, 2000. p.437.
  • DETMANN, E., PAULINO, M.F., ZERVOUDAKIS, J.T. et al. 2001. Suplementaçăo de novilhos mestiços durante o período das águas: parâmetros ingestivos e digestivos. Rev. bras. zootec. (no prelo).
  • DHANOA, M.S., SIDDONS, R.C., FRANCE, J. et al. 1985. A multicompartimental model to describe marker excretion patterns in ruminant faeces. Br. J. Nut., 53(2):663-671.
  • DHANOA, M.S., FRANCE, J., SIDDONS, R.C. et al. 1995. A non-linear compartimental model to describe forage degradation kinetics during incubation in polyester bags in the rumen. Br. J. Nutr., 73(1):3-15.
  • DRAPER, N.R., SMITH, H. 1966. Applied regression analysis. New York: Jonh Willey & Sons. 407p.
  • ELLIS, W.C., MATIS, J.H., HILL, T.M. et al. 1994. Methodology for estimating digestion and passage kinetics of forages. In: FAHEY JR., G.C. (Ed.) Forage quality, evaluation, and utilization Winsconsin: American Society of Agronomy. p.682-756.
  • FRANCE, J., DHANOA, M.S., SIDDONS, R.C. et al. 1988. Estimating the fecal producing by ruminants from faecal marker concentration curves. J. Theor. Biol., 135(2):383-391.
  • GROVUM, W.L., WILLIANS, V.J. 1973. Rate of passage of digesta in sheep. 4. Passage of marker through the alimentary tract and the biological relevance of rate-constants derived from the changes in concentration of markers in faeces. Br. J. Nutr., 30(2):313-329.
  • HOLLEMAN, D.F., WHITE, R.G. 1989. Determination of digesta fill and passage rate from nonabsorbed particulate phase markers using the single dose method. Can. J. Zool., 67(2):488-494.
  • LIRA, V.M.C. Utilizaçăo de diferentes modelos matemáticos e marcadores para simulaçăo de cinética digestiva e de trânsito do capim braquiária (Brachiaria decumbens Stapf.) Viçosa, MG: UFV, 2000, 91p. Dissertaçăo (Mestrado em Zootecnia) - Universidade Federal de Viçosa, 2000.
  • LÓPEZ, S., FRANCE, J., DHANOA, M.S. et al. 1999. Comparison of mathematical models to describe disappearance curves obtained using the polyester bag technique for incubating feeds in the rumen. J. Anim. Sci., 77(5):1875-1888.
  • MOORE, J.A., POND, K.R., POORE, M.H. et al. 1992. Influence of model and marker on digesta kinetic estimate for sheep. J. Anim. Sci., 70(11):3528-3540.
  • MYERS, R.H. 1990. Classical and modern regression with applications 2.ed. Boston: PWS-KENT Publishing Co. 488p.
  • OLIVEIRA, R.L., PEREIRA, J.C., SILVA, P.R.C. et al. 1999. Degradabilidade ruminal da cama de frango e do feno de capim Coast-cross e avaliaçăo de modelos matemáticos para estimativa da taxa de passagem de partículas. R. Bras. Zootec., 28(4):839-849.
  • OWENS, F.N., HANSON, C.F. 1992. External and internal markers for appraising site and extent of digestion in ruminants. J. Dairy Sci., 75(9):2605-2617.
  • REGAZZI, A.J. 1999. Modelos de regressăo năo-linear e aplicaçőes utilizando o SAS (PROC NLIN). Viçosa: UFV. 116p. (Apostila).
  • SNEDECOR,G.W., COCHRAN,W.G. 1976. Statistical methods 6.ed. Iowa: Iowa University Press. 593p.
  • SOARES, J.P.G., AROEIRA, L.J.M., PEREIRA, O.G. et al. 1999. Capim-elefante (Pennisetum purpureum Schum.) sob duas doses de nitrogęnio. Consumo e produçăo de leite. R. Bras. Zootec., 28(4):889-897.
  • SOUZA, G.S. 1998. Introduçăo aos modelos de regressăo linear e năo-linear. Brasília: EMBRAPA-SPI. 505p.
  • SUSMEL, P., STEFANON, B., SPANGHERO, M. et al. 1996. Ability of mathematical models to predict faecal output with a pulse dose of indigestible marker. Brit. J. Nut., 75(4):521-532.
  • UDÉN, P., COLUCCI, P.E, VAN SOEST, P.J. 1980. Investigation of chromium, cerium and cobalt as markers in digesta. Rate of passage studies. J. Sci. Food Agr., 31(7):625-632.
  • VIEIRA, R.A.M., PEREIRA, J.C., MALAFAIA, P.A.M. et al. 1997. Application of non-linear models in the description of in situ degradation profiles of the elephant grass (Pennisetum purpureum Schum., Mineiro variety). Anim. Feed Sci. Tec., 66:197-210.
  • WILLIANS, C.H., DAVID, D.J., IISMA, O. 1962. The determination of chromic oxide in faeces samples by atomic absorption spectrophotometry. J. Agric. Sci., 59(3):381-385.
  • 1
    Zootecnista, MS, estudante de Doutorado, DZO-UFV, bolsista do CNPq E-mail:
  • Datas de Publicação

    • Publicação nesta coleção
      18 Set 2001
    • Data do Fascículo
      Fev 2001

    Histórico

    • Aceito
      19 Nov 2000
    • Recebido
      15 Abr 2000
    Sociedade Brasileira de Zootecnia Universidade Federal de Viçosa / Departamento de Zootecnia, 36570-900 Viçosa MG Brazil, Tel.: +55 31 3612-4602, +55 31 3612-4612 - Viçosa - MG - Brazil
    E-mail: rbz@sbz.org.br