A partially collapsed Gibbs sampler for Bayesian quantile regression

We introduce a set of new Gibbs sampler for Bayesian analysis of quantile re-gression model. The new algorithm, which partially collapsing an ordinary Gibbs sampler, is called Partially Collapsed Gibbs (PCG) sampler. Although the Metropolis-Hastings algorithm has been employed in Bayesian quantile regression, including median regression, PCG has superior convergence properties to an ordinary Gibbs sampler. Moreover, Our PCG sampler algorithm, which is based on a theoretic derivation of an asymmetric Laplace as scale mixtures of normal distributions, requires less computation than the ordinary Gibbs sampler and can significantly reduce the computation involved in approximating the Bayes Factor and marginal likelihood. Like the ordinary Gibbs sampler, the PCG sample can also be used to calculate any associated marginal and predictive distributions. The quantile regression PCG sampler is illustrated by analysing simulated data and the data of length of stay in hospital. The latter provides new insight into hospital perfor-mance. C-code along with an R interface for our algorithms is publicly available on request from the first author. JEL classification: C11, C14, C21, C31, C52, C53

Noise spectra of stochastic pulse sequences: application to large scale magnetization flips in the finite size 2D Ising model

We provide a general scheme to predict and derive the contribution to the noise spectrum of a stochastic sequence of pulses from the distribution of pulse parameters. An example is the magnetization noise spectra of a 2D Ising system near its phase transition. At $T\le T_c$, the low frequency spectra is dominated by magnetization flips of nearly the entire system. We find that both the predicted and the analytically derived spectra fit those produced from simulations. Subtracting this contribution leaves the high frequency spectra which follow a power law set by the critical exponents.Comment: 4 pages, 5 figures. We improved text and included a predicted noise curve in Figure 4. 2 examples from Figure 3 are remove

Evaluation of aerothermal modeling computer programs

Various computer programs based upon the SIMPLE or SIMPLER algorithm were studied and compared for numerical accuracy, efficiency, and grid dependency. Four two-dimensional and one three-dimensional code originally developed by a number of research groups were considered. In general, the accuracy and computational efficieny of these TEACH type programs were improved by modifying the differencing schemes and their solvers. A brief description of each program is given. Error reduction, spline flux and second upwind differencing programs are covered

Robust variable selection in partially varying coefficient single-index model

By combining basis function approximations and smoothly clipped absolute deviation (SCAD) penalty, this paper proposes a robust variable selection procedure for a partially varying coefficient single-index model based on modal regression. The proposed procedure simultaneously selects significant variables in the parametric components and the nonparametric components. With appropriate selection of the tuning parameters, we establish the theoretical properties of our procedure, including consistency in variable selection and the oracle property in estimation. Furthermore, we also discuss the bandwidth selection and propose a modified expectation-maximization (EM)-type algorithm for the proposed estimation procedure. The finite sample properties of the proposed estimators are illustrated by some simulation examples.The research of Zhu is partially supported by National Natural Science Foundation of China (NNSFC) under Grants 71171075, 71221001 and 71031004. The research of Yu is supported by NNSFC under Grant 11261048

The Deuteron Spin Structure Functions in the Bethe-Salpeter Approach and the Extraction of the Neutron Structure Function g1n(x)g_1^n(x)

The nuclear effects in the spin-dependent structure functions $g_1^D$ and $b_2^D$ are calculated in the relativistic approach based on the Bethe-Salpeter equation with a realistic meson-exchange potential. The results of calculations are compared with the non-relativistic calculations. The problem of extraction of the neutron spin structure function, $g_1^n$, from the deuteron data is discussed.Comment: (Talk given at the SPIN'94 International Symposium, September 15-22, 1994, Bloomington, Indiana), 6 pages, 5 figures, Preprint Alberta Thy 29-9

Absolute and differential measurement of water vapor supersaturation using a commercial thin-film sensor

We describe a technique for measuring the water vapor supersaturation of normal air over a temperature range of –40<~T<~0 °C. The measurements use an inexpensive commercial hygrometer which is based on a thin-film capacitive sensor. The time required for the sensor to reach equilibrium was found to increase exponentially with decreasing sensor temperature, exceeding 2 min for T = –30 °C; however, the water vapor sensitivity of the device remained high down to this temperature. After calibrating our measurement procedure, we found residual scatter in the data corresponding to an uncertainty in the absolute water vapor pressure of about ±15%. This scatter was due mainly to long-term drift, which appeared to be intrinsic to the capacitive thin-film sensor. The origin of this drift is not clear, but it effectively limits the applicability of this instrument for absolute measurements. We also found, however, that the high sensitivity of the thin-film sensor makes it rather well suited for differential measurements. By comparing supersaturated and saturated air at the same temperature we obtained a relative measurement uncertainty of about ±1.5%, an order of magnitude better than the absolute measurements

Field-induced structure transformation in electrorheological solids

We have computed the local electric field in a body-centered tetragonal (BCT) lattice of point dipoles via the Ewald-Kornfeld formulation, in an attempt to examine the effects of a structure transformation on the local field strength. For the ground state of an electrorheological solid of hard spheres, we identified a novel structure transformation from the BCT to the face-centered cubic (FCC) lattices by changing the uniaxial lattice constant c under the hard sphere constraint. In contrast to the previous results, the local field exhibits a non-monotonic transition from BCT to FCC. As c increases from the BCT ground state, the local field initially decreases rapidly towards the isotropic value at the body-centered cubic lattice, decreases further, reaching a minimum value and increases, passing through the isotropic value again at an intermediate lattice, reaches a maximum value and finally decreases to the FCC value. An experimental realization of the structure transformation is suggested. Moreover, the change in the local field can lead to a generalized Clausius-Mossotti equation for the BCT lattices.Comment: Submitted to Phys. Rev.

Fluorescent turn-on probes for wash-free mRNA imaging via covalent site-specific enzymatic labeling.

Investigating the many roles RNA plays in cellular regulation and function has increased demand for tools to explore RNA tracking and localization within cells. Our recently reported RNA-TAG (transglycosylation at guanine) approach uses an RNA-modifying enzyme, tRNA-guanine transglycosylase (TGT), to accomplish covalent labeling of an RNA of interest with fluorescent tracking agents in a highly selective and efficient manner. Unfortunately, labeling by this method currently suffers from a high nonspecific fluorescent background and is currently unsuitable for imaging RNA within complex cellular environments. Herein we report the design and synthesis of novel fluorogenic thiazole orange probes that significantly lower nonspecific binding and background fluorescence and, as a result, provide up to a 100-fold fluorescence intensity increase after labeling. Using these fluorogenic labeling agents, we were able to image mRNA expressed in Chinese Hamster Ovary cells in a wash-free manner

Simple choreographies of the planar Newtonian NN-body Problem

In the $N$-body problem, a simple choreography is a periodic solution, where all masses chase each other on a single loop. In this paper we prove that for the planar Newtonian $N$-body problem with equal masses, $N \ge 3$, there are at least $2^{N-3} + 2^{[(N-3)/2]}$ different main simple choreographies. This confirms a conjecture given by Chenciner and etc. in \cite{CGMS02}.Comment: 31pages, 6 figures. Refinements in notations and proof