Genetic parameters estimates for milk, fat and protein yield analyzed by test day models for Murrah buffaloes in Brazil

The objectives of this study were to estimate genetic parameters for test-day milk, fat and protein yields, in Murrah buffaloes. In this study 4,757 complete lactations of Murrah buffaloes were analyzed. The (co) variance components were estimated by restricted maximum likelihood using MTDFREML software. The bi-trait animal test-day models included genetic additive direct and permanent environment effects, as random effects, and the fixed effects of contemporary group (herds-year-month of control) and age of the cow at calving as linear and quadratic covariable. The heritability estimate at first control was 0.19, increased until the third control (0.24), decreasing thereafter, reaching the lowest value at the ninth control (0.09). The highest heritability estimates for fat and protein yield were 0.23 (first control) and 0.33 (third control), respectively. For milk yield, genetic and phenotypic correlation estimates ranged from 0.37 to 0.99 and from 0.52 to 0.94, respectively. Genetic correlations were higher than phenotypic ones. For fat and protein yields, genetic correlation estimates ranged from 0.42 to 0.97.

Genotypic evaluation of Dent composite (Zea mays L) half-sib progenies in late harvest crop conditions

The present work was conducted at FCAV-UNESP, Campus de Jaboticabal, in late harvest crop time of 1991. One hundred interpopulational hybrids obtained from a top-cross between Dent Composite and Flint Composite populations were evaluated. The analyzed characteristics were: Spodoptera frugiperda damages, lodging, damaged ears percentage, and productivity. Estimates of heritability, variance between progenies and genetical gain for each character of available population were determined. Estimates of the progenies variance, like those from heritability for the different analysed characters, indicate that there is an adequate genetical variance for the utilization of that material on subsequent genetic breeding programs, though allowing genetical gains, on the following selection cycles for the characteristics: grain weight, lodging and damages caused by armyworm larvae, whereas for the other characteristics the obtained results were not satisfactory.

Bounds on policy relevant parameters with discrete policy variation

When estimating policy parameters, also known as treatment effects, the assignment to treatment mechanism almost always causes endogeneity and thus bias many of these policy parameters estimates. Additionally, heterogeneity in program impacts is more likely to be the norm than the exception for most social programs. In situations where these issues are present, the Marginal Treatment Effect (MTE) parameter estimation makes use of an instrument to avoid assignment bias and simultaneously to account for heterogeneous effects throughout individuals. Although this parameter is point identified in the literature, the assumptions required for identification may be strong. Given that, we use weaker assumptions in order to partially identify the MTE, i.e. to stablish a methodology for MTE bounds estimation, implementing it computationally and showing results from Monte Carlo simulations. The partial identification we perfom requires the MTE to be a monotone function over the propensity score, which is a reasonable assumption on several economics' examples, and the simulation results shows it is possible to get informative even in restricted cases where point identification is lost. Additionally, in situations where estimated bounds are not informative and the traditional point identification is lost, we suggest a more generic method to point estimate MTE using the Moore-Penrose Pseudo-Invese Matrix, achieving better results than traditional methods.

Genetic parameters of total milk yield and factors describing the shape of lactation curve in dairy buffaloes

The objective of this study was to apply factor analysis to describe lactation curves in dairy buffaloes in order to estimate the phenotypic and genetic association between common latent factors and cumulative milk yield. A total of 31 257 monthly test-day milk yield records from buffaloes belonging to herds located in the state of São Paulo were used to estimate two common latent factors, which were then analysed in a multi-trait animal model for estimating genetic parameters. Estimates of (co)variance components for the two common latent factors and cumulated 270-d milk yield were obtained by Bayesian inference using a multiple trait animal model. Contemporary group, number of milkings per day (two levels) and age of buffalo cow at calving (linear and quadratic) as covariate were included in the model as fixed effects. The additive genetic, permanent environmental and residual effects were included as random effects. The first common latent factor (F1) was associated with persistency of lactation and the second common latent factor (F2) with the level of production in early lactation. Heritability estimates for Fl and F2 were 0.12 and 0.07, respectively. Genetic correlation estimates between El and F2 with cumulative milk yield were positive and moderate (0.63 and 0.52). Multivariate statistics employing factor analysis allowed the extraction of two variables (latent factors) that described the shape of the lactation curve. It is expected that the response to selection to increase lactation persistency is higher than the response obtained from selecting animals to increase lactation peak. Selection for higher total milk yield would result in a favourable correlated response to increase the level of production in early lactation and the lactation persistency.

Topographical body fat distribution links to amino acid and lipid metabolism in healthy non-obese women.

Visceral adiposity is increasingly recognized as a key condition for the development of obesity related disorders, with the ratio between visceral adipose tissue (VAT) and subcutaneous adipose tissue (SAT) reported as the best correlate of cardiometabolic risk. In this study, using a cohort of 40 obese females (age: 25-45 y, BMI: 28-40 kg/m(2)) under healthy clinical conditions and monitored over a 2 weeks period we examined the relationships between different body composition parameters, estimates of visceral adiposity and blood/urine metabolic profiles. Metabonomics and lipidomics analysis of blood plasma and urine were employed in combination with in vivo quantitation of body composition and abdominal fat distribution using iDXA and computerized tomography. Of the various visceral fat estimates, VAT/SAT and VAT/total abdominal fat ratios exhibited significant associations with regio-specific body lean and fat composition. The integration of these visceral fat estimates with metabolic profiles of blood and urine described a distinct amino acid, diacyl and ether phospholipid phenotype in women with higher visceral fat. Metabolites important in predicting visceral fat adiposity as assessed by Random forest analysis highlighted 7 most robust markers, including tyrosine, glutamine, PC-O 44∶6, PC-O 44∶4, PC-O 42∶4, PC-O 40∶4, and PC-O 40∶3 lipid species. Unexpectedly, the visceral fat associated inflammatory profiles were shown to be highly influenced by inter-days and between-subject variations. Nevertheless, the visceral fat associated amino acid and lipid signature is proposed to be further validated for future patient stratification and cardiometabolic health diagnostics.

Dry months in the agricultural region of Ribeirão Preto, state of São Paulo-Brazil: an study based on the extreme value theory

The application of the Extreme Value Theory (EVT) to model the probability of occurrence of extreme low Standardized Precipitation Index (SPI) values leads to an increase of the knowledge related to the occurrence of extreme dry months. This sort of analysis can be carried out by means of two approaches: the block maxima (BM; associated with the General Extreme Value distribution) and the peaks-over-threshold (POT; associated with the Generalized Pareto distribution). Each of these procedures has its own advantages and drawbacks. Thus, the main goal of this study is to compare the performance of BM and POT in characterizing the probability of occurrence of extreme dry SPI values obtained from the weather station of Ribeirão Preto-SP (1937-2012). According to the goodness-of-fit tests, both BM and POT can be used to assess the probability of occurrence of the aforementioned extreme dry SPI monthly values. However, the scalar measures of accuracy and the return level plots indicate that POT provides the best fit distribution. The study also indicated that the uncertainties in the parameters estimates of a probabilistic model should be taken into account when the probability associated with a severe/extreme dry event is under analysis.

Inférence robuste à la présence des valeurs aberrantes dans les enquêtes

Cette thèse comporte trois articles dont un est publié et deux en préparation. Le sujet central de la thèse porte sur le traitement des valeurs aberrantes représentatives dans deux aspects importants des enquêtes que sont : l’estimation des petits domaines et l’imputation en présence de non-réponse partielle. En ce qui concerne les petits domaines, les estimateurs robustes dans le cadre des modèles au niveau des unités ont été étudiés. Sinha & Rao (2009) proposent une version robuste du meilleur prédicteur linéaire sans biais empirique pour la moyenne des petits domaines. Leur estimateur robuste est de type «plugin», et à la lumière des travaux de Chambers (1986), cet estimateur peut être biaisé dans certaines situations. Chambers et al. (2014) proposent un estimateur corrigé du biais. En outre, un estimateur de l’erreur quadratique moyenne a été associé à ces estimateurs ponctuels. Sinha & Rao (2009) proposent une procédure bootstrap paramétrique pour estimer l’erreur quadratique moyenne. Des méthodes analytiques sont proposées dans Chambers et al. (2014). Cependant, leur validité théorique n’a pas été établie et leurs performances empiriques ne sont pas pleinement satisfaisantes. Ici, nous examinons deux nouvelles approches pour obtenir une version robuste du meilleur prédicteur linéaire sans biais empirique : la première est fondée sur les travaux de Chambers (1986), et la deuxième est basée sur le concept de biais conditionnel comme mesure de l’influence d’une unité de la population. Ces deux classes d’estimateurs robustes des petits domaines incluent également un terme de correction pour le biais. Cependant, ils utilisent tous les deux l’information disponible dans tous les domaines contrairement à celui de Chambers et al. (2014) qui utilise uniquement l’information disponible dans le domaine d’intérêt. Dans certaines situations, un biais non négligeable est possible pour l’estimateur de Sinha & Rao (2009), alors que les estimateurs proposés exhibent un faible biais pour un choix approprié de la fonction d’influence et de la constante de robustesse. Les simulations Monte Carlo sont effectuées, et les comparaisons sont faites entre les estimateurs proposés et ceux de Sinha & Rao (2009) et de Chambers et al. (2014). Les résultats montrent que les estimateurs de Sinha & Rao (2009) et de Chambers et al. (2014) peuvent avoir un biais important, alors que les estimateurs proposés ont une meilleure performance en termes de biais et d’erreur quadratique moyenne. En outre, nous proposons une nouvelle procédure bootstrap pour l’estimation de l’erreur quadratique moyenne des estimateurs robustes des petits domaines. Contrairement aux procédures existantes, nous montrons formellement la validité asymptotique de la méthode bootstrap proposée. Par ailleurs, la méthode proposée est semi-paramétrique, c’est-à-dire, elle n’est pas assujettie à une hypothèse sur les distributions des erreurs ou des effets aléatoires. Ainsi, elle est particulièrement attrayante et plus largement applicable. Nous examinons les performances de notre procédure bootstrap avec les simulations Monte Carlo. Les résultats montrent que notre procédure performe bien et surtout performe mieux que tous les compétiteurs étudiés. Une application de la méthode proposée est illustrée en analysant les données réelles contenant des valeurs aberrantes de Battese, Harter & Fuller (1988). S’agissant de l’imputation en présence de non-réponse partielle, certaines formes d’imputation simple ont été étudiées. L’imputation par la régression déterministe entre les classes, qui inclut l’imputation par le ratio et l’imputation par la moyenne sont souvent utilisées dans les enquêtes. Ces méthodes d’imputation peuvent conduire à des estimateurs imputés biaisés si le modèle d’imputation ou le modèle de non-réponse n’est pas correctement spécifié. Des estimateurs doublement robustes ont été développés dans les années récentes. Ces estimateurs sont sans biais si l’un au moins des modèles d’imputation ou de non-réponse est bien spécifié. Cependant, en présence des valeurs aberrantes, les estimateurs imputés doublement robustes peuvent être très instables. En utilisant le concept de biais conditionnel, nous proposons une version robuste aux valeurs aberrantes de l’estimateur doublement robuste. Les résultats des études par simulations montrent que l’estimateur proposé performe bien pour un choix approprié de la constante de robustesse.

Evaluation of water beef buffalo from birth to two years using different growth curves

The buffalo is a domestic animal species of growing world-wide importance. Research to improve genetic improvement programs is important to maintain the productivity of buffalo. The objective this research was to evaluate the growth of Brazilian buffalo to two years of age with different growth curves. Growth curves consolidate the information contained in the weight-age data into three or four biologically meaningful parameters. The data included 31,452 weights at birth and 120, 205, 365, 550 and 730 days of buffalo (n = 5,178) raised on pasture without supplementation. Logistic, Gompertz, quadratic logarithmic, and linear hyperbolic curves (designated L, G, QL, and LH, respectively) were fitted to the data by using proc NUN of SAS (SAS Institute, Inc., Cary, NC, USA). The parameters estimates for L [WT= A * (((1 + exp (-k * AGE)))**-m)] were A = 865.1 +/- 5.42; k= 0.0028 +/- 0.00002; M= 3.808 +/- 0.007; R(2) = 0.95. For G [WT= A * exp (-b * exp (-k * age)] the parameters estimates were A= 967.6 +/- 7.23; k = 0.00217 +/- 0.000015; b = -2.8152 +/- 0.00532. For QL [WT= A + b*age + k*(age*age) + m*log (age)] parameters estimates were A= 37.41 +/- 0.48; k= 0.00019 +/- 6.4E(-6); b= 0.539 +/- 0.006; m= 2.32 +/- 0.23; R(2)=0.96. For LH [WT= A + b*AGE + k*(1/AGE)] the parameters estimates were A= 23.15 +/- 0.44; k=15.16 +/- 0.66; b= 0.707 +/- 0.001; R(2)= 0.96. Each of these curves fit these data equally well and could be used for characterizing growth to two years in beef buffalo.

Comparação do ajuste dos modelos de morgan-mercer-flodin e de johnson-mehl-avrami na reação de precipitação na liga cu-3%al-5%ag

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelos matemáticos utilizados para descrever curvas de crescimento em aves aplicados ao melhoramento genético animal

A utilização de funções matemáticas para descrever o crescimento animal é antiga. Elas permitem resumir informações em alguns pontos estratégicos do desenvolvimento ponderal e descrever a evolução do peso em função da idade do animal. Também é possível comparar taxas de crescimento de diferentes indivíduos em estados fisiológicos equivalentes. Os modelos de curvas de crescimento mais utilizados na avicultura são os derivados da função Richards, pois apresentam parâmetros que possibilitam interpretação biológica e portanto podem fornecer subsídios para seleção de uma determinada forma da curva de crescimento em aves. Também pode-se utilizar polinômios segmentados para descrever as mudanças de tendência da curva de crescimento animal. Entretanto, existem importantes fatores de variação para os parâmetros das curvas, como a espécie, o sistema de criação, o sexo e suas interações. A adequação dos modelos pode ser verificada pelos valores do coeficiente de determinação (R2), do quadrado médio do resíduo (QM res), do erro de predição médio (EPm), da facilidade de convergência dos dados e pela possibilidade de interpretação biológica dos parâmetros. Estudos envolvendo modelagem e descrição da curva de crescimento e seus componentes são amplamente discutidos na literatura. Porém, programas de seleção que visem a progressos genéticos para a forma da curva não são mencionados. A importância da avaliação dos parâmetros dos modelos de curvas de crescimento é ainda mais relevante já que os maiores ganhos genéticos para peso estão relacionados com seleção para pesos em idades próximas ao ponto de inflexão. A seleção para precocidade pode ser auxiliada com base nos parâmetros do modelo associados à variáveis que descrevem esta característica genética dos animais. Esses parâmetros estão relacionados a importantes características produtivas e reprodutivas e apresentam magnitudes diferentes, de acordo com a espécie, o sexo e o modelo utilizados na avaliação. Outra metodologia utilizada são os modelos de regressão aleatória, permitindo mudanças graduais nas covariâncias entre idades ao longo do tempo e predizendo variâncias e covariâncias em pontos contidos ao longo da trajetória estudada. A utilização de modelos de regressões aleatórias traz como vantagem a separação da variação da curva de crescimento fenotípica em seus diferentes efeitos genético aditivo e de ambiente permanente individual, mediante a determinação dos coeficientes de regressão aleatórios para esses diferentes efeitos. Além disto, não há necessidade de utilizar fatores de ajuste para a idade. Esta revisão teve por objetivos levantar os principais modelos matemáticos frequentistas utilizados no estudo de curvas de crescimento de aves, com maior ênfase nos empregados com a finalidade de estimar parâmetros genéticos e fenotípicos.

Estimativas de parâmetros genéticos de características de carcaça, medidas por ultrassonografia, e suas relações com escores visuais na raça Nelore

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Parâmetros genéticos em populações de soja derivadas de cruzamentos simples e múltiplos, conduzidas por três diferentes métodos de avanço de gerações

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Polynomials to model the growth of young bulls in performance tests

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Customer Evolution in Sales Channel Migration

This research has been triggered by an emergent trend in customer behavior: customers have rapidly expanded their channel experiences and preferences beyond traditional channels (such as stores) and they expect the company with which they do business to have a presence on all these channels. This evidence has produced an increasing interest in multichannel customer behavior and it has motivated several researchers to study the customers’ channel choices dynamics in multichannel environment. We study how the consumer decision process for channel choice and response to marketing communications evolves for a cohort of new customers. We assume a newly acquired customer’s decisions are described by a “trial” model, but the customer’s choice process evolves to a “post-trial” model as the customer learns his or her preferences and becomes familiar with the firm’s marketing efforts. The trial and post-trial decision processes are each described by different multinomial logit choice models, and the evolution from the trial to post-trial model is determined by a customer-level geometric distribution that captures the time it takes for the customer to make the transition. We utilize data for a major retailer who sells in three channels – retail store, the Internet, and via catalog. The model is estimated using Bayesian methods that allow for cross-customer heterogeneity. This allows us to have distinct parameters estimates for a trial and an after trial stages and to estimate the quickness of this transit at the individual level. The results show for example that the customer decision process indeed does evolve over time. Customers differ in the duration of the trial period and marketing has a different impact on channel choice in the trial and post-trial stages. Furthermore, we show that some people switch channel decision processes while others don’t and we found that several factors have an impact on the probability to switch decision process. Insights from this study can help managers tailor their marketing communication strategy as customers gain channel choice experience. Managers may also have insights on the timing of the direct marketing communications. They can predict the duration of the trial phase at individual level detecting the customers with a quick, long or even absent trial phase. They can even predict if the customer will change or not his decision process over time, and they can influence the switching process using specific marketing tools

Miocene flora and paleoclimate estimations of Schrotzburg, South Germany

We present a detailed palaeoclimate analysis of the Middle Miocene (uppermost Badenian-lowermost Sarmatian) Schrotzburg locality in S Germany, based on the fossil macro- and micro-flora, using four different methods for the estimation of palaeoclimate parameters: the coexistence approach (CA), leaf margin analysis (LMA), the Climate-Leaf Analysis Multivariate Program (CLAMP), as well as a recently developed multivariate leaf physiognomic approach based on an European calibration dataset (ELPA). Considering results of all methods used, the following palaeoclimate estimates seem to be most likely: mean annual temperature ~15-16°C (MAT), coldest month mean temperature ~7°C (CMMT), warmest month mean temperature between 25 and 26°C, and mean annual precipiation ~1,300 mm, although CMMT values may have been colder as indicated by the disappearance of the crocodile Diplocynodon and the temperature thresholds derived from modern alligators. For most palaeoclimatic parameters, estimates derived by CLAMP significantly differ from those derived by most other methods. With respect to the consistency of the results obtained by CA, LMA and ELPA, it is suggested that for the Schrotzburg locality CLAMP is probably less reliable than most other methods. A possible explanation may be attributed to the correlation between leaf physiognomy and climate as represented by the CLAMP calibration data set which is largely based on extant floras from N America and E Asia and which may be not suitable for application to the European Neogene. All physiognomic methods used here were affected by taphonomic biasses. Especially the number of taxa had a great influence on the reliability of the palaeoclimate estimates. Both multivariate leaf physiognomic approaches are less influenced by such biasses than the univariate LMA. In combination with previously published results from the European and Asian Neogene, our data suggest that during the Neogene in Eurasia CLAMP may produce temperature estimates, which are systematically too cold as compared to other evidence. This pattern, however, has to be further investigated using additional palaeofloras.