Variational approximation for mixtures of linear mixed models

Mixtures of linear mixed models (MLMMs) are useful for clustering grouped data and can be estimated by likelihood maximization through the EM algorithm. The conventional approach to determining a suitable number of components is to compare different mixture models using penalized log-likelihood criteria such as BIC.We propose fitting MLMMs with variational methods which can perform parameter estimation and model selection simultaneously. A variational approximation is described where the variational lower bound and parameter updates are in closed form, allowing fast evaluation. A new variational greedy algorithm is developed for model selection and learning of the mixture components. This approach allows an automatic initialization of the algorithm and returns a plausible number of mixture components automatically. In cases of weak identifiability of certain model parameters, we use hierarchical centering to reparametrize the model and show empirically that there is a gain in efficiency by variational algorithms similar to that in MCMC algorithms. Related to this, we prove that the approximate rate of convergence of variational algorithms by Gaussian approximation is equal to that of the corresponding Gibbs sampler which suggests that reparametrizations can lead to improved convergence in variational algorithms as well.Comment: 36 pages, 5 figures, 2 tables, submitted to JCG

Semiclassical Analysis of Extended Dynamical Mean Field Equations

The extended Dynamical Mean Field Equations (EDMFT) are analyzed using semiclassical methods for a model describing an interacting fermi-bose system. We compare the semiclassical approach with the exact QMC (Quantum Montecarlo) method. We found the transition to an ordered state to be of the first order for any dimension below four.Comment: RevTex, 39 pages, 16 figures; Appendix C added, typos correcte

Mammals, Serra da Concórdia, state of Rio de Janeiro, Brazil.

The Atlantic Forest in the state of Rio de Janeiro is extremely fragmented and reduced in its originalextension. Although there are a great number of studies related to the biome, few lists are available for this region. Thepresent study seeks to get information on the mammals of a remnant forest in Serra da Concórdia, an area from whichno information is currently available. The rapid mammal's survey was accomplished by use of pitfalls, live traps,transects, mist nets, and camera traps. We recorded 37 terrestrial species that represent 22.3 % of Rio de Janeiro statemammal species. Among the recorded species, six are endemic of the Atlantic Forest, one is endemic of the state, andeight are in the state list of threatened or presumably threatened species. The area presents great mammal richness andan increasing effort will probably reveal the presence of additional species

Self-consistent quantal treatment of decay rates within the perturbed static path approximation

The framework of the Perturbed Static Path Approximation (PSPA) is used to calculate the partition function of a finite Fermi system from a Hamiltonian with a separable two body interaction. Therein, the collective degree of freedom is introduced in self-consistent fashion through a Hubbard-Stratonovich transformation. In this way all transport coefficients which dominate the decay of a meta-stable system are defined and calculated microscopically. Otherwise the same formalism is applied as in the Caldeira-Leggett model to deduce the decay rate from the free energy above the so called crossover temperature $T_0$.Comment: 17 pages, LaTex, no figures; final version, accepted for publication in PRE; e-mail: [email protected]

Universal Parametric Correlations of Conductance Peaks in Quantum Dots

We compute the parametric correlation function of the conductance peaks in chaotic and weakly disordered quantum dots in the Coulomb blockade regime and demonstrate its universality upon an appropriate scaling of the parameter. For a symmetric dot we show that this correlation function is affected by breaking time-reversal symmetry but is independent of the details of the channels in the external leads. We derive a new scaling which depends on the eigenfunctions alone and can be extracted directly from the conductance peak heights. Our results are in excellent agreement with model simulations of a disordered quantum dot.Comment: 12 pages, RevTex, 2 Postscript figure

A Uniform Approximation for the Fidelity in Chaotic Systems

In quantum/wave systems with chaotic classical analogs, wavefunctions evolve in highly complex, yet deterministic ways. A slight perturbation of the system, though, will cause the evolution to diverge from its original behavior increasingly with time. This divergence can be measured by the fidelity, which is defined as the squared overlap of the two time evolved states. For chaotic systems, two main decay regimes of either Gaussian or exponential behavior have been identified depending on the strength of the perturbation. For perturbation strengths intermediate between the two regimes, the fidelity displays both forms of decay. By applying a complementary combination of random matrix and semiclassical theory, a uniform approximation can be derived that covers the full range of perturbation strengths. The time dependence is entirely fixed by the density of states and the so-called transition parameter, which can be related to the phase space volume of the system and the classical action diffusion constant, respectively. The accuracy of the approximations are illustrated with the standard map.Comment: 16 pages, 4 figures, accepted in J. Phys. A, special edition on Random Matrix Theor

Acute mesenteric ischaemia in refractory shock on veno-arterial extracorporeal membrane oxygenation

Background: Acute mesenteric ischaemia is a severe complication in critically ill patients, but has never been evaluated in patients on veno-arterial extracorporeal membrane oxygenation (V-A ECMO). This study was designed to determine the prevalence of mesenteric ischaemia in patients supported by V-A ECMO and to evaluate its risk factors, as well as to appreciate therapeutic modalities and outcome. Methods: In a retrospective single centre study (January 2013 to January 2017), all consecutive adult patients who underwent V-A ECMO were included, with exclusion of those dying in the first 24 hours. Diagnosis of mesenteric ischaemia was performed using digestive endoscopy, computed tomography scan or first-line laparotomy. Results: One hundred and fifty V-A ECMOs were implanted (65 for post-cardiotomy shock, 85 for acute cardiogenic shock, including 39 patients after refractory cardiac arrest). Overall, median age was 58 (48-69) years and mortality 56%. Acute mesenteric ischaemia was suspected in 38 patients, with a delay of four (2-7) days after ECMO implantation, and confirmed in 14 patients, that is, a prevalence of 9%. Exploratory laparotomy was performed in six out of 14 patients, the others being too unstable to undergo surgery. All patients with mesenteric ischaemia died. Independent risk factors for developing mesenteric ischaemia were renal replacement therapy (odds ratio (OR) 4.5, 95% confidence interval (CI) 1.3-15.7, p=0.02) and onset of a second shock within the first five days (OR 7.8, 95% CI 1.5-41.3, p=0.02). Conversely, early initiation of enteral nutrition was negatively associated with mesenteric ischaemia (OR 0.15, 95% CI 0.03-0.69, p=0.02). Conclusions: Acute mesenteric ischaemia is a relatively frequent but dramatic complication among patients on V-A ECMO

Marine dipole–dipole controlled source electromagnetic and coincident-loop transient electromagnetic experiments to detect seafloor massive sulphides: effects of three-dimensional bathymetry

Seafloor massive sulphides (SMSs) are regarded as a potential future resource to satisfy the growing global demand of metals including copper, zinc and gold. Aside from mining and retrieving profitable amounts of massive sulphides from the seafloor, the present challenge is to detect and delineate significant SMS accumulations, which are generally located near mid-ocean ridges and along submarine volcanic arc and backarc spreading centres. Currently, several geophysical technologies are being developed to detect and quantify SMS occurrences that often exhibit measurable contrasts in their physical parameters compared to the surrounding host rock. Here, we use a short, fixed-offset controlled source electromagnetic (CSEM) system and a coincident-loop transient electromagnetic (TEM) system, which in theory allow the detection of SMS in the shallow seafloor due to a significant electrical conductivity contrast to their surroundings. In 2016, CSEM and TEM experiments were carried out at several locations near the Trans-Atlantic Geotraverse hydrothermal field to investigate shallow occurrences of massive sulphides below the seafloor. Measurements were conducted in an area that contains distinct SMS sites located several kilometres off-axis from the Mid-Atlantic ridge, some of which are still connected to hydrothermal activity and others where hydrothermal activity has ceased. Based on the quality of the acquired data, both experiments were operationally successful. However, the data analysis indicates bias caused by three-dimensional (3D) effects of the rough bathymetry in the study area and, thus, data interpretation remains challenging. Therefore, we study the influence of 3D bathymetry for marine CSEM and TEM experiments, focusing on shallow 3D conductors located beneath mound-like structures. We analyse synthetic inversion models for attributes associated with 3D distortions of CSEM and TEM data that are not sufficiently accounted for in conventional 1D (TEM) and 2D (CSEM) interpretation schemes. Before an adequate quantification of SMS in the region is feasible, these 3D effects need to be studied to avoid over/underestimation of SMS using the acquired EM data. The sensitivity of CSEM and TEM to bathymetry is investigated by means of 3D forward modelling, followed by 1D (TEM) and 2D (CSEM) inversion of the synthetic data using realistic error conditions. Subsequently, inversion models of the synthetic 3D data are analysed and compared to models derived from the measured data to illustrate that 3D distortions are evident in the recorded data sets