|
|
|
|

*
Dipartimento di Scienze della Produzione Animale, Facoltà di Medicina Veterinaria, Università degli Studi di Udine, 33010 Pagnacco (Ud), Italy and
Federal Research Institute for Agriculture in the Alpine Regions, A 8952, Irdning, Austria
1To whom correspondence should be addressed. E-mail: bruno.stefanon{at}dspa.uniud.it.
| ABSTRACT |
|---|
|
|
|---|
KEY WORDS: modeling neural networks nitrogen excretion purine derivatives cows
| INTRODUCTION |
|---|
|
|
|---|
An introduction to artificial neural networks
An artificial neural network (ANN) is a parallel distributed structure devoted to information processing, consisting of processing elements. As described by Hecht-Nielsen (11
) the processing elements of the network, called nodes, are related to each other by connections; each connection is assigned a weight and each node can receive any number of input connections and can have any number of outgoing connections that fan out the same signal. A generic processing element is provided with a local memory and has a transfer function that can use the local memory and the input signals to generate the output signal. Sarle (12
) has described neural networks as a "bunch of neurons," "simple linear or nonlinear computing elements interconnected in often complex ways and often organized into layers." A particular neural network structure, the multi-layer perceptron (MLP), is discussed in this article. Such a structure is a supervised neural network consisting of three different layers of nodes: input, hidden and output layers. Nodes within the same layer are not connected to each other horizontally but have vertical connections to those of preceeding and following layers (Fig. 1
). The term supervised refers to the type of training scheme adopted during the adaptation phase, when the neural network is exposed to a particular dataset according to a training schedule. In supervised training, the network is supplied with a series of input and output pairs, then the adaptation takes place by comparing the correct output value used for the training process with the model actual output (an estimation of the correct output). As discussed by Sarle (12
), MLP are flexible, nonlinear models able to approximate virtually any function and can be used when there is "little knowledge about the form of the relationship between the independent and the dependent variables"; the complexity of MLP models can be changed "by varying the number of hidden layers and the number of hidden neurons in each hidden layer." If the number of neurons in the hidden layers is small, then the MLP is a parametric model and an alternative to polynomial regression (12
). This new approach can be used for some real-world applications, such as sensor processing, pattern recognition, data analysis and control. Roush and Cravener (13
) used neural networks to predict amino acid contents in feed ingredients, and recently Berg et al. (14
) developed a neural network to predict pork carcass composition from electromagnetic scanning. Neural networks in agriculture would be widely applied if accurate predictions could be obtained using a small number of key variables easy to measure on farms. In this respect, networks with a small number of input variables would be ideal, and their performance, therefore, could challenge those of a traditional multiple regression analysis.
|
The present study is based on the hypothesis that neural network models can be used to study complex physiological aspects like those governing the use of nitrogen in the rumen and the excretion of total nitrogen in the urine of cows.
In this study, a comparison of ANN with traditional linear regression methods has been carried out with the aim of: predicting urine nitrogen excretion using simple parameters measurable at a farm level, predicting the urinary excretion of PD by means of simple parameters affecting microbial synthesis in the rumen, and examining the role of LW and milk yield (MILKY) as possible factors affecting the excretion of PD in cows.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The database used to develop and test the model consisted of 177 nitrogen balances carried out on dry (four experiments) and lactating cows (three experiments) at the Department of Animal Production Science of the University of Udine and at the Federal Institute for Alpine Studies at Irdning (Austria). Methodological and technical details of the trials are reported in the study by Gruber et al. (15
) and in other published articles (16
19
). All of the trials were conducted according to the approved national protocols regarding experiments with animals. The intake of non-structural carbohydrates (CHOINT) was calculated by subtracting crude protein (CP), neutral detergent fiber (NDF), ash and ether extract from dry matter (DM) intake. This gives an estimate of nonstructural carbohydrates (NSC) intake, because the difference between NDF and NDF-bound protein (20
) was not accounted for. In situ rumen degradability trials for CP and non-protein dry matter (NPDM) of the forages and the concentrates fed to the cows were carried out according to the standards described by the Commissione Proteine nella Nutrizione e nellAlimentazione dei Poligastrici (21
), and data were fitted to the model proposed by Ørskov and McDonald (22
). The amount of CP and NPDM effectively degraded in the rumen were computed using the rumen outflow rate for forages and concentrates as predicted by the equations proposed by Owens and Goetsch (23
). The following variables were selected as direct inputs or indirectly used to calculate the input variables: DM intake (DMI; g/d), CP intake (CPINT; g/d), CHOINT (g/d), NDF content of the diet (NDFC; % DM), effective degradability of NPDM (DegNP; g/d) effective degradability of nitrogen (DegCP; g/d), MILKY (kg/d) and LW (kg). The output variables were daily urinary excretion of nitrogen (NURI; g/d) and daily PD nitrogen (PDN) excretion (g/d, calculated as the sum of allantoin and uric acid nitrogen). The original database was split in two for modeling purposes. A first training dataset (n = 130) was used to develop two multiple linear regressions (MLR) aimed at predicting NURI and PDN (24
). The same dataset was then used to obtain two ANN models (25
) devoted to predict NURI and PDN by means of the same variables chosen to develop the regression equations. A second challenge dataset (n = 47) was used to validate the MLR and ANN models. The main dietary, rumen and animal parameters used in the training and in the challenge datasets are presented in Table 1
.
|
The linear model equation to predict NURI was developed using five variables: CPINT, DegNP, NDFC, LW and MILKY. The intercept of the MLR equation did not result significantly different from 0, and, consequently, model with no intercept was recomputed. The second objective of the modeling exercise was to develop an MLR equation to predict PDN with the aim of highlighting the parameters affecting microbial synthesis and the related PDN excretion. A third objective was to investigate the role of LW and milk production level in determining PDN metabolic excretion. For this reason the following variables were adopted as regressors: LW, CHOINT, the ratio between NPDM and CP effectively degraded in the rumen (DNPDCP; g/g), and MILKY per unit of DMI (MILKDM; kg milk/kg DMI). As in the first case, the intercept did not result significantly different from 0, and a model with no intercept was adopted.
Developing an ANN.
Using the same variables as in the linear regression models, the software package NEURAL CONNECTION, Version 2.0 (25
) was used to fit a truly connected MLP model to the training dataset. Fully connected means that every node in each layer is connected to every node in the layer above and below. The structure of the MLP chosen is as follows.
The total input of the jth node of the hidden layer is:
![]() |
where wij is a weight associated with the connection from input layer "i" and hidden layer "j," bj is the bias and xi is the input variable.
The value passed to the single output node from the jth node of the hidden layer is:
![]() |
where
h is the logistic (sigmoidal) transfer function of the hidden layer; note that the use of a logistic transfer function does not imply that there is a direct sigmoidal relationship between input and output layers. Then:
![]() |
The output of the MLP is:
![]() |
where djk is a weight associated with the connection from hidden layer "j" and output layer "k," co is the bias and yj is the passed variable of the hidden layer.
The value passed to the single hidden node to the output layer "k" is:
![]() |
The transfer function
o of the output layer is linear.
The software package NEURAL CONNECTION, Version 2.0 (25
) calculates the bj, wij, co and yjk by means of a back propagation learning rule. This involves repeatedly adjusting the parameters of the model until the error between the predicted values and actual outputs is minimized. The MLP is susceptible to a phenomenon of overtraining, i.e., the model tends to model the noise element of the data as well as the signal. To overcome this, a second training test dataset is used to monitor the predictive capabilities of the model. In the present study the training dataset was randomly split into a training dataset (n = 117, i.e., 90% of the data) and a training test dataset (n = 13, i.e., 10% of the data). The model error was evaluated by the software both on the training and the training test datasets by calculating mean square error between the correct and the actual output values.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
The input variables were included in the linear regression and neural network models because of their potential role in explaining the output variables. CPINT is a direct determinant of NURI as a fraction of dietary protein degraded in the rumen and fermented to produce ammonia escaping in the portal blood (26
). Dietary amino acids absorbed in the gut that are in excess of tissue requirements can be deaminated in the liver and converted to urea. NDFC and DegNP affect rumen environment and nutrient availability for microbial populations. Microbial synthesis and protein fermentation affect ammonia escape from the rumen and, therefore, nitrogen excretion in the urine. The amount of fiber in the ration (NDFC) affects rumen pH (27
). The NDF content of forages has been negatively correlated with CP potential degradability (28
), and a high NDF forage content has been associated with low CP effective degradability in both roughages and concentrates (29
). DegNP can be compared with carbohydrate availability in the rumen, the lack of which can lead to an increased bacterial fermentation of protein (30
) and, consequently, to ammonia production and nitrogen escape in the portal blood. Recent models based on metabolizable protein have assumed that a fraction of the absorbed amino acids supports milk production (31
). For this reason MILKY has been chosen to control the partition of the absorbed protein toward milk production. LW was introduced in the model because metabolic LW affects endogenous urinary nitrogen (32
). The prediction of the urinary excretion of PDN was the second objective of this study. Several researchers have suggested that the metabolic fraction of these compounds is related to animal body mass (6
,8
); therefore, LW was also included in the MLR and ANN models. CHOINT was taken into account because NSC can be rapidly fermented in the rumen (20
). The ratio DNPDCP was considered because microbial synthesis is affected by protein and carbohydrate (or energy) availability in the rumen (33
36
), the latter being represented in the present case by DegNP. Stefanon et al. (10
) suggested that the amount of PD degraded in the tissues can be higher in lactating cows than in dry cows. To isolate the effect of MILKY on metabolic PDN, the ratio between the amount of milk produced per each unit of DMI was computed (MILKDM). A preliminary correlation analysis of the training dataset revealed a strong correlation between MILKY and DMI and that both variables were correlated to PDN excretion. For this reason, the effect of milk production was expressed in terms of a milk:DMI ratio.
Selection of the ANN.
The complexity of a model varies according to the number of hidden layers and the number of neurons in each hidden layer (12
). Increasing the number of nodes in the hidden layer usually improves the performance of the MLP over the training set, however such an advantage is not necessarily reflected on the training test set. Given a certain number of input variables and a certain number of outputs to predict, the neural network architecture is determined by the number of the nodes in the hidden layer. The optimal number of the nodes in the hidden layer is chosen on the basis of the MLP performance measured on the training test set. Adding new nodes in the hidden layer may lead to better results in terms of model prediction ability but on the other hand the model might learn the error (i.e., "the noise") associated with the data. Should this happen, the neural network would lose its predictive power. Some researchers applying Neural Networks have proposed methods to achieve the best model performance; these are concerned either with the best architecture, i.e., the number of nodes to be considered in the hidden layer (37
,38
), or with a strategic interruption of the training phase (39
). In the present study the best model architecture was selected by testing a wide selection of nodes in the single hidden layer of the MLP. A range of one up to fifteen nodes in the hidden layer were studied, monitoring the model performance in terms of the error computed on the training set and the training test set. The error computed on the training test set was used to avoid overtraining as well as to choose the best model architecture. Addition of additional nodes in the hidden layer was evaluated in terms of performance on the training test set. The architecture showing the lowest error on the training test set was then chosen as the best model. The best model architecture for NURI and PDN presented, respectively, four and five nodes in the hidden layer.
Comparison between neural networks and linear equations performances.
The coefficients and statistics of MLR for NURI and PDN are reported in Table 2
. The selected variables always entered in the model with a high level of statistical probability, for both the dependent variables. The parameters of the ANN, reported in Table 3
, refer to the weights of each node of the networks, which can be used to reproduce the ANN model (for details see the ANN: model development paragraph and formulas reported in Fig. 1
).
|
|
|
|
|
The results so far presented suggest some difficulties in the application of neural networks to biological systems. Such systems involve a higher level of complexity than do many physical systems, which can often be successfully modeled by ANN (40
).
Neural networks scenario and model response.
The neural network models obtained for PDN and NURI (training dataset, n = 130) were further evaluated in terms of their response to input variations. The behavior of the output was tested when varying inputs according to sound nutritional and physiological dairy concepts. Such an investigation is available as an option in the NEURAL CONNECTION, Version 2.0 software (25
). Basically, the response of the model is checked by varying an input variable under study against another input variable set to a desired value, while all remaining inputs are kept at their mean value. An arbitrary value between the highest and the lowest value of the training dataset was chosen for each input variable under study. Such a value was increased or decreased by 5% and the corresponding response of the output of the model was recorded. This analysis provides a graphical description of the relationship between a single input variable and the output. Given the large amount of figures needed to show all the input variables effects on NURI and PDN, here only the effect of LW and MILKY on NURI, the effect of CHOINT and DNPDCP on PDN, and the effect on LW and MILKDM on PDN are presented as graphical responses.
The response of the neural network model to MILKY variation is showed in Figure 4A
. The response in terms of NURI to a change in MILKY was studied from a starting value of 17.9 kg/d, against a CPINT fixed at 2291 g/d. A 50% reduction in MILKY generates a 24% increase in NURI excretion (Fig. 4A
), whereas an increase of MILKY from 17.9 kg/d to 26.8 kg/d (+50% MILKY) causes a lower urinary nitrogen excretion of 11.9%. When a fixed amount of protein digested is available, an increased requirement for MILKY causes a lower amount of nitrogen wasted in the urine, a result consistent with current concepts in dairy cow physiology (41
).
|
|
The model also reacted to variations in the ratio of non-protein DM degraded to CP degraded (DNPDCP). A change of -50% and +50% of the basic value of 10.3 of DNPDCP ratio against a fixed value of CHOINT of 4523 g/d, respectively, resulted in a 38.6% increase and in a 24.8% decrease of PD excretion. The ratio between non-protein DM degraded and the CP degraded is a gross index of the balance between energy and nitrogen availability in the rumen. Stokes et al. (43
) reported lower bacterial efficiencies (grams of bacterial N per kilogram of DM digested) as a result of a wider NSC:degradable protein intake ratio in experiments with continuos cultures. This behavior can be confirmed in the present model by the surface response graph in Figure 6
, which shows a reduction of PDN when CHOINT increases if DNPDCP is contemporarily increased.
|
|
| FOOTNOTES |
|---|
3 Abbreviations used: ANN, artificial neural network; CHOINT, non-structural carbohydrates intake; CP, crude protein; CPINT, crude protein intake; DegCP, effective degradability of nitrogen; DegNP, effective degradability of non-protein dry matter; DMI, dry matter intake; DNPDCP, ratio between non-protein dry matter and crude protein effectively degraded in the rumen; LW, live weight; MILKY, milk yield; MILKDM, milk yield per unit of dry matter intake; MLP, multi-layer perceptron; MLR, multiple linear regression; NDF, neutral detergent fiber; NDFC, NDF content of the diet; NPDM, non-protein dry matter; NURI, urinary total nitrogen excretion; PD, purine derivative; PDN, purine derivative nitrogen. ![]()
Manuscript received April 19, 2001. Initial review completed July 23, 2001. Revision accepted September 27, 2001.
| LITERATURE CITED |
|---|
|
|
|---|
1. Wilkerson, V. A., Mertens, D. R. & Casper, D. P. (1997) Prediction of excretion of manure and nitrogen by Holstein dairy cattle. J. Dairy Sci. 80:3193-3204.[Abstract]
2. Jonker, J. S., Kohn, R. A. & Erdman, R. A. (1998) Using milk urea nitrogen to predict nitrogen excretion and utilization efficiency in lactating dairy cows. J. Dairy Sci. 81:2681-2692.[Abstract]
3. Bannink, A., Valk, H. & Van Vuuren, A. M. (1999) Intake and excretion of sodium, potassium, and nitrogen and the effects on urine production by lactating dairy cows. J. Dairy Sci. 82:1008-1018.[Abstract]
4.
Chen, X. B., Subba, D. B.,
rskov, E. R. & Jayasuriya, M.C.N. (1999) Purine Nitrogen Index, potentially a new parameter for rapid feed evaluation in ruminants. Nuclear Based Technologies for Estimating Microbial Protein Supply in Ruminant Livestock 1999Vienna Austria, August 2428, 1998, IAEA-TECDOC-1093 97110..
5. Chen, X. B., DeB Hovell, F. D. & Ørskov, E. R. (1990) Excretion of purine derivatives by ruminants: recycling of allantoin into the rumen via saliva and its fate in the gut. Br. J. Nutr. 63:197-205.[Medline]
6. Chen, X. B., Ørskov, E. R. & DeB Hovell, F. D. (1990) Excretion of purine derivatives by ruminants: endogenous excretion, differences between cattle and sheep. Br. J. Nutr. 63:121-129.[Medline]
7. Zinn, R. A. & Owens, F. N. (1986) A rapid procedure for purine measurements and its use for estimating net ruminal protein synthesis. Can. J. Anim. Sci. 66:157-166.
8. Verbic, J., Chen, X. B., MacLeod, N. A. & Ørskov, E. R. (1990) Excretion of purine derivatives by ruminants: effect of microbial nucleic acid infusion on purine derivative excretion by steers. J. Agri. Sci. Camb. 114:243-248.
9. Chen, X. B., Susmel, P., Stefanon, B. & Ørskov, E. R. (1996) On the use of purine derivatives in spot urine, plasma and milk samples as indicators of microbial protein supply in ruminants. Protein Metabolism and Nutrition 1996Proceedings of the 7th International Symposium on Protein Metabolism and Nutrition, EAAP 81, EZN, Santarem Portugal.
10. Stefanon, B., Mills, C. R., Spanghero, M. & Susmel, P. (1996) An evaluation of purine derivatives as indicators of rumen microbial protein synthesis in dry and lactating cows. Protein Metabolism and Nutrition 1996Proceedings of the 7th International Symposium on Protein Metabolism and Nutrition. EAAP 81, EZN, Santarem, Portugal.
11. Hecht-Nielsen, R. (1991) Neurocomputing 1991 Addison-Wesley Publishing Company, Inc Reading, MA. .
12. Sarle, W. S. (1994) Neural networks and statistical models. Proceedings of the 19th Annual SAS Users Group International Conference 1994 SAS Institute Cary, NC. .
13.
Roush, W. B. & Cravener, T. L. (1997) Artificial neural network prediction of amino acid levels in feed ingredients. Poult. Sci. 76:721-727.
14.
Berg, E. P., Engel, B. A. & Forrest, J. C. (1998) Pork carcass composition derived from a neural network model of electromagnetic scans. J. Anim. Sci. 76:18-22.
15. Gruber, L., Steinwidder, A., Stefanon, B., Steiner, B. & Steinwender, R. (1999) Influence of grassland management in alpine regions and concentrate level on N excretion and milk yield of dairy cows. Livest. Prod. Sci. 61:155-170.
16. Susmel, P., Stefanon, B., Plazzotta, E., Spanghero, M. & Mills, C. R. (1994) The effect of energy and protein intake on the excretion of purine derivatives. J. Agric. Sci. Camb. 123:257-265.
17. Susmel, P., Spanghero, M., Mills, C. R. & Stefanon, B. (1995) Rumen fermentation characteritics and digestibility of diets containing different whey:maize ratios in cows. Anim. Feed Sci. Tech. 53:81-89.
18. Susmel, P., Spanghero, M., Stefanon, B. & Mills, C. R. (1995) Nitrogen balance and partitioning of some nitrogen catabolites in milk and urine of lactaing cows. Liv. Prod. Sci. 44:207-219.
19. Susmel, P., Spanghero, M., Stefanon, B. & Mills, C. R. (1998) Parametri fermentativi ruminali, digeribilità ed escrezione di derivati purinici di bovine alimentate con fonti amilacee diverse. Zoot. Nutr. Anim. 24:3-16.
20. Sniffen, C. J., OConnor, J. D., Van Soest, P. J., Fox, D. G. & Russell, J. B. (1992) A net carbohydrate and protein system for evaluating cattle diets: carbohydrate and protein availability. J. Anim. Sci. 70:3562-3577.[Abstract]
21. Commissione Proteine nella Nutrizione e nellAlimentazione dei Poligastrici (1994) Valutazione degli alimenti di interesse zootecnico: degradabilità e valore proteico degli alimenti per i ruminanti. Zoot. Nutr. Anim. 20:281-291.
22. Oslash;rskov, E. R. & McDonald, I. (1979) The estimation of protein degradability in the rumen from incubation measurements weighted according to rate of passage. J. Agric. Sci. Camb. 92:499-503.
23. Owens, F. N. & Goetsch, A. L. (1986) Digesta passage and microbial protein synthesis. Milligan, L. P. Grovum, W. L. Dobson, A. eds. Control of Digestion and Metabolism in Ruminants 1986:567 Proceedings of the Sixth International Symposium on Ruminant Physiology. Banff, Canada, September 10th-14th, 1984. Reston Book, Prentice-Hall, Englewood Cliffs, NJ..
24. Statistical Package for Social Science. SPSS Advanced Statistics, Version 7.5 1997 SPSS Inc. Chicago, IL .
25. Statistical Package for Social Science. Neural Connection, Version 2.0, Users Guide 1997 SPSS Inc. and Recognition Systems Inc. Chicago, IL .
26. Van Soest, P. J. (1994) Nitrogen metabolism. Nutritional Ecology of the Ruminant 2nd ed. 1994 Cornell University Press Ithaca, NY. .
27. Pitt, R. E., Van Kessel, J. S., Fox, D. G., Pell, A. N., Barry, M. C. & Van Soest, P. J. (1996) Prediction of ruminal volatile fatty acids and pH within the net carbohydrate and protein system. J. Anim. Sci. 74:226-244.[Abstract]
28. Lavrencic, A., Stefanon, B. & Susmel, P. (1997) An evaluation of the Gompertz model in degradability studies of forage chemical components. Anim. Sci. 64:423-431.
29. Susmel, P., Spanghero, M., Stefanon, B., Mills, C. R. & Cargnelutti, C. (1991) Effect of NDF concentration and physical form of fescue hay on rumen degradability, intake and rumen turn-over of cows. Anim. Prod. 53:305-313.
30. Russell, J. B., Sniffen, C. J. & Van Soest, P. J. (1983) Effect of carbohydrate limitation on degradation and utilization of casein by mixed rumen bacteria. J. Dairy Sci. 66:763-775.
31. Baldwin, R. L. (1995) Protein and amino acid models. Modeling Ruminant Digestion and Metabolism 1995:578 Chapman & Hall London, UK. .
32. Oslash;rskov, E. R. & Miller, E. L. (1988) Protein evaluation in ruminants. Ørskov, E. R. eds. Feed Science 1988 Elsevier Publishers B.V Amsterdam. .
33. Institute National de la Recherche Agronomique. (1988) Alimentation des Bovins, Ovins et Caprins (Jarrige, R., ed.), pp. 476. INRA, Paris, France. .
34.
Nocek, J. E. & Russell, J. B. (1988) Protein and energy as an integrated system: relationship of ruminal protein and carbohydrate availability to microbial synthesis and milk production. J. Dairy Sci. 71:2070-2107.
35. Agricultural and Food Research Council Technical Committee on Responses to Nutrients. Report 9: Nutritive Requirements of Ruminant Animals: Protein. Nutrition Abstracts and Reviews, Series B. Livestock Feeds and Feeding 62:787-835.
36. Van Kessel, J. S. & Russell, J. B. (1996) The effect of amino nitrogen on the energetics of ruminal bacteria and its impact on energy spilling. J. Dairy Sci. 79:1237-1243.[Abstract]
37. Anders, U. & Korn, O. (1999) Model selection in neural networks. Neural Networks 12:309-323.[Medline]
38. Lehtokangas, M. (1999) Modeling with constructive backpropagation. Neural Networks 12:707-716.[Medline]
39. Prechelt, L. (1998) Automatic early stopping using cross validation: quantifying the criteria. Neural Networks 11:761-767.[Medline]
40. Wastney, M. E., Patterson, B. H., Linares, O. A., Greif, P. C. & Boston, R. C. (1999) What is modeling?. Investigating Biological Systems Using Modeling: Strategies and Software 1999:382 Academic Press San Diego, CA. .
41. Oldham, J. D. & Emmans, G. C. (1988) Prediction of responses to protein and energy yielding nutrients. Garnsworthy, P. C. eds. Nutrition and Lactation in the Dairy Cow 1988:429 Butterworths London, UK. .
42. Tamminga, S., Van Vuuren, A. M., Van Der Koelen, C. J., Ketelaar, R. S. & Van Der Togt, P. L. (1990) Ruminal behaviour of structural carbohydrates, non-structural carbohydrates and crude protein from concentrate ingredients in dairy cows. Neth. J. Agric. Sci. 38:513-526.
43. Stokes, S. R., Hoover, W. H., Miller, T. K. & Manski, R. P. (1991) Impact of carbohydrate and protein levels on bacterial metabolism in continuous culture. J. Dairy Sci. 74:860-870.[Abstract]
44. Isaacson, H. R., Hinds, F. C., Bryant, M. P. & Owens, F. N. (1975) Efficiency of energy utilization by mixed rumen bacteria in continuous culture. J. Dairy Sci. 58:1645-1659.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||