Estimation of genetic variance for macro- and micro-environmental sensitivity using double hierarchical generalized linear models

Background: Genetic variation for environmental sensitivity indicates that animals are genetically different in their response to environmental factors. Environmental factors are either identifiable (e.g. temperature) and called macro-environmental or unknown and called micro-environmental. The objectives of this study were to develop a statistical method to estimate genetic parameters for macro- and micro-environmental sensitivities simultaneously, to investigate bias and precision of resulting estimates of genetic parameters and to develop and evaluate use of Akaike’s information criterion using h-likelihood to select the best fitting model. Methods: We assumed that genetic variation in macro- and micro-environmental sensitivities is expressed as genetic variance in the slope of a linear reaction norm and environmental variance, respectively. A reaction norm model to estimate genetic variance for macro-environmental sensitivity was combined with a structural model for residual variance to estimate genetic variance for micro-environmental sensitivity using a double hierarchical generalized linear model in ASReml. Akaike’s information criterion was constructed as model selection criterion using approximated h-likelihood. Populations of sires with large half-sib offspring groups were simulated to investigate bias and precision of estimated genetic parameters. Results: Designs with 100 sires, each with at least 100 offspring, are required to have standard deviations of estimated variances lower than 50% of the true value. When the number of offspring increased, standard deviations of estimates across replicates decreased substantially, especially for genetic variances of macro- and micro-environmental sensitivities. Standard deviations of estimated genetic correlations across replicates were quite large (between 0.1 and 0.4), especially when sires had few offspring. Practically, no bias was observed for estimates of any of the parameters. Using Akaike’s information criterion the true genetic model was selected as the best statistical model in at least 90% of 100 replicates when the number of offspring per sire was 100. Application of the model to lactation milk yield in dairy cattle showed that genetic variance for micro- and macro-environmental sensitivities existed. Conclusion: The algorithm and model selection criterion presented here can contribute to better understand genetic control of macro- and micro-environmental sensitivities. Designs or datasets should have at least 100 sires each with 100 offspring.

Hglm: A package for fitting hierarchical generalized linear models

We present the hglm package for fitting hierarchical generalized linear models. It can be used for linear mixed models and generalized linear mixed models with random effects for a variety of links and a variety of distributions for both the outcomes and the random effects. Fixed effects can also be fitted in the dispersion part of the model.

Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Structural damage detection by fuzzy clustering

The development of strategies for structural health monitoring (SHM) has become increasingly important because of the necessity of preventing undesirable damage. This paper describes an approach to this problem using vibration data. It involves a three-stage process: reduction of the time-series data using principle component analysis (PCA), the development of a data-based model using an auto-regressive moving average (ARMA) model using data from an undamaged structure, and the classification of whether or not the structure is damaged using a fuzzy clustering approach. The approach is applied to data from a benchmark structure from Los Alamos National Laboratory, USA. Two fuzzy clustering algorithms are compared: fuzzy c-means (FCM) and Gustafson-Kessel (GK) algorithms. It is shown that while both fuzzy clustering algorithms are effective, the GK algorithm marginally outperforms the FCM algorithm. (C) 2008 Elsevier Ltd. All rights reserved.

New Insight into the Mechanism of Action of Wasp Mastoparan Peptides: Lytic Activity and Clustering Observed with Giant Vesicles

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

Cosmic microwave background and large-scale structure constraints on a simple quintessential inflation model

We derive constraints on a simple quintessential inflation model, based on a spontaneously broken Phi(4) theory, imposed by the Wilkinson Microwave Anisotropy Probe three-year data (WMAP3) and by galaxy clustering results from the Sloan Digital Sky Survey (SDSS). We find that the scale of symmetry breaking must be larger than about 3 Planck masses in order for inflation to generate acceptable values of the scalar spectral index and of the tensor-to-scalar ratio. We also show that the resulting quintessence equation of state can evolve rapidly at recent times and hence can potentially be distinguished from a simple cosmological constant in this parameter regime.

Analysis of beef cattle longitudinal data applying a nonlinear model

The objective of this work was to evaluate the Nelore beef cattle, growth curve parameters using the Von Bertalanffy function in a nested Bayesian procedure that allowed estimation of the joint posterior distribution of growth curve parameters, their (co)variance components, and the environmental and additive genetic components affecting them. A hierarchical model was applied; each individual had a growth trajectory described by the nonlinear function, and each parameter of this function was considered to be affected by genetic and environmental effects that were described by an animal model. Random samples of the posterior distributions were drawn using Gibbs sampling and Metropolis-Hastings algorithms. The data set consisted of a total of 145,961 BW recorded from 15,386 animals. Even though the curve parameters were estimated for animals with few records, given that the information from related animals and the structure of systematic effects were considered in the curve fitting, all mature BW predicted were suitable. A large additive genetic variance for mature BW was observed. The parameter a of growth curves, which represents asymptotic adult BW, could be used as a selection criterion to control increases in adult BW when selecting for growth rate. The effect of maternal environment on growth was carried through to maturity and should be considered when evaluating adult BW. Other growth curve parameters showed small additive genetic and maternal effects. Mature BW and parameter k, related to the slope of the curve, presented a large, positive genetic correlation. The results indicated that selection for growth rate would increase adult BW without substantially changing the shape of the growth curve. Selection to change the slope of the growth curve without modifying adult BW would be inefficient because their genetic correlation is large. However, adult BW could be considered in a selection index with its corresponding economic weight to improve the overall efficiency of beef cattle production.

Soliton clusters inducing insulator-to-metal transition in polyacetylene: a Su-Schrieffer-Heeger model study

The Su-Schrieffer-Heeger (SSH) Hamiltonian has been one of most used models to study the electronic structure of polyacetylene (PA) chains. It has been reported in the literature that in the SSH framework a disordered soliton distribution can not produce a metallic regime. However, in this work (using the same SSH model and parameters) we show that this is possible. The necessary conditions for true metals (non-vanishing density of states and extended wavefunctions around the Fermi level) are obtained for soliton concentration higher than 6% through soliton segregation (clustering). These results are consistent with recent experimental data supporting disorder as an essential mechanism behind the high conductivity of conducting polymers. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Clustering of rare earth in glasses, aluminum effect: experiments and modeling

Luminescent spectra of Eu3+-doped sol-gel glasses have been analyzed during the densification process and compared according to the presence or not of aluminum as a codoping ion. A transition temperature from hydrated to dehydroxyled environments has been found different for doped and codoped samples. However, only slight modifications have been displayed from luminescence measurements beyond this transition. To support the experimental analysis, molecular dynamics simulations have been performed to model the doped and codoped glass structures. Despite no evidence of rare earth clustering reduction due to aluminum has been found, the modeled structures have shown that the luminescent ions are mainly located in aluminum-rich domains. The synthesis of both experimental and numerical analyses has lead us to interpret the aluminum effect as responsible for differences in structure of the luminescent sites rather than for an effective dispersion of the rare earth ions. (C) 2004 Elsevier B.V. All rights reserved.

Branch and bound algorithm for transmission system expansion planning using a transportation model

A method for optimal transmission network expansion planning is presented. The transmission network is modelled as a transportation network. The problem is solved using hierarchical Benders decomposition in which the problem is decomposed into master and slave subproblems. The master subproblem models the investment decisions and is solved using a branch-and-bound algorithm. The slave subproblem models the network operation and is solved using a specialised linear program. Several alternative implementations of the branch-and-bound algorithm have been rested. Special characteristics of the transmission expansion problem have been taken into consideration in these implementations. The methods have been tested on various test systems available in the literature.

NUMERICAL STUDY OF A PERTURBED EINSTEIN-DESITTER COSMOLOGICAL MODEL

Given the ever-increasing scale of structures discovered in the universe, we solve Einstein's equations numerically, under simplifying assumptions, to examine how this lack of uniformity affects the metric of Einstein-de Sitter cosmology. The results confirm the qualitative conclusion of Barrow, that a large density contrast is compatible with much smaller metric perturbations. The contribution of this peculiar gravity to the redshift might complicate studies of peculiar motions of galaxies, although it appears that the distortion is small for nearby clusters.

HIERARCHICAL VERSUS NONHIERARCHICAL PATTERNS OF GENETIC DISTANCES AMONG POPULATIONS - A SIMULATION STUDY

A simulation study was made of the effects of mixing two evolutionary forces (natural selection and random genetic drift), combined in a single data matrix of gene frequencies, on the resulting genetic distances among populations. Twenty-one, kinds of simulated gene frequencies surfaces, for 15 populations linearly distributed over geographic space, were used to construct 21 data matrices, combining different proportions of two types of surfaces (gradients and random surfaces). These matrices were analysed by Unweighted Pair-Group Method - Arithmetic Averages (UPGMA), clustering and Principal Coordinate Analysis. The results obtained show that ordination is more accurate than UPGMA in revealing the spatial patterns in the genetic distances, in comparison with results obtained using the Mantel test comparing directly genetic and geographic distances.

Hierarchical Assembly of CaMoO4 Nano-Octahedrons and Their Photoluminescence Properties

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

Development of a model to measure the efficiency of outsourced companies in the process of aircraft spare parts technical publications

This paper presents a proposal of a model to measure the efficiency of outsourced companies in the aeronautical industry applying the methods DEA and AHP. It also proposes an evaluation in the relation between the variables of the process and the value obtained for the effiiency. The criteria of Quality, Time and Cost were considered the outputs of the process, and those criteria were quantified by AHP for DEA matrix.The number of technical documents received by those outsorced companies were considered the input of the process. The other purpose is to separate the companies in groups considered able to receive an investment to improve their process. Copyright © 2008 SAE International.

Light front fermion model propagation

In this work we consider the propagation of two fermion fields interacting with each other by the exchange of intermediate scalar bosons in the light front. We obtain the corrections up to fourth order in the coupling constant using hierarchical equations in order to obtain the bound state equation (Bethe - Salpeter equation). © 2013 Chinese Physical Society and IOP Publishing Ltd.