Concentrate feeding and milk yield based on field data of milk recorded herds

Field data from 16 051 Finnish milk recorded herds including milk yield (MY), feed consumption, feed analyses, and the herd effect for milk yield (HMILK) obtained from the national breeding value estimation program, were analysed to detect the influence of concentrate feeding on milk production. HMILKs are deviations from the average national level with mean of 45 kg and SD of 722 kg. Mean MY was 6917 kg and mean dry matter intake (DMI) 5679 kg per cow per year. The effect of concentrate feeding on HMILK and MY was studied by using quantitative [amount of energy (FUI) and concentrates (CI) in the diet] and qualitative [proportion of grain (Gc) or compound feed (Mc) in concentrates or CP content (CPc) of concentrates] diet characteristics as dependent variables in multivariate regression analysis. The general linear effect of CI was 1.18 kg MY/kg CI. Production response of CI decreased with increasing CI as indicated by significant interactions between CI and CI classes. Gc showed a negative relationship with HMILK, but CPc proved to be a more important factor affecting HMILK. Feeding grain instead of compound feed was connected with too low protein content in concentrates. Mc was positively correlated with CP content of concentrates. However, the use of compound feed appeared to give a slight increase in HMILK even after accounting for the effect of CP.


Introduction
Typically feeding intensity in dairy cattle is closely related to concentrate feeding.It could be hypothesised that an increase in amount of concentrates would lead into an increase in energy intake and therefore increased milk production.In short term feeding experiments increase in concentrate intake have resulted in smaller responses in milk yield than in longer complete lactation experiments (Wiktorsson 1979).It has been suggested that the effects of energy input are not fully realised short term.The carry over effect of, for example, low level of concentrate feeding, might not influence milk production during short term experiments but will affect subsequent lactation.To fully assess the impact of energy feeding, research has to be focused on whole lactation records, or to data measured across several lactations.Moreover, marginal responses to incremental concentrate feeding should be measured at different levels of input because productive responses usually respond according to the law of diminishing returns when nutrient supply varies around requirements.
In Finland 60% of dairy herds participate in milk recording and about 73% of total milk production comes from milk recorded herds.Most of these herds also estimate their feed consumption.The feed consumption database, which is linked with the milk recording and feed analysis databases, could provide a valuable view to practical on-farm feeding.The ration formulation program used on farms by advisory services is based on Finnish feeding recommendations (Tuori et al. 1995) and feed allowances are based on milk production.The program tends to divide home-grown forages evenly to cows within fill limits and satisfies energy and protein requirements with concentrates.Generally the amounts of forage offered to individual cows does not depend on milk yield and leads to large variations in the forage to concentrate ratio between cows within a herd.
Estimated breeding values of dairy cows are based on 305 day production.They are estimated using a statistical model (animal model) which in Finland includes the most important environmental effects such as herd-year, interaction between calving year and calving month, and calving age by days open effects within lactation number (Lidauer and Mäntysaari 1996).The largest variation in records is due to herd effects.In the evaluation model, the herd-year classification groups animals to subgroups of contemporaries that have produced under equal conditions.Thus, solutions of contemporary groups will express the effect of feeding and management of the herd.The herd-year effects (later on called herd effects) are estimated separately for first lactation and second and third lactations combined.
In the current investigation milk recording data was connected with solutions for herd ef-fect from the national breeding value prediction program.The same approach on a much smaller scale, was used by Agabriel et al. (1993).They found that the herd effects for milk composition, estimated by animal model on seventy-six French dairy farms, was related to feed characteristics such as type of concentrates and forage quality.
The objective of this paper is to describe the relationship between concentrate feeding and milk yield and animal model herd effects of milk yield based on data collected from Finnish dairy farms in 1993.The work is a part of a larger study that targets the development of farm management tools utilising information from monthly milk recording data.

Data
Data of milk production, feed consumption and feed analyses were obtained from milk recording databases.The total number of herds participating milk recording scheme in 1993 was 20 018, and after edits the total number of herds was reduced to 16 051 with all essential data.Average number of cows per herd was 13.7 (SD 5.7).Three quarters of herds had between 8 and 19 cows and only 3% had more than 25 cows.
Milk yield is measured monthly according to the milk recording scheme.Samples for fat, protein and lactose analyses are taken bimonthly (Karjantarkkailutulokset 1993).Feed consumption is monitored on-farm by measuring feed stores in autumn and making feeding plans.Use of different feedstuffs for dairy cows is registered and reported three times per year; at the beginning, middle and at the end of the indoor feeding season.
Crude protein (CP) content was determined using the Kjehldal method from 4014 grain and 2145 hay samples by Viljavuuspalvelu (Mikkeli, Finland).Feed analyses of silages were made by regional laboratories of Valio Ltd. (Helsinki, Vol. 7 (1998): 423-436. Finland).A total of 35 637 silage samples (from 11 027 farms) had been analysed for dry matter (DM) content.Near infrared reflectance (NIR) spectroscopy was used for determining digestible organic matter in the DM (D-value ) (Hellämäki and Moisio 1991) and CP (Hellämäki and Moisio 1983) content of silage samples.Analysed composition was used in the calculation of feed values according to the Finnish feed evaluation system (Tuori et al. 1995).For those farms who had not sent feed samples for analyses feed table values were used.
Feeds were grouped into forages and concentrates.Forages included silage, hay, pasture and other roughages, mainly straw.Concentrates included grain, fibrous by-products (e.g.sugar beet pulp, wheat bran, barley fibre), compound feeds (commercial concentrate mixtures), protein supplements (rape seed meal, soya bean meal, compound feeds with CP content above 170 g/kg), minerals and other concentrates.Total feed consumption of dairy cows in a herd was divided by the number of cows to get the average intake per cow per year.Intake from pasture was estimated by subtracting the registered amount of feeds other than pasture from calculated total energy requirements based on milk production.
Herd effects for the year 1993 for cows in their 2. and 3. lactations were obtained from the national breeding value prediction program run in 1995.Evaluated traits were milk yield, protein and fat content and days open, somatic cell count and live weight.Only herd effects for milk yield (HMILK) are discussed here.

Statistical analysis
For preliminary examination of the data, correlations between milk production and feeding variables were calculated.Next a stepwise regression analysis was made using the selection based on the coefficient of determination (R 2 ).All feeding parameters were used in the stepwise regression analysis.Based on these results, the most important factors were included in analysis of variance and covariance.Factors consid-ered in latter models were total intakes (kg/cow/ year) of dry matter (DMI), crude protein (CPI), forages (FI), concentrates (CI), and energy (feed units/cow/year, FUI); content of CP in the total diet (CPd), and in concentrates (CPc), g/kg; proportion of concentrates in total diet (Cd).Factors above were divided into 5-7 classes, depending on the effect, to allow and detect non-linear relationships.
When the quantitative effects of concentrate feeding were studied, average herd annual milk yield (MY) was used as dependent variable.CI was the most important factor to explain variation in MY (R 2 = 0.36), as could be expected to result from feeding according to recommendations.Inclusion of other factors in the model after CI gave only a small increase in R 2 .The effect of the amount of concentrates offered to cows was studied by calculating the coefficients of regression of MY on CI within different concentrate intake groups.The model was: where µ is intercept, α 1 i is the effect of general level of CI (α 1 1 :<1500, α 1 2 :1500 -1800, α 1 3 :1800-2100, α 1 4 :2100-2400, α 1 5 : >2400 kg/ year per cow), b 1 is the general regression of CI ij across CI classes and b 1 i is the corresponding regression specific to CI class i, where subscript refers to model [1.0].Later CPc was included into [1.0] as a continuous variable (model [1.1]).
The criteria for finding the base model for HMILK was to get the highest possible R 2 and the lowest possible residual SD with a reservation that factors in the model would not be strongly correlated with each other.Bearing these criteria in mind the following factors were selected to the base model: FUI, CPd and CI.Thus, the regression model was where µ is intercept, R i is geographical region, α 2 i is the effect of general level of CI (α 2 1 :<1500, α 2 2 :1500-1800, α 2 3 :1800-2100, α 2 4 :2100-2400, α 2 5 : >2400, kg/year per cow), b 2 i where d 2 i is linear regression for Gc, g 2 i is linear regression for Mc, and h 2i is a regression for CPc while the other effects are the same as in [2.0].
To study the nonlinearity of the effect of CPc the following nested models were used: where α 3 i is a effect of general level of CPc (α 3 1 :100-140, α 3 2 :140-160, α 3 3 :160-180, α 3 4 :180-220 g/kg), b 3 is a general regression of CPc ij across CPc classes and b 3 i is the corresponding regression specific to CPc class i and c 3 i is a general regression of CI ij .

Results
Mean herd MY was 6719 kg (Table 1).In 95% of the herds, MY was between 5231 kg and 8310 kg.The HMILKs had a range of 6189 kg.Mean DMI was 5679 kg per cow per year.The proportion of forage in DMI was 670 g/kg.Contents of CP in the whole diet and in concentrates were 150 and 156 g/kg DM, respectively.Grass silage was the main forage as only a proportion of farms (10%) used hay more than silage.Herds using silage were on average larger than those having hay as the main forage (14.0 vs 10.2 cows per herd).Most of the herds using silage fed hay in amounts recommended by feeding advisers, i.e. 2 kg/day/cow.There were 795 herds which did not use hay.
The parameters most positively correlated with MY and HMILK were CPI, CI and FUI (Table 2).On the other hand these were negatively correlated with the content of grain in concentrates.
The general linear effect of CI on MY was 1.18 kg (P<0.001)(Table 3).When CPc was included in the model [1.0] the response decreased to 1.14 kg MY (P<0.001).The production response of CI decreased with increasing CI as indicated by a significant (P<0.001)interaction between linear effect of CI and CI group effects.
The means of feeding factors in different CI groups (Table 3) show that increasing the amount of concentrates, also increased the intake of compound feed and grain.The proportion of compound feed in the concentrates increased as CI increased whereas grain content slightly decreased.Silage CP content increased slightly, but silage D-value changed very little.Concentra-tion of urea in milk increased as inclusion level of CI increased but the mean values were all within limits (20-30 mg/100 ml) recommended by feeding advisers in Finland.Supply of energy in relation to cow requirements (S/R ratio) increased with CI, but the marginal response to increased supply of MJ of metabolizable energy (ME) did not change.
Coefficients of the base model for HMILK and the solutions for terms for composition of concentrates are shown in Table 4 When Gc was added in model [2.0] as a con-  Responses to increases in CP content of concentrates on HMILK and the corresponding averages of feeding variables in different CPc groups are shown in Table 5.The linear effect of CPc was 11.7 kg HMILK/g CPc (P<0.001) and the response of CPc decreased with increasing CPc as indicated by a significant (P<0.001)interaction between the linear effect of CPc and CPc group effects.Data were restricted to include observations with CPc values only in the range of 100 to 220 g/kg.From low CPc to high CPc, CI increased about 400 kg and a marked change in the composition of concentrates was observed.The content of compound feeds in concentrates increased from 63 to 578 g/kg and grain content decreased from 844 to 413 g/kg.The CP content of silage decreased, but the D-value of silage changed only slightly with increasing level of CPc.Neither the S/R ratio nor the response of MY to increases in ME intake changed across CPc levels.Milk production per kg concentrate intake increased slightly.

Discussion
In nutrition research analyses of field data are rare, thus we could not find any other studies of comparable size to the one presented here.Analysis of field data for nutritional means is difficult because of colinearity of factors, confounding effects due to feeding according to recommendations, and also because of the accuracy of estimation of feed intake is questionable.The present study focused on the effects of amount of concentrates on milk production because it was assumed that the intake of concentrates and their composition can be estimated more accurately than the intake of forages.Management of dairy herd includes various decisions from forage harvesting to choice of protein supplement.All these require expertise from the farmer.Many management practices appear to be associated with concentrate feeding.The effect of quantity of concentrates on milk production The HMILKs are deviations from national average level with a mean near zero (Table 1).A positive HMILK means that the herd is producing more milk than expected on the basis of cows' genetic potential, parity, calving season, calving age and days open.Thus, HMILK expresses the effects of herd feeding and management policies on milk production.As expected the correlation between HMILK and MY was high (Table 2).
Both HMILK and MY were closely associated with the total intake of concentrates (Table 2).This confirms the realisation of feeding recommendations: when cows produce higher yields the amount of concentrates fed will be increased.The majority of herds participating in the milk recording scheme follow feeding recommendations and feeding intensity is most often increased by increasing the amount of concentrates.Usually the amount of forages offered to individual cows does not depend on milk yield which leads into a substantial variation in the forage to concentrate ratio between cows within a herd.
As HMILKs are deviations that do not contain the herds' actual production level, they are not relevant for models using CI or other quantitative variables.For this reason MY was used as dependent variable in these analyses.The linear effect of one kg of concentrate DM was 1.18 kg MY when the general level of herds' concentrate intake varied from under 1500 to over 2400 kg /cow per year (Table 3).This was larger than the response of 0.79 kg milk per kg additional concentrate reported by Agnew et al. (1996) or 0.54 reported in a review of recent Finnish studies (Huhtanen 1998).Rinne et al. (1995) found a response of 0.51 when concentrate intake increased from 6.2 to 8.7 kg per day.All observed responses are much lower than the theoretical value of 2.38 (average energy content of 1 kg concentrates in this data divided by energy requirement for 1 kg milk).
Production responses of milk yield to addi-tional concentrate intake have been smaller in short term feeding experiments than in whole or multiple lactation experiments (Wiktorsson 1979).Residual effects have been variable depending on length of experiment, level of concentrate inclusion and forage feeding (Gordon 1984).At low levels of concentrate feeding, residual responses have been large and positive but at the highest inclusion level they have been negative.Effects might not be fully visible in milk production during an experiment, but they will affect the post-treatment period or subsequent lactation.The present data consists of annual milk production records of whole herds and thus is similar to whole or multiple lactation experiments.Production responses calculated within CI classes diminished from 1.40 to 0.91 when the general level of concentrate intake increased from levels of 1500 kg or less to levels between 2100 and 2400 (Table 3).In the highest CI class (above 2400 kg concentrate per year) this coefficient was slightly higher indicating a curvi-linear response.However, one can question the accuracy of data in the lowest and highest CI groups which include borderline observations.There were outliers at both ends of the distribution although their number was minor compared to the size of the data set.Reduced response with increasing CI has also been shown in feeding experiments.In the literature review of Huhtanen (1998), changes in milk production due to additional concentrate intake (kg/kg DM) were 0.94, 0.80 and 0.64 when concentrate intakes were less than 5.0 kg, between 5.0-7.3kg or over 7.3 kg per day, respectively.Also in their multiple lactation experiment, Spiekers et al. (1991) found a decreasing production response to increases in concentrate intake.
The ratio of supply of energy to requirements (S/R ratio) increased along the intake of concentrates (Table 3).This could, firstly be due to overestimation of energy intake, secondly due to poor utilisation of increased intake of ME or thirdly due to a biased estimation of feed intake.The most apparent reason seems to be an overestimation of the increase in energy intake with increasing CI, due to negative associative effects Kaustell, K. et al.Concentrate feeding and milk yield of herds between forages and concentrates.These effects become important when the intake of high producing cows is limited by physical factors and therefore large amounts of concentrate supplements are needed to meet energy requirements (Huhtanen 1991).Increasing the level of concentrates in the diet results in a reduced rate of cell wall digestion, such that calculated increases in digestibility and metabolisability are not achieved.
Feeding level appears to have great influence upon negative associative effects (Mould 1988), such that these can be large at high levels of feeding.When considerable amounts of highly soluble carbohydrates from concentrates are included in the diet, digestibility of fibre is depressed.At a maintenance level of feeding, prolonged rumen retention time is able to compensate for the reduced rate of digestion, so that fibre digestibility will not decrease.In contrast, at high feeding levels, such as high producing dairy cows, digesta retention decreases more with high concentrate than high forage diets.Therefore a slower rate of digestion with high concentrate diets can not be compensated for by longer rumen retention time.Consequently, the depressive effect of feeding level on diet digestibility increases with the proportion of concentrate in the diet.
According to the Finnish feed tables (Tuori et al. 1995) diet digestibility is calculated to increase by 6 g per kg DM increase in the amount of concentrate whereas the observed change was -1.8 g per kg concentrate DM in Finnish studies (Huhtanen 1998).This difference means that the increase in ME intake was 4 MJ smaller than expected when the amount of concentrate was increased by 1 kg DM/d.From calculated increases of energy intake only about 70% is typically realised.When forage prepared from high quality grass is fed, diet digestibility is not dependent on the proportion of concentrates, i.e. energy content of the diet does not increase by increasing the proportion of concentrates.In the experiment of Agnew et al. (1996) total ration digestibility was unaffected by concentrate inclusion level which varied between 2 and 8 kg/ day.Also Rinne et al. (1995) found no change in the total ration digestibility of organic matter when concentrate allowance was increased from 7 to 10 kg/day.
The second reason for increased S/R ratio along with increased CI is poor utilisation of additional ME intake.This might not be very a obvious reason.In theory utilisation of ME should increase with CI because the proportion of ME in gross energy (GE) increases and efficiency of utilisation of ME in lactation (k L ) should improve (ARC 1980).The third explanation covers the biased estimation of feed intake or forage energy value.This is not likely because the level of concentrate can not affect the intake of forage or feed values in a systematically biased way.Feed intake in the present data was actually the amount of feed given to cows, not the amount of eaten as refusals were not registered.The amount of refusal could have been larger at higher levels of concentrate allowance.
Larger responses to concentrate feeding in the present data compared to those obtained in feeding experiments are mainly due to feeding according to recommendations but there may also be other reasons such as restricted forage feeding.In feeding trials with restricted forage intake, production responses have been greater than in studies with ad libitum feeding (Johnson 1986).In the present study the higher response may partly be due to a limited intake of forage.Differences between responses to concentrates on restricted and ad libitum forage diets are due to substitution of concentrate for forage in the latter.The mean substitution rate (decrease in forage intake, kg DM / increase in concentrate intake, kg DM) in the present data (Table 3) was much smaller than typical values of around 0.4-0.6 reported in feeding experiments and appears to be responsible for the high production responses.This is in agreement with the response of ME which was not changed with increasing concentrate intake (Table 3).The response of ME was surprisingly similar to the results of feeding experiments conducted in our institute (Rinne et al. 1995, Huhtanen 1998).
age feeding, increasing the amount of concentrates reduced the intake of forage with a substitution rate of around 0.50 (Spiekers et al. 1991, Aston et al. 1994band 1995).Rinne et al. (1995) reported a substitution rate of 0.50-0.69 on grass silage based diets cut at different stages of growth.In the experiment of Aston et al. (1994b), substitution rate was smaller with grass silages of high digestibility compared to silage of low digestibility.
In a whole indoor season experiment, Gordon (1984) reported a linear effect of level of concentrate supplementation on the intake of silage DM as described by the equation y=1939-0.26x,where y = intake of silage, kg DM and x = intake of concentrates, kg DM.One possible factor behind a smaller substitution rate of total forage in the present data could be the restricted feeding of silage, which when conducted could cause average silage intakes to increase slightly with increasing CI (Table 3).This shows that feeding plans based on feeding recommendations have been followed in practice, as the ration formulation program used by advisers tends to restrict the allowance of forage for high yielding cows.

Base model and composition of concentrates
The effect of FUI on the HMILK was unexpectedly low in the present study.The theoretical value is as high as 2.27 kg milk per FU according to Finnish feeding recommendations (Tuori et al. 1995).In feeding experiments responses of around 1 kg are often achieved, for example in the experiment of Rinne et al. (1995) it was 1.3 kg milk per FU.The main reason for the small response in the present field data is that HMILK has been corrected for systematic environmental effects but the allowance of feeds is based on actual non-corrected production of milk.Thus quantitative feeding factors do not describe variation in HMILK particularly well.The FUI within CI class consists of both concentrates and forages and changes in forage intake affect the response.
The composition of concentrates changed when CI increased.Based on correlation coefficients, Gc decreased (Table 2) and Mc increased.Also the CP content of concentrates increased with increasing concentrate intake.A high proportion of grain was associated with a low content of CP in concentrates and in the total diet, and with a low CP intake.
The negative effect of Gc on HMILK (Table 4) was mediated by a low content of CP in the concentrate and the total diet, but the grain content as such did not have negative influence.This suggests that there is a shortage of protein in diets where concentrate feeding is based on cereal grains.One reason for insufficient use of protein supplements in practice could be the limits of milk urea content established by feeding advisory authorities.Avoiding high milk urea content might lead to minimal use of protein supplements, although in many cases the negative effect of a high urea content cannot be noticed.Shingfield et al. (1997) criticised, on the basis of evaluation of milk recording data, the current assessment of the upper milk urea limit in Finland (30 mg urea /100 ml milk) because performance of cows in herds with a high urea content was actually increased with little or no adverse effects on reproductive performance.Negative effects of high urea content are more obvious at levels much higher than current Finnish recommendations.
Commercial compound feeds differ from grain by their more complex composition and higher CP content which may explain the positive effect on HMILK.They are composed of a variety of materials besides grain such as byproducts of the food and alcohol industry and they include many kinds of protein and carbohydrate sources.However, although significant, the quantitative effect of proportion of compound feeds in concentrates was relatively small, i.e.HMILK increased 170 kg/year when Mc increased from 0 to 1000 g/kg.In experiments of Huhtanen (1987Huhtanen ( , 1991) ) and Huhtanen et al. (1988) feeding concentrates which were comprised of different types of carbohydrates were compared with feeding the corresponding ingre-Kaustell, K. et al.Concentrate feeding and milk yield of herds dients alone.More complex mixtures of concentrates gave slightly higher milk yields.One effect of replacing starch in concentrates by digestible fibre seems to be an increased silage intake (Aston et al. 1994a).
The overall response of CPc was 11.7 kg HMILK/g CP in kg concentrate DM when the average CP content of concentrates increased from 100 to 220 g/kg DM (Table 5).When this annual response is converted to correspond to daily production, it appears little larger than responses reported by Aston et al. (1994a).In addition to different carbohydrate sources they also fed four different levels of concentrate CP from 120 to 240 g/kg, which resulted in a mean production response of 0.028 kg milk per additional g/kg CP in concentrates.They concluded that when cows were given silage ad libitum and 9 kg concentrates per day, milk production was more affected by CP content than source of carbohydrate.In the present data the importance of sufficient CP intake can also be seen.
Production response of CPc decreased when CPc level increased (Table 5).In the highest CPc level group the response was even negative.The response to CPc was lower when CI was included in the model, compared to that of CP content of concentrates alone.While cows with higher milk yield are fed with higher amounts of concentrates, CP content of concentrates also increased.The more a cow produces the more concentrates and concentrates richer in protein are used.Average silage CP content was lower when CPc was higher (Table 5).This indicates that the ration formulation scheme which in 1993 was based on digestible crude protein, is widely adopted on Finnish farms.
An increase in concentrate CP content has usually caused an increase in silage intake and milk yield (Aston et al. 1994a, Sutton et al. 1994).Replacement of 1.15 kg concentrate by rape seed meal induced a 0.69 kg increase in daily silage DM intake in the experiment of Rinne et al. (1995).In the present data the average intake of silage was increased by up to concentrate CP content of 160-180 g CP per kg DM.

Conclusions
The response to increased amount of concentrates in the diet was greater than generally observed in feeding experiments but much smaller than that derived theoretically.A larger response in field data is obviously related to the strict application of feeding recommendations indicated by a much smaller substitution rate than observed in feeding experiments.As a result of only small decreases in forage DMI with increased amounts of concentrates in the diet, the marginal responses to incremental ME were similar to values reported in feeding experiments.
The proportion of grain in concentrates showed a negative relationship to HMILK, but CP content of concentrates proved to be more important.Feeding grain instead of compound feed was typically associated with too low a protein content in concentrates.The proportion of compound feed was positively correlated with CP content of concentrates.However, the use of compound feed seemed to give slight increase in HMILK even after taking into account the effect of CP.
. The overall response to FUI was 0.34 kg HMILK per FU (model [2.0]).It was significantly different (P<0.001) in different CI classes, being higher in the lowest and highest classes than in the middle classes.Neither the linear response of FUI nor the interaction between FUI and CI classes changed when CPc, Gc or Mc were added in the model.The overall response of CPd was on average 15.2 kg HMILK per g CPd (model [2.0]).It differed significantly (P<0.01) between CI classes with a positive trend along increasing CI classes.The interaction between CPd and CI classes lost its significance when Mc was added in the model although the linear effect of CPd remained significant.

Table 1 .
Description of the data.

Table 3 .
Milk yield (MY) response to concentrate intake (CI) and characterisation of feeding factors in different concentrate intake groups.
1 Model with CI intake group and linear regression of CI within CI intake groups. 2 General regression across all CI groups. 3Within CI group regression. 4Model [1.1] plus regression effect on concentrate CP.CP = crude protein; S/R =supply of energy in relation to the requirements of the cow; D-value = content of digestible organic matter in DM;

Table 4 .
Regression coefficients from the base model for herd effect of milk by animal model (HMILK) and the solutions for refined models with the composition of concentrates.

Table 5 .
Effect of concentrate crude protein content (CPc) on herd milk solutions (HMILK) and characterisation of feeding factors in different CPc groups.
1 HMILK ij = µ+ CPc-group + b CPc ij + b i CPc ij . 2 Estimate for b. 3 Within CPc-group estimates of b+b ij . 4Model [3.0] + c i CI ij .CP = crude protein; D-value = content of digestible organic matter in DM; S/R =supply of energy in relation to cow requirements; MY = milk yield, kg/cow per year