Acessibilidade / Reportar erro

Prediction of weight and percentage of salable meat from Brazilian market lambs by subjective conformation and fatness scores

ABSTRACT

This study assessed the use of conformation and fatness scores of the EUROP sheep carcass grading system to predict weight and percentage of salable meat from Brazilian market lambs. Data were collected from in vivo, carcass, and retail production from 252 uncastrated lambs. Evaluated models included single regressions, two multivariate models, and one determined by the stepwise procedure. Conformation was moderately correlated with weight of salable meat. Fatness scores were correlated with rump perimeter, carcass width, and thoracic depth with coefficients of −0.33, −0.32, and −0.23, respectively. Body weight was the best single predictor for weight of salable meat and cold carcass yield for percentage of salable meat. All multivariate models for weight of salable meat prediction were significant. Stepwise regression with body weight, leg perimeter, thoracic depth, rump perimeter, and fatness scores predicted 98% of weight of salable meat variation. For percentage of salable meat prediction, stepwise regression with cold carcass yield, leg perimeter, and conformation score was significant. The EUROP conformation and fatness scores can be used in Brazil for the prediction of lamb meat production.

Key Words:
regression models; retail cuts; sheep

Introduction

Carcass grading system facilitates trade by the description of the most important commercial traits, directs products for appropriate niche markets, helps marketing, and is used as a public policy tool to regulate the sector (Price, 1995Price, M. A. 1995. Development of carcass grading and classification systems. p.173-199. In: Quality and grading of carcasses of meat animals. Jones, S. D. M., ed. CRC Press, Boca Raton.).

Muscularity and fatness are carcass traits used as basis for grading systems, evaluated subjectively or objectively, which have great impact on the carcass value. Muscularity indicates the amount of muscle tissue, determined by conformation or muscle to bone ratio and fatness describes the external and internal deposition of fat in the carcass.

Despite showing weak relationship with the weight and percentage of salable meat, subjective carcass assessment, for conformation and fatness, is considered important for both farmers and industry (Nsoso et al., 2000Nsoso, S. J.; Young, M. J. and Beatson, P. R. 2000. A review of carcass conformation in sheep: assessment, genetic control and development. Small Ruminant Research 35:89-96.) and is used in various grading systems.

Australia and New Zealand, accounting for 13% of world production and 68% of the exported volume of sheep meat (FAO, 2015), use automated methods, described by Stanford et al. (1998Stanford, K.; Jones, S. D. M. and Price, M. A. 1998. Methods of predicting lamb carcass composition: A review. Small Ruminant Research 29:241-254.), which predict the ratio of muscle and fat in the carcass. However, similar to many EU countries, the United Kingdom, a traditional sheep producer, uses a subjective grading system based on conformation and fatness of the carcasses (AHDB Beef & Lamb, 2012).

Automated methods for carcass grading are more elaborate, have higher accuracy, but are more expensive, and the use of the method is justified mainly by the volume produced. Subjective methods of grading are cheaper, quicker, and easier to apply, despite the lower accuracy (Stanford et al., 1998Stanford, K.; Jones, S. D. M. and Price, M. A. 1998. Methods of predicting lamb carcass composition: A review. Small Ruminant Research 29:241-254.). Our objective was to assess if EUROP conformation and fatness scores can be used to predict retail cut production of Brazilian market lambs.

Material and Methods

The study was conducted in a commercial slaughterhouse registered in the Federal Inspection Service, located in São Paulo State. All procedures for handling and slaughter of animals in the company were conducted in accordance with the Regulation of Industrial and Sanitary Inspection of Animal Products in the country (Brasil, 1952) and Normative Instruction No. 56, of November 06, 2008, on the general procedures of Recommendations of Good Practices of Welfare for Farm Animals with Economic Interest, covering the production and transportation systems (Brasil, 2008).

The animals were evaluated pre-slaughter in the fasting area, held in the corral of the slaughterhouse, with 252 non-castrated male lambs chosen at random. Body weight (BW) was recorded with the use of a mobile mechanical scale with a capacity of 300 kg and body condition score (BCS) was determined according to the methodology described by Russell et al. (1969Russel, A. J. F.; Doney, J. M. and Gunn, R. G. 1969. Subjective assessment of body fat in live sheep. Journal Agriculture Science 72:451-454.), scoring 0.5 of a score.

The slaughter procedure was performed following the steps of electrical stunning, bleeding, and removal of the head, feet, skin (and fleece), blood, digestive system, heart, lungs, liver, kidneys, and internal fat depots. After the slaughter, each carcass had its hot carcass weight (HCW) recorded, which was used to calculate the hot carcass yield (HCY = [HCW ? BW] × 100). Cold carcass weight (CCW) was recorded after chilling for 24 h at a temperature of 4 °C and was used to calculate the cold carcass yield (CCY = [CCW ? BW] × 100) and chilling losses (CL = {[HCW - CCW] ? HCW} × 100).

In the cold carcass, linear measurements of carcass width (Wr), thoracic depth (Th), rump width (G), and rump perimeter (D) were determined as described by Carrasco et al. (2009Carrasco, S.; Ripoll, G.; Sanz, A.; Álvarez-Rodríguez, J.; Panea, B.; Revilla, R. and Joy, M. 2009. Effect of feeding system on growth and carcass characteristics of Churra Tensina light lambs. Livestock Science 121:56-63.), plus external carcass length (K, distance between the base of the neck and base of the tail from the back with a straight line) and leg perimeter (LP, maximum length around the leg).

Under the EU lamb carcass classification (EUROP) scheme (Commission Regulation (EC) No 823/98, 1998), each carcass was assessed and classified one hour after the slaughter. The classification of all carcasses was done subjectively by a trained technician. The carcass conformation scores (CS) describe the development of carcass profiles, particularly the essential components of the round, back, and shoulder. Under the EUROP system, five CS are defined, represented by the letters E, U, R, O, and P. The letters represent an incremental scale ranging from P, which denotes poor conformation, to E, representing excellent conformation.

The carcass fatness scores (FS) describe the amount of fat on the outside of the carcass and in the thoracic cavity. Five scores are defined, represented by the numbers 1, 2, 3, 4, and 5. This incremental scale ranges from 1, which denotes low fatness, to 5, denoting very high fatness. To be consistent with the statistical methods employed, CS scales were transformed to numerical classification units of 5 to 1, with 5 representing score E and 1 representing score P.

Retail cuts were based on New Zealand lamb carcass break-out (Beef + Lamb New Zealand Ltd., 2013): neck (NK, neck bone-in); shoulder (SH, shoulder banjo cut, bone-in); breast and flap (BF); loin (LO, loin with musculus longissimus lumborum bone-in between first and last lumbar vertebra); leg (LG, short cut leg, bone-in, chump-off); frenched rack (FR, rack, fully frenched); shoulder rack (SR), and boneless chump (BL). Each trimmed cut was weighed and the data was used to determine the weight of salable meat (WSM) by: WSM, in kg = (NK + SH + BF + LO + LG + FR + SR + BL). With WSM, we determined the percentage of salable meat (PSM, in % = [WSM / BW] × 100).

Data analysis was performed with the aid of the statistical package R (R Development Core Team, 2012). Descriptive analysis was performed and the Spearman correlation coefficients (cor function, stats package, R Development Core Team, 2012), among all variables were determined. Simple linear regressions (y = a + bx) and multiple regression analysis, among the variables collected in vivo and in the carcass (independent variables), and WSM and PSM (dependent variables) were assessed. For multiple regressions, three models were evaluated, one composed of CS + FS, another of HCW + CS + FS, and the third obtained by the stepwise command (lm and step functions, stats package, R Development Core Team, 2012). For all regressions, the coefficient of determination (R2), the root mean square error (RMSE, rmse function, hydroGOF package, Zambrano-Bigarini, 2012), and the Akaike information criterion (AIC, AIC function, stats package, R Development Core Team, 2012) were determined.

R statistical coding for model fitting

#Data input#

ad <- read.table("C:/Estatisticas/R/dados2.txt", header=TRUE, dec=",")

#Single regression#

sing_reg <- lm(dependent variable~independent variable, data=ad)

anova(sing_reg)

summary.lm(sing_reg)

coef <- coefficients(sing_reg); coef

R2 <- (cor(sing_reg$fitted.values, ad$dependent variable, method="pearson"))^2; R2

rmse <- rmse(sing_reg$fitted.value, ad$dependent variable,na.rm = TRUE); rmse

CIA <- step(sing_reg, direction="both")

#Multiple regression#

mult_reg1<- lm(dependent variable~conformation + fatness, data=ad)

anova(mult_reg1)

summary.lm(mult_reg1)

coef <- coefficients(mult_reg1); coef

R2 <- (cor(mult_reg1$fitted.values, ad$dependent variable, method="pearson"))^2; R2

rmse <- rmse(mult_reg1$fitted.value, ad$dependent variable,na.rm = TRUE); rmse

CIA <- step(mult_reg1, direction="both")

mult_reg2 <- lm(dependent variable~HCW + conformation + fatness, data=ad)

anova(mult_reg2)

summary.lm(mult_reg2)

coef <- coefficients(mult_reg2); coef

R2 <- (cor(mult_reg2$fitted.values, ad$dependent variable, method="pearson"))^2; R2

rmse <- rmse(mult_reg2$fitted.value, ad$dependent variable,na.rm = TRUE); rmse

CIA <- step(mult_reg2, direction="both")

mult_stpew <- lm(dependente variable~all independente variables, data=ad)

teste <- step(mult_stpew, direction="both")

stepw_sign <- lm(dependent variable ~ significant independente variables , data=ad)

anova(stepw_sign)

summary.lm(stepw_sign)

Results

The average BW of the animals was 36.1 kg with a mean of 2.69 BCS. Carcasses had an average HCW of 17.77 kg with yield of 49.19%. By the subjective assessment of the carcasses, mean scores of 2.41 and 1.59 were obtained for conformation and fatness, respectively. At the end of the manufacture of the cuts, an average WSM of 12.47 kg was obtained with 71.27% of PSM (Table 1).

Table 1
Means, standard deviation (SD), and ranges of assessed traits

The conformation scores were positively correlated with WSM, while there was no significant correlation between FS and WSM or PSM (Table 2). The carcass conformation showed moderate positive correlation coefficients with BW, HCW, CCW, CCY, and external carcass length (K) and weak negative coefficients with rump width (G) and rump perimeter (D).

Table 2
Significant Spearman correlation coefficients (P<0.05) among animal, carcass traits, and meat production

The fatness scores showed only negative coefficients with D, carcass width (Wr), and thoracic depth (Th). There was a strong correlation between WSM and BW, HCW, CCW, K, D, Th, and leg perimeter (LP) and moderate coefficients with BCS and Wr. The yield of retail cuts had only significant negative correlations with HCY and CCY and positive with LP.

For WSM, significant regressions were obtained (P<0.01) with BW, HCW, CS, CCW, K, D, Wr, Th, and LP. The body weight was the best single predictor of WSM, with coefficient of determination of 0.95, while CS had a value of 0.22. Among the linear measures of the carcass, Th was the best predictor (Table 3). Only HCY, CCY, and LP showed significant results (P<0.05) for PSM prediction, with coefficients of 0.42, 0.45, and 0.19, respectively (Table 4).

Table 3
Single regression equations to predict weight of salable meat (WSM)
Table 4
Single regression equations to predict percentage of salable meat (PSM)

The three multivariate models for WSM prediction were significant (P<0.01), in which the stepwise analysis produced the model WSM = BW + LP + Th + D + FS. The model CS + FS accounted for 22% of WSM variation, with better results for HCW + CS + FS and the stepwise model, shown above, with 92 and 98% of variation explained, respectively (Table 5). For PSM, only the stepwise model PSM = CCY + LP + CS explained a significant proportion of the variation (P<0.01).

Table 5
Multiple regression equations to predict weight and percentage of salable meat

Discussion

Body and carcass weight are the strongest single predictors and the most commonly used traits in prediction models of meat production for total production and yield or to predict the tissue composition of the carcass, or even meat quality traits, such as tenderness (Smith et al., 1969Smith, G. C.; Carpenter, Z. L. and King, G. T. 1969. Ovine carcass cutability. Journal of Animal Science 29:272-282.; Safari et al., 2001Safari, E.; Hopkins, D. L. and Fogarty, N. M. 2001. Diverse lamb genotypes 4. Predicting the yield of saleable meat and high value trimmed cuts from carcass measurements. Meat Science 58:207-214.; Brady et al., 2003Brady, A. S.; Belk, K. E.; LeValley, S. B.; Dalsted, N. L.; Scanga, J. A.; Tatum, J. D. and Smith, G. C. 2003. An evaluation of the lamb vision system as a predictor of lamb carcass red meat yield percentage. Journal of Animal Science 81:1488-1498.; Díaz et al., 2004Díaz, M. T.; Cañeque, V.; Lauzuruca, S.; Velasco, S.; Ruíz de Huidobro, F. and Pérez, C. 2004. Prediction of suckling lamb carcass composition from objective and subjective carcass measurements. Meat Science 66:895-902.; Lambe et al., 2009Lambe, N. R.; Navajas, E. A.; Bünger, L.; Fisher, A. V.; Roehe, R. and Simm, G. 2009. Prediction of lamb carcass composition and meat quality using combinations of post-mortem measurements. Meat Science 81:711-719.).

The results showed the importance of weight for lamb carcass grading systems, in which the main lamb-producing countries, such as Australia, New Zealand, and the European Union, use hot or cold carcass weight to determine the price, which helps marketing. According to Price (1995Price, M. A. 1995. Development of carcass grading and classification systems. p.173-199. In: Quality and grading of carcasses of meat animals. Jones, S. D. M., ed. CRC Press, Boca Raton.), the reason for using weight as a grading criterion is the standardization of the production, which facilitates the processing of carcasses, manufacture, and packaging of the cuts.

Although Brazil has an important global market share of animal protein in the areas of beef, pork, and poultry (Ferraz and Felício, 2010Ferraz, J. B. S. and Felício, P. E. 2010. Production systems - An example from Brazil. Meat Science 84:238-243.), the Brazilian sheep industry has low international market share, accounting for only 1.0% of total production. The insignificance of the sector, even to the domestic market, means that there is little, if any, investment in technologies for carcass grading or prediction of the carcass yield.

Officially, Brazil has a system of sheep carcass grading regulated by the Ordinance number 307 of December 26, 1990 (Brasil, 1990). Despite using muscularity and fatness traits, the national system is complex, and there are no official records of grading. In this scenario, subjective grading methods can be used to assess, guide the production, and predict meat production.

However, for WSM and PSM, CS and FS were not good predictors, either individually, or in combination. Jones et al. (1993Jones, S. D. M.; Robertson, W. M. and Price, M. A. 1993. The assessment of saleable meat yield in lamb carcasses. p.190. In: International Congress on Meat Science and Technology, Calgary.) and Stanford et al. (1997Stanford, K.; Woloschuk, C. M.; McClelland, L. A.; Jones, S. D. M. and Price, M. A. 1997. Comparison of objective external carcass measurements and subjective conformation scores for prediction of lamb carcass quality. Canadian Journal of Animal Science 77:217-223.) found that subjective traits can show good prediction results when used in heterogeneous lots of animals. Einarsson et al. (2014Einarsson, E.; Eythórsdóttir, E.; Smith, C. R. and Jónmundsson, J. V. 2014. The ability of video image analysis to predict lean meat yield and EUROP score of lamb carcasses. Animal 8:1170-1177.) and Rius-Vilarrasa et al. (2009Rius-Vilarrasa, E.; Bünger, L.; Maltin, C.; Matthews, K. R. and Roehe, R. 2009. Evaluation of Video Image Analysis (VIA) technology to predict meat yield of sheep carcasses on-line under UK abattoir conditions. Meat Science 82:94-100.) obtained R2 values of 0.49 and 0.41, respectively, with the EUROP system for prediction of meat yield and primal cut yield, the same way as in the work of Johansen et al. (2006Johansen, J.; Aastveit, A. H.; Egelandsdal, B.; Kvaal, K. and Røe, M. 2006. Validation of the EUROP system for lamb classification in Norway; repeatability and accuracy of visual assessment and prediction of lamb carcass composition. Meat Science 74:497-509.). As commercial market lambs were used in the current study directly from the slaughterhouse, the demand made by the industry for uniformity in purchased lambs provided a certain homogeneity, which can be seen by the BW variation observed.

Although the results for CS and FS were not satisfactory, these traits were included in the multivariate models generated by the stepwise regression procedure. In the case of WSM, this inclusion indicates certain importance of FS among other traits that were not included in the model, but had relationships with WSM, such as K and Wr (Table 2). Díaz et al. (2004Díaz, M. T.; Cañeque, V.; Lauzuruca, S.; Velasco, S.; Ruíz de Huidobro, F. and Pérez, C. 2004. Prediction of suckling lamb carcass composition from objective and subjective carcass measurements. Meat Science 66:895-902.) used the stepwise procedure for multivariate prediction of weight and yield of lean, fat, and bone from lamb carcasses. Likewise, FS was included in the models, in this case, for weight and yield of fat.

With the exception of the measure G, the linear measurements of the carcass showed good prediction of WSM, both in the simple regression and as part of the model generated by stepwise regression. These results support the potential of using video image analysis (VIA), which uses prediction equations based on external carcass measurements, for the prediction of meat production (Cunha et al., 2004Cunha, B. C. N.; Belk, K. E.; Scanga, J. A.; LeValley, S. B.; Tatum, J. D. and Smith, G. C. 2004. Development and validation of equations utilizing lamb vision system output to predict lamb carcass fabrication yields. Journal of Animal Science 82:2069-2076.; Hopkins et al., 2004Hopkins D. L.; Safari, E.; Thompson, J. M. and Smith, C. R. 2004. Video image analysis in the Australian meat industry - precision and accuracy of predicting lean meat yield in lamb carcasses. Meat Science 67:269-274.; Rius-Vilarrasa et al., 2009Rius-Vilarrasa, E.; Bünger, L.; Maltin, C.; Matthews, K. R. and Roehe, R. 2009. Evaluation of Video Image Analysis (VIA) technology to predict meat yield of sheep carcasses on-line under UK abattoir conditions. Meat Science 82:94-100.; Einarsson et al., 2014Einarsson, E.; Eythórsdóttir, E.; Smith, C. R. and Jónmundsson, J. V. 2014. The ability of video image analysis to predict lean meat yield and EUROP score of lamb carcasses. Animal 8:1170-1177.). Although the perimeter measures have showed relationship with WSM and PSM, the use of these online would be impractical by the method of obtainment, unlike the measurements obtained in the dorsal and lateral images of carcasses, captured in VIA systems.

Conclusions

The EUROP conformation and fatness scores can be used to predict the weight and percentage of salable meat from market lambs in Brazil. The combination of conformation and fatness with carcass linear measurements increases the accuracy to predict the meat production.

Acknowledgments

The authors thank Fundação de Amparo à Pesquisa do Estado de São Paulo for the financial support (Project no. 09/15600-8).

References

  • AHDB Beef & Lamb. 2012. Understanding lambs & carcases for better returns. Agriculture and Horticulture Development Board, AHDB, Warwickshire, UK.
  • Beef + Lamb New Zealand Ltd. 2013. Available at: <http://beeflambnz.com/Global/Microsite/index.html> Accessed on: Mar. 2, 2010.
    » http://beeflambnz.com/Global/Microsite/index.html
  • Brady, A. S.; Belk, K. E.; LeValley, S. B.; Dalsted, N. L.; Scanga, J. A.; Tatum, J. D. and Smith, G. C. 2003. An evaluation of the lamb vision system as a predictor of lamb carcass red meat yield percentage. Journal of Animal Science 81:1488-1498.
  • Brasil. 1952. Decreto nº 30.691, de 29 de março de 1952. Aprova o novo Regulamento da Inspeção Industrial e Sanitária de Produtos de Origem Animal. Diário Oficial da União, Brasília, DF, 7 jul. 1952.
  • Brasil. 1990. Portaria nº 307, de 26 de dezembro de 1990. Aprova o Sistema Nacional de Tipificação de Carcaças Ovinas. Diário Oficial da União, Brasília, DF, 27 dez. 1990.
  • Brasil. 2008. Instrução Normativa nº 56, de 6 de novembro de 2008. Estabelecer os procedimentos gerais de Recomendações de Boas Práticas de Bem-Estar para Animais de Produção e de Interesse Econômico - REBEM, abrangendo os sistemas de produção e o transporte. Diário Oficial da União, Brasília, DF, 7 nov. 2008.
  • Carrasco, S.; Ripoll, G.; Sanz, A.; Álvarez-Rodríguez, J.; Panea, B.; Revilla, R. and Joy, M. 2009. Effect of feeding system on growth and carcass characteristics of Churra Tensina light lambs. Livestock Science 121:56-63.
  • Commission Regulation (EC) No 823/98. 1998. Commission Regulation (EC) No 823/98 of 20 April 1998 amending Regulation (EEC) No 461/93 laying down detailed rules for the Community scale for the classification of carcases of ovine animals.
  • Cunha, B. C. N.; Belk, K. E.; Scanga, J. A.; LeValley, S. B.; Tatum, J. D. and Smith, G. C. 2004. Development and validation of equations utilizing lamb vision system output to predict lamb carcass fabrication yields. Journal of Animal Science 82:2069-2076.
  • Díaz, M. T.; Cañeque, V.; Lauzuruca, S.; Velasco, S.; Ruíz de Huidobro, F. and Pérez, C. 2004. Prediction of suckling lamb carcass composition from objective and subjective carcass measurements. Meat Science 66:895-902.
  • Einarsson, E.; Eythórsdóttir, E.; Smith, C. R. and Jónmundsson, J. V. 2014. The ability of video image analysis to predict lean meat yield and EUROP score of lamb carcasses. Animal 8:1170-1177.
  • FAO - Food and Agriculture Organization of the United Nations. 2015. FAOSTAT database. Available at: <http://faostat3.fao.org/browse/Q/QA/E>. Accessed on: Feb. 15, 2016.
    » http://faostat3.fao.org/browse/Q/QA/E
  • Ferraz, J. B. S. and Felício, P. E. 2010. Production systems - An example from Brazil. Meat Science 84:238-243.
  • Hopkins D. L.; Safari, E.; Thompson, J. M. and Smith, C. R. 2004. Video image analysis in the Australian meat industry - precision and accuracy of predicting lean meat yield in lamb carcasses. Meat Science 67:269-274.
  • Johansen, J.; Aastveit, A. H.; Egelandsdal, B.; Kvaal, K. and Røe, M. 2006. Validation of the EUROP system for lamb classification in Norway; repeatability and accuracy of visual assessment and prediction of lamb carcass composition. Meat Science 74:497-509.
  • Jones, S. D. M.; Robertson, W. M. and Price, M. A. 1993. The assessment of saleable meat yield in lamb carcasses. p.190. In: International Congress on Meat Science and Technology, Calgary.
  • Lambe, N. R.; Navajas, E. A.; Bünger, L.; Fisher, A. V.; Roehe, R. and Simm, G. 2009. Prediction of lamb carcass composition and meat quality using combinations of post-mortem measurements. Meat Science 81:711-719.
  • Nsoso, S. J.; Young, M. J. and Beatson, P. R. 2000. A review of carcass conformation in sheep: assessment, genetic control and development. Small Ruminant Research 35:89-96.
  • Price, M. A. 1995. Development of carcass grading and classification systems. p.173-199. In: Quality and grading of carcasses of meat animals. Jones, S. D. M., ed. CRC Press, Boca Raton.
  • Rius-Vilarrasa, E.; Bünger, L.; Maltin, C.; Matthews, K. R. and Roehe, R. 2009. Evaluation of Video Image Analysis (VIA) technology to predict meat yield of sheep carcasses on-line under UK abattoir conditions. Meat Science 82:94-100.
  • Russel, A. J. F.; Doney, J. M. and Gunn, R. G. 1969. Subjective assessment of body fat in live sheep. Journal Agriculture Science 72:451-454.
  • Safari, E.; Hopkins, D. L. and Fogarty, N. M. 2001. Diverse lamb genotypes 4. Predicting the yield of saleable meat and high value trimmed cuts from carcass measurements. Meat Science 58:207-214.
  • Smith, G. C.; Carpenter, Z. L. and King, G. T. 1969. Ovine carcass cutability. Journal of Animal Science 29:272-282.
  • Stanford, K.; Woloschuk, C. M.; McClelland, L. A.; Jones, S. D. M. and Price, M. A. 1997. Comparison of objective external carcass measurements and subjective conformation scores for prediction of lamb carcass quality. Canadian Journal of Animal Science 77:217-223.
  • Stanford, K.; Jones, S. D. M. and Price, M. A. 1998. Methods of predicting lamb carcass composition: A review. Small Ruminant Research 29:241-254.
  • Zambrano-Bigiarini, M. 2012. hydroGOF: Goodness-of-fit functions for comparison of simulated and observed hydrological time series. R package version 0.3-5. Available at: <http://CRAN.R-project.org/package=hydroGOF>. Accessed on: Feb. 17, 2016.
    » http://CRAN.R-project.org/package=hydroGOF

Publication Dates

  • Publication in this collection
    Oct 2016

History

  • Received
    15 Mar 2016
  • Accepted
    04 Aug 2016
Sociedade Brasileira de Zootecnia Universidade Federal de Viçosa / Departamento de Zootecnia, 36570-900 Viçosa MG Brazil, Tel.: +55 31 3612-4602, +55 31 3612-4612 - Viçosa - MG - Brazil
E-mail: rbz@sbz.org.br