Parameters estimation in the model for in situ degradability of mertens and loften

Estimação dos parâmetros no modelo para degradabilidade in situ de Mertens e Loften

ZOOTECNIA E MEDICINA VETERINÁRIA

Parameters estimation in the model for in situ degradability of mertens and loften¹

Estimação dos parâmetros no modelo para degradabilidade in situ de Mertens e Loften

Joel Augusto Muniz^I; Taciana Villela Savian^II; João Domingos Scalon^III

^IEngenheiro Agrônomo, Doutor, Professor Titular - Departamento de Ciências Exatas/DEX - Universidade Federal de Lavras/UFLA - Cx. P. 3037 - 37200-000 - Lavras, MG - [e-mail]

^IIZootecnista, Mestre em Estatística e Experimentação Agropecuária - Departamento de Ciências Exatas/DEX - Universidade Federal de Lavras/UFLA - Cx. P. 3037 - 37200-000 - Lavras, MG - [e-mail]

^IIIEstatístico, PhD, Professor Associado - Departamento de Ciências Exatas/DEX - Universidade Federal de Lavras/UFLA - Cx. P. 3037 - 37200-000 - Lavras, MG - [e-mail]

ABSTRACT

The aim of this study is to evaluate the behavior of the parameters of the degradation model proposed by Mertens & Loften (1980) fitted to the results of a rehearsal of degradability in situ. In the experiment, we evaluate the potentially degradable residue of neutral detergent fiber (NDF) of coastcross grass (Cynodon dactylon x Cynodon nlemfunensis) cut at 60 days, with three replications. The potentially degradable residue of NDF is studied using fifteen incubation times (0; 0.5; 1; 3; 6; 9; 12; 18; 24; 36; 48; 56; 72; 96 and 120 hours). The experimental plot is comprised of a non-lactating cow with a permanent ruminal fistula. It is obtained mean and individual fits for the animals. Variances of parameter estimators is also obtained through both the covariance matrix of the parameters and the jackknife method, then resulting expressions for the estimate of the confidence interval for the parameters of the model. The results shows that the jackknife method presents larger variance estimate for the parameters of the model of Mertens & Loften (1980), resulting in confidence intervals of greater amplitude and less precise parameter estimates for both the individual and mean fits.

Index terms: Nonlinear regression, ruminal degradability, jackknife method, coast-cross grass.

RESUMO

Objetiva-se avaliar o comportamento dos parâmetros do modelo de degradação proposto por Mertens & Loften (1980) ajustado aos resultados de um ensaio de degradabilidade in situ. No experimento é avaliado o resíduo potencialmente degradável da fibra em detergente neutro (FDN) da gramínea coastcross (Cynodon dactylon x Cynodon nlemfunensis) cortada aos 60 dias, com três repetições. O resíduo potencialmente degradável da FDN é estudado utilizando quinze tempos de incubação (0; 0,5; 1; 3; 6; 9; 12; 18; 24; 36; 48; 56; 72; 96 e 120 horas). A parcela experimental é constituída por uma vaca não lactante, com fístula ruminal permanente. São obtidos ajustes médios e individuais para os animais. Obtem-se também as variâncias dos estimadores dos parâmetros por meio da matriz de variância e covariância dos parâmetros e pelo método jackknife, propondo-se expressões para a estimação do intervalo de confiança para os parâmetros do modelo. Os resultados mostram que o método de jackknife apresenta maior estimativa de variância para os parâmetros do modelo de Mertens & Loften (1980), resultando em intervalos de confiança de maior amplitude e estimativas dos parâmetros menos precisas, nos ajustes individuais e médios.

Termos para indexação: Regressão não-linear, degradabilidade ruminal, método jackknife, gramínea coastcross.

INTRODUCTION

Quantitatively, carbohydrates are the most important nutrients in the diet of ruminants. Vegetables contain approximately 75% of these elements that are the primary source of energy for the microorganisms of the rumen and for the host animal. The main carbohydrates found in the vegetables are cellulose, hemicelulose, pectins, fructosans, and starch. The microorganisms of the rumen present high capacity of utilization of these carbohydrates, in that the ruminants use the final products of the ruminal fermentation, more specifically the fatty acids that are absorbed, supplying most of the energy demands of the animals.

For a long time, the feeding of ruminants was inadequate, because it was just based on the amount, and not on the quality of the feeds. Today, the characterization of feeds according to their chemical composition, and the constitution of their different fractions, degraded or not, in the rumen, is one of the principal objectives of nutritionists when balancing rations that provide nutrients for the growth and development of the microorganisms in the rumen and, consequently, of the animal.

The in situ technique, used to determine the utilization of the different fractions of the feeds, is used today as reference technique and consists of determining the disappearance of the nutritious component analyzed that remains packed in nylon bag, in different periods of time. This technique was already known in 1930, when Quinn et al. (1938) used this method to evaluate the digestion of the feedstuff in the rumen of fistulated sheeps.

According to Mertens (1993), the first evaluations of digestion processes and conversion, those depended on retention times, were qualitative and based on the visual interpretation of digestion curves. The description of the process was difficult because the digestion curves showed nonlinear behavior, and seemed to not fit to the kinetics of typical chemical reactions. According to Mertens (1993), Waldo, in 1970, was the first to report that degradation profiles were combinations of digestible and indigestible material, and that the potentially digestible fraction would follow kinetics of first order. Mertens (1977) verified that the digestion of the neutral detergent fiber (NDF) presented a period in that degradation of the component indeed didn't happen. Mertens & Loften (1980) suggested the inclusion of a parameter, for the estimates of the parameters of the model of first order of Waldo et al. (1972) that contemplated this period for degradability in vitro and in situ of NDF, dry matter (DM) and nitrogen (N).

Souza (1998) considers the nonlinear regression model in the form y_t,= f(x_t, θ⁰) + ε_t, where t = 1, ..., n; f(x_t, θ⁰) has known functional form; x_t is a k dimensional vector formed by observations in exogenous variables, θ⁰∈ Θ is a p dimensional parameter belonging to the parametric space Θ e ε_t is an experimental error, not directly visible. The same author mentions that, in a similar way to linear models, the process of parameter estimates in a nonlinear model can be obtained by the minimization of the error sum of squares, obtaining the system of normal nonlinear equations, which does not present a closed-form solution for the parameter estimate, and therefore, it must be obtained by iterative processes. The success in the use of Gauss-Newton's algorithm, as an iterative method, will depend on the appropriate choice of both the function response and starting values. Although some general orientations exist for determination of starting values, the choice process is a procedure resolved by the researcher. Several alternatives for the determination of those values are presented in Draper & Smith (1998) and Gallant (1987).

Some important statistical considerations, usually ignored in the studies of the curves of ruminal degradation, are the heterogeneity of variances of the percentage of degradation of the food in time, and the existence of autocorrelation of the regression residuals. Violating the assumptions of homogeneity of variances and no autocorrelation, then the estimates are likely to be biased and the variances of the parameters are subestimated, respectively. (SOUZA, 1998).

In agreement with Hoffman & Vieira (1998), in presence of heterogeneity of variances, the weighted least squares procedure is more appropriate for supplying unbiased estimators with minimum variance, while in the presence of both heterogeneity of variances and autocorrelation of the residuals, the generalized least squares method is more efficient than both weighted and ordinary least squares methods.

Usually, in basic regression models, it is assumed that the errors are not correlated, that is, that they are independent, which is not appropriate when one works with chronological series of data. In such a case, the error corresponding to an observation is correlated with the error of the previous observation (HOFFMANN & VIEIRA, 1998). According to Morettin & Toloi (2004), the general characteristic of the autocorrelation of the residuals is the one of a systematic variation of the values among successive observations. When this happens, it is said that the residuals present autocorrelation.

According to Mood et al. (1974), the deficiency observed in the point estimator is that by itself conveys no information about its closeness to the parameter value and, therefore, supply limited information regarding the parameters. In practice, only rarely the point estimate is exactly equal to the parameter. This situation makes possible the inference to be complemented, whenever possible, with presuppositions concerning the probabilities of θ being close or not to its point estimate. That can be done by the construction of confidence intervals so that, with a chosen degree of confidence (probability) the value of the parameter will be captured inside the interval. Several methods may exist for obtaining those intervals, though, the question is to determine which is the best interval. It is considered the more appropriate the one that presents minimum amplitude.

Vieira (1995) compared three statistical models, among them the one of Mertens & Loften (1980), for estimate of the kinetics of ruminal degradation in vitro and in situ of elephant grass (Pennisetum purpureum Schum. cv. Miner), cut at 61, 82, 103, 124, and 145 days after planting. Also, Feitosa (1999) used the same model to describe the degradation of the dry matter, crude protein, and NDF of coastcross grass hay.

Lira (2000), studying different models and markers for simulation of the digestive kinetics and of passage of braquiaria grass (Brachiaria decumbens Stapf.) made use of the model of Mertens & Loften (1980) to predict the degradation of NDF of the mentioned grass in two seasons (dry and rainy).

The main aim of this work is to estimate the parameters of the degradability model of Mertens & Loften (1980) and to calculate the confidence intervals for the parameters by using both the covariance matrix of the parameters and the jackknife resampling method using data of an experiment of degradability in situ with fistulated cows.

MATERIAL AND METHODS

To illustrate the methodology, experimental data was used (REIS, 2000) of the potentially degradable residue, of the neutral detergent fiber (NDF), of coastcross grass (Cynodon dactylon x Cynodon nlemfunensis), cut at 60 days of age. The degradation profile was evaluated at fifteen incubation times (0; 0.5; 1; 3; 6; 9; 12; 18; 24; 36; 48; 56; 72; 96 and 120 hours), the experimental portion been constituted by a non-lactating cow, with a permanent ruminal fistula.

The model used to describe the potentially degradable residue was of Mertens & Loften (1980), given as:

where R(t) is the residue after the incubation in the rumen in the time t (%); D is the degradable fraction (%); c is the degradation rate (%/h); t is the time of incubation (hours); I is the insoluble fraction and non-degradable (%) and L is the colonization time of the particles or lag time (hours).

This model is considered nonlinear and segmented, therefore the partial derivatives in relation to the parameters, for the larger times than the colonization time (t>L), continue in function of the parameters themselves, as shown by the following equations:

Several iterative methods have been proposed in the literature for obtaining the least squares estimates of the parameters of a nonlinear regression model. One of them is Gauss-Newton's method that makes use of the partial derivatives, and consists of the development of Taylor series to the term of first order of the function f(Xi, θ) around the point θ⁰.

The iterative formula known with Gauss-Newton's method is given by θ¹ = θ0 + (X ' X)^-1X', where X is the matrix of partial derivatives in relation to the parameters. This process is repeated placing θ¹ in the place of θ⁰(vector of starting estimates) until some convergence criterion is accepted. The speed in the convergence depends on the complexity of the model in study and, mainly, on the quality of the starting values, necessary in any iterative method.

The matrix of partial derivatives (X) is of the type _nX_p, being n the number of incubation times of the samples in the rumen (n=15) and p the number of parameters of the model (p=4) and assumes the form X = (a_m1a_m2a_m3a_m4) where the elements of this matrix are vectors column of type 15x1, and m assumes the values of the incubation times, in other words, 0; 0,5; 1; 3; 6; 9; 12; 18; 24; 36; 48; 56; 72; 96, and 120.

If in the vectors column m>L then, the partial derivatives are:

If in the vectors column m>L then, the partial derivatives are:

In this case, X'X is a symmetrical matrix (4×4) and its elements (b_ij) can be written in the following way,

where is the greatest integer function, in other words, = is the greatest integer not superior to x.

Draper & Smith (1998) present the estimates of the asymptotic matrix of covariances in the following manner:

with j = 1,..., p, where p is the number of parameters.

In this manner, the confidence interval is defined for the parameter as , where is the estimate of the j-th parameter, ^t is the superior percentile of the Student's t-distribution with n-p degrees of freedom and is the estimate of the standard error of the estimate of the j-th parameter.

Since there are no closed-form expressions for the variance estimators of the parameters of the model in study, the jackknife resampling method is also used to build the confidence intervals for the parameter estimates of the model of Mertens & Loften (1980), Those intervals are obtained in the following way: first, the estimates of the parameters is obtained, through the SAS program (1995), by using the data set without the first incubation time. Soon after, the process is repeated considering the original data set without the second incubation time and, successively, until all the partial jackknife estimates have been calculated. This process allowed to obtain a group of 15 partial estimates of the parameters. The estimates were applied to the following formula: E^*_j = nE^{* - (n - 1)} , with j = 1,..., 15, where E^*_j is the estimate of the pseudo-value for the parameter when the j-th incubation time was excluded, n is the number of observations, E^* is the estimate of the parameters with all the incubation times and is the value of the partial parameter estimate obtained by the fitted model by using the SAS program (1995) when the j-th incubation time was excluded. Then, the mean was obtained, which represents the jackknife estimator and the variance of this estimator in a conventional way, making possible the construction of the confidence interval as follows: where is the superior percentile of the Student's t-distribution with n-1 degrees of freedom, and S is the standard deviation of the pseudo-values E^*_j.

Parameter estimates were obtained for both the individual curves and the average curve of the model in question. To verify the heterogeneity of variances, the variance of the potentially degradable residues was calculated for each incubation time, making possible the calculation of maximum F of Hartley (PEARSON & HARTLEY, 1970), obtained by the quotient between the maximum and minimum residual variance, in that the non-significance shows that there is no heterogeneity of variances. To verify the presence of residual auto-correlation it was used both the macro %AR (y, p), implemented in the proc model (SAS, 1995) and the correlogram analysis.

RESULTS AND DISCUSSION

Just as in linear regression analysis, it is extremely important in nonlinear regression to have tools that allow goodness-of-fit of the model in a given application. The typical representations of informal diagnosis in nonlinear regression involve the plots of residuals against the independent variables (or ordered sequence) to check for nonconstant variance (heterocedasticity), influential observations and specification errors in the response function.

Figure 1 shows the plot of the residuals against the ordered sequence of the independent variable (incubation times), for the individual fits of the animals. It is noticed that the plot of residuals, for grass cut at 60 days, does not demonstrate evidences of any heterocedasticity pattern. A typical heterocedasticity pattern that can be found is mentioned by Vieira (1995), whereas the evaluation of the model of Mertens & Loften (1980) linearized for logarithmic transformation, tended to present the band of residuals widening to the right (larger times) showing nonconstant variance.

Table 1 shows the values of mean degradation of NDF (in percentage) and their variances. It is observed that, with the increase of the incubation time of the samples in the rumen, there is a decrease in the residual variances of the degradation of NDF of coastcross grass. The relationship between the maximum and the minimum variance is of 1.20 for the cut of the grass at 60 days, not showing the presence of heterogeneity. In other words, this relationship is not significant for 3 experimental groups (3 animals) and 14 degrees of freedom, with a significance level of 5%, according to the Hartley's maximum F test.

The ratio between the maximum and minimum. variance was also used by Mazzini (2001) when studying growth curves of bovines of the Hereford breed. This author tells us that, as the age of the animals increased, there is an increment in the variances of the body weights. This may be also a possible source of heterocedasticity pattern.

It is observed in Table 2 that the autocorrelation parameters (Ø₁ and Ø₂) are not significant at a significance level of 5%, then not being necessary the fit of autoregressive structure errors. Mazzini (2001), Mazzini et al. (2003) and Mazzini et al. (2005), when studying growth curves of bovine through several functions considering autocorrelation of the residues, found that some functions didn't fit to an autoregressive model of 1st or 2nd order when fitting mean curves; however, when fitting individual curves there were animals that presented these error structures.

The parameter estimate values obtained by the jackknife method didn't suffer great alterations. Whereas, the variance estimates of the parameter estimates, obtained by this method presented significant increases in its values, resulting in confidence intervals of greater amplitude, and less precise parameter estimates.

This increase in the variance estimate of the parameter estimates, when compared with the estimates obtained by the covariance matrix of the parameters, can be explained by the fact that when excluding the j-th sample observation, the program converged for values differing of the partial estimates, affecting directly the estimate of the pseudo-value of the parameter directly and its variance.

This methodology was also used by Pereira et al. (2005) when studying the prediction of the nitrogen mineralized in Latossolo through nonlinear models. When comparing confidence intervals obtained by using the covariance matrix of the parameters and by using the jackknife method, the authors also verified that this second method provided greater variance estimates of the parameters.

Table 3

CONCLUSIONS

The model of Mertens & Loften (1980) was adapted to describe the degradability in situ of neutral detergent fiber of coastcross grass (Cynodon dactylon x Cynodon nlemfunensis).

The jackknife method provided greater variance for the parameter estimates of the model of Mertens & Loften (1980), resulting in confidence interval of greater amplitude and less precise parameter estimates in both individual and mean fits.

(Received in april 27, 2006 and approved in april 2, 2007)

• DRAPER, N. R.; SMITH, H. Applied regression analysis. New York: J. Wiley & Sons, 1998. 706p.
• FEITOSA, J. V. Ensaios de degradabilidade in situ: uma abordagem estatística. 1999. 117p. Dissertação (Mestrado em Zootecnia)-Universidade de São Paulo, Jaboticabal, 1999.
• GALLANT, A. R. Nonlinear statistical models. New York: J. Wiley, 1987. 610 p.
• HOFFMANN, R.; VIEIRA, S. Análise de regressão: uma introdução à econometria. São Paulo: HUCITEC, 1998. 379p.
• LIRA, V. M. C. Utilização de diferentes modelos matemáticos e marcadores para simulação da cinética digestiva e de trânsito do capim braquiária (Brachiaria decumbens Stapf.). 2000. 90p. Dissertação (Mestrado em Zootecnia)-Universidade Federal de Viçosa, Viçosa, 2000.
• MAZZINI, A. R. de A. Análise da curva de crescimento de machos Hereford considerando heterogeneidade de variâncias e autocorrelação dos erros. 2001. 94p. Dissertação (Mestrado em Estatística e Experimentação Agropecuária)-Universidade Federal de Lavras, Lavras, 2001.
• MAZZINI, A. R. de A. et al. Análise da curva de crescimento de machos Hereford. Ciência e Agrotecnologia. Lavras, v.27, n.5, p.1105-1112. 2003.
• MAZZINI, A. R. de A. et al. Curva de crescimento de novilhos Hereford: heterocedasticidade e resíduos autoregressivos. Ciência Rural. Santa Maria, v.35, n.2, p.422-427. 2005.
• McDONALD, I. A revised model for the estimation of protein degradability in the rumen. Journal of Agricultural Science, v.96, p.251-252.1981.
• MERTENS, D. R. Dietary fiber components: relatioship to the rate and extent of ruminal digestion. Federation Proc., v.36, n.2, p.187-192, 1977.
• MERTENS, D. R. Rate and extent of digestion. In: FORBES, J. M.; FRANCE, J. (Ed.) Qualitative aspects of ruminant digestion and metabolism Wallingford, UK: Cambridge University, 1993. Cap. 2, p.13-51.
• MERTENS, D. R.; LOFTEN, J. R. the effects of starch on forage fiber digestion kinetics in vitro Journal of Dairy Science, v.63, p.1437-46. 1980.
• MOOD, A. M; GRAYBILL, F. A; BOES, D. C. Introduction to the theory of statistics. Tokyo: McGraw-Hill Kogakusha, 1974. 564 p
• MORETTIN, P. A.; TOLOI, C. M. de C. Previsão de séries temporais. São Paulo: Atual, 2004. 436p.
• PEARSON, E. S.; HARTLEY, H. O. Biometrika tables for statisticians. Cambridge: Cambridge University, 1970. v.1.
• PEREIRA, J. M. et al. Non linear models to predict nitrogen mineralization in an oxisol. Scientia Agrícola, v.62, n.4, p.395-400, 2005.
• QUINN, J. I.; VAN DER WATH, J. G.; MYBURGH, S. Studies on the alimentary tract of merino sheep in South Africa. Description of experimental technique. Journal of Veterinary Science Animal Indian, Ondersteport, v.11, n.2, p.341-360, 1938.
• REIS, S. T. dos. Valor nutricional de gramíneas tropicais em diferentes idades de corte 2000. 99p. Dissertação (Mestrado em Zootecnia)-Universidade Federal de Lavras, Lavras, 2000.
• SAS INSTITUTE Inc. SAS/ETS^® User's Guide. Version 6. 2.ed. Cary, 1995.
• SOUZA, G. da S. Introdução aos modelos de regressão linear e não-linear. Brasília: Embrapa-SPI/Embrapa-SEA, 1998. 489p.
• VIEIRA, R.A.M. Modelos matemáticos para estimativa de parâmetros da cinética de degradação do capim elefante (Pennisetum purpureum, Schum., cv. Mineiro) em diferentes idades de corte 1995. 88p. Dissertação (Mestrado em Zootecnia)-Universidade Federal de Viçosa, Viçosa, 1995.
• WALDO, D. R.; SMITH, L. W.; COX, L. E. Model of cellulose disappearance from the rumen. Journal of Dairy Science, v.55, n.1, 1972.

1

Trabalho conduzido com apoio do CNPq.

Publication Dates

• Publication in this collection
  13 Nov 2008
• Date of issue
  Oct 2008

History

• Accepted
  02 Apr 2007
• Received
  27 Apr 2006