Squared Extrapolation Methods (SQUAREM): A New Class of Simple and Efficient Numerical Schemes for Accelerating the Convergence of the EM Algorithm

We derive a new class of iterative schemes for accelerating the convergence of the EM algorithm, by exploiting the connection between fixed point iterations and extrapolation methods. First, we present a general formulation of one-step iterative schemes, which are obtained by cycling with the extrapolation methods. We, then square the one-step schemes to obtain the new class of methods, which we call SQUAREM. Squaring a one-step iterative scheme is simply applying it twice within each cycle of the extrapolation method. Here we focus on the first order or rank-one extrapolation methods for two reasons, (1) simplicity, and (2) computational efficiency. In particular, we study two first order extrapolation methods, the reduced rank extrapolation (RRE1) and minimal polynomial extrapolation (MPE1). The convergence of the new schemes, both one-step and squared, is non-monotonic with respect to the residual norm. The first order one-step and SQUAREM schemes are linearly convergent, like the EM algorithm but they have a faster rate of convergence. We demonstrate, through five different examples, the effectiveness of the first order SQUAREM schemes, SqRRE1 and SqMPE1, in accelerating the EM algorithm. The SQUAREM schemes are also shown to be vastly superior to their one-step counterparts, RRE1 and MPE1, in terms of computational efficiency. The proposed extrapolation schemes can fail due to the numerical problems of stagnation and near breakdown. We have developed a new hybrid iterative scheme that combines the RRE1 and MPE1 schemes in such a manner that it overcomes both stagnation and near breakdown. The squared first order hybrid scheme, SqHyb1, emerges as the iterative scheme of choice based on our numerical experiments. It combines the fast convergence of the SqMPE1, while avoiding near breakdowns, with the stability of SqRRE1, while avoiding stagnations. The SQUAREM methods can be incorporated very easily into an existing EM algorithm. They only require the basic EM step for their implementation and do not require any other auxiliary quantities such as the complete data log likelihood, and its gradient or hessian. They are an attractive option in problems with a very large number of parameters, and in problems where the statistical model is complex, the EM algorithm is slow and each EM step is computationally demanding.

Ant community richness and composition across a gradient from Eucalyptus plantations to secondary Atlantic Forest

Secondary forests and exotic tree plantations are expanding across tropical landscapes. However, our current understanding of the value of these human-dominated forest landscapes for invertebrate biodiversity conservation is still very poor. In this paper, we use the leaf-litter ant fauna to assess invertebrate diversity in one commercially managed Eucalyptus plantation (four years old), two abandoned plantations of different regeneration ages (16 and 31 years), and one neighboring secondary Atlantic Forest in Southeastern Brazil. There was a clear gradient in species richness from the secondary forest to the managed Eucalyptus plantation; richness and diversity peaked in secondary forest and in the older regenerating Eucalyptus plantation. Significantly more species were recorded in secondary forest samples than in Eucalyptus plantations, but Eucalyptus plantations had a similar level of richness. Furthermore, a non-metric multidimensional scaling analysis revealed clear differences in species composition between the younger managed Eucalyptus plantation (understory absent) and habitats with sub-developed or developed understory. Eucalyptus plantations were characterized by an assemblage of widespread, generalist species very different from those known to occur in core forest habitats of southeastern Brazil. Our results indicate that while older regenerating Eucalyptus plantations can provide habitat to facilitate the persistence of generalist ant species, it is unlikely to conserve most of the primary forest species, such as specialized predators, Dacetini predators, and nomadic species.

Epilithic community development on artificial reefs deployed along a cross-shelf environmental gradient off Paraná state, southern Brazil

Concrete modules were deployed on the bottom of the 11, 18 and 30 meters isobaths along a cross-shelf hydrographic gradient off Paraná State, Southern Brazil, with the purpose of studying the colonization of sessile epilithic macroinvertebrates on artificial surfaces. After one year of submersion a total of 63 species of epilithic organisms were identified, dominated by Ostrea puelchana, Chthamalus bisinuatus, Balanus cf spongicola, Astrangia cf rathbuni, Didemnum spp, poryphers and bryozoans. Diversity index and percent cover at reef stations placed at 11, 18 and 30 meters isobaths were respectively 2.28 and 66.7%, 2.79 and 96.6% and 1.66 and 77.4%. Differences of general community structure among the three assemblages were not clearly related to the general environmental conditions at the bottom layers near the reef stations. Turbidity and larval abundance are discussed as important factors affecting colonization processes. Results indicate that depths between 15-20 meters are more suitable for the implementation of large scale artificial reef systems in the inner shelf off Paraná and, possibly, throughout the inner shelves off southern Brazil with similar hydrographic conditions.

Accelerating aging chemotherapy on aged animals: Protective effect from nutraceutical modulation

The aim of this study was to test a novel phytocompound in an experimental model of antitumor-induced immunosuppression. Five groups of mice were considered: young (Y) and aged (A) that were given intraperitoneally 10 doses of cyclophosphamide (CPX, 25mg/kg/bw) or CPX plus (150 mg/kg/bw) of the nutraceutical DTS (Denshichi-Tochiu-Sen), and control. After sacrifice, macrophage chemotaxis and serum levels of IFN-gamma, IL-2, and GM-CSF were determined. Liver and urinary bladder were examined histologically, as were the liver and kidney for redox enzymes. CPX significantly decreased macrophage chemotaxis and all cytokines (p < 0.05, A >> Y). DTS restored macrophage function and cytokine concentration (p < 0.001) and partly improved the necro-inflammatory score and substance P receptor expression in the bladder and the redox status in liver and kidney (p < 0.05). Such data suggest that DTS effectively prevents CPX-induced immune suppression and oxidative-inflammatory damage, which are particularly enhanced in aged organisms.

Transmission electron microscopy for characterization of acrosomal damage after Percoll gradient centrifugation of cryopreserved bovine spermatozoa

The objective of this study was to characterize acrosomal ultrastructure following discontinuous Percoll gradient centrifugation of cryopreserved bovine sperm. Semen was collected from six bulls of different breeds and three ejaculates per bull were evaluated. Frozen semen samples were thawed and the acrosomal region of sperm cells was evaluated by transmission electron microscopy (TEM) before (n = 18) and after (n = 18) Percoll centrifugation. The evaluation of 20 sperm heads from each of the 36 samples analyzed ensured that a large number of cells were investigated. The data were subjected to analysis of variance at a level of significance of 5%. Percoll centrifugation reduced the percentage of sperm exhibiting normal acrosomes (from 61.77 to 30.24%), reduced the percentage of sperm presenting atypical acrosome reactions (from 28.38 to 4.84%) and increased the percentage of sperm exhibiting damage in the acrosome (from 6.14 to 64.26%). The percentage of sperm with typical acrosome reactions was not significantly different before (3.70%) and after (0.67%) centrifugation. TEM distinguished four different types of acrosomal status and enabled ultrastructural characterization of acrosomal injuries. The percentage of sperm exhibiting normal acrosomes decreased and damage in the acrosome was the most frequent acrosomal injury with the Percoll gradient centrifugation protocol utilized.

Constraints on cold dark matter accelerating cosmologies and cluster formation

We discuss the properties of homogeneous and isotropic flat cosmologies in which the present accelerating stage is powered only by the gravitationally induced creation of cold dark matter (CCDM) particles (Omega(m) = 1). For some matter creation rates proposed in the literature, we show that the main cosmological functions such as the scale factor of the universe, the Hubble expansion rate, the growth factor, and the cluster formation rate are analytically defined. The best CCDM scenario has only one free parameter and our joint analysis involving baryonic acoustic oscillations + cosmic microwave background (CMB) + SNe Ia data yields (Omega) over tilde = 0.28 +/- 0.01 (1 sigma), where (Omega) over tilde (m) is the observed matter density parameter. In particular, this implies that the model has no dark energy but the part of the matter that is effectively clustering is in good agreement with the latest determinations from the large- scale structure. The growth of perturbation and the formation of galaxy clusters in such scenarios are also investigated. Despite the fact that both scenarios may share the same Hubble expansion, we find that matter creation cosmologies predict stronger small scale dynamics which implies a faster growth rate of perturbations with respect to the usual Lambda CDM cosmology. Such results point to the possibility of a crucial observational test confronting CCDM with Lambda CDM scenarios through a more detailed analysis involving CMB, weak lensing, as well as the large-scale structure.

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.

Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids

One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory [Perdew-Burke-Ernzerhof (PBE) GGA] and a recently proposed modification designed specifically for solids (PBEsol) are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules, and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints stemming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.

Nonempirical hyper-generalized-gradient functionals constructed from the Lieb-Oxford bound

A simple and completely general representation of the exact exchange-correlation functional of density-functional theory is derived from the universal Lieb-Oxford bound, which holds for any Coulomb-interacting system. This representation leads to an alternative point of view on popular hybrid functionals, providing a rationale for why they work and how they can be constructed. A similar representation of the exact correlation functional allows to construct fully nonempirical hyper-generalized-gradient approximations (HGGAs), radically departing from established paradigms of functional construction. Numerical tests of these HGGAs for atomic and molecular correlation energies and molecular atomization energies show that even simple HGGAs match or outperform state-of-the-art correlation functionals currently used in solid-state physics and quantum chemistry.

Abundance of sedentary consumers and sessile organisms along the wave exposure gradient of subtropical rocky shores of the south-west Atlantic

Sedentary consumers play an important role on populations of prey and, hence, their patterns of abundance, distribution and coexistence on shores are important to evaluate their potential influence on ecosystem dynamics. Here, we aimed to describe their spatio-temporal distribution and abundance in relation to wave exposure in the intertidal rocky shores of the south-west Atlantic to provide a basis for further understanding of ecological processes in this system. The abundance and composition of the functional groups of sessile organisms and sedentary consumers were taken by sampling the intertidal of sheltered and moderately exposed shores during a period of one year. The sublittoral fringe of sheltered areas was dominated by macroalgae, while the low midlittoral was dominated by bare rock and barnacles. In contrast, filter-feeding animals prevailed at exposed shores, probably explaining the higher abundance of the predator Stramonita haemastoma at these locations. Limpets were more abundant at the midlittoral zone of all shores while sea urchins were exclusively found at the sublittoral fringe of moderately exposed shores, therefore, adding grazing pressure on these areas. The results showed patterns of coexistence, distribution and abundance of those organisms in this subtropical area, presumably as a result of wave action, competition and prey availability. It also brought insights on the influence of top-down and bottom-up processes in this area.

Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets

Outgassing of carbon dioxide (CO(2)) from rivers and streams to the atmosphere is a major loss term in the coupled terrestrial-aquatic carbon cycle of major low-gradient river systems (the term ""river system"" encompasses the rivers and streams of all sizes that compose the drainage network in a river basin). However, the magnitude and controls on this important carbon flux are not well quantified. We measured carbon dioxide flux rates (F(CO2)), gas transfer velocity (k), and partial pressures (p(CO2)) in rivers and streams of the Amazon and Mekong river systems in South America and Southeast Asia, respectively. F(CO2) and k values were significantly higher in small rivers and streams (channels <100 m wide) than in large rivers (channels >100 m wide). Small rivers and streams also had substantially higher variability in k values than large rivers. Observed F(CO2) and k values suggest that previous estimates of basinwide CO(2) evasion from tropical rivers and wetlands have been conservative and are likely to be revised upward substantially in the future. Data from the present study combined with data compiled from the literature collectively suggest that the physical control of gas exchange velocities and fluxes in low-gradient river systems makes a transition from the dominance of wind control at the largest spatial scales (in estuaries and river mainstems) toward increasing importance of water current velocity and depth at progressively smaller channel dimensions upstream. These results highlight the importance of incorporating scale-appropriate k values into basinwide models of whole ecosystem carbon balance.

Forest structure and live aboveground biomass variation along an elevational gradient of tropical Atlantic moist forest (Brazil)

Live aboveground biomass (AGB) is an important source of uncertainty in the carbon balance from the tropical regions in part due scarcity of reliable estimates of live AGB and its variation across landscapes and forest types. Studies of forest structure and biomass stocks of Neotropical forests are biased toward Amazonian and Central American sites. In particular, standardized estimates of aboveground biomass stocks for the Brazilian Atlantic forest are rarely available. Notwithstanding the role of environmental variables that control the distribution and abundance of biomass in tropical lowland forests has been the subject of considerable research, the effect of short, steep elevational gradients on tropical forest structure and carbon dynamics is not well known. In order to evaluate forest structure and live AGB variation along an elevational gradient (0-1100 m a.s.l.) of coastal Atlantic Forest in SE Brazil, we carried out a standard census of woody stems >= 4.8 cm dbh in 13 1-ha permanent plots established on four different sites in 2006-2007. Live AGB ranged from 166.3 Mg ha(-1) (bootstrapped 95% CI: 1444,187.0) to 283.2 Mg ha(-1) (bootstrapped 95% CI: 253.0,325.2) and increased with elevation. We found that local-scale topographic variation associated with elevation influences the distribution of trees >50 cm dbh and total live AGB. Across all elevations, we found more stems (64-75%) with limited crown illumination but the largest proportion of the live AGB (68-85%) was stored in stems with highly illuminated or fully exposed crowns. Topography, disturbance and associated changes in light and nutrient supply probably control biomass distribution along this short but representative elevational gradient. Our findings also showed that intact Atlantic forest sites stored substantial amounts of carbon aboveground. The live tree AGB of the stands was found to be lower than Central Amazonian forests, but within the range of Neotropical forests, in particular when compared to Central American forests. Our comparative data suggests that differences in live tree AGB among Neotropical forests are probably related to the heterogeneous distribution of large and medium-sized diameter trees within forests and how the live biomass is partitioned among those size classes, in accordance with general trends found by previous studies. In addition, the elevational variation in live AGB stocks suggests a large spatial variability over coastal Atlantic forests in Brazil, clearly indicating that it is important to consider regional differences in biomass stocks for evaluating the role of this threatened tropical biome in the global carbon cycle. (C) 2010 Elsevier B.V. All rights reserved.

Sequential injections as an alternative to gradient exploitation for implementing differential kinetic analysis in a flow injection system

A novel flow-based strategy for implementing simultaneous determinations of different chemical species reacting with the same reagent(s) at different rates is proposed and applied to the spectrophotometric catalytic determination of iron and vanadium in Fe-V alloys. The method relies on the influence of Fe(II) and V(IV) on the rate of the iodide oxidation by Cr(VI) under acidic conditions, the Jones reducing agent is then needed Three different plugs of the sample are sequentially inserted into an acidic KI reagent carrier stream, and a confluent Cr(VI) solution is added downstream Overlap between the inserted plugs leads to a complex sample zone with several regions of maximal and minimal absorbance values. Measurements performed on these regions reveal the different degrees of reaction development and tend to be more precise Data are treated by multivariate calibration involving the PLS algorithm The proposed system is very simple and rugged Two latent variables carried out ca 95% of the analytical information and the results are in agreement with ICP-OES. (C) 2010 Elsevier B V. All rights reserved.

Accelerating autonomous learning by using heuristic selection of actions

This paper investigates how to make improved action selection for online policy learning in robotic scenarios using reinforcement learning (RL) algorithms. Since finding control policies using any RL algorithm can be very time consuming, we propose to combine RL algorithms with heuristic functions for selecting promising actions during the learning process. With this aim, we investigate the use of heuristics for increasing the rate of convergence of RL algorithms and contribute with a new learning algorithm, Heuristically Accelerated Q-learning (HAQL), which incorporates heuristics for action selection to the Q-Learning algorithm. Experimental results on robot navigation show that the use of even very simple heuristic functions results in significant performance enhancement of the learning rate.

Ultra-Fast Gradient LC Method for Omeprazole Analysis Using a Monolithic Column: Assay Development, Validation, and Application to the Quality Control of Omeprazole Enteric-Coated Pellets

A method was optimized for the analysis of omeprazole (OMZ) by ultra-high speed LC with diode array detection using a monolithic Chromolith Fast Gradient RP 18 endcapped column (50 x 2.0 mm id). The analyses were performed at 30 degrees C using a mobile phase consisting of 0.15% (v/v) trifluoroacetic acid (TFA) in water (solvent A) and 0.15% (v/v) TFA in acetonitrile (solvent B) under a linear gradient of 5 to 90% B in 1 min at a flow rate of 1.0 mL/min and detection at 220 nm. Under these conditions, OMZ retention time was approximately 0.74 min. Validation parameters, such as selectivity, linearity, precision, accuracy, and robustness, showed results within the acceptable criteria. The method developed was successfully applied to OMZ enteric-coated pellets, showing that this assay can be used in the pharmaceutical industry for routine QC analysis. Moreover, the analytical conditions established allow for the simultaneous analysis of OMZ metabolites, 5-hydroxyomeprazole and omeprazole sulfone, in the same run, showing that this method can be extended to other matrixes with adequate procedures for sample preparation.