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.

On the Generalized Mittag-Leffler Function and its Application in a Fractional Telegraph Equation

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

On anomalous diffusion and the fractional generalized Langevin equation for a harmonic oscillator

The fractional generalized Langevin equation (FGLE) is proposed to discuss the anomalous diffusive behavior of a harmonic oscillator driven by a two-parameter Mittag-Leffler noise. The solution of this FGLE is discussed by means of the Laplace transform methodology and the kernels are presented in terms of the three-parameter Mittag-Leffler functions. Recent results associated with a generalized Langevin equation are recovered.

Modeling energy utilization and growth parameter description for broiler chickens

Two experiments were conducted to develop and evaluate a model to estimate ME requirements and determine Gompertz growth parameters for broilers. The first experiment was conducted to determine maintenance energy requirements and the efficiencies of energy utilization for fat and protein deposition. Maintenance ME (ME m) requirements were estimated to be 157.8, 112.1, and 127.2 kcal of ME/kg 0.75 per day for broilers at 13, 23, and 32°C, respectively. Environmental temperature (T) had a quadratic effect on maintenance requirements (ME m = 307.87 - 15.63T + 0.3105T 2; r 2= 0.93). Energy requirements for fat and protein deposition were estimated to be 13.52 and 12.59 kcal of ME/g, respectively. Based on these coefficients, a model was developed to calculate daily ME requirements: ME = BW 0.75 (307.87 - 15.63T + 0.3105 T 2) + 13.52 G f + 12.59 G p. This model considers live BW, the effects of environmental temperature, and fractional fat (G f) and protein (G p) deposition. The second experiment was carried out to estimate the growth parameters of Ross broilers and to collect data to evaluate the ME requirement model proposed. Live BW, empty feather-free carcass, weight of the feathers, and carcass chemical compositions were analyzed until 16 wk of age. Parameters of Gompertz curves for each component were estimated. Males had higher growth potential and higher capacity to deposit nutrients than females, except for fat deposition. Data of BW and body composition collected in this experiment were fitted into the energy model proposed herein and the equations described by Emmans (1989) and Chwalibog (1991). The daily ME requirements estimated by the model determined in this study were closer to the ME intake observed in this trial compared with other models. ©2005 Poultry Science Association, Inc.

The fractional and nonlinear magneto-flexible rod

This paper deals with a system involving a flexible rod subjected to magnetic forces that can bend it while simultaneously subjected to external excitations produces complex and nonlinear dynamic behavior, which may present different types of solutions for its different movement-related responses. This fact motivated us to analyze such a mechanical system based on modeling and numerical simulation involving both, integer order calculus (IOC) and fractional order calculus (FOC) approaches. The time responses, pseudo phase portraits and Fourier spectra have been presented. The results obtained can be used as a source for conduct experiments in order to obtain more realistic and more accurate results about fractional-order models when compared to the integer-order models. © Published under licence by IOP Publishing Ltd.

The fractional-nonlinear robotic manipulator: Modeling and dynamic simulations

In this paper, we applied the Riemann-Liouville approach and the fractional Euler-Lagrange equations in order to obtain the fractional-order nonlinear dynamics equations of a two link robotic manipulator. The aformentioned equations have been simulated for several cases involving: integer and non-integer order analysis, with and without external forcing acting and some different initial conditions. The fractional nonlinear governing equations of motion are coupled and the time evolution of the angular positions and the phase diagrams have been plotted to visualize the effect of fractional order approach. The new contribution of this work arises from the fact that the dynamics equations of a two link robotic manipulator have been modeled with the fractional Euler-Lagrange dynamics approach. The results reveal that the fractional-nonlinear robotic manipulator can exhibit different and curious behavior from those obtained with the standard dynamical system and can be useful for a better understanding and control of such nonlinear systems. © 2012 American Institute of Physics.

Estruturas shrimp e propriedades dinâmicas no modelo dissipativo do acelerador de Fermi

Pós-graduação em Física - IGCE

TYPE-ZERO SADDLE-NODE BIFURCATIONS AND STABILITY REGION ESTIMATION OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEMS

A complete characterization of the stability boundary of a class of nonlinear dynamical systems that admit energy functions is developed in this paper. This characterization generalizes the existing results by allowing the type-zero saddle-node nonhyperbolic equilibrium points on the stability boundary. Conceptual algorithms to obtain optimal estimates of the stability region (basin of attraction) in the form of level sets of a given family of energy functions are derived. The behavior of the stability region and the corresponding estimates are investigated for parameter variation in the neighborhood of a type-zero saddle-node bifurcation value.

Coarse-grained peptide models: conformational sampling, peptide association and dynamical properties for peptides

In dieser Arbeit wird ein vergröbertes (engl. coarse-grained, CG) Simulationsmodell für Peptide in wässriger Lösung entwickelt. In einem CG Verfahren reduziert man die Anzahl der Freiheitsgrade des Systems, so dass manrngrössere Systeme auf längeren Zeitskalen untersuchen kann. Die Wechselwirkungspotentiale des CG Modells sind so aufgebaut, dass die Peptid Konformationen eines höher aufgelösten (atomistischen) Modells reproduziert werden.rnIn dieser Arbeit wird der Einfluss unterschiedlicher bindender Wechsel-rnwirkungspotentiale in der CG Simulation untersucht, insbesondere daraufhin,rnin wie weit das Konformationsgleichgewicht der atomistischen Simulation reproduziert werden kann. Im CG Verfahren verliert man per Konstruktionrnmikroskopische strukturelle Details des Peptids, zum Beispiel, Korrelationen zwischen Freiheitsgraden entlang der Peptidkette. In der Dissertationrnwird gezeigt, dass diese “verlorenen” Eigenschaften in einem Rückabbildungsverfahren wiederhergestellt werden können, in dem die atomistischen Freiheitsgrade wieder in die CG-Strukturen eingefügt werden. Dies gelingt, solange die Konformationen des CG Modells grundsätzlich gut mit der atomistischen Ebene übereinstimmen. Die erwähnten Korrelationen spielen einerngrosse Rolle bei der Bildung von Sekundärstrukturen und sind somit vonrnentscheidender Bedeutung für ein realistisches Ensemble von Peptidkonformationen. Es wird gezeigt, dass für eine gute Übereinstimmung zwischen CG und atomistischen Kettenkonformationen spezielle bindende Wechselwirkungen wie zum Beispiel 1-5 Bindungs- und 1,3,5-Winkelpotentiale erforderlich sind. Die intramolekularen Parameter (d.h. Bindungen, Winkel, Torsionen), die für kurze Oligopeptide parametrisiert wurden, sind übertragbarrnauf längere Peptidsequenzen. Allerdings können diese gebundenen Wechselwirkungen nur in Kombination mit solchen nichtbindenden Wechselwirkungspotentialen kombiniert werden, die bei der Parametrisierung verwendet werden, sind also zum Beispiel nicht ohne weiteres mit einem andere Wasser-Modell kombinierbar. Da die Energielandschaft in CG-Simulationen glatter ist als im atomistischen Modell, gibt es eine Beschleunigung in der Dynamik. Diese Beschleunigung ist unterschiedlich für verschiedene dynamische Prozesse, zum Beispiel für verschiedene Arten von Bewegungen (Rotation und Translation). Dies ist ein wichtiger Aspekt bei der Untersuchung der Kinetik von Strukturbildungsprozessen, zum Beispiel Peptid Aggregation.rn

Non-exponential return time distributions for vorticity extremes explained by fractional Poisson processes

Serial correlation of extreme midlatitude cyclones observed at the storm track exits is explained by deviations from a Poisson process. To model these deviations, we apply fractional Poisson processes (FPPs) to extreme midlatitude cyclones, which are defined by the 850 hPa relative vorticity of the ERA interim reanalysis during boreal winter (DJF) and summer (JJA) seasons. Extremes are defined by a 99% quantile threshold in the grid-point time series. In general, FPPs are based on long-term memory and lead to non-exponential return time distributions. The return times are described by a Weibull distribution to approximate the Mittag–Leffler function in the FPPs. The Weibull shape parameter yields a dispersion parameter that agrees with results found for midlatitude cyclones. The memory of the FPP, which is determined by detrended fluctuation analysis, provides an independent estimate for the shape parameter. Thus, the analysis exhibits a concise framework of the deviation from Poisson statistics (by a dispersion parameter), non-exponential return times and memory (correlation) on the basis of a single parameter. The results have potential implications for the predictability of extreme cyclones.

On the Equivalence of the Riemann-Liouville and the Caputo Fractional Order Derivatives in Modeling of Linear Viscoelastic Materials

Mathematics Subject Classification: 26A33

Nonlinear Time-Fractional Differential Equations in Combustion Science

MSC 2010: 34A08 (main), 34G20, 80A25

Matrix-Variate Statistical Distributions and Fractional Calculus

MSC 2010: 15A15, 15A52, 33C60, 33E12, 44A20, 62E15 Dedicated to Professor R. Gorenflo on the occasion of his 80th birthday

Table 1 Core site, fractional abundance of individual isoprenoid GDGTs and TEX86 values

TEXL86 and TEXH86 are organic palaeothermometers based on the lipids of Group 1 Crenarchaeota, recently proposed as a modified version of the original TEX86 index, but with significantly improved geographical coverage. Since few data from the global core top calibration are from the Pacific, this study was carried out to assess whether the global core top calibration is regionally biased or not. The result of principal components analysis of the fractional abundance of GDGTs, an analysis of variance (ANOVA) and the comparison of the residuals of TEXH 86 derived sea surface temperature (SST) estimates of the Pacific subset with that of the global data set suggest that the Pacific subset has a similar TEXH 86-SST relationship with the global data set. However, the regression line through the Pacific data and an ANOVA on the residuals of TEXL 86 derived SST estimates suggest otherwise. The contradictory findings are likely to stem from the large scatter in the Pacific TEXL 86 values in the mid temperature range. While regionality does not seem to exert a strong bias on TEXL 86 and TEXH 86 calibration, it appears that there is a strong need to resolve the large scatter in the global data set, especially in the mid and high latitudes, in order to improve the calibration for a better SST estimation.

Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators: the Riemann-Liouville case

In this paper we study eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator $\Delta_+^{(\alpha,\beta,\gamma)}:= D_{x_0^+}^{1+\alpha} +D_{y_0^+}^{1+\beta} +D_{z_0^+}^{1+\gamma},$ where $(\alpha, \beta, \gamma) \in \,]0,1]^3$, and the fractional derivatives $D_{x_0^+}^{1+\alpha}$, $D_{y_0^+}^{1+\beta}$, $D_{z_0^+}^{1+\gamma}$ are in the Riemann-Liouville sense. Applying operational techniques via two-dimensional Laplace transform we describe a complete family of eigenfunctions and fundamental solutions of the operator $\Delta_+^{(\alpha,\beta,\gamma)}$ in classes of functions admitting a summable fractional derivative. Making use of the Mittag-Leffler function, a symbolic operational form of the solutions is presented. From the obtained family of fundamental solutions we deduce a family of fundamental solutions of the fractional Dirac operator, which factorizes the fractional Laplace operator. We apply also the method of separation of variables to obtain eigenfunctions and fundamental solutions.