Expectation Propagation with Factorizing Distributions: A Gaussian Approximation and Performance Results for Simple Models

We discuss the expectation propagation (EP) algorithm for approximate Bayesian inference using a factorizing posterior approximation. For neural network models, we use a central limit theorem argument to make EP tractable when the number of parameters is large. For two types of models, we show that EP can achieve optimal generalization performance when data are drawn from a simple distribution.

Modelling of species in hood fires by conditional moment closure

Carbon monoxide, the chief killer in fires, and other species are modelled for a series of enclosure fires. The conditions emulate building fires where CO is formed in the rich, turbulent, nonpremixed flame and is transported frozen to lean mixtures by the ceiling jet which is cooled by radiation and dilution. Conditional moment closure modelling is used and computational domain minimisation criteria are developed which reduce the computational cost of this method. The predictions give good agreement for CO and other species in the lean, quenched-gas stream, holding promise that this method may provide a practical means of modelling real, three-dimensional fire situations. (c) 2005 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Conditional variance reduction by measurements on correlated field modes

We consider the case of two cavity modes of the electromagnetic field, which are coupled via the action of a parametric amplifier. The fields are allowed to leak from the cavity and homodyne measurement is performed on one of the modes. Because of the correlations between the modes, this leads to a reduction of the variance in a quadrature of the other mode, although no measurement is performed on it directly. We discuss how this relates to the Einstein-Podolky-Rosen Gedankenexperiment.

Limiting conditional distributions for birth-death processes

In a recent paper [16], one of us identified all of the quasi-stationary distributions for a non-explosive, evanescent birth-death process for which absorption is certain, and established conditions for the existence of the corresponding limiting conditional distributions. Our purpose is to extend these results in a number of directions. We shall consider separately two cases depending on whether or not the process is evanescent. In the former case we shall relax the condition that absorption is certain. Furthermore, we shall allow for the possibility that the minimal process might be explosive, so that the transition rates alone will not necessarily determine the birth-death process uniquely. Although we shall be concerned mainly with the minimal process, our most general results hold for any birth-death process whose transition probabilities satisfy both the backward and the forward Kolmogorov differential equations.

Evaluation of computational methods for the reconstruction of HLA haplotypes

Human leukocyte antigen (HLA) haplotypes are frequently evaluated for population history inferences and association studies. However, the available typing techniques for the main HLA loci usually do not allow the determination of the allele phase and the constitution of a haplotype, which may be obtained by a very time-consuming and expensive family-based segregation study. Without the family-based study, computational inference by probabilistic models is necessary to obtain haplotypes. Several authors have used the expectation-maximization (EM) algorithm to determine HLA haplotypes, but high levels of erroneous inferences are expected because of the genetic distance among the main HLA loci and the presence of several recombination hotspots. In order to evaluate the efficiency of computational inference methods, 763 unrelated individuals stratified into three different datasets had their haplotypes manually defined in a family-based study of HLA-A, -B, -DRB1 and -DQB1 segregation, and these haplotypes were compared with the data obtained by the following three methods: the Expectation-Maximization (EM) and Excoffier-Laval-Balding (ELB) algorithms using the arlequin 3.11 software, and the PHASE method. When comparing the methods, we observed that all algorithms showed a poor performance for haplotype reconstruction with distant loci, estimating incorrect haplotypes for 38%-57% of the samples considering all algorithms and datasets. We suggest that computational haplotype inferences involving low-resolution HLA-A, HLA-B, HLA-DRB1 and HLA-DQB1 haplotypes should be considered with caution.

Serial conditional discrimination and temporal bisection in rats selectively lesioned in the dentate gyrus

In positive serial conditional discrimination, animals respond during a target stimulus when it is preceded by a feature stimulus, but they do not respond when the same target stimulus is presented alone. Moreover, the feature and target stimuli are separated from each other by an empty interval. The present work aimed to investigate if two durations (4 or 16s) of the same feature stimulus (light) could modulate the operant responses of rats to different levers (A and B) during a 5-s target stimulus (tone). In the present study, lever A was associated with the 4-s light, and lever B was associated with the 16-s light. A 5-s empty interval was included between the light and the tone. In the same training procedure, the rats were also presented with the 5-s tone without the preceding light stimuli. In these trials, the responses were not reinforced. We evaluated the hippocampal involvement of these behavioral processes by selectively lesioning the dentate gyrus with colchicine. Once trained, the rats were submitted to a test using probe trials without reinforcement. They were presented with intermediate durations of the feature stimulus (light) to obtain a temporal bisection curve recorded during the exposure to the target stimuli. The rats from both groups learned to respond with high rates during tones preceded by light and with low rates during tones presented alone, which indicated acquisition of the serial conditional discrimination. The rats were able to discriminate between the 4- and 16-s lights by correctly choosing lever A or B. In the test, the temporal bisection curves from both experimental groups showed a bisection point at the arithmetic mean between 4 and 16s. Such processes were not impaired by the dentate gyrus lesion. Thus, our results showed that different durations of a feature stimulus could result in conditional properties. However, this processing did not appear to depend on the dentate gyrus alone. (C) 2011 Published by Elsevier B.V.

Rat forming incisor requires a rigorous ECM remodeling modulated by MMP/RECK balance

Reversion-inducing-cysteine-rich protein with Kazal motifs (RECK) is a single membrane-anchored MMP-regulator and regulates matrix metalloproteinases (MMP) 2, 9 and 14. In turn, MMPs are endopeptidases that play a pivotal role in remodeling ECM. In this work, we decided to evaluate expression pattern of RECK in growing rat incisor during, specifically focusing out amelogenesis process. Based on different kinds of ameloblasts, our results showed that RECK expression was conducted by secretory and post-secretory ameloblasts. At the secretory phase, RECK was localized in the infra-nuclear region of the ameloblast, outer epithelium, near blood vessels, and in the stellate reticulum. From the transition to the maturation phases, RECK was strongly expressed by non-epithelial immuno-competent cells (macrophages and/or dendritic-like cells) in the papillary layer. From the transition to the maturation stage, RECK expression was increased. RECK mRNA was amplified by RT-PCR from whole enamel organ. Here, we verified the presence of RECK mRNA during all stages of amelogenesis. These events were governed by ameloblasts and by non-epithelial cells residents in the enamel organ. Concluding, we found differential expression of MMPs-2, -9 and RECK in the different phases of amelogenesis, suggesting that the tissue remodeling is rigorously controlled during dental mineralization.

On the relation between the conditional moment closure and unsteady flamelets

We consider the relation between the conditional moment closure (CMC) and the unsteady flamelet model (FM). The CMC equations were originally constructed as global equations, while FM was derived asymptotically for a thin reaction zone. The recent tendency is to use FM-type equations as global equations. We investigate the possible consequences and suggest a new version of FM: coordinate-invariant FM (CIFM). Unlike FM, CIFM complies with conditional properties of the exact transport equations which are used effectively in CMC. We analyse the assumptions needed to obtain another global version of FM: representative interactive flamelets (RIF), from original FM and demonstrate that, in homogeneous turbulence, one of these assumptions is equivalent to the main CMC hypothesis.

Fitting mixtures of Kent distributions to aid in joint set identification

When examining a rock mass, joint sets and their orientations can play a significant role with regard to how the rock mass will behave. To identify joint sets present in the rock mass, the orientation of individual fracture planer can be measured on exposed rock faces and the resulting data can be examined for heterogeneity. In this article, the expectation-maximization algorithm is used to lit mixtures of Kent component distributions to the fracture data to aid in the identification of joint sets. An additional uniform component is also included in the model to accommodate the noise present in the data.

Teleportation improvement by conditional measurements on the two-mode squeezed vacuum

We show that by making conditional measurements on the Einstein-Podolsky-Rosen (EPR) squeezed vacuum [T. Opatrny, G. Kurizki, and D.-G. Welsch, Phys. Rev. A 61, 032302 (2000)], one can improve the efficacy of teleportation for both the position-difference, momentum-sum, and number-difference, phase-sum continuous variable teleportation protocols. We investigate the relative abilities of the standard and conditional EPR states, and show that by conditioning we can improve the fidelity of teleportation of coherent states from below to above the (F) over bar =2/3 boundary, thereby achieving unambiguously quantum teleportation.

Growth of subtropical ECM fungi with different nitrogen sources using a new floating culture technique

Eight species of ectomycorrhizal (ECM) fungi in the genera Amanita. Gymnoboletus, Lactarius, and Russula were isolated from subtropical plant communities in eastem Australia. Two species were isolated from each of rainforest, Nothofagus forest, Eucalyptus forest, and Eucalyptus dominated wallum (heath) forest. These communities differ strongly in their soluble soil nitrogen (N) composition. The ability of the fungi to use inorganic (nitrate, ammonium) and organic (amide, peptide, protein) nitrogen sources was determined. As the fungi did not grow in liquid culture, a 'floating culture' technique was devised that allows hyphal growth on a screen floating on liquid medium. With some exceptions, fungal biomass production in floating culture closely reflected fungal growth on solid media assessed by total colony glucosamine content. Most isolates grown in floating culture had similar glucosamine concentrations on all N sources, with isolate specific concentrations ranging from 6 to 12 mug glucosamine g(-1) DW. However, Russula spp. had up to 1.7-fold higher glucosamine concentrations when growing with glutamine or ammonium compared to nitrate, glutathione or protein. Floating cultures supplied with 0.5, 1.5. 4.5, or 10 mm N mostly produced greatest biomass with 4.5 mM N. In vitro nitrate reductase activity (NRA) ranged from very low (0.03 mumol NO2- g(-1) fw h(-1)) in Russula sp. (wallum) to high (2.16 mumol NO2- g(-1) fw h(-1)) in Gymnoboletus sp. (rainforest) and mirrored the fungi's ability to use nitrate as a N source. All Russula spp. (wallum, Nothofagus and Eucalyptus forests), Lactarills sp, (rainforest) and.4manita sp. (wallum) utilized ammonium and glutamine but had little ability to use other N sources. In contrast,Amanita species (Nothofagus and Eucalyptus forests) grew on all N sources but produced most biomass with ammonium and glutamine. Only Gymnoboletus sp. (rainforest) showed similar growth with nitrate and ammonium as N sources. Fungal N source use was not associated with taxonomic groups, but is discussed in the context of soil N sources in the different habitats.

Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.

Effect of instructed extinction on verbal and autonomic indices of Pavlovian learning with fear-relevant and fear-irrelevant conditional stimuli

We investigated the effects of conditional stimulus fear-relevance and of instructed extinction on human Pavlovian conditioning as indexed by electrodermal responses and verbal ratings of conditional stimulus unpleasantness. Half of the participants (n = 64) were trained with pictures of snakes and spiders (fear-relevant) as conditional stimuli, whereas the others were trained with pictures of flowers and mushrooms (fear-irrelevant) in a differential aversive Pavlovian conditioning procedure. Half of the participants in each group were instructed after the completion of acquisition that no more unconditional stimuli were to be presented. Extinction of differential electrodermal responses required more trials after training with fear-relevant pictures. Moreover, there was some evidence that verbal instructions did not affect extinction of second interval electrodermal responses to fear-relevant pictures. However, neither fear-relevance nor instructions affected the changes in rated conditional stimulus pleasantness. This dissociation across measures is interpreted as reflecting renewal of Pavlovian learning.