Genomic analysis of feed efficiency traits in beef cattle using random regression models.

Feed efficiency plays a major role in the overall profitability and sustainability of the beef cattle industry, as it is directly related to the reduction of the animal demand for input and methane emissions. Traditionally, the average daily feed intake and weight gain are used to calculate feed efficiency traits. However, feed efficiency traits can be analysed longitudinally using random regression models (RRMs), which allow fitting random genetic and environmental effects over time by considering the covariance pattern between the daily records. Therefore, the objectives of this study were to: (1) propose genomic evaluations for dry matter intake (DMI), body weight gain (BWG), residual feed intake (RFI) and residual weight gain (RWG) data collected during an 84-day feedlot test period via RRMs; (2) compare the goodness-of-fit of RRM using Legendre polynomials (LP) and B-spline functions; (3) evaluate the genetic parameters behaviour for feed efficiency traits and their implication for new selection strategies. The datasets were provided by the EMBRAPA-GENEPLUS beef cattle breeding program and included 2920 records for DMI, 2696 records for BWG and 4675 genotyped animals. Genetic parameters and genomic breeding values (GEBVs) were estimated by RRMs under ssGBLUP for Nellore cattle using orthogonal LPs and B-spline. Models were compared based on the deviance information criterion (DIC). The ranking of the average GEBV of each test week and the overall GEBV average were compared by the percentage of individuals in common and the Spearman correlation coefficient (top 1%, 5%, 10% and 100%). The highest goodness-of-fit was obtained with linear B-Spline function considering heterogeneous residual variance. The heritability estimates across the test period for DMI, BWG, RFI and RWG ranged from 0.06 to 0.21, 0.11 to 0.30, 0.03 to 0.26 and 0.07 to 0.27, respectively. DMI and RFI presented within-trait genetic correlations ranging from low to high magnitude across different performance test-day. In contrast, BWG and RWG presented negative genetic correlations between the first 3 weeks and the other days of performance tests. DMI and RFI presented a high-ranking similarity between the GEBV average of week eight and the overall GEBV average, with Spearman correlations and percentages of individuals selected in common ranging from 0.95 to 1.00 and 93 to 100, respectively. Week 11 presented the highest Spearman correlations (ranging from 0.94 to 0.98) and percentages of individuals selected in common (ranging from 85 to 94) of BWG and RWG with the average GEBV of the entire period of the test. In conclusion, the RRM using linear B-splines is a feasible alternative for the genomic evaluation of feed efficiency. Heritability estimates of DMI, RFI, BWG and RWG indicate enough additive genetic variance to achieve a moderate response to selection. A new selection strategy can be adopted by reducing the performance test to 56 days for DMI and RFI selection and 77 days for BWG and RWG selection.

Autoregressive and random regression test-day models for multiple lactations in genetic evaluation of Brazilian Holstein cattle.

Autoregressive (AR) and random regression (RR) models were fitted to test-day records from the first three lactations of Brazilian Holstein cattle with the objective of comparing their efficiency for national genetic evaluations. The data comprised 4,142,740 records of milk yield (MY) and somatic cell score (SCS) from 274,335 cows belonging to 2,322 herds. Although heritabilities were similar between models and traits, additive genetic variance estimates using AR were 7.0 (MY) and 22.2% (SCS) higher than those obtained from RR model. On the other hand, residual variances were lower in both traits when estimated through AR model. The rank correlation between EBV obtained from AR and RR models was 0.96 and 0.94 (MY) and 0.97 and 0.95 (SCS), respectively, for bulls (with 10 or more daughters) and cows. Estimated annual genetic gains for bulls (cows) obtained using AR were 46.11 (49.50) kg for MY and -0.019 (-0.025) score for SCS; whereas using RR these values were 47.70 (55.56) kg and -0.022 (-0.028) score. Akaike information criterion was lower for AR in both traits. Although AR model is more parsimonious, RR model assumes genetic correlations different from the unity within and across lactations. Thus, when these correlations are relatively high, these models tend to yield to similar predictions; otherwise, they will differ more and RR model would be theoretically sounder.

Genetic analysis of milk production traits of Tunisian Holsteins using random regression test-day model with Legendre polynomials.

OBJECTIVE: The objective of this study was to estimate genetic parameters of milk, fat, and protein yields within and across lactations in Tunisian Holsteins using a random regression test-day (TD) model. METHODS: A random regression multiple trait multiple lactation TD model was used to estimate genetic parameters in the Tunisian dairy cattle population. Data were TD yields of milk, fat, and protein from the first three lactations. Random regressions were modeled with third-order Legendre polynomials for the additive genetic, and permanent environment effects. Heritabilities, and genetic correlations were estimated by Bayesian techniques using the Gibbs sampler. RESULTS: All variance components tended to be high in the beginning and the end of lactations. Additive genetic variances for milk, fat, and protein yields were the lowest and were the least variable compared to permanent variances. Heritability values tended to increase with parity. Estimates of heritabilities for 305-d yield-traits were low to moderate, 0.14 to 0.2, 0.12 to 0.17, and 0.13 to 0.18 for milk, fat, and protein yields, respectively. Within-parity, genetic correlations among traits were up to 0.74. Genetic correlations among lactations for the yield traits were relatively high and ranged from 0.78±0.01 to 0.82±0.03, between the first and second parities, from 0.73±0.03 to 0.8±0.04 between the first and third parities, and from 0.82±0.02 to 0.84±0.04 between the second and third parities. CONCLUSION: These results are comparable to previously reported estimates on the same population, indicating that the adoption of a random regression TD model as the official genetic evaluation for production traits in Tunisia, as developed by most Interbull countries, is possible in the Tunisian Holsteins.

Partitioning of functional gene expression data using principal points.

BACKGROUND: DNA microarrays offer motivation and hope for the simultaneous study of variations in multiple genes. Gene expression is a temporal process that allows variations in expression levels with a characterized gene function over a period of time. Temporal gene expression curves can be treated as functional data since they are considered as independent realizations of a stochastic process. This process requires appropriate models to identify patterns of gene functions. The partitioning of the functional data can find homogeneous subgroups of entities for the massive genes within the inherent biological networks. Therefor it can be a useful technique for the analysis of time-course gene expression data. We propose a new self-consistent partitioning method of functional coefficients for individual expression profiles based on the orthonormal basis system. RESULTS: A principal points based functional partitioning method is proposed for time-course gene expression data. The method explores the relationship between genes using Legendre coefficients as principal points to extract the features of gene functions. Our proposed method provides high connectivity in connectedness after clustering for simulated data and finds a significant subsets of genes with the increased connectivity. Our approach has comparative advantages that fewer coefficients are used from the functional data and self-consistency of principal points for partitioning. As real data applications, we are able to find partitioned genes through the gene expressions found in budding yeast data and Escherichia coli data. CONCLUSIONS: The proposed method benefitted from the use of principal points, dimension reduction, and choice of orthogonal basis system as well as provides appropriately connected genes in the resulting subsets. We illustrate our method by applying with each set of cell-cycle-regulated time-course yeast genes and E. coli genes. The proposed method is able to identify highly connected genes and to explore the complex dynamics of biological systems in functional genomics.

Random Regression Models Are Suitable to Substitute the Traditional 305-Day Lactation Model in Genetic Evaluations of Holstein Cattle in Brazil.

The aim of this study was to compare two random regression models (RRM) fitted by fourth (RRM4) and fifth-order Legendre polynomials (RRM5) with a lactation model (LM) for evaluating Holstein cattle in Brazil. Two datasets with the same animals were prepared for this study. To apply test-day RRM and LMs, 262,426 test day records and 30,228 lactation records covering 305 days were prepared, respectively. The lowest values of Akaike's information criterion, Bayesian information criterion, and estimates of the maximum of the likelihood function (-2LogL) were for RRM4. Heritability for 305-day milk yield (305MY) was 0.23 (RRM4), 0.24 (RRM5), and 0.21 (LM). Heritability, additive genetic and permanent environmental variances of test days on days in milk was from 0.16 to 0.27, from 3.76 to 6.88 and from 11.12 to 20.21, respectively. Additive genetic correlations between test days ranged from 0.20 to 0.99. Permanent environmental correlations between test days were between 0.07 and 0.99. Standard deviations of average estimated breeding values (EBVs) for 305MY from RRM4 and RRM5 were from 11% to 30% higher for bulls and around 28% higher for cows than that in LM. Rank correlations between RRM EBVs and LM EBVs were between 0.86 to 0.96 for bulls and 0.80 to 0.87 for cows. Average percentage of gain in reliability of EBVs for 305-day yield increased from 4% to 17% for bulls and from 23% to 24% for cows when reliability of EBVs from RRM models was compared to those from LM model. Random regression model fitted by fourth order Legendre polynomials is recommended for genetic evaluations of Brazilian Holstein cattle because of the higher reliability in the estimation of breeding values.

Genetic Parameters for Milk Yield and Lactation Persistency Using Random Regression Models in Girolando Cattle.

A total of 32,817 test-day milk yield (TDMY) records of the first lactation of 4,056 Girolando cows daughters of 276 sires, collected from 118 herds between 2000 and 2011 were utilized to estimate the genetic parameters for TDMY via random regression models (RRM) using Legendre's polynomial functions whose orders varied from 3 to 5. In addition, nine measures of persistency in milk yield (PSi) and the genetic trend of 305-day milk yield (305MY) were evaluated. The fit quality criteria used indicated RRM employing the Legendre's polynomial of orders 3 and 5 for fitting the genetic additive and permanent environment effects, respectively, as the best model. The heritability and genetic correlation for TDMY throughout the lactation, obtained with the best model, varied from 0.18 to 0.23 and from -0.03 to 1.00, respectively. The heritability and genetic correlation for persistency and 305MY varied from 0.10 to 0.33 and from -0.98 to 1.00, respectively. The use of PS7 would be the most suitable option for the evaluation of Girolando cattle. The estimated breeding values for 305MY of sires and cows showed significant and positive genetic trends. Thus, the use of selection indices would be indicated in the genetic evaluation of Girolando cattle for both traits.

Computational investigation of stochastic Zika virus optimal control model using Legendre spectral method.

This study presents a computational investigation of a stochastic Zika virus along with optimal control model using the Legendre spectral collocation method (LSCM). By accumulation of stochasticity into the model through the proposed stochastic differential equations, we appropriating the random fluctuations essential in the progression and disease transmission. The stability, convergence and accuracy properties of the LSCM are conscientiously analyzed and also demonstrating its strength for solving the complex epidemiological models. Moreover, the study evaluates the various control strategies, such as treatment, prevention and treatment pesticide control, and identifies optimal combinations that the intervention costs and also minimize the proposed infection rates. The basic properties of the given model, such as the reproduction number, were determined with and without the presence of the control strategies. For R 0 < 0 , the model satisfies the disease-free equilibrium, in this case the disease die out after some time, while for R 0 > 1 , then endemic equilibrium is satisfied, in this case the disease spread in the population at higher scale. The fundamental findings acknowledge the significant impact of stochastic phonemes on the robustness and effectiveness of control strategies that accelerating the need for cost-effective and multi-faceted approaches. In last the results provide the valuable insights for public health department to enabling more impressive mitigation of Zika virus outbreaks and management in real-world scenarios.

Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes.

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21-0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal.

Optimal control of nonlinear fractional order delay systems governed by Fredholm integral equations based on a new fractional derivative operator.

We aim to develop a direct transcription approach for solving a notable category of optimal control problems governed by nonlinear fractional Fredholm integral equations having delays in both input and output signals. The foundation of the new methodology is based on a multi domains decomposition scheme by utilizing the fractional-order Legendre functions. A new fractional derivative operator associated with the fractional basis is introduced by using the Caputo fractional derivative operator. With the use of derivative and delay operators, one can transform the dynamical system related to the fractional control problem into a new system containing algebraic equations. A wide variety of challenging test problems are studied to provide a detailed explanation of the designed approach.

Deconvolution of Tc-99m-Mercaptoacetyltriglycine Renograms with the Concomitant Use of a Sparse Legendre Polynomial Representation and the Moore-Penrose Pseudo-inverse.

OBJECTIVES: This study aimed to introduce an improved deconvolution technique for Tc-99m-mercaptoacetyltriglycine renograms based on the combination of a sparse Legendre polynomial representation and the Moore-Penrose inversion matrix (LG). This method reduces the effect of noise on the measurement of renal retention function transit time (TT). METHODS: The stability and accuracy of the proposed method were tested using a renal database containing Monte Carlo-simulated studies and real adult patient data. Two clinical parameters, namely, split function (SF) and mean TT (meanTT), obtained with LG were compared with values calculated with the established method that combines matrix deconvolution and a three-point linear smoothing (F121) as recommended by the 2008 International Scientific Committee of Radionuclides in Nephrourology consensus on renal TT measurements. RESULTS: For simulated data, the root mean square error (RMSE) between the theoretical non-noisy renal retention curve (RRC) and the results of the deconvolution methods applied to the noisy RRC were up to two times lower with LG (p<0.001). The RMSE of the reconvoluted renogram and the theoretical one was also lower for LG (p<0.001) and showed better preservation of the original signal. The SF was neither improved nor degraded by the proposed method. For patient data, no statistically significant difference was found between the SF for the LG method compared with the database values, and the meanTT better agreed with the physician's diagnosis than the matrix or clinical software (Hermes) outputs. A visual improvement of the RRC was also observed. CONCLUSION: By combining the sparse Legendre representation of the renogram curves and the Moore-Penrose matrix inverse techniques, we obtained improved noise reduction in the deconvoluted data, leading to better elimination of non-physiological signals -as negative values- and the avoidance of the smear effect of conventional smoothing on the vascular peak, which both influenced the meanTT measurement.

Design and experimental evaluation of block-pulse functions and Legendre polynomials observer for attitude-heading reference system.

The main purpose of this paper is design and implementation of a new linear observer for an attitude and heading reference system (AHRS), which includes three-axis accelerometers, gyroscopes, and magnetometers in the presence of sensors and modeling uncertainties. Since the increase of errors over time is the main difficulty of low-cost micro electro mechanical systems (MEMS) sensors producing instable on-off bias, scale factor (SF), nonlinearity and random walk errors, development of a high-precision observer to improve the accuracy of MEMS-based navigation systems is considered. First, the duality between controller and estimator in a linear system is presented as the base of design method. Next, Legendre polynomials together with block-pulse functions are applied for the solution of a common linear time-varying control problem. Through the duality theory, the obtained control solution results in the block-pulse functions and Legendre polynomials observer (BPLPO). According to product properties of the hybrid functions in addition to the operational matrices of integration, the optimal control problem is simplified to some algebraic equations which particularly fit with low-cost implementations. The improved performance of the MEMS AHRS owing to implementation of BPLPO has been assessed through vehicle field tests in urban area compared with the extended Kalman filter (EKF).

Legendre Polynomials: A Fully Automatic Method for Noise Reduction in 99mTc-Mercaptoacetyltriglycine Renogram Analysis.

Our purpose was to develop a fully automatic method to deal with the presence of high levels of noise interfering with quantitative analysis of fast, dynamic mercaptoacetyltriglycine renogram images. Methods: A method based on Legendre polynomials to fit and filter time-activity curves was proposed. The method was applied to a renal database that contains Monte Carlo (MC)-simulated studies and real adult patient data. Clinically relevant parameters such as relative function, time to maximum uptake (Tmax), and half-emptying time (T1/2) were obtained with the proposed method, the 1-2-1 filter (F121) method recommended in the 2018 guidelines of the European Association of Nuclear Medicine, and a state-of-the-art commercial software program (Hermes) currently used in routine nuclear medicine. Results: The root mean squared error between reference time-activity curves and the same curves with Poisson noise added was about 2 times lower for the Legendre method than for F121. The left relative function for MC and patient data was statistically equivalent for Hermes, Legendre, and F121 (P < 0.001). For MC studies, the Legendre technique performed better that the Hermes method regarding the known values of Tmax (P < 0.05), and the T1/2 determination was significantly improved (P < 0.05). For patient data, the Legendre and F121 methods were less influenced by noise in the data than the Hermes method, particularly for T1/2. Conclusion: In dynamic nuclear medicine imaging, Legendre polynomials appear to be a promising, fully automatic noise-removal tool that is routinely applicable, accurate, and robust.

A numerical solution by alternative Legendre polynomials on a model for novel coronavirus (COVID-19).

Coronavirus disease (COVID-19) is an infectious disease caused by a newly discovered coronavirus. This paper provides a numerical solution for the mathematical model of the novel coronavirus by the application of alternative Legendre polynomials to find the transmissibility of COVID-19. The mathematical model of the present problem is a system of differential equations. The goal is to convert this system to an algebraic system by use of the useful property of alternative Legendre polynomials and collocation method that can be solved easily. We compare the results of this method with those of the Runge-Kutta method to show the efficiency of the proposed method.

Impact of the Order of Legendre Polynomials in Random Regression Model on Genetic Evaluation for Milk Yield in Dairy Cattle Population.

The random regression test-day model has become the most commonly adopted model for routine genetic evaluations in dairy populations, which allows accurately accounting for genetic and environmental effects over lactation. The objective of this study was to explore appropriate random regression test-day models for genetic evaluation of milk yield in a Holstein population with a relatively small size, which is the common situation in regional or independent breeding companies to preform genetic evaluation. Data included 419,567 test-day records from 54,417 cows from the first lactation. Variance components and breeding values were estimated using a random regression test-day model with different orders (from first to fifth) of Legendre polynomials (LP) and accounted for homogeneous or heterogeneous residual variance across the lactation. Models were compared based on Akaike information criterion (AIC), Bayesian information criterion (BIC), and predictive ability. In general, models with a higher order of LP showed better goodness of fit based on AIC and BIC values. However, models with third, fourth, and fifth order of LP led to similar estimates of genetic parameters and predictive ability. Models with assumption of heterogeneous residual variances achieved better goodness of fit but yielded similar predictive ability compared with those with assumption of homogeneous residual variances. Therefore, a random regression model with third order of LP is suggested to be an appropriate model for genetic evaluation of milk yield in local Chinese Holstein populations.

Genetic parameters for faecal egg count, packed-cell volume and body-weight in Santa Inês lambs.

Worm infection is one of the main factors responsible for economic losses in sheep breeding in Brazil. Random regression analysis was used to estimate genetic parameters for the factors faecal egg-count (FEC), packed-cell volume (PCV) and body weight (BW) in Santa Inês lambs. Data from 119 female, offspring of nine rams, were collected between December, 2005 and December, 2006, from the experimental flock of Embrapa Tabuleiros Costeiros, the Brazilian Agricultural Research Corporation located in Frei Paulo, SE, Brazil. After weaning, females were drenched until the faecal egg count had dropped to zero. Two natural challenges were undertaken. FEC heritability was extremely variable, this increasing from 0.04 to 0.27 in the first challenge and from 0.01 to 0.52 during the second. PCV heritability peaks were 0.31 and 0.12 in the first and second challenges, respectively. In the second challenge, BW heritability was close to 0.90. The genetic correlations among these traits did not differ from zero. There is the possibility of increasing parasite resistance in Santa Inês by selecting those animals with lower FEC. Selection to increase resistance will not adversely affect lamb-growth, although lambs with a slow growth-rate may be more susceptible to infection.

Genetic parameters for stayability to consecutive calvings in Zebu cattle.

Longer-lived cows tend to be more profitable and the stayability trait is a selection criterion correlated to longevity. An alternative to the traditional approach to evaluate stayability is its definition based on consecutive calvings, whose main advantage is the more accurate evaluation of young bulls. However, no study using this alternative approach has been conducted for Zebu breeds. Therefore, the objective of this study was to compare linear random regression models to fit stayability to consecutive calvings of Guzerá, Nelore and Tabapuã cows and to estimate genetic parameters for this trait in the respective breeds. Data up to the eighth calving were used. The models included the fixed effects of age at first calving and year-season of birth of the cow and the random effects of contemporary group, additive genetic, permanent environmental and residual. Random regressions were modeled by orthogonal Legendre polynomials of order 1 to 4 (2 to 5 coefficients) for contemporary group, additive genetic and permanent environmental effects. Using Deviance Information Criterion as the selection criterion, the model with 4 regression coefficients for each effect was the most adequate for the Nelore and Tabapuã breeds and the model with 5 coefficients is recommended for the Guzerá breed. For Guzerá, heritabilities ranged from 0.05 to 0.08, showing a quadratic trend with a peak between the fourth and sixth calving. For the Nelore and Tabapuã breeds, the estimates ranged from 0.03 to 0.07 and from 0.03 to 0.08, respectively, and increased with increasing calving number. The additive genetic correlations exhibited a similar trend among breeds and were higher for stayability between closer calvings. Even between more distant calvings (second v. eighth), stayability showed a moderate to high genetic correlation, which was 0.77, 0.57 and 0.79 for the Guzerá, Nelore and Tabapuã breeds, respectively. For Guzerá, when the models with 4 or 5 regression coefficients were compared, the rank correlations between predicted breeding values for the intercept were always higher than 0.99, indicating the possibility of practical application of the least parameterized model. In conclusion, the model with 4 random regression coefficients is recommended for the genetic evaluation of stayability to consecutive calvings in Zebu cattle.

A New Computational Model of High-Order Stochastic Simulation Based on Spatial Legendre Moments.

Multiple-point simulations have been introduced over the past decade to overcome the limitations of second-order stochastic simulations in dealing with geologic complexity, curvilinear patterns, and non-Gaussianity. However, a limitation is that they sometimes fail to generate results that comply with the statistics of the available data while maintaining the consistency of high-order spatial statistics. As an alternative, high-order stochastic simulations based on spatial cumulants or spatial moments have been proposed; however, they are also computationally demanding, which limits their applicability. The present work derives a new computational model to numerically approximate the conditional probability density function (cpdf) as a multivariate Legendre polynomial series based on the concept of spatial Legendre moments. The advantage of this method is that no explicit computations of moments (or cumulants) are needed in the model. The approximation of the cpdf is simplified to the computation of a unified empirical function. Moreover, the new computational model computes the cpdfs within a local neighborhood without storing the high-order spatial statistics through a predefined template. With this computational model, the algorithm for the estimation of the cpdf is developed in such a way that the conditional cumulative distribution function (ccdf) can be computed conveniently through another recursive algorithm. In addition to the significant reduction of computational cost, the new algorithm maintains higher numerical precision compared to the original version of the high-order simulation. A new method is also proposed to deal with the replicates in the simulation algorithm, reducing the impacts of conflicting statistics between the sample data and the training image (TI). A brief description of implementation is provided and, for comparison and verification, a set of case studies is conducted and compared with the results of the well-established multi-point simulation algorithm, filtersim. This comparison demonstrates that the proposed high-order simulation algorithm can generate spatially complex geological patterns while also reproducing the high-order spatial statistics from the sample data.

Recurrence relations for orthogonal polynomials for PDEs in polar and cylindrical geometries.

This paper introduces two families of orthogonal polynomials on the interval (-1,1), with weight function [Formula: see text]. The first family satisfies the boundary condition [Formula: see text], and the second one satisfies the boundary conditions [Formula: see text]. These boundary conditions arise naturally from PDEs defined on a disk with Dirichlet boundary conditions and the requirement of regularity in Cartesian coordinates. The families of orthogonal polynomials are obtained by orthogonalizing short linear combinations of Legendre polynomials that satisfy the same boundary conditions. Then, the three-term recurrence relations are derived. Finally, it is shown that from these recurrence relations, one can efficiently compute the corresponding recurrences for generalized Jacobi polynomials that satisfy the same boundary conditions.

Genetic parameters for growth characteristics of free-range chickens under univariate random regression models.

Repeated measures from the same individual have been analyzed by using repeatability and finite dimension models under univariate or multivariate analyses. However, in the last decade, the use of random regression models for genetic studies with longitudinal data have become more common. Thus, the aim of this research was to estimate genetic parameters for body weight of four experimental chicken lines by using univariate random regression models. Body weight data from hatching to 84 days of age (n = 34,730) from four experimental free-range chicken lines (7P, Caipirão da ESALQ, Caipirinha da ESALQ and Carijó Barbado) were used. The analysis model included the fixed effects of contemporary group (gender and rearing system), fixed regression coefficients for age at measurement, and random regression coefficients for permanent environmental effects and additive genetic effects. Heterogeneous variances for residual effects were considered, and one residual variance was assigned for each of six subclasses of age at measurement. Random regression curves were modeled by using Legendre polynomials of the second and third orders, with the best model chosen based on the Akaike Information Criterion, Bayesian Information Criterion, and restricted maximum likelihood. Multivariate analyses under the same animal mixed model were also performed for the validation of the random regression models. The Legendre polynomials of second order were better for describing the growth curves of the lines studied. Moderate to high heritabilities (h(2) = 0.15 to 0.98) were estimated for body weight between one and 84 days of age, suggesting that selection for body weight at all ages can be used as a selection criteria. Genetic correlations among body weight records obtained through multivariate analyses ranged from 0.18 to 0.96, 0.12 to 0.89, 0.06 to 0.96, and 0.28 to 0.96 in 7P, Caipirão da ESALQ, Caipirinha da ESALQ, and Carijó Barbado chicken lines, respectively. Results indicate that genetic gain for body weight can be achieved by selection. Also, selection for body weight at 42 days of age can be maintained as a selection criterion.

Estimation of genetic parameters for test-day milk yield in Girolando cows using a random regression model / [Estimação de parâmetros genéticos para produção de leite no dia do teste em vacas da raça Girolando utilizando modelo de regressão aleatória]

The objective of this study was to estimate the components of variance and genetic parameters of test-day milk yield in first lactation Girolando cows, using a random regression model. A total of 126,892 test-day milk yield (TDMY) records of 15,351 first-parity Holstein, Gyr, and Girolando breed cows were used, obtained from the Associação Brasileira dos Criadores de Girolando. To estimate the components of (co) variance, the additive genetic functions and permanent environmental covariance were estimated by random regression in three functions: Wilmink, Legendre Polynomials (third order) and Linear spline Polynomials (three knots). The Legendre polynomial function showed better fit quality. The genetic and permanent environment variances for TDMY ranged from 2.67 to 5.14 and from 9.31 to 12.04, respectively. Heritability estimates gradually increased from the beginning (0.13) to mid-lactation (0.19). The genetic correlations between the days of the control ranged from 0.37 to 1.00. The correlations of permanent environment followed the same trend as genetic correlations. The use of Legendre polynomials via random regression model can be considered as a good tool for estimating genetic parameters for test-day milk yield records.(AU)

O objetivo deste estudo foi estimar os componentes de variância e os parâmetros genéticos da produção de leite no dia do teste (TDMY) em vacas Girolando de primeira lactação, usando modelo de regressão aleatória. Foram utilizados 126.892 registros de produção de leite no dia controle de 15.351 vacas primíparas das raças Holandesa, Gir e Girolando, obtidas na Associação Brasileira dos Criadores de Girolando. Para estimar os componentes de (co) variância, as funções genéticas aditivas e de covariância ambiental permanente foram estimadas por regressão aleatória em três funções: Wilmink, polinômios de Legendre (terceira ordem) e polinômios splines lineares (três nós). A função polinomial de Legendre apresentou melhor qualidade de ajuste. As variâncias genéticas e de ambiente permanente para produção de leite no dia do controle variaram de 2,67 a 5,14 e de 9,31 a 12,04, respectivamente. As estimativas de herdabilidade aumentaram gradativamente do início (0,13) para o meio da lactação (0,19). As correlações genéticas entre os dias do controle variaram de 0,37 a 1,00. As correlações de ambiente permanente seguiram a mesma tendência das correlações genéticas. A utilização dos polinômios de Legendre via modelos de regressão aleatória pode ser considerada como uma boa ferramenta para estimação de parâmetros genéticos da produção de leite no dia do teste.(AU)