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.

Spontaneous Interblink Time Distributions in Patients with Graves` Orbitopathy and Normal Subjects

PURPOSE. To determine the shape of spontaneous interblink time interval distributions obtained in a long observation period in normal subjects and patients with Graves` orbitopathy. METHODS. The magnetic search coil technique was used to register the spontaneous blinking activity during 1 hour of video observation of two groups of 10 subjects each (normal controls aged 27-61 years, mean +/- SD = 46.0 +/- 13.6; patients with Graves` orbitopathy aged 33-61 years, mean +/- SD +/- 46.7 +/- 8.9). The spontaneous blink rate of each subject was calculated for the entire period of observation and for 56 five-minute bins. Histograms of the interblink time interval were plotted for each measurement of blink rate. RESULTS. Neither the overall mean blink rate (controls, 19.8 +/- 4.9; Graves`, 17.6 +/- 5.4) nor the interblink time (controls, 5.2 +/- 3.1, Graves`, 7.9 +/- 3.5) differed between the two groups. There was a large variation of both measurements when the 5-minute bins were considered. The interblink time distribution of all subjects was highly positively skewed when the 1-hour period was measured. A significant number of the 5-minute bin distributions deviated from the overall pattern and became symmetric. CONCLUSIONS. The normal blinking process is characterized by highly positively skewed interblink time distributions. This result means that most blinks have a short time interval, and occasionally a small number of blinks have long time intervals. The different patterns of distribution described in the early literature probably represent artifacts because of the small samples analyzed. (Invest Ophthalmol Vis Sci. 2011;52:3419-3424) DOI:10.1167/iovs.10-7060

Drosophila melanogaster lipins are tissue-regulated and developmentally regulated and present specific subcellular distributions

Lipins constitute a novel family of Mg2+-dependent phosphatidate phosphatases that catalyze the dephosphorylation of phosphatidic acid to yield diacylglycerol, an important intermediate in lipid metabolism and cell signaling. Whereas a single lipin is detected in less complex organisms, in mammals there are distinct lipin isoforms and paralogs that are differentially expressed among tissues. Compatible with organism tissue complexity, we show that the single Drosophila Lpin1 ortholog (CG8709, here named DmLpin) expresses at least three isoforms (DmLpinA, DmLpinK and DmLpinJ) in a temporal and spatially regulated manner. The highest levels of lipin in the fat body, where DmLpinA and DmLpinK are expressed, correlate with the highest levels of triacylglycerol (TAG) measured in this tissue. DmLpinK is the most abundant isoform in the central nervous system, where TAG levels are significantly lower than in the fat body. In the testis, where TAG levels are even lower, DmLpinJ is the predominant isoform. Together, these data suggest that DmLpinA might be the isoform that is mainly involved in TAG production, and that DmLpinK and DmLpinJ could perform other cellular functions. In addition, we demonstrate by immunofluorescence that lipins are most strongly labeled in the perinuclear region of the fat body and ventral ganglion cells. In visceral muscles of the larval midgut and adult testis, lipins present a sarcomeric distribution. In the ovary chamber, the lipin signal is concentrated in the internal rim of the ring canal. These specific subcellular localizations of the Drosophila lipins provide the basis for future investigations on putative novel cellular functions of this protein family.

Bayesian nonlinear regression models with scale mixtures of skew-normal distributions: Estimation and case influence diagnostics

The purpose of this paper is to develop a Bayesian analysis for nonlinear regression models under scale mixtures of skew-normal distributions. This novel class of models provides a useful generalization of the symmetrical nonlinear regression models since the error distributions cover both skewness and heavy-tailed distributions such as the skew-t, skew-slash and the skew-contaminated normal distributions. The main advantage of these class of distributions is that they have a nice hierarchical representation that allows the implementation of Markov chain Monte Carlo (MCMC) methods to simulate samples from the joint posterior distribution. In order to examine the robust aspects of this flexible class, against outlying and influential observations, we present a Bayesian case deletion influence diagnostics based on the Kullback-Leibler divergence. Further, some discussions on the model selection criteria are given. The newly developed procedures are illustrated considering two simulations study, and a real data previously analyzed under normal and skew-normal nonlinear regression models. (C) 2010 Elsevier B.V. All rights reserved.

Fragment mass distributions in the fission of heavy nuclei by intermediate- and high-energy probes

Recent experiments have shown that the multimode approach for describing the fission process is compatible with the observed results. Asystematic analysis of the parameters obtained by fitting the fission-fragment mass distribution to the spontaneous and low-energy data has shown that the values for those parameters present a smooth dependence upon the nuclear mass number. In this work, a new methodology is introduced for studying fragment mass distributions through the multimode approach. It is shown that for fission induced by energetic probes (E > 30 MeV) the mass distribution of the fissioning nuclei produced during the intranuclear cascade and evaporation processes must be considered in order to have a realistic description of the fission process. The method is applied to study (208)Pb, (238)U, (239)Np and (241)Am fission induced by protons or photons.

Spallation product distributions and neutron multiplicities for accelerator-driven system using the CRISP code

Neutron multiplicities for several targets and spallation products of proton-induced reactions in thin targets of interest to an accelerator-driven system obtained with the CRISP code have been reported. This code is a Monte Carlo calculation that simulates the intranuclear cascade and evaporationl fission competition processes. Results are compared with experimental data, and agreement between each other can be considered quite satisfactory in a very broad energy range of incitant particles and different targets.

Influence of the halo upon angular distributions for elastic scattering and breakup

The angular distributions for elastic scattering and breakup of halo nuclei are analysed using a near-side/far-side decomposition within the framework of the dynamical eikonal approximation. This analysis is performed for (11)Be impinging on Pb at 69 MeV/nucleon. These distributions exhibit very similar features. In particular they are both near-side dominated, as expected from Coulomb-dominated reactions. The general shape of these distributions is sensitive mostly to the projectile-target interactions, but is also affected by the extension of the halo. This suggests the elastic scattering not to be affected by a loss of flux towards the breakup channel. (C) 2010 Elsevier B.V. All rights reserved.

Challenging measurement of the (16)O+(27)Al elastic and inelastic angular distributions up to large angles

The (16)O+(27)Al elastic and inelastic angular distributions have been measured in a broad angular range (13 degrees < theta(lab) < 52 degrees) at about 100 MeV incident energy. The use of the MAGNEX large acceptance magnetic spectrometer and of the ray-reconstruction analysis technique has been crucial in order to provide, in the same experiment, high-resolution energy spectra and cross-section measurements distributed over more than seven orders of magnitude down to hundreds of nb/sr. (C) 2011 Elsevier B.V. All rights reserved.

Properties and Inference on the Skew-Curved-Symmetric Family of Distributions

In this article, we study some results related to a specific class of distributions, called skew-curved-symmetric family of distributions that depends on a parameter controlling the skewness and kurtosis at the same time. Special elements of this family which are studied include symmetric and well-known asymmetric distributions. General results are given for the score function and the observed information matrix. It is shown that the observed information matrix is always singular for some special cases. We illustrate the flexibility of this class of distributions with an application to a real dataset on characteristics of Australian athletes.

On estimation and influence diagnostics for the Grubbs` model under heavy-tailed distributions

The Grubbs` measurement model is frequently used to compare several measuring devices. It is common to assume that the random terms have a normal distribution. However, such assumption makes the inference vulnerable to outlying observations, whereas scale mixtures of normal distributions have been an interesting alternative to produce robust estimates, keeping the elegancy and simplicity of the maximum likelihood theory. The aim of this paper is to develop an EM-type algorithm for the parameter estimation, and to use the local influence method to assess the robustness aspects of these parameter estimates under some usual perturbation schemes, In order to identify outliers and to criticize the model building we use the local influence procedure in a Study to compare the precision of several thermocouples. (C) 2008 Elsevier B.V. All rights reserved.

QUASISTATIONARY DISTRIBUTIONS AND FLEMING-VIOT PROCESSES IN FINITE SPACES

Consider a continuous-time Markov process with transition rates matrix Q in the state space Lambda boolean OR {0}. In In the associated Fleming-Viot process N particles evolve independently in A with transition rates matrix Q until one of them attempts to jump to state 0. At this moment the particle jumps to one of the positions of the other particles, chosen uniformly at random. When Lambda is finite, we show that the empirical distribution of the particles at a fixed time converges as N -> infinity to the distribution of a single particle at the same time conditioned on not touching {0}. Furthermore, the empirical profile of the unique invariant measure for the Fleming-Viot process with N particles converges as N -> infinity to the unique quasistationary distribution of the one-particle motion. A key element of the approach is to show that the two-particle correlations are of order 1/N.

A compound class of Weibull and power series distributions

In this paper we introduce the Weibull power series (WPS) class of distributions which is obtained by compounding Weibull and power series distributions where the compounding procedure follows same way that was previously carried out by Adamidis and Loukas (1998) This new class of distributions has as a particular case the two-parameter exponential power series (EPS) class of distributions (Chahkandi and Gawk 2009) which contains several lifetime models such as exponential geometric (Adamidis and Loukas 1998) exponential Poisson (Kus 2007) and exponential logarithmic (Tahmasbi and Rezaei 2008) distributions The hazard function of our class can be increasing decreasing and upside down bathtub shaped among others while the hazard function of an EPS distribution is only decreasing We obtain several properties of the WPS distributions such as moments order statistics estimation by maximum likelihood and inference for a large sample Furthermore the EM algorithm is also used to determine the maximum likelihood estimates of the parameters and we discuss maximum entropy characterizations under suitable constraints Special distributions are studied in some detail Applications to two real data sets are given to show the flexibility and potentiality of the new class of distributions (C) 2010 Elsevier B V All rights reserved

Axial ligand influence on geometries, charge distributions and electronic structures of iron tetraazamacrocycle [Fe((II))TIM(X)(Y)](2+) complexes assessed by Density Functional Theory

We employed the Density Functional Theory along with small basis sets, B3LYP/LANL2DZ, for the study of FeTIM complexes with different pairs of axial ligands (CO, H(2)O, NH(3), imidazole and CH(3)CN). These calculations did not result in relevant changes of molecular quantities as bond lengths, vibrational frequencies and electronic populations supporting any significant back-donation to the carbonyl or acetonitrile axial ligands. Moreover, a back-donation mechanism to the macrocycle cannot be used to explain the observed changes in molecular properties along these complexes with CO or CH(3)CN. This work also indicates that complexes with CO show smaller binding energies and are less stable than complexes with CH(3)CN. Further, the electronic band with the largest intensity in the visible region (or close to this region) is associated to the transition from an occupied 3d orbital on iron to an empty pi* orbital located at the macrocycle. The energy of this Metal-to-Ligand Charge Transfer (MLCT) transition shows a linear relation to the total charge of the macrocycle in these complexes as given by Mulliken or Natural Population Analysis (NPA) formalisms. Finally, the macrocycle total charge seems to be influenced by the field induced by the axial ligands. (C) 2011 Elsevier Ltd. All rights reserved.

Influence of specimens' design and manufacturing process on microtensile bond strength to enamel: laboratory and FEA comparison

This study evaluated the effect of specimens' design and manufacturing process on microtensile bond strength, internal stress distributions (Finite Element Analysis - FEA) and specimens' integrity by means of Scanning Electron Microscopy (SEM) and Laser Scanning Confocal Microscopy (LCM). Excite was applied to flat enamel surface and a resin composite build-ups were made incrementally with 1-mm increments of Tetric Ceram. Teeth were cut using a diamond disc or a diamond wire, obtaining 0.8 mm² stick-shaped specimens, or were shaped with a Micro Specimen Former, obtaining dumbbell-shaped specimens (n = 10). Samples were randomly selected for SEM and LCM analysis. Remaining samples underwent microtensile test, and results were analyzed with ANOVA and Tukey test. FEA dumbbell-shaped model resulted in a more homogeneous stress distribution. Nonetheless, they failed under lower bond strengths (21.83 ± 5.44 MPa)c than stick-shaped specimens (sectioned with wire: 42.93 ± 4.77 MPaª; sectioned with disc: 36.62 ± 3.63 MPa b), due to geometric irregularities related to manufacturing process, as noted in microscopic analyzes. It could be concluded that stick-shaped, nontrimmed specimens, sectioned with diamond wire, are preferred for enamel specimens as they can be prepared in a less destructive, easier, and more precise way.

Revision of the species of the genus Cathorops (Siluriformes: Ariidae) from Mesoamerica and the Central American Caribbean, with description of three new species

The ariid genus Cathorops includes species that occur mainly in estuarine and freshwater habitats of the eastern and western coasts of southern Mexico, Central and South America. The species of Cathorops from the Mesoamerica (Atlantic slope) and Caribbean Central America are revised, and three new species are described: C. belizensis from mangrove areas in Belize; C. higuchii from shallow coastal areas and coastal rivers in the Central American Caribbean, from Honduras to Panama; and C. kailolae from río Usumacinta and lago Izabal basins in Mexico and Guatemala. Additionally, C. aguadulce, from the río Papaloapan basin in Mexico, and C. melanopus from the río Motagua basin in Guatemala and Honduras, are redescribed and their geographic distributions are revised.