Evaluation of equations to predict dry matter intake of dairy heifers

Daily pen dry matter intakes (DMI, n = 9,275) were collected over a 28-mo period at the University of Wisconsin's Integrated Dairy Research Facility. Heifers were housed in pens containing 8 Holstein or Holstein × Jersey crossbred heifers/pen. Heifer diets were formulated to energy and protein requirement twice monthly, with feed intake, dietary nutrient density, and ambient temperature recorded daily. Heifers were weighed at 60-d intervals, and mean pen body weights (BW) were estimated for each day between the weigh dates using the interval average daily gain as a regression coefficient. Prediction of heifer DMI was evaluated using the equations of NRC (2001), Quigley et al. (1986), or alternative random effects mixed models or nonlinear exponential models. The effects of breed, BW, temperature and neutral detergent fiber deviation (NDFdv) were considered as independent variables. Holstein and crossbred heifer DMI was predicted with reasonable precision standard error (SE) < 0.86 kg/d], by the NRC (2001) or Quigley et al. (1986) equations, but heifer DMI was over- or underpredicted for heifers >500 kg, respectively. Improved heifer DMI prediction equations were achieved with exponential models. For Holsteins (SE = 0.71 kg/d), the prediction equation was: DMI (kg/d) = 15.79 × 1 - e^{(−0.00210 × BW)}] − 0.0820 × NDFdv, where NDFdv = (dietary neutral detergent fiber as a % of dry matter) - {22.07 + 0.08714 × BW] - 0.00007383 × (BW)²]}. For crossbred heifers (SE = 0.60 kg/d), the prediction equation was: DMI (kg/d) = 13.48 × 1 - e^{(−0.00271 × BW)}] - 0.0824 × NDFdv where NDFdv = (dietary neutral detergent fiber as a % of dry matter) - {23.11 + 0.07968 × BW] - 0.00006252 × (BW)²]}. Alternative exponential DMI model equations when dietary neutral detergent fiber is unknown were also developed. The Holstein DMI equation (SE = 0.73 kg/d) was: DMI (kg/d) = 15.36 × 1 - e^{(−0.00220 × BW)}], and the crossbred DMI equation (SE = 0.81 kg/d) was: DMI (kg/d) = 12.91 × 1 - e^{(−0.00295 × BW)}].