Robust linear mixed models with normal/independent distributions and Bayesian MCMC implementation

Linear mixed effects models have been widely used in analysis of data where responses are clustered around some random effects, so it is not reasonable to assume independence between observations in the same cluster. In most biological applications, it is assumed that the distributions of the random effects and of the residuals are Gaussian. This makes inferences vulnerable to the presence of outliers. Here, linear mixed effects models with normal/independent residual distributions for robust inferences are described. Specific distributions examined include univariate and multivariate versions of the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted and Markov chain Monte Carlo is used to carry out the posterior analysis. The procedures are illustrated using birth weight data on rats in a texicological experiment. Results from the Gaussian and robust models are contrasted, and it is shown how the implementation can be used for outlier detection. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process in linear mixed models, and they are easily implemented using data augmentation and MCMC techniques.

Stable two-dimensional dispersion-managed soliton

The existence of a dispersion-managed soliton in two-dimensional nonlinear Schrodinger equation with periodically varying dispersion has been explored. The averaged equations for the soliton width and chirp are obtained which successfully describe the long time evolution of the soliton. The slow dynamics of the soliton around the fixed points for the width and chirp are investigated and the corresponding frequencies are calculated. Analytical predictions are confirmed by direct partial differential equation (PDE) and ordinary differential equation (ODE) simulations. Application to a Bose-Einstein condensate in optical lattice is discussed. The existence of a dispersion-managed matter-wave soliton in such system is shown.

Cytokines and dietary energy restriction in stable chronic obstructive pulmonary disease patients

Tumour necrosis factor (TNF)-α has been found to be increased in malnourished chronic obstructive pulmonary disease (COPD) patients; however, the main cause of this phenomenon remains undetermined. In normal subjects, TNF-α production may be induced by dietary energy deprivation. The aim of this study was to investigate if stable COPD patients present alterations of inflammatory mediators after 48 h of dietary energy restriction. Fourteen COPD patients were admitted to the hospital while receiving an experimental diet with an energy content of approximately one-third of their energy needs. Clinical evaluation, nutritional assessment and serum levels of interleukin (IL)-6, TNF-α and C-reactive protein, and secretion of TNF-α by peripheral blood monocytes were assessed on admission and after the experimental diet. For reference values of the laboratory parameters, blood was collected from 10 healthy, elderly subjects. COPD patients showed significantly higher serum concentrations of IL-6 than control subjects, however, the experimental diet was not associated with statistically significant changes in the inflammatory mediators. The findings of this study, although preliminary because of the limited degree and duration of the energy restriction, suggest that the elevated levels of tumour necrosis factor-α, previously described in undernourished or weight-losing chronic obstructive pulmonary disease patients, may not be linked to a decrease of dietary energy intake.

Tracing metabolic routes of feed ingredients in tissues of broiler chickens using stable isotopes

The present study aimed to quantify the proportion of 13C from energy and protein feed ingredients that follow the metabolic routing of the liver and muscle in broiler chickens. A stable isotope of carbon technique was used that is based on the isotopic discrimination that occurs in the plants during the photosynthesis process. One-day-old male chicks were subjected to treatments based on free choice of energy and protein sources. Rice bran (R) and soybean meal (S), C3 plants, have higher isotopic ratios than corn (C), a C4 plant, and corn gluten meal (G). Choices were R+S, C+G, R+G, C+S, or R+C+G+S. A complete feed (CF) was a sixth treatment. Feed intake and BW were measured at 30 d of age, when liver and breast muscle were collected for isotopic analysis. Treatments affected the amount of feed intake and the choices of energy or protein sources. Complete feed had the largest intake, differing from the other treatments that had free-choice feeding. Final BW was a direct reflection of consumption by these birds in all treatments. The isotopic results indicated that the 13C/12C ratio was generally higher in breast muscle than in liver, probably because of higher protein content. Moreover, in the liver, the proportion of 13C retained from the energy ingredient was greater than the proportion from the protein ingredient. That is in contrast to muscle, where the proportion of 13C retained from the protein ingredient was greater than from the energy ingredient that was self-selected.

Template-based fluoropolymer ferroelectrets with multiple layers of tubular channels

Significant efforts are devoted to developing new ferroelectrets with well-controlled void distributions or uniform voids and with good long-term and thermal stability of the piezoelectricity. Here, we describe the concept, the fabrication, and the most relevant properties of fluoropolymer ferroelectret systems with three separate films of fluoroethylenepropylene (FEP), alternating with two polytetrafluoroethylene (PTFE) templates. The FEP films are selectively fused by means of a lamination process. Two practically identical PTFE templates are used, which have parallel rectangular openings (1.5×30 mm 2) separated by PTFE ridges of 1.5 mm width. After removing the PTFE templates, a three-layer FEP-film sandwich with tubular channels is obtained. We demonstrate that such FEP-film systems exhibit significant and stable piezoelectricity after charging under a high DC voltage. The resulting piezoelectric effect may be further improved by carefully assembling and arranging the PTFE templates during preparation. ©2010 IEEE.

A maximum principle for infinite time asymptotically stable impulsive dynamic control systems

We consider an infinite horizon optimal impulsive control problems for which a given cost function is minimized by choosing control strategies driving the state to a point in a given closed set C ∞. We present necessary conditions of optimality in the form of a maximum principle for which the boundary condition of the adjoint variable is such that non-degeneracy due to the fact that the time horizon is infinite is ensured. These conditions are given for conventional systems in a first instance and then for impulsive control problems. They are proved by considering a family of approximating auxiliary interval conventional (without impulses) optimal control problems defined on an increasing sequence of finite time intervals. As far as we know, results of this kind have not been derived previously. © 2010 IFAC.

Andean exchange-rate regimes, 1994-2003: a brief for stable but flexible" regimes"

Includes bibliography

Sensitivity, specificity and predictive values of linear and nonlinear indices of heart rate variability in stable angina patients

Background: Decreased heart rate variability (HRV) is related to higher morbidity and mortality. In this study we evaluated the linear and nonlinear indices of the HRV in stable angina patients submitted to coronary angiography. Methods. We studied 77 unselected patients for elective coronary angiography, which were divided into two groups: coronary artery disease (CAD) and non-CAD groups. For analysis of HRV indices, HRV was recorded beat by beat with the volunteers in the supine position for 40 minutes. We analyzed the linear indices in the time (SDNN [standard deviation of normal to normal], NN50 [total number of adjacent RR intervals with a difference of duration greater than 50ms] and RMSSD [root-mean square of differences]) and frequency domains ultra-low frequency (ULF) ≤ 0,003 Hz, very low frequency (VLF) 0,003 - 0,04 Hz, low frequency (LF) (0.04-0.15 Hz), and high frequency (HF) (0.15-0.40 Hz) as well as the ratio between LF and HF components (LF/HF). In relation to the nonlinear indices we evaluated SD1, SD2, SD1/SD2, approximate entropy (-ApEn), α1, α2, Lyapunov Exponent, Hurst Exponent, autocorrelation and dimension correlation. The definition of the cutoff point of the variables for predictive tests was obtained by the Receiver Operating Characteristic curve (ROC). The area under the ROC curve was calculated by the extended trapezoidal rule, assuming as relevant areas under the curve ≥ 0.650. Results: Coronary arterial disease patients presented reduced values of SDNN, RMSSD, NN50, HF, SD1, SD2 and -ApEn. HF ≤ 66 ms§ssup§2§esup§, RMSSD ≤ 23.9 ms, ApEn ≤-0.296 and NN50 ≤ 16 presented the best discriminatory power for the presence of significant coronary obstruction. Conclusion: We suggest the use of Heart Rate Variability Analysis in linear and nonlinear domains, for prognostic purposes in patients with stable angina pectoris, in view of their overall impairment. © 2012 Pivatelli et al.; licensee BioMed Central Ltd.

Improving methods in gap ecology: Revisiting size and shape distributions using a model selection approach

Questions: We assess gap size and shape distributions, two important descriptors of the forest disturbance regime, by asking: which statistical model best describes gap size distribution; can simple geometric forms adequately describe gap shape; does gap size or shape vary with forest type, gap age or the method used for gap delimitation; and how similar are the studied forests and other tropical and temperate forests? Location: Southeastern Atlantic Forest, Brazil. Methods: Analysing over 150 gaps in two distinct forest types (seasonal and rain forests), a model selection framework was used to select appropriate probability distributions and functions to describe gap size and gap shape. The first was described using univariate probability distributions, whereas the latter was assessed based on the gap area-perimeter relationship. Comparisons of gap size and shape between sites, as well as size and age classes were then made based on the likelihood of models having different assumptions for the values of their parameters. Results: The log-normal distribution was the best descriptor of gap size distribution, independently of the forest type or gap delimitation method. Because gaps became more irregular as they increased in size, all geometric forms (triangle, rectangle and ellipse) were poor descriptors of gap shape. Only when small and large gaps (> 100 or 400m2 depending on the delimitation method) were treated separately did the rectangle and isosceles triangle become accurate predictors of gap shape. Ellipsoidal shapes were poor descriptors. At both sites, gaps were at least 50% longer than they were wide, a finding with important implications for gap microclimate (e.g. light entrance regime) and, consequently, for gap regeneration. Conclusions: In addition to more appropriate descriptions of gap size and shape, the model selection framework used here efficiently provided a means by which to compare the patterns of two different types of forest. With this framework we were able to recommend the log-normal parameters μ and σ for future comparisons of gap size distribution, and to propose possible mechanisms related to random rates of gap expansion and closure. We also showed that gap shape varied highly and that no single geometric form was able to predict the shape of all gaps, the ellipse in particular should no longer be used as a standard gap shape. © 2012 International Association for Vegetation Science.

Stable singularities of co-rank one quasi homogeneous map germs from (ℂn+1, 0) to (ℂn, 0), n = 2, 3

In this article, we investigate the geometry of quasi homogeneous corank one finitely determined map germs from (ℂn+1, 0) to (ℂn, 0) with n = 2, 3. We give a complete description, in terms of the weights and degrees, of the invariants that are associated to all stable singularities which appear in the discriminant of such map germs. The first class of invariants which we study are the isolated singularities, called 0-stable singularities because they are the 0-dimensional singularities. First, we give a formula to compute the number of An points which appear in any stable deformation of a quasi homogeneous co-rank one map germ from (ℂn+1, 0) to (ℂn, 0) with n = 2, 3. To get such a formula, we apply the Hilbert's syzygy theorem to determine the graded free resolution given by the syzygy modules of the associated iterated Jacobian ideal. Then we show how to obtain the other 0-stable singularities, these isolated singularities are formed by multiple points and here we use the relation among them and the Fitting ideals of the discriminant. For n = 2, there exists only the germ of double points set and for n = 3 there are the triple points, named points A1,1,1 and the normal crossing between a germ of a cuspidal edge and a germ of a plane, named A2,1. For n = 3, there appear also the one-dimensional singularities, which are of two types: germs of cuspidal edges or germs of double points curves. For these singularities, we show how to compute the polar multiplicities and also the local Euler obstruction at the origin in terms of the weights and degrees. © 2013 Pushpa Publishing House.

Stable regions around Pluto

In a previous work, GiuliattiWinter et al. found several stable regions for test particles in orbit around Pluto associated with families of periodic orbits obtained in the circular, restricted three-body problem. They have shown that a possible eccentricity of the Pluto-Charon binary slightly reduces but does not destroy any of these stable regions. In thiswork, we extended their results by analysing the cases with the orbital inclination (I) equal to zero and considering the argument of pericentre (w) equal to 90°, 180° and 270°. We explore the influence of the orbital inclination of the particles in these stable regions. In this case, the initial inclination varies from 10° to 170° in steps of 10°. We also present a sample of results for the longitude of the ascending node Ω = 90°, considering the cases I = 20°, 50°, 130° and 180°. Our results show that stable regions are present in all of the inclined cases, except when the initial inclination of the particles is equal to 110°. A sample of 3D trajectories of quasi-periodic orbits were found related to the periodic orbits obtained in the planar case by Giuliatti Winter et al. © 2013 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Single-particle momentum distributions of Efimov states in mixed-species systems

We solve the three-body bound-state problem in three dimensions for mass imbalanced systems of two identical bosons and a third particle in the universal limit where the interactions are assumed to be of zero range. The system displays the Efimov effect and we use the momentum-space wave equation to derive formulas for the scaling factor of the Efimov spectrum for any mass ratio assuming either that two or three of the two-body subsystems have a bound state at zero energy. We consider the single-particle momentum distribution analytically and numerically and analyze the tail of the momentum distribution to obtain the three-body contact parameter. Our findings demonstrate that the functional form of the three-body contact term depends on the mass ratio, and we obtain an analytic expression for this behavior. To exemplify our results, we consider mixtures of lithium with either two caesium or rubidium atoms which are systems of current experimental interest. © 2013 American Physical Society.

Division of labor in stable social hierarchy of the independent-founding wasp Mischocyttarus (Monocyttarus) cassununga, von Ihering (Hymenoptera, Vespidae)

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

Instantaneously measured traits may detect non-plastic ecophysiological performances in response to drought, explaining distributions of Styrax species in the Cerrado

We analyzed the differences between irrigated and non-irrigated plants of three congeneric Styrax species that present distinct distribution patterns in the physiognomies of the Cerrado vegetation in Brazil. Styrax ferrugineus showed a stomatal conductance (gs) unresponsive to soil water deficit in potted plants. This may explain the high gas exchange and photochemical efficiency found in this species, which is well adapted to the Cerrado sensu stricto (s. str.), a savanna-type vegetation. S. camporum, which is widely distributed in the Cerrado sensu lato (s. l.) areas, was the only species that exhibited increased intrinsic water use efficiency on the days of maximum water deficit. This result distinguishes S. camporum from S. pohlii, which is a forest species, since the gs of both species decreased during the days of maximum water stress. In contrast to other studies, we propose that instantaneously measured traits, such as leaf gas exchange rates and chlorophyll fluorescence, may be used to detect non-plastic performances in response to environmental stress, helping explain distinct geographical distributions of congeneric species in the Cerrado vegetation. © 2013 Springer-Verlag Berlin Heidelberg.

Diastolic function is associated with quality of life and exercise capacity in stable heart failure patients with reduced ejection fraction

Exercise capacity and quality of life (QOL) are important outcome predictors in patients with systolic heart failure (HF), independent of left ventricular (LV) ejection fraction (LVEF). LV diastolic function has been shown to be a better predictor of aerobic exercise capacity in patients with systolic dysfunction and a New York Heart Association (NYHA) classification >II. We hypothesized that the currently used index of diastolic function E/e' is associated with exercise capacity and QOL, even in optimally treated HF patients with reduced LVEF. This prospective study included 44 consecutive patients aged 55±11 years (27 men and 17 women), with LVEF,0.50 and NYHA functional class I-III, receiving optimal pharmacological treatment and in a stable clinical condition, as shown by the absence of dyspnea exacerbation for at least 3 months. All patients had conventional transthoracic echocardiography and answered the Minnesota Living with HF Questionnaire, followed by the 6-min walk test (6MWT). In a multivariable model with 6MWT as the dependent variable, age and E/e' explained 27% of the walked distance in 6MWT (P=0.002; multivariate regression analysis). No association was found between walk distance and LVEF or mitral annulus systolic velocity. Only normalized left atrium volume, a sensitive index of diastolic function, was associated with decreased QOL. Despite the small number of patients included, this study offers evidence that diastolic function is associated with physical capacity and QOL and should be considered along with ejection fraction in patients with compensated systolic HF.