Bayesian longitudinal data analysis with mixed models and thick-tailed distributions using MCMC

Linear mixed effects models are frequently used to analyse longitudinal data, due to their flexibility in modelling the covariance structure between and within observations. Further, it is easy to deal with unbalanced data, either with respect to the number of observations per subject or per time period, and with varying time intervals between observations. In most applications of mixed models to biological sciences, a normal distribution is assumed both for the random effects and for the residuals. This, however, makes inferences vulnerable to the presence of outliers. Here, linear mixed models employing thick-tailed distributions for robust inferences in longitudinal data analysis are described. Specific distributions discussed include the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted, and the Gibbs sampler and the Metropolis-Hastings algorithms are used to carry out the posterior analyses. An example with data on orthodontic distance growth in children is discussed to illustrate the methodology. Analyses based on either the Student-t distribution or on the usual Gaussian assumption are contrasted. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process for modelling distributions of the random effects and of residuals in linear mixed models, and the MCMC implementation allows the computations to be performed in a flexible manner.

Nitric oxide production in blowfly hemolymph after yeast inoculation

Although insects lack the adaptive immune response of the mammalians, they manifest effective innate immune responses that include both cellular and humoral components. Cellular responses are mediated by hemocytes and Immoral responses include the activation of proteolytic cascades that initiate many events, including NO production. In this work, we determined NO production in Chrysomya megaccphala hemolymph and hemocytes after yeast inoculation. Assays were performed with non-infected controls (NIL), saline-injected larvae (SIL) or larvae injected with Saccharomyces cerevisiae (YIL). The hemolymph of injected groups was collected 0.5, 1, 2, 4, 12, 24 or 48 h post-injection. NO levels in SIL were comparable to those measured in NIL until 12 h, which might be considered the basal production, increasing at 24 and 48 h post-injection, probably in response to the increased larval fragility after cuticle rupture. YIL exhibited significantly higher levels of NO than were found in other groups, peaking at 24 h. L-NAME and EDTA caused a significant reduction of NO production in YIL at this time, suggesting the activity of a Ca2+ -dependent NOS. Plasmatocytes and granular cells phagocytosed the yeasts. Plasmatocytes initiated the nodule formation and granular cells were the only hemocyte type to produce NO. These results permit us to conclude that yeasts induced augmented NO production in C. megacephala hemolymph and granular cells are the hemocyte type involved with the generation of this molecule. (c) 2005 Elsevier B.V. All rights reserved.

Studies on the respiration rate of free and immobilized cells of Cephalosporium acremonium in cephalosporin C production

Bioprocesses using filamentous fungi immobilized in inert supports present many advantages when compared to conventional free cell processes. However, assessment of the real advantages of the unconventional process demands a rigorous study of the limitations to diffusional mass transfer of the reagents, especially concerning oxygen. In this work, a comparative study was carried out on the cephalosporin C production process in defined medium containing glucose and sucrose as main carbon and energy sources, by free and immobilized cells of Cephalosporium acremonium ATCC 48272 in calcium alginate gel beads containing alumina. The effective diffusivity of oxygen through the gel beads and the effectiveness factors related to the respiration rate of the microorganism were determined experimentally. By applying Monod kinetics, the respiration kinetics parameters were experimentally determined in independent experiments in a complete production medium. The effectiveness factor experimental values presented good agreement with the theoretical values of the approximated zero-order effectiveness factor, considering the dead core model. Furthermore, experimental results obtained with immobilized cells in a 1.7-L tower bioreactor were compared with those obtained in 5-L conventional fermenter with free cells. It could be concluded that it is possible to attain rather high production rates working with relatively large diameter gel beads (ca. 2.5 mm) and sucrose consumption-based productivity was remarkably higher with immobilized cells, i.e., 0.33 gCPC/kg sucrose/h against 0.24 gCPC/kg sucrose/h in the aerated stirred tank bioreactor process. (C) 1999 John Wiley & Sons, Inc.

Controlling Excitation Inversion of a Cooper Pair Box Interacting with a Nanomechanical Resonator

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

Interactions between TLR2, TLR4, and mannose receptors with gp43 from Paracoccidioides brasiliensis induce cytokine production by human monocytes

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

Production and evaluation of lettuce seedlings and plants in hydroponic floating system

The investigation was conducted under greenhouse conditions at UNESP, Brazil, to evaluate the effects of the commercial substrates Plantmax and Sunshine and its mixtures with charred rice hulls, and the effects of floating and conventional seedling production systems in lettuce cultivation. The experiment was a randomized block experimental design, in 4 x 4 factorial treatment arrangement with two replications. Seedlings produced in the floating system were more precocious in relation to the conventional system. From transplant to harvest, the seedling production systems had similar mean time of production.

Estimates of covariance functions for growth of Nelore cattle applying a parametric correlation structure to model within-animal correlations

A total of 20,065 weights recorded on 3016 Nelore animals were used to estimate covariance functions for growth from birth to 630 days of age, assuming a parametric correlation structure to model within-animal correlations. The model of analysis included fixed effects of contemporary groups and age of dam as quadratic covariable. Mean trends were taken into account by a cubic regression on orthogonal polynomials of animal age. Genetic effects of the animal and its dam and maternal permanent environmental effects were modelled by random regressions on Legendre polynomials of age at recording. Changes in direct permanent environmental effect variances were modelled by a polynomial variance function, together with a parametric correlation function to account for correlations between ages. Stationary and nonstationary models were used to model within-animal correlations between different ages. Residual variances were considered homogeneous or heterogeneous, with changes modelled by a step or polynomial function of age at recording. Based on Bayesian information criterion, a model with a cubic variance function combined with a nonstationary correlation function for permanent environmental effects, with 49 parameters to be estimated, fitted best. Modelling within-animal correlations through a parametric correlation structure can describe the variation pattern adequately. Moreover, the number of parameters to be estimated can be decreased substantially compared to a model fitting random regression on Legendre polynomial of age. © 2004 Elsevier B.V. All rights reserved.

Double Higgs production at the LHC as a robust test of little Higgs models

We analyze double Higgs boson production at the LHC in the context of Little Higgs models. In double Higgs production, the diagrams involved are directly related to those that cause the cancellation of the quadratic divergence of the Higgs self-energy, so this mode provides a robust prediction for this class of models. We find that in extensions of this model with the inclusion of a so-called T-parity, there is a significant enhancement in the cross sections as compared to the Standard Model. © 2006 American Institute of Physics.

A study of the effect of overshooting deep convection on the water content of the TTL and lower stratosphere from Cloud Resolving Model simulations

Simulations of overshooting, tropical deep convection using a Cloud Resolving Model with bulk microphysics are presented in order to examine the effect on the water content of the TTL (Tropical Tropopause Layer) and lower stratosphere. This case study is a subproject of the HIBISCUS (Impact of tropical convection on the upper troposphere and lower stratosphere at global scale) campaign, which took place in Bauru, Brazil (22° S, 49° W), from the end of January to early March 2004. Comparisons between 2-D and 3-D simulations suggest that the use of 3-D dynamics is vital in order to capture the mixing between the overshoot and the stratospheric air, which caused evaporation of ice and resulted in an overall moistening of the lower stratosphere. In contrast, a dehydrating effect was predicted by the 2-D simulation due to the extra time, allowed by the lack of mixing, for the ice transported to the region to precipitate out of the overshoot air. Three different strengths of convection are simulated in 3-D by applying successively lower heating rates (used to initiate the convection) in the boundary layer. Moistening is produced in all cases, indicating that convective vigour is not a factor in whether moistening or dehydration is produced by clouds that penetrate the tropopause, since the weakest case only just did so. An estimate of the moistening effect of these clouds on an air parcel traversing a convective region is made based on the domain mean simulated moistening and the frequency of convective events observed by the IPMet (Instituto de Pesquisas Meteorológicas, Universidade Estadual Paulista) radar (S-band type at 2.8 Ghz) to have the same 10 dBZ echo top height as those simulated. These suggest a fairly significant mean moistening of 0.26, 0.13 and 0.05 ppmv in the strongest, medium and weakest cases, respectively, for heights between 16 and 17 km. Since the cold point and WMO (World Meteorological Organization) tropopause in this region lies at ∼ 15.9 km, this is likely to represent direct stratospheric moistening. Much more moistening is predicted for the 15-16 km height range with increases of 0.85-2.8 ppmv predicted. However, it would be required that this air is lofted through the tropopause via the Brewer Dobson circulation in order for it to have a stratospheric effect. Whether this is likely is uncertain and, in addition, the dehydration of air as it passes through the cold trap and the number of times that trajectories sample convective regions needs to be taken into account to gauge the overall stratospheric effect. Nevertheless, the results suggest a potentially significant role for convection in determining the stratospheric water content. Sensitivity tests exploring the impact of increased aerosol numbers in the boundary layer suggest that a corresponding rise in cloud droplet numbers at cloud base would increase the number concentrations of the ice crystals transported to the TTL, which had the effect of reducing the fall speeds of the ice and causing a ∼13% rise in the mean vapour increase in both the 15-16 and 16-17 km height ranges, respectively, when compared to the control case. Increases in the total water were much larger, being 34% and 132% higher for the same height ranges, but it is unclear whether the extra ice will be able to evaporate before precipitating from the region. These results suggest a possible impact of natural and anthropogenic aerosols on how convective clouds affect stratospheric moisture levels.

A nonlinear model for the long-term hydro-thermal generation scheduling problem over multiple areas with transmission constraints

This paper presents a nonlinear model with individual representation of plants for the centralized long-term hydrothermal scheduling problem over multiple areas. In addition to common aspects of long-term scheduling, this model takes transmission constraints into account. The ability to optimize hydropower exchange among multiple areas is important because it enables further minimization of complementary thermal generation costs. Also, by considering transmission constraints for long-term scheduling, a more precise coupling with shorter horizon schedules can be expected. This is an important characteristic from both operational and economic viewpoints. The proposed model is solved by a sequential quadratic programming approach in the form of a prototype system for different case studies. An analysis of the benefits provided by the model is also presented. ©2009 IEEE.

Poverty impacts of trade integration with the European Union: lessons for Ecuador

Includes bibliography

Dynamics of a continuous flash fermentation for butanol production

The ABE (acetone, butanol, ethanol) fermentation is characterized by its low productivity. In this paper, this issue is overcome with an innovative industrial process that employs the flash fermentation technology. The process consists of three interconnected units, as follows: fermentor, cell retention system (tangential microfiltration) and vacuum flash vessel (responsible for the continuous recovery of butanol from the broth). The dynamic behaviour of the process is described by a nonlinear mathematical model with kinetic parameters determined experimentally. From simulations of the mathematical model the dynamic characteristics of the process were investigated. Analyzes of the open-loop dynamic behavior of the process, after step perturbations in the manipulated variables, determined the best control structures for the process. Copyright © 2010, AIDIC Servizi S.r.l.

Buffalo milk production in brazil and colombia: Genotype by environment interaction

The objective of this study was to determine whether there is a genotype by environment interaction (GxE) for dairy buffaloes in Brazil and Colombia. The (co)variance components were estimated by using a bi-trait repeatability animal model with the REML method. Each trait consisted in the milk yield obtained in both countries. Contemporary group (herd, year and season of parity) and age at parity (linear and quadratic covariate) fixed effects, along with the additive genetic, permanent environment, and the residual random effects were included in the model. Genetic, permanent environmental and residual variance and heritabilities were different for both countries. The genetic correlations for milk yield between Brazil and Colombia were low (between 0.10 and 0.13), indicating a GxE interaction between both countries. Knowing that this interaction influences the genetic progress of buffalo populations in Brazil and Colombia, we recommend choosing sires tested in the country they will be used, along with conducting joint genetic evaluations that consider GxE interaction effects.

Observation of a diffractive contribution to dijet production in proton-proton collisions at √s=7 TeV

The cross section for dijet production in proton-proton collisions at √s=7 TeV is presented as a function of a variable that approximates the fractional momentum loss of the scattered proton in single-diffractive events. The analysis is based on an integrated luminosity of 2.7 nb-1 collected with the CMS detector at the LHC at low instantaneous luminosities, and uses events with jet transverse momentum of at least 20 GeV. The dijet cross section results are compared to the predictions of diffractive and nondiffractive models. The low- data show a significant contribution from diffractive dijet production, observed for the first time at the LHC. The associated rapidity gap survival probability is estimated. © 2013 CERN. Published by the American Physical Society.

Search for supersymmetry in events with photons and low missing transverse energy in pp collisions at √s=7TeV

Many models of new physics, including versions of supersymmetry (SUSY), predict production of events with low missing transverse energy, electroweak gauge bosons, and many energetic final-state particles. The stealth SUSY model yields this signature while conserving R-parity by means of a new hidden sector in which SUSY is approximately conserved. The results of a general search for new physics, with no requirement on missing transverse energy, in events with two photons and four or more hadronic jets are reported. The study is based on a sample of proton-proton collisions at s=7TeV corresponding to 4.96fb-1 of integrated luminosity collected with the CMS detector in 2011. Based on good agreement between the data and the standard model expectation, the data are used to determine model-independent cross-section limits and a limit on the squark mass in the framework of stealth SUSY. With this first study of its kind, squark masses less than 1430 GeV are excluded at the 95% confidence level. © 2012 CERN.