Next Article in Journal
In Vitro Evaluation of the Adsorption Efficacy of Biochar Materials on Aflatoxin B1, Ochratoxin A, and Zearalenone
Previous Article in Journal
Canine Pyometra: A Short Review of Current Advances
Previous Article in Special Issue
Data Analytics of Broiler Growth Dynamics and Feed Conversion Ratio of Broilers Raised to 35 d under Commercial Tropical Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Environmental and Farm-Associated Factors on Live Performance Parameters of Broilers Raised under Commercial Tropical Conditions

by
Gustavo A. Quintana-Ospina
1,2,
Maria C. Alfaro-Wisaquillo
1,
Edgar O. Oviedo-Rondon
1,*,
Juan R. Ruiz-Ramirez
2,
Luis C. Bernal-Arango
2 and
Gustavo D. Martinez-Bernal
2
1
Prestage Department of Poultry Science, North Carolina State University, Raleigh, NC 27695-7608, USA
2
Grupo BIOS Inc., Envigado 055420, Antioquia, Colombia
*
Author to whom correspondence should be addressed.
Animals 2023, 13(21), 3312; https://doi.org/10.3390/ani13213312
Submission received: 28 July 2023 / Revised: 2 October 2023 / Accepted: 16 October 2023 / Published: 25 October 2023
(This article belongs to the Special Issue The Application of Biostatistics in Animals)

Abstract

:

Simple Summary

In commercial poultry operations, chickens are subjected to different housing, environmental, and management conditions, which can drive the performance results at slaughter age. Although modern poultry housing might incorporate environmental electronic sensors, the data collected are not often analyzed to determine the impact of the environment on broiler live performance results. On the other hand, the variability observed in performance can depend on complex interactions between housing characteristics and management that are difficult to elucidate to understand the potential effects on chicken growth. Gathering and analyzing commercial records could be challenging since the datasets produced are large, without structure, and variable. Consequently, the objective of this study was to evaluate the effects of environmental, management, and housing factors on broiler live performance using statistical and machine learning techniques. Several datasets of different sizes were evaluated. Results indicated the vital importance of temperature control during the first three weeks of live performance and livability. Additionally, the results of random forest and artificial neural network analyses indicated the greater relevance of sex, transportation distance of day-old chicks from the hatchery to the farm, and the farm altitude on the BW and feed conversion ratio of broiler flock live performance parameters.

Abstract

Although temperature, relative humidity, and farm-associated factors are known to affect broiler live performance, data about the impact of these variables under commercial operations are still scarce. This study aimed to evaluate the effect of temperature, relative humidity, a thermal humidity index, management, and farm-associated factors on BW, BW gain, feed conversion ratio (FCR), and mortality of broilers raised to 35 d under commercial tropical conditions. The data analyzed included performance records of Ross 308 AP broiler flocks placed between 2018 and 2020. Environmental monitoring information was obtained from electronic sensors that captured data hourly from 80 flocks in 29 farms. Farm-associated factors were gathered using a survey of 86 farms. Three data analyses were conducted in parallel. Correlation analyses, one-way ANOVA, and machine learning techniques were employed. Results indicated that BW and BW gain were reduced, and FCR worsened (p < 0.001) up to 21 d when chickens were mainly exposed to temperatures 2.5 °C lower than the recommended optimums for each age period. At the same time, mortality at 28 and 35 d increased. In conclusion, all farm-associated factors affected chicken live performance. Variable importance analysis indicated that performance results at 14 and 21 d were significant to predict BW at 35. At the same time, sex, distance between the hatchery and farm, and farm altitude accounted for the most significant contributions from the farm-associated factors.

1. Introduction

During the last five decades, broiler chickens have been subjected to continuous genetic selection [1] to improve BW and feed conversion ratio (FCR). However, these changes have resulted in energy metabolism and expenditure alterations [2,3,4,5,6], making them very sensitive to environmental conditions [7,8]. Extensive research [9,10,11,12,13,14] has demonstrated the optimum T (thermoneutral zone) for ideal performance at all ages and the physiological consequences of keeping birds under higher (heat stress) or lower (cold stress) than the recommended T.
Under heat stress, Goo et al. [15] indicated that chickens subject to 27.8 °C from 21 to 35 d of age resulted in a final BW of 171 g less than their counterparts raised at 20 °C (1676 vs. 1847 g). Similarly, Awad et al. [16] reported a significant increase in the FCR of 0.06 g:g when chickens were exposed to 34 °C, 6 h daily from 22 to 35 d. A recent meta-analysis [17], including seven trials for BW gain and eleven experiments for FCR, demonstrated that heat-stressed chickens gained, on average, 151.4 g less. At the same time, FCR worsened by 0.17 g/g. On the other hand, several authors [5,6,18] have shown the effects of a low T on broiler performance. Ipek and Sahan [18] indicated that broilers subject to cold stress were significantly 111 g lighter on average from 1 to 5 wk of age, while FCR was higher by 0.36 g/g on average from 1 to 3 wk. Zhou et al. [6] demonstrated that chickens exposed to 16 °C between 14 and 17 d of age presented a BW of 18.84 g less and an FCR of 0.23 g/g more at the end of the trial. Similarly, Mohammadalipour et al. [5] reported an FCR 0.09 g/g poorer between 28 and 42 d (2.12 vs. 2.03 g/g) in chickens raised under a cold T from 21 to 42 d. In the same study, higher mortality due to ascites was detected in chickens raised in a cold T compared to normal conditions.
Studies evaluating the effects of relative humidity (RH) on broiler live performance have shown conflicting results across the years depending on the ambient T [19,20,21]. Thus, different T-humidity indices (THIs) have been developed to account for the variation of both environmental parameters [22,23,24,25,26]. Purswell et al. [27] observed that chickens raised between 49 and 63 d of age under a THI (THI = 0.85 × Tdb + 0.15 × Twb) value of 15.3 °C resulted in the best BW compared to chickens subjected to a THI of 26.2 °C (4547 vs. 3873 g). Concomitantly, FCR was substantially increased during this period (6.00 vs. 3.09 g/g) at this level, compared to treatments from 14.8 to 25.8 °C.
In addition to farm environmental conditions, other factors like the transportation of day-old chicks from hatchery to farm, farm location, and equipment, among others, have been related to suboptimal live performance and high mortality in broilers [28,29,30,31]. In a study conducted by Bergoug et al. [28], broilers were lighter at placement (38.4 vs. 40.0 g), 7 (153.3 vs. 158.8 g), 14 (409.5 vs. 423.0 g), and 21 d (872.7 vs. 890.8 g) when day-old chick transportation time from hatchery to farm increased from 0 to 4, and even 10 h. In other reports, Rachmawati et al. [29] showed that when farm altitude was greater than 700 m above sea level (m.a.s.l.), BW at slaughter age was reduced by 184 g (1697 vs. 1881 g). On the other hand, Garces et al. [32] and Toledo et al. [30] have indicated that litter type may affect broilers’ survival rate, BW, BW gain, and FCR.
In commercial operations, only a few reports [31,33] have described the effects of the house, season, equipment, and factors not conventionally evaluated in broiler production. Even though T and RH are usually monitored daily or hourly in broiler farms, and management factors are highly related to flock performance, data about the impact of these variables under commercial operations are still scarce.
The use of commercial data is always challenging [34,35] due to the large volumes of data, the high variability and multidimensionality, repeated measurements over time (time series data), high correlation and multicollinearity, and the lack of structure. The classical objective of statistical models is to identify the exact effect of a set of variables on another. The unstructured, correlated, and multicollinear data from commercial farms include many uncontrolled and unknown effects. If traditional statistics methodologies do not include those confounding variables, the effects measured for the variables included in the models will be biased. Consequently, the prediction ability is reduced.
In contrast, using modern data analytics and machine learning (ML) tools, it is possible to explore environmental sensor data, cluster using farm traits and categories, and identify causality and variable importance that cannot be determined with traditional statistics [34,35,36]. The ML modeling methods tolerate non-linearity, non-normality, and multicollinearity because little to no assumptions are being made. The ML models are designed to make the most accurate predictions, and methodologies include cross-validation for future utilization of models [34,35]. Confounding effects may not have the same impact in ML because the objective is not to identify the exact effect of a variable on another. The ML modeling aims to identify the most likely value of a dependent variable given a set of predictors. ML methodologies like cross-validations [37] tend to identify if the variables included in the model are affected by unidentified confounder variables. The model selection with the best predictability performance should be the least affected by the intrinsic confounding variables in these datasets [37].
This study evaluated the effect of the T, RH, a THI, growth, management, and farm-associated factors on BW, BW gain, FCR, and mortality of chickens raised to 35 d under commercial tropical conditions.

2. Materials and Methods

2.1. Database and Statistical Software

Ten datasets were obtained from a broiler integration located in Colombia. One dataset contained live performance parameters from 1347 male and 1353 female Ross 308 AP broiler flocks collected weekly from the placement up to slaughter age in 86 farms related to one poultry integrator company. House sizes vary between 10 × 50 m and 16 × 150 m, and some had two or three floors. However, since not all flocks were processed at the same age between 36 and 45 d, the cutoff point on the data was 35 d. Performance records included BW, BW gain, FCR, and mortality of all flocks placed between 2018 and 2020. The weekly performance was obtained from the averages of random samples of the individual BW of 1% of the flocks and feed intake from the records per house. All flocks were subjected to control feeding across three major poultry production regions in the country. Consequently, feed amounts offered were predetermined and recorded. The regions were labeled as Northwest, Midwest, and East. Groups of farms were fed diets of similar composition since feed formulation details are common to the whole company. Farm-associated conditions were collected using a survey for all 86 farms and are detailed in Table 1 and Table 2.
Only 80 flocks had complete hourly data of continuous environmental monitoring with all electronic sensors between January and July 2020 in 29 farms. Electronic environmental monitoring only started in this company in 2019, and clean data were processed in 2020. Some flocks were eliminated from the dataset because at least one of the environmental parameters recorded did not have all the data or presented several outliers due to faulty electronic sensors. These data corresponded to 38 male and 42 female Ross 308 AP broiler flocks. The other nine datasets mentioned contained the broiler house T and RH records data logged automatically every hour over the entire growing cycle. The electronic sensors and datalogger equipment used were provided by Asimetrix Inc. (Asimetrix Inc., Durham, NC, USA). The T and RH were automatically captured using two sensors located within each broiler house and sent via the internet to a central database administered by the provider. Environmental records were shared for analysis as raw data. Then, data obtained from the two sensors were averaged for each hour and identified for each farm and production region to be later used for analysis.
Data were cleaned, organized, and analyzed using R (R Core Team, 2021) in RStudio (RStudio Team, Boston, MA, USA) and JMP Pro 15 (SAS Institute, Cary, NC, USA). Three parallel analyses were carried out with the data collected.

2.2. Data Analysis 1: Performance and Environment

The hourly average environmental records captured in 29 farms during the five weeks of each of the 80 flocks were imported to statistical software and merged by farm, date, and time. This corresponds to 67,200 lines of environmental data for all 80 flocks. Then, the placement date and the date at 35 d were added to the environmental dataset to estimate the age of chickens within each flock. These data were highly variable, and other authors have reported that traditional statistical procedures did not find good relationships with performance [36]. Consequently, we calculated the amount of time (hours) chickens were exposed to non-recommended environmental conditions from hatch up to 35 d or within 840 h. The Aviagen management recommendations [38] for house T and RH were used as optimum values per age. The ideal range of the T was ±1.5 °C of the recommended average values. For each flock, exposure time was calculated to be lower or greater than 1.5, 2.5, 5.0, and 7.0 °C from the recommendations with an RH below 50 and above 75%. Additionally, a THI was computed for each hour based on the equation described by Berman et al. [39] as follows:
T H I = 3.43 + 1.058 × T 0.293 × R H + 0.0164 × T × R H + 35.7
where
  • T = Temperature;
  • RH = Relative humidity.
Index boundaries for each day of age were calculated based on the recommended T and RH combination. Exposure time to THI levels out of limits was also considered for analysis. Since feed allowance was controlled during the growth of chickens, this variable was not included in the data analysis.
Correlation analyses were conducted between the BW, BW gain, FCR, and mortality parameters of the 80 flocks with environmental monitoring and hours of exposure to non-recommended T, RH, and THI levels. Since these flocks were subjected to control feeding, feed intake was not considered for these analyses. Correlations were fit using either week or cumulative periods of exposure from placement to each week of age up to 35 d. Linear regressions were carried out on correlated variables with live performance parameters as a response and environmental exposure time as a predictor.

2.3. Data Analysis 2: Farm Management Factors and Performance

Farm-associated and management factors of 86 farms were joined to the broiler performance records of 2700 flocks (1347 male and 1353 female flocks) using the farm name as a link. Continuous variables from the farm management survey conducted in 86 farms, like altitude, litter reuse cycle number, and downtime between flocks, were clustered by FCR at 35 d to determine groups to compare. The clusterization method used was the k-means procedure [40] in JMP 15 (SAS Institute, Cary, NC, USA). Several clusters were developed based on an iterative fitting process. Briefly, the k-means procedure algorithm randomly took many starting points as cluster seeds or centers and computed the distance between each pair of datapoints. Then, it iteratively moved the centers to minimize the total within-cluster variance [40]. The k-means cluster with the lowest AIC and BIC was chosen for further analysis. Broiler flocks with missing data related to management conditions or those excluded within the clusterization procedure were not considered for analysis.
After clustering all possible variables, data were analyzed in a one-way ANOVA with variable cluster (continuous variables) or factors (categorical variables) as main effects and BW, FCR, or cumulative mortality at 35 d as responses. Mean separation was performed using the LS means method using Tukey’s or student’s t-test at a significance level of alpha = 0.05.
The management survey indicated that the company was divided into three major poultry production regions: Northwest, Midwest, and East, which received chicks from two company-owned hatcheries (Hatchery 1 and Hatchery 2). Farms also had two technification levels, a high level that accounted for farms with greater capacity, more houses, and more retrofitted houses and a low technification level for farms with less capacity and few houses with open-sided structures. Farms were located up to 531 km or 11 h away from the hatchery. These two variables were clustered into three groups, near (73 ± 71 km), intermediate (356 ± 169 km), and far (454 ± 64 km) for distance in km. Three groups accounted for the distance measured in hours, near (2.1 ± 1.6 h), intermediate (7.0 ± 2.9 h), and far (8.2 ± 1.8 h). Also, farms had an average altitude of 1381 ± 449 m.a.s.l., which were grouped into low (1041 ± 349 m.a.s.l.), middle (1446 ± 498 m.a.s.l.), and high (1670 ± 434 m.a.s.l.) altitude. Poultry houses were equipped with rice hulls or wood shavings as litter material, and the frequency at which the farm recycled the litter was clustered into 0.8 ± 0.7, 1.3 ± 1.2, 5.8 ± 2.2, and 12.0 ± 0.01 times. The downtime between flocks resulted in three ranks, 12.7 ± 1.0, 13.5 ± 1.2, and 15.0 ± 1.6 d.
Other characteristics linked to house infrastructure were also detailed. It was observed that some farms had poultry houses with 1 or 2 stories, which allowed them to place chicks on both floors, and the house floor was composed of two materials—soil or concrete. Additionally, different systems of water storage were reported by farm managers. Water tanks that had any protection against sun and rain (covered), tanks outside the house with no environmental protection (uncovered), and water tanks located underground.

2.4. Data Analysis 3: Prediction of Performance with ML

Three supervised ML methods were evaluated in this analysis [35,37]. These methodologies were multiple linear regression (MLR), random forest (RF), and artificial neural networks (ANN). An MLR model was fitted using the lm function from the R stats package. Additionally, an RF and ANN were trained and validated using the caret package [37].
Farm-associated factors were joined to the performance records from all flocks to the model chickens’ BW and FCR at 35 d. BW and FCR at 35 d were used as responses, while performance records and farm-associated conditions were used as predictors to fit these models. Only broiler performance records up to 21 d of age were included to obtain a prediction two weeks before slaughter age.
Data were automatically divided into random training and testing datasets using a five-fold cross-validation procedure to fit the models. Models were then trained using four-folds, and the predictability performance was assessed with the hold-out fold, repeating this step for each of the five folds [37,40,41]. This k-fold cross-validation inevitably results in a lower coefficient of determination (R2) but more confidence for predictability [42]. Categorical predictors such as farm region or litter type were transformed into dummy or indicator variables that denote the absence or presence of a factor with 0 and 1, respectively, and data were normalized in a range from 0 to 1 [37,43].
The R2 and root mean square error (RMSE) were calculated [37] and averaged across five folds for each model for model evaluation [42]. In the MLR, collinearity was evaluated with the variance inflation factor (VIF). Variable selection was carried out based on the p-value. Thus, collinear and non-significant predictors were manually removed to obtain the final model.
A variable importance analysis [44,45] was conducted in all models. This analysis indicates how much a specific variable contributes to the predictive response in a given model to obtain accurate predictions. In general, variable importance is calculated with the relative influence of each variable, and the specific calculation method varies among the three methods evaluated in the analyses described here. In the MLR, the variable contribution was estimated as the predictor percentage sum of squares over the sum of squares of the model. For RF, the permutation importance method [46] was employed to determine variable importance and was expressed as the increase in mean square error (MSE). The ANN utilized the connection weights approach described by Olden et al. [47] to calculate the order of variable importance.

3. Results

3.1. Data Analysis: Performance and Environment

The house environmental data corresponded to 29 farms and 80 flocks analyzed. Correlation coefficients of live performance parameters and exposure time to the T below or above recommended values were determined. The dataset did not contain flocks exposed to temperatures above 9 °C from the recommended values of the standard genetic line.
The flocks with higher r and R2 are shown in Figure 1 and Figure 2. Overall, live performance results were mainly correlated (p < 0.05) to exposure to the non-recommended T up to 21 d, except for cumulative mortality, which included cumulative exposure values up to 28 d. Only 10 flocks out of the 80 evaluated with continuous environmental monitoring had a T above the target during the first 4 wk of life. No correlations between T above target and BW, BW gain, and FCR (p > 0.05) were observed throughout the growing cycle.
Negative correlations of BW (p < 0.05) were detected when chickens were exposed to a T 5 °C below the recommended values at 7, 14, and 21 d. Linear regression estimates from the same week of age indicated (p < 0.05) that chickens lost 1.51, 3.50, and 4.91 g, respectively, each hour they were exposed to these conditions. Also, cooler environmental T affected FCR at the second and third week of age (p < 0.05). Chickens exposed to a T 2.5 °C lower than the recommended from 7 to 14 d or 14 to 21 were less efficient with 0.011 and 0.004 g/g each additional hour under these conditions at 14 at 21 d, respectively. No significant correlations (p > 0.05) between a low T and mortality were observed during the first three weeks of age. However, a non-recommended T up to four weeks of age affected cumulative mortality at four and five weeks of age. As time under a T below 2.5 °C from 0 to 28 d increased (p < 0.05), mortality augmented on average 0.018% at 28 and 35 d per each additional hour of exposure.
No significant correlations (p > 0.05) between BW or BW gain with environmental T was observed at 28 d. FCR significant correlations (p < 0.05) with time were exposed to a T below the recommended values during the last week. Positive correlations (p < 0.05) and regression analyses indicated that each hour exposed to a T below 2.5 °C up to the fourth week of age reduced the FCR by 0.004 g/g at 35 d.
On the other hand, exposure to a T above the recommended values affected mortality (p < 0.05) from 28 to 35 d. Regression analyses resulted in the following equations. Mortality = 0.208 + 0.028 × Hours of exposure (p < 0.001; R2 = 0.88). It indicates that chicken flocks exposed to 7 °C above the target T from 7 to 14 d of age resulted in a mortality 0.028% greater in the last week of growth for each hour of exposition. Consequently, the same variable similarly affected cumulative mortality at 35 d (0.208% per hour of exposure to T 7 °C above target from 7 to 14 d). No other correlations (p < 0.05) or relationships were observed among flock performance variables with the T above target.
Results from correlations and regression analyses between live performance parameters and times of exposure to non-recommended values of RH are detailed in Table 3. Significant correlations (p < 0.05) between FCR or mortality and time exposed to non-recommended RH values were detected at 7, 21, and 28 d of age. However, the correlation coefficients estimated were between 0.41 and −0.51 (Table 3). The regression analyses exhibited poor goodness-of-fit with R2 lower than 0.26 (Table 3). Similarly, correlation coefficients between live performance parameters and THI levels indicated that only FCR at 21 d was positively associated (r = 0.69; p = 0.025) with exposure at lower levels of THI. The equation fitted at 21 d was FCR = 1.284 + 0.007 × Hours of exposure (p = 0.025; R2 = 0.49). This linear regression indicated that the FCR worsened by 0.007 g/g for each hour that chickens were exposed to THI levels lower than the recommended from 14 to 21 d.

3.2. Data Analysis 2: Farm Management and Infrastructure Factors and Performance

Significant effects of farm-associated factors (p < 0.05) were detected in male and female BW, FCR, and cumulative mortality at 35 d (Table 4). In this case, 86 farms and 2700 flocks were analyzed. Male chickens from the Northwest region of the country resulted in the lowest BW (p < 0.001) but, at the same time, the best FCR (p < 0.001). In contrast, females from the East region were observed to show the worst live performance with the lightest chickens and the highest (p < 0.001) FCR. Cumulative mortality at 35 d in farms from the East region was, on average, 0.39% higher in both males and females (p < 0.001) than in farms located in the Midwest (2.78 vs. 2.39%). On the technification level, farms with high levels produced heavier male chickens (p < 0.001) but were less efficient (p < 0.001) than their counterparts in the low category. Contrarily, low-level farms produced female chickens 18 g heavier (p < 0.001) and 0.017 g/g more efficiently (p < 0.001) than those raised in high-technification farms.
The live performance of chickens raised at different geographic altitudes also varied (p < 0.001). Chickens raised at higher elevation (1670 ± 434 m.a.s.l.) resulted in the best performance with males and females that were 106 and 114 g (p < 0.001) heavier, 0.120 and 0.098 g/g (p < 0.001) more efficient, with 1.11 and 0.63% (p < 0.001) less mortality at 35 d compared to chickens placed on middle-altitude farms (1446 ± 498 m.a.s.l.). Flocks raised in lower altitudes (1041 ± 349 m.a.s.l.) resulted in intermediate responses.
Conflicting results were observed in chickens from different hatcheries (Table 5). Male chickens from Hatchery 1 displayed less FCR (p < 0.001) and the lowest BW than chickens from Hatchery 2. Different from males, females exhibited the highest BW (p < 0.001) and lowest FCR (p < 0.001) when chicks were obtained from Hatchery 1. At the same time, females from Hatchery 2 were 42 g lighter and 0.038 g/g less efficient than their counterparts from Hatchery 1 at 35 d. No significant effects of chick source (p > 0.05) were identified in cumulative mortality at 35 d for males or females. Similar results were observed on chick transportation measured in the distance between the hatchery and the farm in km or time (h) needed to reach the farm. Male chicks subject to long trips (454 ± 64 km) were 97 g heavier (p < 0.001) at 35 d (1909 vs. 1812 g) than chickens subject to intermediate trips (356 ± 169 km). Still, the male and female chickens raised in farms near (73 ± 71 km) to the hatcheries were significantly (p < 0.001) 0.098 and 0.025 g/g more efficient than those chickens traveling 356 ± 169 km (1.473 vs. 1.571 g/g) or 454 ± 64 km (1.473 vs. 1.498 g/g), respectively, to reach the farm. Chickens traveling intermediate distances (356 ± 169 km) presented on average 0.86% and 0.74% more (p < 0.001) male and female mortality, respectively, (p < 0.001) than chickens from the short (73 ± 71 km) and long 454 ± 64 km clusters.
Differences due to litter type (p < 0.05) indicated that male chickens raised in rice hulls were heavier (p < 0.001) and had worse FCR (p < 0.001) compared to those raised in wood shavings. In contrast, females accounted for lower BW (p < 0.001) and worse FCR (p < 0.001) than their counterparts placed in wood shavings. Additionally, the number of times the litter was recycled affected (p < 0.05) BW, FCR, and mortality at 35 d. Flocks subject to non-recycled litters or with the minimal number of times (0.8 ± 0.7) resulted in males and females on average 80 g and 65 g lighter (p < 0.001), respectively, than flocks from all other clusters. Flocks in this cluster of low litter recycling were also (p < 0.001) the least efficient (1.556 g/g). They had the highest mortality (p < 0.001) for both males (0.63% more) and females (0.74% more) compared to flocks raised in litter recycled among all other groups. Finally, flocks with an average downtime of 13.5 d resulted (p < 0.001) in the heaviest (1906 g for males and 1778 g for females) and most efficient (1.457 g/g for males and 1.466 g/g for females) flocks at 35 d with the lowest mortality (2.51%) compared to flocks whose downtime was shorter (12.7 ± 1.0) or longer (15.0 ± 1.6).

3.3. Data Analysis 3: Prediction of Performance with ML

BW Prediction: The results of the five-fold cross-validation process presented in Table 6 indicated that the RF had the better fit (R2 = 0.78), followed by the ANN (R2 = 0.73) and MLR (R2= 0.59). Hyperparameter tuning [37,42,43] indicated that the best fit of RF was achieved with 37 parameters, while the best fit of ANN included three nodes. On fitting the FCR, goodness-of-fit metrics for five-fold cross-validation resulted in an R2 of 0.18, 0.48, and 0.45 for MLR, RF, and ANN, respectively. Variable importance analyses [44,45,46,47] demonstrated that the MLR for BW included only eight non-collinear and significant predictors from which sex accounted for 89% of the total response, feed intake at 14 d with 4.41%, and BW gain at 7 d with 1.82% (Table 6). The RF model obtained for BW indicated that BW at 21 d, sex, and distance from hatchery to farm in km was ranked as the first three most important to predict the final response at 35 d (Figure 3a). Other farm factors were also ranked as critical factors for prediction. These factors included farm altitude (21.98%) and the region where the farm was located (14.05%). Among farm-associated factors, the downtime between flocks presented a 14.72% increase in MSE, becoming the most important factor. Similarly, the ANN showed that mortality at 14 d, BW at 21 d, and BW gain at 21 d were vital in the model (Figure 3b). However, this model detected predictors related to farm infrastructure like water tank location, type of house, and farm area within the top 20 factors to predict BW at 35 d.
FCR prediction: In MLR, the FCR at 21 d was the most critical variable, with 42.41% of the total response, followed by downtime between flocks with 18.67%, feed intake at 14 d with 17.63%, and percentage of retrofitted houses (infrastructure) with 7.34%. The other predictors in the model contributed less than 4% each. A total of 10 variables were included in this model. In the variable importance analysis from the RF of FCR at 35 d, the FCR at 21 d was the most important variable (63.26%) compared to other predictors, followed by sex (24.07%) and farm elevation (18.59%) (Figure 3c). Like the BW models, BW, BW gain, and FCR at 14 and 21 d were within the 20 most important factors, and distance between farm and hatchery (16.86%) was listed within the ten first ones. The ANN demonstrated that the response strongly depended on the feed intake, FCR, and BW results at 21 d, followed by the percentage of retrofitted and open-sided houses in the farm (Figure 3d). Again, the BW and BW gain at 14 d and the downtime between flocks are essential factors for predicting FCR at 35 d.

4. Discussion

4.1. Broiler Performance and Environment

In the dataset explored in this study, no effects of the T above the target were observed, except for mortality between 28 and 35 d. These farms are located in the Andean tropical conditions, which are mountains with microclimates where the T at night and early in the morning always reduces below brooding targets. Considering the low housing insulation and open-sided structures, a low T could be a more prevalent problem than heat stress in these mountain farms.
During the first three weeks of age, cumulative exposure to 2.5 °C below the recommended T correlated to poorer BW, BW gain and FCR, and high mortality during the production cycle. In a recent study, Su et al. [48] demonstrated that chickens exposed to 3 °C lower than the control group from 8 d onwards resulted in similar FCR and BW up to three and six weeks, respectively, compared to chickens kept under thermoneutral conditions. In contrast, chickens exposed to 12 °C lower than the control group from 8 d onwards were significantly lighter and less efficient than their counterparts from two up to six weeks of age.
In another report, Candido et al. [13] indicated that a mild cold T during the early stages (27 °C first week and 24 °C in the second week) did not affect BW and FCR from 0 to 21 d, while the most significant BW gain was observed in chickens subjected to these Ts. According to those studies, 3 °C below the ideal T was still within the versatile T range for chicks, preventing a reduction in the live performance. Although in the current study, no significant correlations of T lower than 1.5 nor 2.5 °C with live performance parameters were detected within the first week of age, the sensitivity of chickens to lower T seemed to increase with age since more exposure to mild cold T worsened FCR at 14 d and BW and FCR at 21 d. Other authors that have reported similar results [2,3,5,6] have indicated that an increase in the basal metabolic rate and energy metabolism due to cold stress led to the rise in chicken energy requirements, which explained the diminished BW or increased FCR observed here. Zhou et al. [6] also suggested that these effects were related to a redistribution of nutrients during the growth toward thermoregulatory responses.
On the other hand, chickens mainly exposed to intermittently T below 2.5 °C from recommendations up to four weeks of age were correlated with high mortality and worsened FCR at 28 and 35 d. Similarly, Bruzual et al. [49] demonstrated that chicks raised to 12 d under cooler brooding T (26 ± 4 °C first week and 24 ± 4 °C second week) presented 1.71% more mortality compared to chicks subject to warmer T (32 ± 4 °C first weel and 29 ± 4 °C second week). Another report indicated that the highest mortality at 42 d resulted from chicks exposed to 26.7 °C in the first week, compared to other treatments (4.79 vs. 2.36%) that ranged between 29.4 and 35 °C [50]. The information obtained from the analyzed dataset suggested that a low T in commercial operations might lead to mortality at the end of the production cycle. It is also necessary to consider that previous reports from the literature described in the present study have been conducted under experimental conditions, achieving desired environmental T and RH through environmentally controlled chambers, which significantly differed from environmental conditions under commercial operations that mainly reflect variability in house environmental conditions.
Additionally, environmental conditions in the tropics may remarkably influence the poultry house environment due to the open-sided structure, representing 90.14% of all houses evaluated in this company. Thus, cycling T between day and night or hourly, precipitations, region, geographic altitude, cold wind drafts from mountains, ventilation and heating systems, and other variables could drive the environmental response observed in this poultry company. It is assumed that mainly the highly variable indoor environmental conditions (Figure 4) resulted in more cumulative hours of exposure to the house T and RH below recommendations. These house variable conditions affected broiler flocks during the first two weeks since it is more difficult for chickens to adapt to variable environments [51], and consequently, live performance at slaughter age was negatively affected.
On fitting the RH data ranged from 30.32 to 92%, only a few mild correlations were observed up to 28 d, while linear models among correlated variables had low R2 (0.17–0.26). Similarly, Zhou et al. [21] reported that average daily feed intake, FCR, and mortality were not affected when chickens were exposed to RH of 35, 60, or 85%. In contrast, Weaver and Meijerhof [19] detected broilers to be on average 32 g lighter at 42 d when they were either subjected to weekly increments of 8% of RH from 40 to 80% or raised at a constant 75% RH. Yahav [20] indicated that BW at 35 d was reduced when chickens were reared in environments with RH less or greater than 60–65% compared to other treatments (40–45, 50–55, and 70–75%), while no differences among treatments were observed on FCR or BW at 28 d. However, it was described [20] that the response varied at different ambient Ts (28 or 30 °C). Then, RH could rely on an ambient T to affect live performance. At 28 °C, Yahav [20] determined that the heaviest broilers were observed between 60 and 65% of RH compared to other treatments. Chickens within RH treatments of 50–55% and 70–75% but at 30 °C improved their BW (3.74%), while the 40–45% RH treatment obtained the lightest chickens in both scenarios at 28 (1438 g) and 30 °C (1398 g). Zhou et al. [21] also indicated that although average daily feed intake was not affected by RH treatments at 26 °C, this parameter was reduced by 35 and 85% RH at 31 °C.
THI analyses indicated that as hours of exposure from 14 and 21 d under THI levels below the combined Aviagen recommended T and RH increased, the FCR at 21 d worsened by 0.007 g/g. In contrast, Purswell et al. [27] did not detect effects on live performance parameters between 49 and 63 d when chickens were exposed to a THI between 14.8 and 20.7 °C. The BW, BW gain, feed intake, and FCR diminished in that study with higher THI levels. It is assumed that the well-developed feather cover that those animals could have at 49 d helped to resist lower THI levels and resulted in similar responses at slaughter age. Contrarily, chickens at 21 d from the current study may not have a good development of feathers, which causes exposure to slightly lower THI to increase mortality. According to the results presented herein, BW, BW gain, and FCR of chickens were more affected by a low T during the first three weeks of age. At the same time, mortality was associated with an increase of up to 125 h between 0 and 28 d in the exposure to a T 2.5 °C lower than recommended. The RH between 30.32 and 92% seemed to be a parameter that did not significantly affect the live performance of chickens possible to detect in this data analysis. The THI employed here [38] did not show a significant relationship with performance variables.
Finally, the methodology described herein based on hours of exposure to different environmental conditions can improve the understanding of the effects of T and RH since mean environmental values, maximum and minimum T per day, could not reflect the actual fluctuation within the poultry house. In conclusion, a lower T than recommended affected the live performance of broilers during the whole production cycle. On the contrary, only a few correlations were observed with a T above the target and RH, while the THI did not depict the effects observed with a lower T.

4.2. Farm Management and Infrastructure Factors and Performance

On farm-associated factors, evaluating 86 farms and 2700 flocks, chickens that were raised at the highest altitude (1670 ± 434 m.a.s.l.) were the heaviest (1836 g) and the most efficient (1.466 g/g) in both males and females compared to those chickens raised at a medium elevation (1446 ± 498 m.a.s.l.). Rachmawati et al. [29] indicated that the heaviest chickens at 35 d were observed when farms were located above 700 m.a.s.l. In comparison, no differences were detected between chickens raised either in lowlands (<400 m) or middle lands (400–700 m). In that study, the authors suggested that a higher BW could result from a greater feed intake associated with an increment in the maintenance requirements due to a cooler ambient T [29]. In contrast, it is assumed that higher altitudes may represent lower environmental Ts, reducing the effects of heat stress observed in low-altitude regions.
In the present study, the rice hulls as litter type also worsened the FCR of both males and females compared to wood shaving. The litter material like wood shavings, sawdust, sand, rice hulls, wheat straw, and grass in which chickens are raised is a factor that has been widely evaluated in poultry [32,52,53,54,55]. Garcês et al. [32] reported that chickens raised to 35 d in rice hulls did not significantly differ in BW, FI, or FCR compared to those raised in wood shavings. Nevertheless, the survival rate of chickens reared in rice hulls was 4.8% lower than their counterparts reared in wood shavings [32]. Still, the number of replicates in that study was limited to only three pens per treatment, which could not be enough to reproduce the results, while the current data analysis utilized, on average, 601 flocks for each litter type.
Similarly, using two replicates per treatment, Ramadan and El-Khloya [56] showed that live performance parameters or carcass traits were not affected by the litter type when using five different types of litter, including wood shavings and rice hulls. Toledo et al. [30] conducted a comprehensive meta-analysis indicating that broilers raised in wood shavings presented higher BW and better FCR than those reared on rice hulls or alternative litter, respectively. It has been suggested that rice hulls account for the more significant proliferation of fungi and, consequently, mycotoxins [57], worsening live performance in broilers.
In addition, the results presented herein indicated that BW in males and females improved when the litter was recycled more than once. Abougabal [58] observed no significant differences in the BW from 1 to 6 wk of age nor in the FCR at the end of the experiment when evaluating four litter treatments (new litter, 50% new litter–50% reused litter, 100% reused litter, and 100% reused treated litter). Still, the author suggested that recycled litter was not hazardous for broiler production performance. In contrast, an experiment conducted by Garcés-Gudino et al. [59] in tropical conditions indicated that BW (1922 vs. 1753 g), BW gain (53.7 vs. 48.8 g/d), FCR (1.588 vs. 1.633 g/g), and mortality (1.61 vs. 3.2%) at 35 d were improved when the same litter was used for two or three production cycles. It has been mentioned that the recycling process may have beneficial impacts on the performance of broilers since the litter bacterial environment could work as a probiotic or a direct-fed microbial supplementation in chickens that contributes to improving the immunity and the response against coccidia challenges, which usually affect the live performance of broilers [59,60,61,62].
The distance from the hatchery to the farm measured in both kilometers and time spent to reach the farm showed that chicks traveling the longest distances were heavier at 35 d. In contrast, Bergoug et al. [28] indicated that chicks traveling 4 and 10 h had less BW from the placement up to 21 d, where the effect of transportation time disappeared, compared to chicks that traveled less than 5 min to the farm (0 h, control). However, it was further investigated with hatchery shipping reports that the results presented herein could vary because chicks intended to travel long distances are usually scheduled to start the journey at night, which could reduce stress generated by heat, dehydration, and feed and water deprival, in comparison to chicks delivered during the day. Those aspects related to stress during transportation could be more detrimental to chick development and health than the duration of the trip per se [63,64,65,66]. Still, males and females were more efficient when post-hatch transportation was the shortest.
When evaluating these factors together with the growth parameters, the variable importance analyses from the supervised ML models demonstrated that the BW and FCR at 35 d are highly dependent on the live performance that chickens may exhibit during the second and third week of age. Other parameters like sex, the distance between the hatchery and the farm, and farm elevation were important to predict BW.

4.3. Prediction of Performance with ML

The MLR presented the lowest predictability for BW and FCR at 35 d. In contrast, the RF and ANN had the best model fit to predict both the BW and FCR at 35 d, with RF being the method with the greatest predictability capacity. It suggests that both RF and ANN could be used to predict performance results in broiler operations. Nevertheless, it is important to consider the pros and cons of each ML technique before scaling to commercial settings. Some studies have demonstrated the use of ANN to predict responses in broiler performance [67], behavior [68,69], and health [70,71]. Although those studies exhibited excellent model predictability performance (R2 > 0.99), those studies were conducted using small datasets from experimental units, including a few individuals, and under controlled conditions with specific treatments and pre-established data structures. Those characteristics do not represent the normal variability observed in data from commercial broiler operations. A few studies have explored using RF to predict responses in poultry, mainly for broiler breeders under precision feeding [72,73,74] or using sound data [75]. Only a few reports have used broiler commercial data to predict relations among the variables [34,76], indicating that RF was also the best predictive methodology.
The variation of the FCR was explained in less than 50% with these ML techniques under five-fold cross-validation. Based on the present data, it was concluded that the contribution of management, environment, and infrastructure to the FCR at 35 d is less than 50%. Most likely, its variability is associated with other factors like nutrient intake, which was not included in those models. Finally, BW at 21 d, sex, distance from hatchery to farm, feed intake, and farm altitude were the five most important variables to predict the BW at 35 d. In comparison, FCR at 21 d, sex, farm altitude, feed intake at 21, and FCR at 14 d accounted for the most essential factors in predicting FCR at 35.

5. Conclusions

In conclusion, temperatures lower than 2.5 °C added up to the variable environmental conditions reduced the BW and BW gain at 7, 14, and 21 d, while the FCR was worsened at 14 (0.011 g/g), 21 (0.004 g/g), and 35 (0.004 g/g) d of age. Not enough flocks were exposed to temperatures above 1.5 °C in these datasets to detect the effect of heat stress early in life. In addition, mortality increased by 0.018% on average at 28 and 35 d when chickens were exposed to temperatures lower than 2.5 °C. All farm-associated and chick source factors affected the performance of broilers, and sex, distance from the hatchery to the farm, and farm altitude were within the five most important variables to determine broiler BW at 35 d. The intermediate chick transportation distances (356 ± 169 km or 7.0 ± 2.9 h) worsened the final BW, FCR, and increased mortality due to environmental conditions during the hot hours of the day. The best performance was observed in high altitude (1670 ± 434 m.a.s.l.). Wood shavings, instead of rice hulls as the litter type, and litter reuse are positive for broiler flock performance.

Author Contributions

Conceptualization, E.O.O.-R. and G.A.Q.-O.; methodology, G.A.Q.-O. and E.O.O.-R.; validation, J.R.R.-R., L.C.B.-A., and G.D.M.-B.; formal analysis, G.A.Q.-O.; data curation, M.C.A.-W.; writing—original draft preparation, G.A.Q.-O.; review and editing, G.A.Q.-O. and E.O.O.-R. All co-authors contributed to the final manuscript revision. All authors have read and agreed to the published version of the manuscript.

Funding

The NC Agricultural Foundation, grant number 20-42, supported the research/extension assistantship of G.A.Q.-O.

Institutional Review Board Statement

These data correspond to records of commercial broiler production, and no chickens were subjected to experiments for this publication.

Informed Consent Statement

Informed consent was obtained from all participants involved in the study.

Data Availability Statement

The data used in this study are the property of Grupo BIOS Inc., Colombia, and are, therefore, not publicly available.

Acknowledgments

The authors thank Grupo BIOS Inc. (Envigado, Colombia) for providing the data for evaluation and constant support while conducting all the analyses. Furthermore, the authors thank all staff and students of the Prestage Department of Poultry Science at North Carolina State University.

Conflicts of Interest

The co-authors (G.A.Q.-O., J.R.R.-R., L.C.B.-A. and G.D.M.-B.) are currently employees of Grupo BIOS Inc., the organization that provided the data for evaluation. M.C.A.-W. is an employee of the Elanco Animal Health company. All of the authors declare no conflict of interest.

References

  1. Zuidhof, M.J.; Schneider, B.L.; Carney, V.L.; Korver, D.R.; Robinson, F.E. Growth, efficiency, and yield of commercial broilers from 1957, 1978, and 2005. Poult. Sci. 2014, 93, 2970–2982. [Google Scholar] [CrossRef]
  2. Chen, X.Y.; Li, R.; Wang, M.; Geng, Z.Y. Identification of differentially expressed genes in hypothalamus of chicken during cold stress. Mol. Biol. Rep. 2014, 41, 2243–2248. [Google Scholar] [CrossRef]
  3. Yang, X.; Luo, Y.H.; Zeng, Q.F.; Zhang, K.Y.; Ding, X.M.; Bai, S.P.; Wang, J.P. Effects of low ambient temperatures and dietary vitamin C supplement on growth performance, blood parameters, and antioxidant capacity of 21-day-old broilers. Poult. Sci. 2014, 93, 898–905. [Google Scholar] [CrossRef] [PubMed]
  4. Tallentire, C.W.; Leinonen, I.; Kyriazakis, I. Breeding for efficiency in the broiler chicken: A review. Agron. Sustain. Dev. 2016, 36, 66. [Google Scholar] [CrossRef]
  5. Mohammadalipour, R.; Rahmani, H.R.; Jahanian, R.; Riasi, A.; Mohammadalipour, M.; Nili, N. Effect of early feed restriction on physiological responses, performance and ascites incidence in broiler chickens raised in normal or cold environment. Animal 2017, 11, 219–226. [Google Scholar] [CrossRef]
  6. Zhou, H.J.; Kong, L.L.; Zhu, L.X.; Hu, X.Y.; Busye, J.; Song, Z.G. Effects of cold stress on growth performance, serum biochemistry, intestinal barrier molecules, and adenosine monophosphate-activated protein kinase in broilers. Animal 2021, 15, 100138. [Google Scholar] [CrossRef]
  7. Yousaf, A.; Jabbar, A.; Rajput, N.; Memon, A.; Shahnawaz, R.; Mukhtar, N.; Farooq, F.; Abbas, M.; Khalil, R. Effect of environmental heat stress on performance and carcass yield of broiler chicks. World Vet. J. 2019, 9, 26–30. [Google Scholar] [CrossRef]
  8. Kang, D.; Shim, K. Early heat exposure effect on the heat shock proteins in broilers under acute heat stress. Poult. Sci. 2021, 100, 100964. [Google Scholar] [CrossRef] [PubMed]
  9. Siegel, H.V.; Van Kampen, M. Energy relationships in growing chickens given daily injections of corticosterone. Br. Poult. Sci. 1984, 25, 477–485. [Google Scholar] [CrossRef]
  10. Faria Filho, D.E.; Rosa, P.S.; Vieira, B.S.; Macari, M.; Furlan, R.L. Protein levels and environmental temperature effects on carcass characteristics, performance, and nitrogen excretion of broiler chickens from 7 to 21 days of age. Braz. J. Poult. Sci. 2005, 7, 247–253. [Google Scholar] [CrossRef]
  11. Smith, S.M.; Vale, W.W. The role of the hypothalamic-pituitary-adrenal axis in neuroendocrine responses to stress. Dialogues Clin. Neurosci. 2006, 8, 383. [Google Scholar] [CrossRef]
  12. Cassuce, D.C.; Tinôco, I.D.; Baêta, F.C.; Zolnier, S.; Cecon, P.R.; Vieira, M.D. Thermal comfort temperature update for broiler chickens up to 21 days of age. Eng. Agríc. Jaboticabal. 2013, 33, 28–36. [Google Scholar] [CrossRef]
  13. Candido, M.G.L.; Tinoco, I.d.F.F.; Pinto, F.d.A.d.C.; Santos, N.T.; Roberti, R.P. Determination of thermal comfort zone for early-stage broilers. Eng. Agric. Jaboticabal. 2016, 36, 760–767. [Google Scholar] [CrossRef]
  14. Naga Raja Kumari, K.; Narendra Nath, D. Ameliorative measures to counter heat stress in poultry. World’s Poult. Sci. J. 2018, 74, 117–130. [Google Scholar] [CrossRef]
  15. Goo, D.; Kim, J.H.; Park, G.H.; de los Reyes, J.B.; Kil, D.Y. Effect of heat stress and stocking density on growth performance, breast meat quality, and intestinal barrier function in broiler chickens. Animal 2019, 9, 107. [Google Scholar] [CrossRef]
  16. Awad, E.A.; Najaa, M.; Zulaikha, Z.A.; Zulkifli, I.; Soleimani, A.F. Effects of heat stress on growth performance, selected physiological and immunological parameters, caecal microflora, and meat quality in two broiler strains. Asian-Australas J. Anim. Sci. 2020, 33, 778–787. [Google Scholar] [CrossRef]
  17. Liu, L.; Ren, M.; Ren, K.; Jin, Y.; Yan, M. Heat stress impacts on broiler performance: A systematic review and meta-analysis. Poult. Sci. 2020, 96, 6205–6211. [Google Scholar] [CrossRef]
  18. Ipek, A.; Sahan, U. Effects of cold stress on broiler performance and ascites susceptibility. Asian-Aust. J. Anim. Sci. 2006, 19, 734–738. [Google Scholar] [CrossRef]
  19. Weaver, W.; Meijerhof, R. The effect of different levels of relative humidity and air movement on litter conditions, ammonia levels, growth, and carcass quality for broiler chickens. Poult. Sci. 1990, 70, 746–755. [Google Scholar] [CrossRef]
  20. Yahav, S. Relative humidity at moderate ambient temperatures: Its effect on male broiler chickens and turkeys. Br. Poult. Sci. 2000, 41, 94–100. [Google Scholar] [CrossRef]
  21. Zhou, Y.; Li, X.; Zhang, M.; Feng, J. Effect of relative humidity at either acute or chronic moderate temperature on growth performance and droppings’ corticosterone metabolites of broilers. J. Integr. Agric. 2019, 18, 152–159. [Google Scholar] [CrossRef]
  22. Zulovich, J.M.; DeShazer, J.A. Estimating egg production declines at high environmental temperatures and humidities. ASAE Pap. 1990, 904021, 15. [Google Scholar]
  23. Xin, H.; DeShazer, J.A.; Beck, M.M. Responses of prefasted growing turkeys to acute heat exposure. Trans. ASAE 1992, 35, 315–318. [Google Scholar] [CrossRef]
  24. Brown-Brand, T.M.; Beck, M.M.; Schulte, D.D.; Parkhurst, A.M.; DeShazer, J.A. Temperature humidity index for growing tom turkeys. Proc. Trans. ASAE 1997, 40, 203–209. [Google Scholar] [CrossRef]
  25. Tao, X.; Xin, H. Acute synergistic effects of air temperature, humidity, and velocity on homeostasis of market–size broilers. Trans. ASAE 2003, 46, 491–497. [Google Scholar] [CrossRef]
  26. Chepete, H.J.; Chimbombi, E.; Tsheko, R. Production performance and temperature-humidity index of Cobb 500 broilers reared in open-sided naturally ventilated houses in Botswana. In Proceedings of the 7th International Symposium, Beijing, China, 18 May 2005. [Google Scholar] [CrossRef]
  27. Purswell, J.L.; Dozier III, W.A.; Olanrewaju, H.A.; Davis, J.D.; Xin, H.; Gates, R.S. Effect of temperature-humidity index on live performance in broiler chickens grown from 49 to 63 days of age. In Proceedings of the 9th International Livestock Environment Symposium, Valencia, Spain, 8–12 July 2012. [Google Scholar] [CrossRef]
  28. Bergoug, H.; Guinebretière, M.; Tong, Q.; Roulston, N.; Romanini, C.E.B.; Exadaktylos, V.; Berckmans, D.; Garain, P.; Demmers, T.G.M.; McGonnell, I.M.; et al. Effect of transportation duration of 1-day-old chicks on postplacement production performances and pododermatitis of broilers up to slaughter age. Poult. Sci. 2013, 92, 3300–3309. [Google Scholar] [CrossRef]
  29. Rachmawati, L.S.; Sjofjan, O.; Natsir, D.M.H. Effect of elevation altitude rearing and population on carcass quality of the broilers. J. Agric. Vet. Sci. 2016, 9, 36–41. [Google Scholar] [CrossRef]
  30. Toledo, T.D.S.D.; Pich, C.S.; Roll, A.A.; Dai Prá, M.A.; Leivas Leite, F.; Gonçalves Xavier, E.; Roll, V.F.B. The effect of litter materials on broiler performance: A systematic review and meta-analysis. Br. Poult. Sci. 2019, 60, 605–616. [Google Scholar] [CrossRef]
  31. Yerpes, M.; Llonch, P.; Manteca, X. Factors associated with cumulative first-week mortality in broiler chicks. Animal 2020, 10, 310. [Google Scholar] [CrossRef]
  32. Garcês, A.; Afonso, S.M.S.; Chilundo, A.; Jairoce, C.T.S. Evaluation of different litter materials for broiler production in a hot and humid environment: 2. Productive performance and carcass characteristics. Trop. Anim. Health Prod. 2017, 49, 369–374. [Google Scholar] [CrossRef]
  33. van Limbergen, T.; Sarrazin, S.; Chantziaras, I.; Dewulf, J.; Ducatelle, R.; Kyriazakis, I.; McMullin, P.; Mendez, J.; Niemi, J.K.; Papasolomontos, S.; et al. Risk factors for poor health and performance in European broiler production systems. BMC Vet. Res. 2020, 16, 287. [Google Scholar] [CrossRef]
  34. Pitesky, M.; Gendreau, J.; Bond, T.; Carrasco-Medanic, R. Data challenges and practical aspects of machine learning-based statistical methods for the analyses of poultry data to improve food safety and production efficiency. Cab. Rev. 2020, 15, 1–11. [Google Scholar] [CrossRef]
  35. Neethirajan, S. The role of sensors, big data and machine learning in modern animal farming. Sens. Bio-Sens. 2020, 29, 100367. [Google Scholar] [CrossRef]
  36. Wen, P.; Li, L.; Xue, H.; Jia, Y.; Gao, L.; Li, R.; Huo, L. Comprehensive evaluation method of the poultry house indoor environment based on grey relation analysis and analytic hierarchy process. Poult. Sci. 2021, 101, 101587. [Google Scholar] [CrossRef] [PubMed]
  37. Kuhn, M. Building predictive models in R using the caret package. J. Stat. Softw. 2008, 28, 1–26. [Google Scholar] [CrossRef]
  38. Aviagen. Broiler Management Handbook; Aviagen, Ltd.: Hunstville, AL, USA, 2018. [Google Scholar]
  39. Berman, A.; Horovitz, T.; Kaim, M.; Gacitua, H. A comparison of THI indices leads to a sensible heat-based heat stress index for shaded cattle that aligns temperature and humidity stress. Int. J. Biometeorol. 2016, 60, 1453–1462. [Google Scholar] [CrossRef]
  40. Hastie, T.; Tibshirani, R.; Friedman, J.H. Prototype methods and nearest-neighbors. In The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd ed.; Springer: New York, NY, USA, 2009; pp. 459–485. [Google Scholar]
  41. Rahimi, I.; Behmanesh, R. Improve poultry farm efficiency in Iran: Using combination neural networks, decision trees, and data envelopment analysis (DEA). Int. J. Oper. Res. 2012, 2, 69–84. [Google Scholar] [CrossRef]
  42. Bischl, B.; Lang, M.; Kotthoff, L.; Schiffner, J.; Richter, J.; Studerus, E.; Casalicchio, G.; Jones, Z.M. MLR: Machine Learning in R. J. Mach. Learn. Res. 2016, 17, 5938–5942. [Google Scholar]
  43. Kuhn, M.; Johnson, K. Data pre-processing. In Applied Predictive Modeling; Springer: New York, NY, USA, 2013; pp. 27–59. [Google Scholar]
  44. Wei, P.; Lu, Z.; Song, J. Variable importance analysis: A comprehensive review. Reliab. Eng. Syst. 2015, 142, 399–432. [Google Scholar] [CrossRef]
  45. Kazemitabar, S.J.; Amini, A.A.; Bloniarz, A.; Talwalkar, A. Variable importance using decision trees. In Proceedings of the 31st International Conference on Neural Information Processing Systems, Long Beach, CA, USA, 4 December 2017. [Google Scholar]
  46. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  47. Olden, J.D.; Joy, M.K.; Death, R.G. An accurate comparison of methods for quantifying variable importance in artificial neural networks using simulated data. Ecol. Modell. 2004, 178, 389–397. [Google Scholar] [CrossRef]
  48. Su, Y.; Li, S.; Xin, H.; Li, J.; Li, X.; Zhang, R.; Li, J.; Bao, J. Proper cold stimulation starting at an earlier age can enhance immunity and improve adaptability to cold stress in broilers. Poult. Sci. 2020, 99, 129–141. [Google Scholar] [CrossRef] [PubMed]
  49. Bruzual, J.J.; Peak, S.; Brake, J.; Peebles, E. Effects of relative humidity during the last five days of incubation and brooding temperature on performance of broiler chicks from young broiler breeders. Poult. Sci. 2000, 79, 1385–1391. [Google Scholar] [CrossRef]
  50. Deaton, J.W.; Branton, S.L.; Simmons, J.D.; Lott, B.D. The effect of brooding temperature on broiler performance. Poult. Sci. 1996, 75, 1217–1220. [Google Scholar] [CrossRef] [PubMed]
  51. Yahav, S.; Straschnow, A.; Plavnik, I.; Hurwitz, S. Effects of diurnally cycling versus constant temperatures on chicken growth and food intake. Br. Poult. Sci. 1996, 37, 43–54. [Google Scholar] [CrossRef]
  52. Angelo, J.C.; Gonzalez, E.; Kondo, N.; Anzai, N.H.; Cabral, M.M. Material de cama: Qualidade, quantidade e efeito sobre o desempenho de frangos de corte. R. Bras. Zootec. 1997, 26, 121–130. [Google Scholar]
  53. Araújo, J.S.; Oliveira, V.D.; Braga, G.C. Desempenho de frangos de corte criados em diferentes tipos de cama e taxa de lotação. Cienc. Anim. Bras. 2007, 8, 59–64. [Google Scholar]
  54. Ramadan, S.G.A.; Mahboub, H.D.; Helal, M.A.; Gaafar, K.M. Behaviour, welfare and performance of broiler chicks reared on different litter materials. Assiut Vet. Med. J. 2013, 59, 9–18. [Google Scholar] [CrossRef]
  55. Brito, D.A.P.; Brito, D.R.B.; Gomes, A.M.N.; Cunha, A.S.; Filho, U.A.S.; Pinheiro, A.A. Desempenho produtivo e rendimento de carcaça de frangos criados em diferentes materiais de cama aviária. Cienc. Anim. Bras. 2016, 17, 192–197. [Google Scholar] [CrossRef]
  56. Ramadan, S.G.A.; El-Khloya, S.Z. Do alternative litter materials affect performance, welfare and immune response of broiler chicks? Alex. J. Vet. Sci. 2017, 52, 133–141. [Google Scholar] [CrossRef]
  57. Llewellyn, G.C.; Sherertz, P.C.; Armstrong, C.W.; Miller, G.B., Jr.; Reynolds, J.D.; Kimbrough, T.D.; Bean, G.A.; Hagler, W.M., Jr.; Haney, C.A.; Trempus, C.S.; et al. Mycotoxigenic isolates and toxin production on buckwheat and rice hulls used as bedding materials. J. Ind. Microbiol. Biotechnol. 1988, 3, 351–356. [Google Scholar] [CrossRef]
  58. Abougabal, M.S. Possibility of broiler production on reused litter. Egypt. Poult. Sci. 2019, 39, 405–421. [Google Scholar] [CrossRef]
  59. Garcés-Gudino, J.; Merino, R.; Cevallos-Gordón, A.L. Litter reuse reduces Eimeria spp. oocyst counts and improves the performance in broiler chickens reared in a tropical zone in Ecuador. Europ. Poult. Sci. 2017, 82, 1–9. [Google Scholar] [CrossRef]
  60. Chapman, H.D. Anticoccidial drugs and their effects upon the development of immunity to Eimeria infections in poultry. Avian Pathol. 1999, 28, 521–535. [Google Scholar] [CrossRef] [PubMed]
  61. Cressman, M.D.; Yu, Z.; Nelson, M.C.; Moeller, S.J.; Lilburn, M.S.; Zerby, H.N. Interrelations between the microbiotas in the litter and in the intestines of commercial broiler chickens. Appl. Environ. Microbiol. 2010, 76, 6572–6582. [Google Scholar] [CrossRef] [PubMed]
  62. Smith, J. Broiler production without antibiotics: United States field perspectives. Anim. Feed Sci. Technol. 2019, 250, 93–98. [Google Scholar] [CrossRef]
  63. Decuypere, E.; Tona, K.; Bruggeman, V.; Bamelis, F. The day-old chick: A crucial hinge between breeders and broilers. World’s Poult. Sci. J. 2001, 57, 127–138. [Google Scholar] [CrossRef]
  64. Oviedo-Rondón, E.O.; Wineland, M.J.; Small, J.; Cutchin, H.; McElroy, A.; Barri, A.; Martin, S. Effect of incubation temperatures and chick transportation conditions on bone development and leg health. J. Appl. Poult. Res. 2009, 18, 671–678. [Google Scholar] [CrossRef]
  65. Abreu, L.H.P.; Junior, T.Y.; Campos, A.T.; Bahuti, M.; Fassani, E. Cloacal and surface temperatures of broilers subject to thermal stress. Eng. Agric. 2017, 37, 877–886. [Google Scholar] [CrossRef]
  66. Yerpes, M.; Llonch, P.; Manteca, X. Effect of environmental conditions during transport on chick weight loss and mortality. Poult. Sci. 2021, 100, 129–137. [Google Scholar] [CrossRef]
  67. Moharrery, A.; Kargar, A. Artificial neural network for prediction of plasma hormones, liver enzymes and performance in broilers. J. Anim. Feed Sci. 2007, 16, 293–304. [Google Scholar] [CrossRef]
  68. Fang, C.; Zhang, T.; Zheng, H.; Huang, J.; Cuan, K. Pose estimation and behavior classification of broiler chickens based on deep neural networks. Comput. Electron. Agric. 2021, 180, 105863. [Google Scholar] [CrossRef]
  69. Li, G.; Zhao, Y.; Porter, Z.; Purswell, J. Automated measurement of broiler stretching beahviors under four stocking densities via faster region-based convolutional neural network. Animal 2021, 15, 100059. [Google Scholar] [CrossRef]
  70. Işık, Y.; Kayabaşı, A. Automatic classification of healthy and sick broilers in terms of avian influenza by using neural networks. Müh. Bil. Ve Araş. Derg. 2022, 24, 212–226. [Google Scholar] [CrossRef]
  71. Cuan, K.; Zhang, T.; Huang, J.; Fang, C.; Guan, Y. Detection of avian influenza-infected chickens based on a chicken sound convolutional neural network. Comput. Electron. Agric. 2020, 178, 105688. [Google Scholar] [CrossRef]
  72. You, J.; van der Klein, S.A.; Lou, E.; Zuidhof, M.J. Application of random forest classification to predict daily oviposition events in broiler breeders fed by precision feeding system. Comput. Electron. Agric. 2020, 175, 105526. [Google Scholar] [CrossRef]
  73. You, J.; Lou, E.; Afrouziyeh, M.; Zukiwsky, N.M.; Zuidhof, M.J. A supervised machine learning method to detect anomalous real-time broiler breeder body weight data recorded by a precision feeding system. Comput. Electron. Agric. 2021, 185, 106171. [Google Scholar] [CrossRef]
  74. You, J.; Lou, E.; Afrouziyeh, M.; Zukiwsky, N.M.; Zuidhof, M.J. Using an artificial neural network to predict the probability of oviposition events of precision-fed broiler breeder hens. Poult. Sci. 2021, 100, 101187. [Google Scholar] [CrossRef] [PubMed]
  75. Sun, Z.; Gao, M.; Zhang, M.; Lv, M.; Wang, G. Research on recognition method of broiler overlapping sounds based on random forest and confidence interval. Comput. Electron. Agric. 2023, 209, 107801. [Google Scholar] [CrossRef]
  76. Cahyaningtyas, C.; Manongga, D.; Sembiring, I. Algorithm comparison and feature selection for classification of broiler chicken harvest. J. Tek. Inform. (JUTIF) 2022, 3, 1717–1727. [Google Scholar] [CrossRef]
Figure 1. Relationship between the live performance parameters of Ross 308 AP broilers from 7 to 14 d and the hours of exposure to temperatures below the recommended values by Aviagen [38]. n = 80 flocks.
Figure 1. Relationship between the live performance parameters of Ross 308 AP broilers from 7 to 14 d and the hours of exposure to temperatures below the recommended values by Aviagen [38]. n = 80 flocks.
Animals 13 03312 g001
Figure 2. Relationship between live performance parameters of Ross 308 AP broilers from 21 to 35 d and the hours of exposure to temperatures below the recommended values by Aviagen [38]. n = 80 flocks.
Figure 2. Relationship between live performance parameters of Ross 308 AP broilers from 21 to 35 d and the hours of exposure to temperatures below the recommended values by Aviagen [38]. n = 80 flocks.
Animals 13 03312 g002
Figure 3. Variable importance of (a) BW with RF, (b) BW with ANN, (c) FCR with RF, and (d) FCR with ANN at 35 d on Ross 308 AP chickens raised under commercial tropical conditions.
Figure 3. Variable importance of (a) BW with RF, (b) BW with ANN, (c) FCR with RF, and (d) FCR with ANN at 35 d on Ross 308 AP chickens raised under commercial tropical conditions.
Animals 13 03312 g003
Figure 4. Mean ± SD of house temperature from flocks that presented (a) more than 450 h (53.6% of the time) within the temperatures recommended by Aviagen [38] for chickens between 1 and 35 d and (b) flocks affected by cumulative exposure to temperatures 2.5 °C lower than recommended. These graphics demonstrate the higher variability or wider confidence interval observed in houses (b) with flocks with more cumulative hours below target.
Figure 4. Mean ± SD of house temperature from flocks that presented (a) more than 450 h (53.6% of the time) within the temperatures recommended by Aviagen [38] for chickens between 1 and 35 d and (b) flocks affected by cumulative exposure to temperatures 2.5 °C lower than recommended. These graphics demonstrate the higher variability or wider confidence interval observed in houses (b) with flocks with more cumulative hours below target.
Animals 13 03312 g004
Table 1. Descriptive statistics for farm-associated categorical variables.
Table 1. Descriptive statistics for farm-associated categorical variables.
VariableCategoryNorthwestMidwestEast
n%n%n%
--------------------------------------- (farms 1) ---------------------------------------
Farms 4248.841719.772731.40
Technification level 2High78.14910.471011.63
Low3540.7089.301719.77
Chick sourceHatchery 14248.8400.0000.00
Hatchery 200.001719.772731.40
Litter typeRice hulls00.0011.162731.40
Wood shavings4248.841618.6000.00
--------------------------------------- (houses 3) ---------------------------------------
Houses 36749.8015320.7621729.44
Type of houseOpen-sided35848.3111815.9219426.18
Retrofitted152.02212.83233.10
Controlled00.00121.6200.00
House stories 4126235.7414319.5121729.60
210113.78101.3600.00
House floor typeSoil27337.24618.329613.10
Concrete9012.289212.5512116.51
Water storage systemCovered20227.4111816.0114820.08
Uncovered15721.30415.56699.36
Underground20.2700.0000.00
1 Variable percentage based on the number of farms by factor level over the total number of farms (n = 86); 2 technification levels. High level = farms with greater capacity, more houses, and more retrofitted houses; low = farms with less capacity and few houses with open-sided structures; 3 variable percentage based on the number of houses by factor level over the total number of houses (n = 737). Number of houses in each farm vary between 1 and 22. 4 Number of houses with one or two stories or levels enabled to place chicks.
Table 2. Descriptive statistics for farm-associated continuous variables.
Table 2. Descriptive statistics for farm-associated continuous variables.
VariableNorthwestMidwestEast
MeanSDMinMaxMeanSDMinMaxMeanSDMinMax
Altitude (m.a.s.l.)158839810502353118430091519551184466992253
Litter reuse cycles (# times)1.91.850125.584.150120.440.5702
Distance from the hatchery (km)39.118.6580188.829.4132245444.640358531
Distance from the hatchery (h)1.50.60.135.40.746.58.31.1611
Downtime between flocks (d)14.10.4141616.31.4142113.01.61120
Farm area (m2)64354607162024,99212,6558761144031,68077959238124048,956
SD = standard deviation; Min = minimum value; and Max = maximum value. # = average number of times that litter was used.
Table 3. Correlation coefficients (r) and regression estimates of live performance parameters of Ross 308 AP broilers from 0 to 28 d and hours of exposure to a relative humidity below and above the recommended values by Aviagen [38].
Table 3. Correlation coefficients (r) and regression estimates of live performance parameters of Ross 308 AP broilers from 0 to 28 d and hours of exposure to a relative humidity below and above the recommended values by Aviagen [38].
AgeLive PerformanceRelative Humidity (%)Age of Exposure (d)rp-Value CorrelationnInterceptEstimateR2R2 AdjRMSEp-Value LinearRegression
7FCR<500 to 7−0.51<0.001460.936−0.0020.260.240.06<0.001
21Week mortality<5014 to 21−0.460.015280.155−0.0040.210.180.070.015
21Cumulative mortality<5014 to 21−0.450.015290.209−0.0030.200.170.060.015
21FCR<50 or >7514 to 210.41<0.001691.2770.0010.170.160.06<0.001
28FCR<5021 to 280.410.035271.3600.0030.170.130.060.035
28FCR>750 to 280.410.001681.3840.00020.170.160.060.001
28FCR<50 or >750 to 280.41<0.001701.3720.00020.170.160.06<0.001
FCR = feed conversion ratio at each age. n = 80.
Table 4. Effect of farm-associated factors on BW, FCR, and cumulative mortality at 35 d of male and female Ross 308 AP broilers raised under commercial tropical conditions.
Table 4. Effect of farm-associated factors on BW, FCR, and cumulative mortality at 35 d of male and female Ross 308 AP broilers raised under commercial tropical conditions.
ItemCategory/Cluster MeanSDMalesFemales
nBWFCRMortalitynBWFCRMortality
Flocks--(g)---(g/g)--- (%) --flocks-- (g) ---(g/g)---(%)--
RegionNorthwest 3811836 c1.461 c2.98 a3851729 b1.514 c2.68 ab
Midwest 1481947 a1.493 b2.52 b1541751 a1.533 b2.27 b
East 6691882 b1.515 a2.85 a6701673 c1.557 a2.71 a
SEM ± 40.0030.10 40.0030.09
CV % 4.13.326.6 4.32.927.3
Technification levelHigh 4971889 a1.507 a2.945021690 b1.550 a2.72 a
Low 6871866 b1.487 b2.796941708 a1.533 b2.54 b
SEM ± 30.0020.07 30.0020.07
CV % 4.43.626.8 4.73.128.2
Altitude (m.a.s.l.)10413493761873 b1.522 b2.69 b3591730 b1.519 b1.96 c
14464981721792 c1.577 a3.82 a6451660 c1.573 a3.06 a
16704346501898 a1.457 c2.71 b2051774 a1.475 c2.43 b
SEM ± 40.0030.09 40.0020.09
CV % 4.02.226.1 3.92.027.4
Source of variation -------------------------------------- p-value --------------------------------------
Region <0.001<0.0010.019 <0.001<0.0010.003
Type of administration <0.001<0.0010.087 <0.001<0.0010.008
Altitude <0.001<0.001<0.001 <0.001<0.001<0.001
Litter typeRice hulls 6721883 a1.514 a2.856761674 b1.556 a2.69 a
Wood shaving 5261867 b1.470 b2.875331733 a1.519 b2.54 b
SEM ± 30.0020.07 30.0020.07
CV % 4.43.426.8 4.32.928.9
Litter reuse cycles number0.80.73481819 b1.556 a3.47 a7971671 c1.561 a2.89 a
1.31.27671898 a1.465 c2.59 c2971769 a1.481 c2.10 b
5.82.2321898 a1.564 a3.41 ab691722 b1.562 a2.10 b
12.00.0511898 a1.479 b2.53 bc461718 b1.522 b2.24 b
SEM ± 80.0040.17 60.0030.15
CV % 4.02.426.0 4.02.328.1
Downtime between flocks (d)12.71.03141822 c1.562 a3.44 a6781662 c1.569 a2.94 a
13.51.26621906 a1.457 c2.51 b1501778 a1.466 c2.05 b
15.01.62221859 b1.510 b3.11 a3811735 b1.517 b2.31 b
SEM ± 40.0020.09 40.0020.09
CV % 4.02.225.9 3.92.128.1
Source of variation ---------------------------------------- p-value -------------------------------------
Litter type 0.001<0.0010.657 <0.001<0.0010.011
Litter reuse cycle number <0.001<0.001<0.001 <0.001<0.001<0.001
Downtime between flocks <0.001<0.001<0.001 <0.001<0.001<0.001
a–c Means in columns followed by different superscript letters are statistically different with Tukey’s or student’s t-test (p < 0.05). BW = body weight; FCR = feed conversion ratio. m.a.s.l. = meters above sea level.
Table 5. Effect of chick source and hatchery distance on BW, FCR, and cumulative mortality at 35 d of male and female Ross 308 AP broilers raised under commercial tropical conditions.
Table 5. Effect of chick source and hatchery distance on BW, FCR, and cumulative mortality at 35 d of male and female Ross 308 AP broilers raised under commercial tropical conditions.
ItemCategory/Cluster MeanSDMalesFemales
nBWFCRMortalitynBWFCRMortality
Flocks--(g)---(g/g)--- (%) --Flocks-- (g) ---(g/g)---(%)--
Chick sourceHatchery 1 3811836 b1.461 b3.003851729 a1.514 b2.65
Hatchery 2 8171894 a1.511 a2.798241687 b1.552 a2.62
SEM ± 30.0030.07 30.0020.07
CV % 4.23.426.8 4.53.00.1
Distance from hatchery (km)73714611875 b1.455 c2.71 b3221758 a1.491 c2.18 b
3561692571812 c1.572 a3.54 a7221665 c1.570 a2.92 a
454644801909 a1.492 b2.65 b1651740 b1.503 b2.21 b
SEM ± 40.0020.09 40.0020.09
CV % 4.02.326.2 3.92.028.1
Distance from hatchery (h)2.11.64291873 a1.449 c2.73 b2581755 a1.486 c2.21 b
7.02.92471808 b1.573 a3.58 a7141663 b1.571 a2.93 a
8.21.85221909 a1.496 b2.63 b2371750 a1.507 b2.16 b
SEM ± 40.0023.85 40.0023.85
CV % 4.02.226.2 3.92.028.1
Source of variation ---------------------------------------- p-value ----------------------------------------
Chick source <0.001<0.0010.146 <0.001<0.0010.349
Distance from the hatchery (km) <0.001<0.001<0.001 <0.001<0.001<0.001
Distance from the hatchery (h) <0.001<0.001<0.001 <0.001<0.001<0.001
a–c Means in columns followed by different superscript letters are statistically different with Tukey’s or student’s t-test (p < 0.05).
Table 6. Variable contribution of MLR fit with growth and farm-associated conditions on BW and FCR at 35 d of Ross 308 AP chickens reared under commercial tropical conditions.
Table 6. Variable contribution of MLR fit with growth and farm-associated conditions on BW and FCR at 35 d of Ross 308 AP chickens reared under commercial tropical conditions.
VariableSum SqContribution (%) 1F-Valuep-Value
BW
Females13,108,737.089.272325.7<0.001
Feed Intake 14 d646,918.04.41114.8<0.001
BW Gain at 7 d267,502.81.8247.5<0.001
FCR at 21 d209,950.51.4337.2<0.001
Percentage of Retrofitted200,095.71.3635.5<0.001
Percentage of Open-Sided119,749.40.8221.2<0.001
Region79,587.10.5414.1<0.001
Percentage of Underground Tanks51,800.20.359.20.002
FCR
FCR at 21 d0.4142.41162.21<0.001
Downtime Between Flocks0.1818.6771.41<0.001
Feed Intake 14 d0.1717.6367.43<0.001
Percentage of Retrofitted0.077.3428.09<0.001
Percentage of Underground Tanks0.043.9415.06<0.001
Percentage of Open-Sided0.032.6910.290.001
Percentage of Concrete Floor0.022.569.800.002
Feed Intake 21 d0.022.228.510.004
Month at Placement0.011.475.620.018
Mortality at 21 d0.011.064.040.045
1 Variable contribution of MLR models expressed as the percentage of the variable sum of squares over the sum of squares of the model. Sum sq = Sum of squares; BW = body weight; FCR = feed conversion ratio.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Quintana-Ospina, G.A.; Alfaro-Wisaquillo, M.C.; Oviedo-Rondon, E.O.; Ruiz-Ramirez, J.R.; Bernal-Arango, L.C.; Martinez-Bernal, G.D. Effect of Environmental and Farm-Associated Factors on Live Performance Parameters of Broilers Raised under Commercial Tropical Conditions. Animals 2023, 13, 3312. https://doi.org/10.3390/ani13213312

AMA Style

Quintana-Ospina GA, Alfaro-Wisaquillo MC, Oviedo-Rondon EO, Ruiz-Ramirez JR, Bernal-Arango LC, Martinez-Bernal GD. Effect of Environmental and Farm-Associated Factors on Live Performance Parameters of Broilers Raised under Commercial Tropical Conditions. Animals. 2023; 13(21):3312. https://doi.org/10.3390/ani13213312

Chicago/Turabian Style

Quintana-Ospina, Gustavo A., Maria C. Alfaro-Wisaquillo, Edgar O. Oviedo-Rondon, Juan R. Ruiz-Ramirez, Luis C. Bernal-Arango, and Gustavo D. Martinez-Bernal. 2023. "Effect of Environmental and Farm-Associated Factors on Live Performance Parameters of Broilers Raised under Commercial Tropical Conditions" Animals 13, no. 21: 3312. https://doi.org/10.3390/ani13213312

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop