Likelihood-free estimation of model evidence

Statistical methods of inference typically require the likelihood function to be computable in a reasonable amount of time. The class of “likelihood-free” methods termed Approximate Bayesian Computation (ABC) is able to eliminate this requirement, replacing the evaluation of the likelihood with simulation from it. Likelihood-free methods have gained in efficiency and popularity in the past few years, following their integration with Markov Chain Monte Carlo (MCMC) and Sequential Monte Carlo (SMC) in order to better explore the parameter space. They have been applied primarily to estimating the parameters of a given model, but can also be used to compare models. Here we present novel likelihood-free approaches to model comparison, based upon the independent estimation of the evidence of each model under study. Key advantages of these approaches over previous techniques are that they allow the exploitation of MCMC or SMC algorithms for exploring the parameter space, and that they do not require a sampler able to mix between models. We validate the proposed methods using a simple exponential family problem before providing a realistic problem from human population genetics: the comparison of different demographic models based upon genetic data from the Y chromosome.

Bayesian parameter estimation for latent Markov random fields and social networks

Undirected graphical models are widely used in statistics, physics and machine vision. However Bayesian parameter estimation for undirected models is extremely challenging, since evaluation of the posterior typically involves the calculation of an intractable normalising constant. This problem has received much attention, but very little of this has focussed on the important practical case where the data consists of noisy or incomplete observations of the underlying hidden structure. This paper specifically addresses this problem, comparing two alternative methodologies. In the first of these approaches particle Markov chain Monte Carlo (Andrieu et al., 2010) is used to efficiently explore the parameter space, combined with the exchange algorithm (Murray et al., 2006) for avoiding the calculation of the intractable normalising constant (a proof showing that this combination targets the correct distribution in found in a supplementary appendix online). This approach is compared with approximate Bayesian computation (Pritchard et al., 1999). Applications to estimating the parameters of Ising models and exponential random graphs from noisy data are presented. Each algorithm used in the paper targets an approximation to the true posterior due to the use of MCMC to simulate from the latent graphical model, in lieu of being able to do this exactly in general. The supplementary appendix also describes the nature of the resulting approximation.

Characterization of a Listeria monocytogenes Scott A isolate with high tolerance towards high hydrostatic pressure

An isolate of L. monocytogenes Scott A that is tolerant to high hydrostatic pressure (HHP), named AK01, was isolated upon a single pressurization treatment of 400 MPa for 20 min and was further characterized. The survival of exponential- and stationary-phase cells of AK01 in ACES [N-(2-acetamido)-2-aminoethanesulfonic acid] buffer was at least 2 log units higher than that of the wild type over a broad range of pressures (150 to 500 MPa), while both strains showed higher HHP tolerance (piezotolerance) in the stationary than in the exponential phase of growth. In semiskim milk, exponential-phase cells of both strains showed lower reductions upon pressurization than in buffer, but again, AK01 was more piezotolerant than the wild type. The piezotolerance of AK01 was retained for at least 40 generations in rich medium, suggesting a stable phenotype. Interestingly, cells of AK01 lacked flagella, were elongated, and showed slightly lower maximum specific growth rates than the wild type at 8, 22, and 30°C. Moreover, the piezotolerant strain AK01 showed increased resistance to heat, acid, and H2O2 compared with the wild type. The difference in HHP tolerance between the piezotolerant strain and the wild-type strain could not be attributed to differences in membrane fluidity, since strain AK01 and the wild type had identical in situ lipid melting curves as determined by Fourier transform infrared spectroscopy. The demonstrated occurrence of a piezotolerant isolate of L. monocytogenes underscores the need to further investigate the mechanisms underlying HHP resistance of food-borne microorganisms, which in turn will contribute to the appropriate design of safe, accurate, and feasible HHP treatments.

Combined action of S-carvone and mild heat treatment on Listeria monocytogenes Scott A

The combined action of the plant-derived volatile, S-carvone, and mild heat treatment on the food-borne pathogen, Listeria monocytogenes, was evaluated. The viability of exponential phase cultures grown at 8 °C could be reduced by 1·3 log units after exposure to S-carvone (5 mmol l−1) for 30 min at 45 °C, while individual treatment with S-carvone or exposure to 45 °C for 30 min did not result in a loss in viability. Other plant-derived volatiles, namely carvacrol, cinnamaldehyde, thymol and decanal, were also found to reduce the viability of L. monocytogenes in combination with the same mild heat treatment at concentrations of 1·75 mmol l−1, 2·5 mmol l−1, 1·5 mmol l−1 and 2 mmol l−1, respectively. These findings show that essential oil compounds can play an important role in minimally processed foods, and can be used in the concept of Hurdle Technology to reduce the intensity of heat treatment or other individual hurdles.

On the time-dependent analysis of Gamow decay

Gamow's explanation of the exponential decay law uses complex 'eigenvalues' and exponentially growing 'eigenfunctions'. This raises the question, how Gamow's description fits into the quantum mechanical description of nature, which is based on real eigenvalues and square integrable wavefunctions. Observing that the time evolution of any wavefunction is given by its expansion in generalized eigenfunctions, we shall answer this question in the most straightforward manner, which at the same time is accessible to graduate students and specialists. Moreover, the presentation can well be used in physics lectures to students.

Heterogeneous surface expression of EspA translocon filaments by Escherichia coli O157:H7 is controlled at the posttranscriptional level

Type III secretion systems of enteric bacteria enable translocation of effector proteins into host cells. Secreted proteins of verotoxigenic Escherichia coli O157 strains include components of a translocation apparatus, EspA, -B, and -D, as well as "effectors" such as the translocated intimin receptor (Tir) and the mitochondrion-associated protein (Map). This research has investigated the regulation of LEE4 translocon proteins, in particular EspA. EspA filaments could not be detected on the bacterial cell surface when E. coli O157:H7 was cultured in M9 minimal medium but were expressed from only a proportion of the bacterial population when cultured in minimal essential medium modified with 25 mM HEPES. The highest proportions of EspA-filamented bacteria were detected in late exponential phase, after which filaments were lost rapidly from the bacterial cell surface. Our previous research had shown that human and bovine E. coli O157:H7 strains exhibit marked differences in EspD secretion levels. Here it is demonstrated that the proportion of the bacterial population expressing EspA filaments was associated with the level of EspD secretion. The ability of individual bacteria to express EspA filaments was not controlled at the level of LEE1-4 operon transcription, as demonstrated by using both beta-galactosidase and green fluorescent protein (GFP) promoter fusions. All bacteria, whether expressing EspA filaments or not, showed equivalent levels of GFP expression when LEEI-4 translational fusions were used. Despite this, the LEE4-espADB mRNA was more abundant from populations with a high proportion of nonsecreting bacteria (low secretors) than from populations with a high proportion of secreting and therefore filamented bacteria (high secretors). This research demonstrates that while specific environmental conditions are required to induce LEEI-4 expression, a further checkpoint exists before EspA filaments are produced on the bacterial surface and secretion of effector proteins occurs. This checkpoint in E. coli O157:H7 translocon expression is controlled by a posttranscriptional mechanism acting on LEE4-espADB mRNA. The heterogeneity in EspA filamentation could arise from phase-variable expression of regulators that control this posttranscriptional mechanism.

Environmental regulation and colonization attributes of the long polar fimbriae (LPF) of Escherichia coli O157:H7

The ability of Escherichia coli O157:H7 to colonize the intestinal epithelia is dependent on the expression of intimin and other adhesins. The chromosome of E. coli O157:H7 carries two loci encoding long polar fimbriae (LPF). These fimbriae mediate adherence to epithelial cells and are associated with colonization of the intestine. In order to increase our knowledge about the conditions controlling their expression and their role in colonization of an animal model, the environmental cues that promote expression of lpf genes and the role of E. coli O157:H7 LPF in intestinal colonization of lambs were investigated. We found that expression of lpf1 was regulated in response to growth phase, osmolarity, and pH; that lpf2 transcription was stimulated during late exponential growth and iron depletion; and that LPF impacts the ability of E. coli O157:H7 to persist in the intestine of infected 6-week-old lambs.

Experimental and modeled UV erthemal irradiance under overcast conditions: the role of cloud optical depth

This paper evaluates the relationship between the cloud modification factor (CMF) in the ultraviolet erythe- mal range and the cloud optical depth (COD) retrieved from the Aerosol Robotic Network (AERONET) "cloud mode" algorithm under overcast cloudy conditions (confirmed with sky images) at Granada, Spain, mainly for non-precipitating, overcast and relatively homogenous water clouds. Empirical CMF showed a clear exponential dependence on experimental COD values, decreasing approximately from 0.7 for COD=10 to 0.25 for COD=50. In addition, these COD measurements were used as input in the LibRadtran radia tive transfer code allowing the simulation of CMF values for the selected overcast cases. The modeled CMF exhibited a dependence on COD similar to the empirical CMF, but modeled values present a strong underestimation with respect to the empirical factors (mean bias of 22 %). To explain this high bias, an exhaustive comparison between modeled and experimental UV erythemal irradiance (UVER) data was performed. The comparison revealed that the radiative transfer simulations were 8 % higher than the observations for clear-sky conditions. The rest of the bias (~14 %) may be attributed to the substantial underestimation of modeled UVER with respect to experimental UVER under overcast conditions, although the correlation between both dataset was high (R2 ~ 0.93). A sensitive test showed that the main reason responsible for that underestimation is the experimental AERONET COD used as input in the simulations, which has been retrieved from zenith radiances in the visible range. In this sense, effective COD in the erythemal interval were derived from an iteration procedure based on searching the best match between modeled and experimental UVER values for each selected overcast case. These effective COD values were smaller than AERONET COD data in about 80 % of the overcast cases with a mean relative difference of 22 %.

Numerical simulation of the Filchner overflow

The plume of Ice Shelf Water (ISW) flowing into the Weddell Sea over the Filchner sill contributes to the formation of Antarctic Bottom Water. The Filchner overflow is simulated using a hydrostatic, primitive equation three-dimensional ocean model with a 0.5–2 Sv ISW influx above the Filchner sill. The best fit to mooring temperature observations is found with influxes of 0.5 and 1 Sv, below a previous estimate of 1.6 ± 0.5 Sv based on sparse mooring velocities. The plume first moves north over the continental shelf, and then turns west, along slope of the continental shelf break where it breaks up into subplumes and domes, some of which then move downslope. Other subplumes run into the eastern submarine ridge and propagate along the ridge downslope in a chaotic manner. The next, western ridge is crossed by the plume through several paths. Despite a number of discrepancies with observational data, the model reproduces many attributes of the flow. In particular, we argue that the temporal variability shown by the observations can largely be attributed to the unstable structure of the flow, where the temperature fluctuations are determined by the motion of the domes past the moorings. Our sensitivity studies show that while thermobaricity plays a role, its effect is small for the flows considered. Smoothing the ridges out demonstrate that their presence strongly affects the plume shape around the ridges. An increase in the bottom drag or viscosity leads to slowing down, and hence thickening and widening of the plume

The maximum rate of mammal evolution

How fast can a mammal evolve from the size of a mouse to the size of an elephant? Achieving such a large transformation calls for major biological reorganization. Thus, the speed at which this occurs has important implications for extensive faunal changes, including adaptive radiations and recovery from mass extinctions. To quantify the pace of large-scale evolution we developed a metric, clade maximum rate, which represents the maximum evolutionary rate of a trait within a clade. We applied this metric to body mass evolution in mammals over the last 70 million years, during which multiple large evolutionary transitions occurred in oceans and on continents and islands. Our computations suggest that it took a minimum of 1.6, 5.1, and 10 million generations for terrestrial mammal mass to increase 100-, and 1,000-, and 5,000- fold, respectively. Values for whales were down to half the length (i.e., 1.1, 3, and 5 million generations), perhaps due to the reduced mechanical constraints of living in an aquatic environment. When differences in generation time are considered, we find an exponential increase in maximum mammal body mass during the 35 million years following the Cretaceous–Paleogene (K–Pg) extinction event. Our results also indicate a basic asymmetry in macroevolution: very large decreases (such as extreme insular dwarfism) can happen at more than 10 times the rate of increases. Our findings allow more rigorous comparisons of microevolutionary and macroevolutionary patterns and processes. Keywords: haldanes, biological time, scaling, pedomorphosis

Sparse shape representation using the Laplace-Beltrami eigenfunctions and its application to modeling subcortical structures

We present a new sparse shape modeling framework on the Laplace-Beltrami (LB) eigenfunctions. Traditionally, the LB-eigenfunctions are used as a basis for intrinsically representing surface shapes by forming a Fourier series expansion. To reduce high frequency noise, only the first few terms are used in the expansion and higher frequency terms are simply thrown away. However, some lower frequency terms may not necessarily contribute significantly in reconstructing the surfaces. Motivated by this idea, we propose to filter out only the significant eigenfunctions by imposing l1-penalty. The new sparse framework can further avoid additional surface-based smoothing often used in the field. The proposed approach is applied in investigating the influence of age (38-79 years) and gender on amygdala and hippocampus shapes in the normal population. In addition, we show how the emotional response is related to the anatomy of the subcortical structures.

Femtosecond pulse shaping using differential evolutionary algorithm and wavelet operators

Evolutionary meta-algorithms for pulse shaping of broadband femtosecond duration laser pulses are proposed. The genetic algorithm searching the evolutionary landscape for desired pulse shapes consists of a population of waveforms (genes), each made from two concatenated vectors, specifying phases and magnitudes, respectively, over a range of frequencies. Frequency domain operators such as mutation, two-point crossover average crossover, polynomial phase mutation, creep and three-point smoothing as well as a time-domain crossover are combined to produce fitter offsprings at each iteration step. The algorithm applies roulette wheel selection; elitists and linear fitness scaling to the gene population. A differential evolution (DE) operator that provides a source of directed mutation and new wavelet operators are proposed. Using properly tuned parameters for DE, the meta-algorithm is used to solve a waveform matching problem. Tuning allows either a greedy directed search near the best known solution or a robust search across the entire parameter space.

Classification and fault detection methods for fuel cell monitoring and quality control

In this paper, various types of fault detection methods for fuel cells are compared. For example, those that use a model based approach or a data driven approach or a combination of the two. The potential advantages and drawbacks of each method are discussed and comparisons between methods are made. In particular, classification algorithms are investigated, which separate a data set into classes or clusters based on some prior knowledge or measure of similarity. In particular, the application of classification methods to vectors of reconstructed currents by magnetic tomography or to vectors of magnetic field measurements directly is explored. Bases are simulated using the finite integration technique (FIT) and regularization techniques are employed to overcome ill-posedness. Fisher's linear discriminant is used to illustrate these concepts. Numerical experiments show that the ill-posedness of the magnetic tomography problem is a part of the classification problem on magnetic field measurements as well. This is independent of the particular working mode of the cell but influenced by the type of faulty behavior that is studied. The numerical results demonstrate the ill-posedness by the exponential decay behavior of the singular values for three examples of fault classes.

High yield production of a soluble bifidobacterial β-galactosidase (BbgIV) inE. coliDH5α with improved catalytic efficiency for the synthesis of prebiotic galactooligosaccharides

The bifidobacterial β-galactosidase (BbgIV) was produced in E. coli DH5α at 37 and 30 °C in a 5 L bioreactor under varied conditions of dissolved oxygen (dO2) and pH. The yield of soluble BbgIV was significantly (P < 0.05) increased once the dO2 dropped to 0–2% and remained at such low values during the exponential phase. Limited dO2 significantly (P < 0.05) increased the plasmid copy number and decreased the cells growth rate. Consequently, the BbgIV yield increased to its maximum (71–75 mg per g dry cell weight), which represented 20–25% of the total soluble proteins in the cells. In addition, the specific activity and catalytic efficiency of BbgIV were significantly (P < 0.05) enhanced under limited dO2 conditions. This was concomitant with a change in the enzyme secondary structure, suggesting a link between the enzyme structure and function. The knowledge generated from this work is very important for producing BbgIV as a biocatalyst for the development of a cost-effective process for the synthesis of prebiotic galactooligosaccharides from lactose.

The structure of paranoia in the general population

Background: Psychotic phenomena appear to form a continuum with normal experience and beliefs, and may build on common emotional interpersonal concerns. Aims: We tested predictions that paranoid ideation is exponentially distributed and hierarchically arranged in the general population, and that persecutory ideas build on more common cognitions of mistrust, interpersonal sensitivity and ideas of reference. Method: Items were chosen from the Structured Clinical Interview for DSM-IV Axis II Disorders (SCID-II) questionnaire and the Psychosis Screening Questionnaire in the second British National Survey of Psychiatric Morbidity (n = 8580), to test a putative hierarchy of paranoid development using confirmatory factor analysis, latent class analysis and factor mixture modelling analysis. Results: Different types of paranoid ideation ranged in frequency from less than 2% to nearly 30%. Total scores on these items followed an almost perfect exponential distribution (r = 0.99). Our four a priori first-order factors were corroborated (interpersonal sensitivity; mistrust; ideas of reference; ideas of persecution). These mapped onto four classes of individual respondents: a rare, severe, persecutory class with high endorsement of all item factors, including persecutory ideation; a quasi-normal class with infrequent endorsement of interpersonal sensitivity, mistrust and ideas of reference, and no ideas of persecution; and two intermediate classes, characterised respectively by relatively high endorsement of items relating to mistrust and to ideas of reference. Conclusions: The paranoia continuum has implications for the aetiology, mechanisms and treatment of psychotic disorders, while confirming the lack of a clear distinction from normal experiences and processes.