Background
Gene expression levels in a given cell can be influenced by different factors, namely pharmacological or medical treatments. The response to a given stimulus is usually different for different genes and may depend on time. One of the goals of modern molecular biology is the high-throughput identification of genes associated with a particular treatment or a biological process of interest. From methodological and computational point of view, analyzing high-dimensional time course microarray data requires very specific set of tools which are usually not included in standard software packages. Recently, the authors of this paper developed a fully Bayesian approach which allows one to identify differentially expressed genes in a `one-sample' time-course microarray experiment, to rank them and to estimate their expression profiles. The method is based on explicit expressions for calculations and, hence, very computationally efficient.
Results
The software package BATS (Bayesian Analysis of Time Series) presented here implements the methodology described above. It allows an user to automatically identify and rank differentially expressed genes and to estimate their expression profiles when at least 5-6 time points are available. The package has a user-friendly interface. BATS successfully manages various technical difficulties which arise in time-course microarray experiments, such as a small number of observations, non-uniform sampling intervals and replicated or missing data.
Conclusions
BATS is a free user-friendly software for the analysis of both simulated and real microarray time course experiments. The software, the user manual and a brief illustrative example are freely available online at the BATS website: http://www.na.iac.cnr.it/bats
Bayesian approaches
Time course microarray
Software
We introduce some numerical approximations to a quasilinear problem proposed by G. I. Barenblatt to describe non-equilibrium two-phase fluid flows in permeable porous media, which apply to secondary oil recovery from natural reservoirs. Taking into account the theoretical results of global existence and uniqueness, we approximate the solutions by three numerical schemes, namely, the Diagonal First Order schemes (DFO and DFO2) and the Diagonal Second Order scheme (DSO). For DFO schemes convergence is proved. The schemes' behaviour is analysed and discussed through some numerical experiments.
In this paper we introduce a computation algorithm to trace car paths on road networks, whose load evolution is modeled by conservation laws. This algorithm is composed of two parts: computation of solutions to conservation equations on each road and localization of car position resulting by interactions with waves produced on roads. Some applications and examples to describe the behavior of a driver traveling in a road network are shown. Moreover, a convergence result for wave front tracking approximate solutions, with BV initial data on a single road, is established.
We consider existence and qualitative properties of traveling wave solutions of a new free boundary problem which describes fluid flow in diatomite rocks and in particular the phenomenon of hydraulic fracturing. For certain parameter values discontinuities of the traveling waves may appear near their free boundaries.
We provide Monte Carlo estimators for the derivatives of functionals of Poisson-driven systems, and apply the theoretical results to models of interest in stochastic geometry and insurance
The paper tackles the problem of approximately reconstructing a real function defined on the surface of the unit sphere in the Euclidean q-dimensional space, with q>1, starting from function's samples at scattered sites. Two new operators are introduced for continuous and discrete approximation at scattered sites. Moreover precise error estimates as well as Marcinkiewicz-Zygmund inequalities are derived in every Lp space, giving concrete bounds for all the involved constants.
The problem of estimating a complex measure made up by a linear combination of Dirac distributions centered on points of the complex plane from a finite number of its complex moments affected by additive i.i.d. Gaussian noise is considered. A random measure is defined whose expectation approximates the unknown measure under suitable conditions. An estimator of the approximating measure is then proposed as well as a new discrete transform of the noisy moments that allows computing in estimate of the Unknown measure. A small simulation study is also performed to experimentally check, the goodness of the approximations.