Variational Markov Chain Monte Carlo for Bayesian smoothing of non-linear niffusions

In this paper we develop set of novel Markov chain Monte Carlo algorithms for Bayesian smoothing of partially observed non-linear diffusion processes. The sampling algorithms developed herein use a deterministic approximation to the posterior distribution over paths as the proposal distribution for a mixture of an independence and a random walk sampler. The approximating distribution is sampled by simulating an optimized time-dependent linear diffusion process derived from the recently developed variational Gaussian process approximation method. Flexible blocking strategies are introduced to further improve mixing, and thus the efficiency, of the sampling algorithms. The algorithms are tested on two diffusion processes: one with double-well potential drift and another with SINE drift. The new algorithm's accuracy and efficiency is compared with state-of-the-art hybrid Monte Carlo based path sampling. It is shown that in practical, finite sample, applications the algorithm is accurate except in the presence of large observation errors and low observation densities, which lead to a multi-modal structure in the posterior distribution over paths. More importantly, the variational approximation assisted sampling algorithm outperforms hybrid Monte Carlo in terms of computational efficiency, except when the diffusion process is densely observed with small errors in which case both algorithms are equally efficient.

Variational Markov chain Monte Carlo for Bayesian smoothing of non-linear diffusions

In this paper we develop set of novel Markov Chain Monte Carlo algorithms for Bayesian smoothing of partially observed non-linear diffusion processes. The sampling algorithms developed herein use a deterministic approximation to the posterior distribution over paths as the proposal distribution for a mixture of an independence and a random walk sampler. The approximating distribution is sampled by simulating an optimized time-dependent linear diffusion process derived from the recently developed variational Gaussian process approximation method. The novel diffusion bridge proposal derived from the variational approximation allows the use of a flexible blocking strategy that further improves mixing, and thus the efficiency, of the sampling algorithms. The algorithms are tested on two diffusion processes: one with double-well potential drift and another with SINE drift. The new algorithm's accuracy and efficiency is compared with state-of-the-art hybrid Monte Carlo based path sampling. It is shown that in practical, finite sample applications the algorithm is accurate except in the presence of large observation errors and low to a multi-modal structure in the posterior distribution over paths. More importantly, the variational approximation assisted sampling algorithm outperforms hybrid Monte Carlo in terms of computational efficiency, except when the diffusion process is densely observed with small errors in which case both algorithms are equally efficient. © 2011 Springer-Verlag.

Coactivation of the shoulder and arm muscles during closed kinetic chain exercises on an unstable surface

Introduction: The purpose of this study was to compare the electromyography index of muscle coactivation of the following muscle pairs: posterior deltoid and pectoralis major (PD/PM); triceps brachii and biceps brachii (TB/BB); and serratus anterior and upper trapezius (SA/UT) during three different closed kinetic chain exercises (wall-press, bench-press and push-up) on an unstable surface at the maximal load. Methods: A total of 20 healthy sedentary men participated in the study. Integral linear values were obtained from three sustained contractions of six seconds each for the three proposed exercises. Mean coactivation index values were compared using the mixed-effects linear model, with a five percent significance level. Results: Electromyography indexes of muscle coactivation showed significant differences for the PD/PM and TB/BB muscle pairs. No differences were found between exercises for the SA/UT muscle pair. Conclusion: Our results seem to differ from those of previous studies, which reported that the similarity in exercises performed is responsible for the comparable muscle activation levels.

Evaluating Patellar Kinematics Through Magnetic Resonance Imaging During Open- and Closed-Kinetic-Chain Exercises

Purpose: To evaluate patellar kinematics of volunteers Without knee pain at rest and during isometric contraction in open- and closed-kinetic-chain exercises. Methods: Twenty individuals took part in this study. All were submitted to magnetic resonance imaging (MRI) during rest and voluntary isometric contraction (VIC) in the open anti closed kinetic chain at 15 degrees, 30 degrees, and 45 degrees of knee flexion. Through MRI and using medical e-film software, the following measurements were evaluated: sulcus angle, patellar-tilt angle, and bisect offset. The mixed-effects linear model was used for comparison between knee positions, between rest and isometric contractions, and between (he exercises. Results: Data analysis revealed that the sulcus angle decreased as knee flexion increased and revealed increases with isometric contractions in both the open and closed kinetic chain for all knee-flexion angles. The patellar-tilt angle decreased with isometric contractions in both the open and closed kinetic chain for every knee position. However, in the closed kinetic chain, patellar tilt increased significantly with the knee flexed at 15 degrees. The bisect offset increased with the knee flexed at 15 degrees during isometric contractions and decreased as knee flexion increased during both exercises. Conclusion: VIC in the last degrees of knee extension may compromise patellar dynamics. On the other hand, it is possible to favor patellar stability by performing muscle contractions with the knee flexed at 30 degrees and 45 degrees in either the open or closed kinetic chain.

Degradation of detergent (linear alkylbenzene sulfonate) in an anaerobic stirred sequencing-batch reactor containing granular biomass

This study aimed to determine the efficiency of an anaerobic stirred sequencing-batch reactor containing granular biomass for the degradation of linear alkylbenzene sulfonate (LAS), a surfactant present in household detergent. The bioreactor was monitored for LAS concentrations in the influent, effluent and sludge, pH, chemical oxygen demand, bicarbonate alkalinity, total solids, and volatile solids. The degradation of LAS was found to be higher in the absence of co-substrates (53%) than in their presence (24-37%). Using the polymerase chain reaction and denaturing gradient gel electrophoresis (PCR/DGGE), we identified populations of microorganisms from the Bacteria and Archaea domains. Among the bacteria, we identified uncultivated populations of Arcanobacterium spp. (94%) and Opitutus spp. (96%). Among the Archaea, we identified Methanospirillum spp. (90%), Methanosaeta spp. (98%), and Methanobacterium spp. (96%). The presence of methanogenic microorganisms shows that LAS did not inhibit anaerobic digestion. Sampling at the last stage of reactor operation recovered 61 clones belonging to the domain bacteria. These represented a variety of phyla: 34% shared significant homology with Bacteroidetes, 18% with Proteobacteria, 11% with Verrucomicrobia, 8% with Fibrobacteres, 2% with Acidobacteria, 3% with Chlorobi and Firmicutes, and 1% with Acidobacteres and Chloroflexi. A small fraction of the clones (13%) were not related to any phylum. Published by Elsevier Ltd.

Model for facilities or vendors location in a global scale considering several echelons in the Chain

The facilities location problem for companies with global operations is very complex and not well explored in the literature. This work proposes a MILP model that solves the problem through minimization of the total logistic cost. Main contributions of the model are the pioneer carrying cost calculation, the treatment given to the take-or-pay costs and to the international tax benefits such as drawback and added value taxes in Brazil. The model was successfully applied to a real case of a chemical industry with industrial plants and sales all over the world. The model application recommended a totally new sourcing model for the company.

Non-linear modal analysis for beams subjected to axial loads: Analytical and finite-element solutions

A rigorous derivation of non-linear equations governing the dynamics of an axially loaded beam is given with a clear focus to develop robust low-dimensional models. Two important loading scenarios were considered, where a structure is subjected to a uniformly distributed axial and a thrust force. These loads are to mimic the main forces acting on an offshore riser, for which an analytical methodology has been developed and applied. In particular, non-linear normal modes (NNMs) and non-linear multi-modes (NMMs) have been constructed by using the method of multiple scales. This is to effectively analyse the transversal vibration responses by monitoring the modal responses and mode interactions. The developed analytical models have been crosschecked against the results from FEM simulation. The FEM model having 26 elements and 77 degrees-of-freedom gave similar results as the low-dimensional (one degree-of-freedom) non-linear oscillator, which was developed by constructing a so-called invariant manifold. The comparisons of the dynamical responses were made in terms of time histories, phase portraits and mode shapes. (C) 2008 Elsevier Ltd. All rights reserved.

Linear minimum mean square filter for discrete-time linear systems with Markov jumps and multiplicative noises

In this paper we obtain the linear minimum mean square estimator (LMMSE) for discrete-time linear systems subject to state and measurement multiplicative noises and Markov jumps on the parameters. It is assumed that the Markov chain is not available. By using geometric arguments we obtain a Kalman type filter conveniently implementable in a recurrence form. The stationary case is also studied and a proof for the convergence of the error covariance matrix of the LMMSE to a stationary value under the assumption of mean square stability of the system and ergodicity of the associated Markov chain is obtained. It is shown that there exists a unique positive semi-definite solution for the stationary Riccati-like filter equation and, moreover, this solution is the limit of the error covariance matrix of the LMMSE. The advantage of this scheme is that it is very easy to implement and all calculations can be performed offline. (c) 2011 Elsevier Ltd. All rights reserved.

Double generalized linear model for tissue culture proportion data: a Bayesian perspective

Joint generalized linear models and double generalized linear models (DGLMs) were designed to model outcomes for which the variability can be explained using factors and/or covariates. When such factors operate, the usual normal regression models, which inherently exhibit constant variance, will under-represent variation in the data and hence may lead to erroneous inferences. For count and proportion data, such noise factors can generate a so-called overdispersion effect, and the use of binomial and Poisson models underestimates the variability and, consequently, incorrectly indicate significant effects. In this manuscript, we propose a DGLM from a Bayesian perspective, focusing on the case of proportion data, where the overdispersion can be modeled using a random effect that depends on some noise factors. The posterior joint density function was sampled using Monte Carlo Markov Chain algorithms, allowing inferences over the model parameters. An application to a data set on apple tissue culture is presented, for which it is shown that the Bayesian approach is quite feasible, even when limited prior information is available, thereby generating valuable insight for the researcher about its experimental results.

Interaction of RTD-1, a cyclic antimicrobial defensin from Rhesus macaque leukocytes, and its open chain analogue with membrane mimics

In contrast to other mammalian defensins, rhesus theta defensin-1 (RTD-1) is composed of just 18 amino acids with the backbone cyclized through peptide bonds. Antibacterial activities of both the native cyclic peptide and a linear form were examined, showing that the cyclic form was 3-fold more active than the open chain analogue, oRTD-1, although both peptides adopt very similar structures in water. It was suggested that the additional charges at the termini of oRTD-1 are the cause for its lower antimicrobial activity. Therefore, we studied the interaction of both peptides with membrane mimics composed of zwitterionic (PC) and negatively charged (PG) phospholipids, major lipid components of erythrocyte and bacterial cell membranes, respectively. Microcalorimetry showed that RTD-1 and oRTD-1 did not affect the phase behavior of PC liposomes, while in PG liposomes both peptides induced new phase transitions above the chain melting transition of the lipid. The shape and fraction differed between both peptides, depending also on their concentration, which will be discussed in terms of their antimicrobial activity.

The effect of long-chain polyunsaturated fatty acids intake during pregnancy on adiposity of healthy full-term offspring at birth

Objective - The adjusted effect of long-chain polyunsaturated fatty acid (LCPUFA) intake during pregnancy on adiposity at birth of healthy full-term appropriate-for-gestational age neonates was evaluated. Study Design - In a cross-sectional convenience sample of 100 mother and infant dyads, LCPUFA intake during pregnancy was assessed by food frequency questionnaire with nutrient intake calculated using Food Processor Plus. Linear regression models for neonatal body composition measurements, assessed by air displacement plethysmography and anthropometry, were adjusted for maternal LCPUFA intakes, energy and macronutrient intakes, prepregnancy body mass index and gestational weight gain. Result - Positive associations between maternal docosahexaenoic acid intake and ponderal index in male offspring (β=0.165; 95% confidence interval (CI): 0.031–0.299; P=0.017), and between n-6:n-3 LCPUFA ratio intake and fat mass (β=0.021; 95% CI: 0.002–0.041; P=0.034) and percentage of fat mass (β=0.636; 95% CI: 0.125–1.147; P=0.016) in female offspring were found. Conclusion - Using a reliable validated method to assess body composition, adjusted positive associations between maternal docosahexaenoic acid intake and birth size in male offspring and between n-6:n-3 LCPUFA ratio intake and adiposity in female offspring were found, suggesting that maternal LCPUFA intake strongly influences fetal body composition.

Dynamics of a backlash chain

This paper studies the dynamical properties of a system with distributed backlash and impact phenomena. This means that it is considered a chain of masses that interact with each other solely by means of backlash and impact phenomena. It is observed the emergence of non-linear phenomena resembling those encountered in the Fermi-Pasta-Ulam problem.

Tetranuclear COPPER(II) complexes with macrocyclic and open-chain disiloxane ligands as catalyst precursors for hydrocarboxylation and oxidation of alkanes and 1-phenylethanol

Two new tetranuclear complexes [Cu-4(mu-O)(L-1)-Cl-4] and [Cu-4(mu(4)-O)(L-2)(2)Cl-4] (2), where H2L1 is a macrocyclic ligand resulting from [2+2] condensation of 2,6-diformy1-4-methylphanol (DFF) and 1,3-bis(aminopropy1)tetramethyldisiloxane, and HL2 is a 1:2 condensation product: of DFF with trimethylsilyl p-aminobenzoate, have been prepared. The structures of the products were established by Xray diffraction. The complexes have been characterised by FTIR, UV/Vis spectroscopy, ES1 mass-spectrometry and magnetic susceptibility measurements. The latter revealed that the letrftriuclear complexes can be descr bed as two ferromagnetically coupled dinuclear units, in which the two copper(II) ions interact antiferromacinetically. The ccimpi.iunds act as homogeneous catalyst precursors for a number of single-pot reactions, including (I) hydrocarbaxylation, with CO, H2O and K2S2O8, of a variety of linear and cyclic (n = 5-8) alkanes into the corresponding Cn+1 carboxylic acids, (ii) peroxidative oxidation of cyclohexane, and (iii) solvent-free microwave-assisted oxidation of 1-phenyletha.nol.

Supply chain optimization of residual forestry biomass for bioenergy production: the case study of Portugal

Within a large set of renewable energies being explored to tackle energy sourcing problems, bioenergy can represent an attractive solution if effectively managed. The supply chain design supported by mathematical programming can be used as a decision support tool to the successful bioenergy production systems establishment. This strategic decision problem is addressed in this paper where we intent to study the design of the residual forestry biomass to bioelectricity production in the Portuguese context. In order to contribute to attain better solutions a mixed integer linear programming (MILP) model is developed and applied in order to optimize the design and planning of the bioenergy supply chain. While minimizing the total supply chain cost the production energy facilities capacity and location are defined. The model also includes the optimal selection of biomass amounts and sources, the transportation modes selection, and links that must be established for biomass transportation and products delivers to markets. Results illustrate the positive contribution of the mathematical programming approach to achieve viable economic solutions. Sensitivity analysis on the most uncertain parameters was performed: biomass availability, transportation costs, fixed operating costs and investment costs. (C) 2015 Elsevier Ltd. All rights reserved.

The Effect of Long-Chain Polyunsaturated Fatty Acids Intake During Pregnancy on Adiposity of Healthy Full-Term Offspring at Birth

OBJECTIVE: The adjusted effect of long-chain polyunsaturated fatty acid (LCPUFA) intake during pregnancy on adiposity at birth of healthy full-term appropriate-for-gestational age neonates was evaluated. STUDY DESIGN: In a cross-sectional convenience sample of 100 mother and infant dyads, LCPUFA intake during pregnancy was assessed by food frequency questionnaire with nutrient intake calculated using Food Processor Plus. Linear regression models for neonatal body composition measurements, assessed by air displacement plethysmography and anthropometry, were adjusted for maternal LCPUFA intakes, energy and macronutrient intakes, prepregnancy body mass index and gestational weight gain. RESULT: Positive associations between maternal docosahexaenoic acid intake and ponderal index in male offspring (β=0.165; 95% confidence interval (CI): 0.031-0.299; P=0.017), and between n-6:n-3 LCPUFA ratio intake and fat mass (β=0.021; 95% CI: 0.002-0.041; P=0.034) and percentage of fat mass (β=0.636; 95% CI: 0.125-1.147; P=0.016) in female offspring were found. CONCLUSION: Using a reliable validated method to assess body composition, adjusted positive associations between maternal docosahexaenoic acid intake and birth size in male offspring and between n-6:n-3 LCPUFA ratio intake and adiposity in female offspring were found, suggesting that maternal LCPUFA intake strongly influences fetal body composition.