An Adaptive Interacting Wang-Landau Algorithm for Automatic Density Exploration

While statisticians are well-accustomed to performing exploratory analysis in the modeling stage of an analysis, the notion of conducting preliminary general-purpose exploratory analysis in the Monte Carlo stage (or more generally, the model-fitting stage) of an analysis is an area which we feel deserves much further attention. Towards this aim, this paper proposes a general-purpose algorithm for automatic density exploration. The proposed exploration algorithm combines and expands upon components from various adaptive Markov chain Monte Carlo methods, with the Wang-Landau algorithm at its heart. Additionally, the algorithm is run on interacting parallel chains -- a feature which both decreases computational cost as well as stabilizes the algorithm, improving its ability to explore the density. Performance is studied in several applications. Through a Bayesian variable selection example, the authors demonstrate the convergence gains obtained with interacting chains. The ability of the algorithm's adaptive proposal to induce mode-jumping is illustrated through a trimodal density and a Bayesian mixture modeling application. Lastly, through a 2D Ising model, the authors demonstrate the ability of the algorithm to overcome the high correlations encountered in spatial models.Comment: 33 pages, 20 figures (the supplementary materials are included as appendices

Developmental changes in mesenteric artery reactivity in embryonic and newly hatched chicks

At birth, the intestine becomes the sole site for nutrient absorption requiring a dramatic increase in blood flow. The vascular changes accompanying this transition have been partly characterized in mammals. We investigated, using wire myography, the developmental changes in chick mesenteric artery (MA) reactivity. Rings of the MA from 15-day (E15) and 19-day (E19) chicken embryos (total incubation 21 days) as well as non-fed 0–3-h-old (NH3h) and first-fed 1-day-old (NH1d) newly hatched chicks contracted in response to KCl, norepinephrine (NE), U46619, and endothelin (ET)-1 and relaxed in response to acetylcholine (ACh), sodium nitroprusside (SNP), and forskolin indicating the presence of electro- and pharmaco-mechanical coupling as well as cGMP- and cAMP-mediated relaxation. In ovo development and transition to ex ovo life was accompanied by alterations in the response of the MAs, but a different developmental trajectory was observed for each reactivity pathway tested. Thus, the contractile efficacy of KCl underwent a linear increase (E15 < E19 < NH3h < NH1d). The efficacy of NE and U46619 increased in ovo, but not ex ovo (E15 < E19 = NH3h = NH1d) and the efficacy of ET-1 peaked at E19 (E15 < E19 > NH3h = NH1d). The relaxations elicited by ACh (endothelium-dependent), SNP, and forskolin did not undergo significant developmental changes. In conclusion, the ability of chick MAs to constrict in response to pharmacological stimuli increases during the embryonic period, but no dramatic changes are induced by hatching or the first feeding. Maturation of vasodilator mechanisms precedes that of vasoconstrictor mechanisms. Alterations of the delicate balance between vasoconstrictors and vasodilators may play an important role in perinatal intestinal diseases

Plasmodium vivax Sub-Patent Infections after Radical Treatment Are Common in Peruvian Patients: Results of a 1-Year Prospective Cohort Study

BACKGROUND: There is an increasing body of literature reporting treatment failure of the currently recommended radical treatment of Plasmodium vivax infections. As P. vivax is the main malaria species outside the African continent, emerging tolerance to its radical treatment regime could have major consequences in countries like Peru, where 80% of malaria cases are due to P. vivax. Here we describe the results of a 1-year longitudinal follow up of 51 confirmed P. vivax patients living around Iquitos, Peruvian Amazon, and treated according to the Peruvian national guidelines. METHODOLOGY: Each month a blood sample for microscopy and later genotyping was systematically collected. Recent exposure to infection was estimated by detecting antibodies against the P. vivax circumsporozoite protein (CSP) and all PCR confirmed P. vivax infections were genotyped with 16 polymorphic microsatellites. RESULTS: During a 1-year period, 84 recurrent infections, 22 positive also by microscopy, were identified, with a median survival time to first recurrent infection of 203 days. Most of them (71%) were asymptomatic; in 13 patients the infection persisted undetected by microscopy for several consecutive months. The genotype of mostly recurrent infections differed from that at day 0 while fewer differences were seen between the recurrent infections. The average expected heterozygosity was 0.56. There was strong linkage disequilibrium (I(A) (s) = 0.29, p<1.10(-4)) that remained also when analyzing only the unique haplotypes, suggesting common inbreeding. CONCLUSION: In Peru, the P. vivax recurrent infections were common and displayed a high turnover of parasite genotypes compared to day 0. Plasmodium vivax patients, even when treated according to the national guidelines, may still represent an important parasite reservoir that can maintain transmission. Any elimination effort should consider such a hidden reservoir

Cosmic ray oriented performance studies for the JEM-EUSO first level trigger

JEM-EUSO is a space mission designed to investigate Ultra-High Energy Cosmic Rays and Neutrinos (E > 5 ⋅ 1019 eV) from the International Space Station (ISS). Looking down from above its wide angle telescope is able to observe their air showers and collect such data from a very wide area. Highly specific trigger algorithms are needed to drastically reduce the data load in the presence of both atmospheric and human activity related background light, yet retain the rare cosmic ray events recorded in the telescope. We report the performance in offline testing of the first level trigger algorithm on data from JEM-EUSO prototypes and laboratory measurements observing different light sources: data taken during a high altitude balloon flight over Canada, laser pulses observed from the ground traversing the real atmosphere, and model landscapes reproducing realistic aspect ratios and light conditions as would be seen from the ISS itself. The first level trigger logic successfully kept the trigger rate within the permissible bounds when challenged with artificially produced as well as naturally encountered night sky background fluctuations and while retaining events with general air-shower characteristics

Experimental damage detection of cracked beams by using nonlinear characteristics of forced response

Experimental evaluation of the flexural forced vibrations of a steel cantilever beam having a transverse surface crack extending uniformly along the width of the beam was performed, where an actual fatigue crack was introduced instead - as usual - of a narrow slot. The nonlinear aspects of the dynamic response of the beam under harmonic excitation were considered and the relevant quantitative parameters were evaluated, in order to relate the nonlinear resonances to the presence and size of the crack. To this end, the existence of sub- and super-harmonic components in the Fourier spectra of the acceleration signals was evidenced, and their amplitudes were quantified. In particular, the acceleration signals were measured in different positions along the beam axis and under different forcing levels at the beam tip. The remarkable relevance of the above mentioned nonlinear characteristics, and their substantial independence on force magnitude and measurement point were worthily noted in comparison with the behavior of the intact beam. Thus, a reliable method of damage detection was proposed which was based on simple tests requiring only harmonically forcing and acceleration measuring in any point non-necessarily near the crack Then, the time-history of the acceleration recorded at the beam tip was numerically processed in order to obtain the time-histories of velocity and displacement. The nonlinear features of the forced response were described and given a physical interpretation in order to define parameters suitable for damage detection. The efficiency of such parameters was discussed with respect to the their capability of detecting damage and a procedure for damage detection was proposed which was able to detect even small cracks by using simple instruments. A finite element model of the cantilever beam was finally assembled and tuned in order to numerically simulate the results of the experimental tests. (C) 2012 Elsevier Ltd. All rights reserved