60-XX PROBABILITY THEORY AND STOCHASTIC PROCESSES (For additional applications, see 11Kxx, 62-XX, 90-XX, 91-XX, 92-XX, 93-XX, 94-XX)
Refine
Year of publication
Document Type
- ZIB-Report (13)
- Article (1)
- Doctoral Thesis (1)
- Master's Thesis (1)
Language
- English (16)
Is part of the Bibliography
- no (16)
Keywords
- Markov State Models (3)
- Mixed-Integer Programming (3)
- NESS (3)
- Non-reversible Markov Processes (3)
- 2-photon microscopy (1)
- Adaptive Importance Sampling (1)
- Bayesian inverse problems (1)
- Girsanov (1)
- Karhunen–Loève expansion (1)
- Metadynamics (1)
Institute
- Numerical Mathematics (9)
- Mathematical Optimization (7)
- Computational Molecular Design (3)
- Mathematical Optimization Methods (2)
- Modeling and Simulation of Complex Processes (2)
- Visual Data Analysis (2)
- Visual and Data-centric Computing (2)
- Image Analysis in Biology and Materials Science (1)
- Uncertainty Quantification (1)
The paper deals with the workload and busy period for the M/GI/1
system under FCFS discipline, where the customers may become
impatient during their waiting for service with generally
distributed maximal waiting times and also during their service
with generally distributed maximal service times depending
on the time waited for service. This general impatience mechanism,
originally introduced by Kovalenko (1961) and considered
by Daley (1965), too, covers the special cases of impatience
on waiting times as well as impatience on sojourn times,
for which Boxma et al. (2010), (2011) gave new results
and outlined special cases recently. Our unified approach
bases on the vector process of workload and busy time.
Explicit representations for the LSTs of workload and
busy period are given in case of phase-type distributed
impatience.
This article extends the framework of Bayesian inverse problems in infinite-dimensional parameter spaces, as advocated by Stuart (Acta Numer. 19:451–559, 2010) and others, to the case of a heavy-tailed prior measure in the family of stable distributions, such as an infinite-dimensional Cauchy distribution, for which polynomial moments are infinite or undefined. It is shown that analogues of the Karhunen–Loève expansion for square-integrable random variables can be used to sample such measures. Furthermore, under weaker regularity assumptions than those used to date, the Bayesian posterior measure is shown to depend Lipschitz continuously in the Hellinger metric upon perturbations of the misfit function and observed data.
Neurotransmission at chemical synapses relies on the calcium-induced fusion of synaptic vesicles with the presynaptic membrane. The distance to the calcium channels determines the release probability and thereby the postsynaptic signal. Suitable models of the process need to capture both the mean and the variance observed in electrophysiological measurements of the postsynaptic current. In this work, we propose a method to directly compute the exact first- and second-order moments for signals generated by a linear reaction network under convolution with an impulse response function, rendering computationally expensive numerical simulations of the underlying stochastic counting process obsolete. We show that the autocorrelation of the process is central for the calculation of the filtered signal’s second-order moments, and derive a system of PDEs for the cross-correlation functions (including the autocorrelations) of linear reaction networks with time-dependent rates. Finally, we employ our method to efficiently compare different spatial coarse graining approaches for a specific model of synaptic vesicle fusion. Beyond the application to neurotransmission processes, the developed theory can be applied to any linear reaction system that produces a filtered stochastic signal.
This paper proposes a new method for probabilistic analysis of online algorithms. It is based on the notion of stochastic dominance. We develop the method for
the online bin coloring problem introduced by Krumke et al (2008). Using methods for the stochastic
comparison of Markov chains we establish the result that the performance of the online algorithm GreedyFit is stochastically better than the performance of the algorithm OneBin for any number of items processed. This result gives a more realistic
picture than competitive analysis and explains the behavior observed in simulations.
Some connections between importance sampling and enhanced sampling methods in molecular dynamics
(2017)
Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Schütte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling.
Following axon pathfinding, growth cones transition from stochastic filopodial exploration to the formation of a limited number of synapses. How the interplay of filopodia and synapse assembly ensures robust connectivity in the brain has remained a challenging problem. Here, we developed a new 4D analysis method for filopodial dynamics and a data-driven computational model of synapse formation for R7 photoreceptor axons in developing Drosophila brains. Our live data support a 'serial synapse formation' model, where at any time point only a single 'synaptogenic' filopodium suppresses the synaptic competence of other filopodia through competition for synaptic seeding factors. Loss of the synaptic seeding factors Syd-1 and Liprin-α leads to a loss of this suppression, filopodial destabilization and reduced synapse formation, which is sufficient to cause the destabilization of entire axon terminals. Our model provides a filopodial 'winner-takes-all' mechanism that ensures the formation of an appropriate number of synapses.
Despite the success of constraint programming (CP) for scheduling, the much wider penetration of mixed integer programming (MIP) technology into business applications means that many practical scheduling problems are being addressed with MIP, at least as an initial approach. Furthermore, there has been impressive and well-documented improvements in the power of generic MIP solvers over the past decade.
We empirically demonstrate that on an existing set of resource allocation and scheduling problems standard MIP and CP models are now competitive with the state-of-the-art manual decomposition approach. Motivated by this result, we formulate two tightly coupled hybrid models based on constraint integer programming (CIP) and demonstrate that these models,
which embody advances in CP and MIP, are able to out-perform the CP, MIP, and decomposition models. We conclude that both MIP and CIP
are technologies that should be considered along with CP for solving scheduling problems.
We deal with an infinite-server system where the
service speed is governed by a stationary and ergodic
process with countably many states. Applying a random
time transformation such that the service speed
becomes one, the sojourn time of a class of virtual
requests with given required service time is equal
in distribution to an additive functional defined
via a stationary version of the time-changed process.
Thus bounds for the expectation of functions of additive
functionals yield bounds for the expectation
of functions of virtual sojourn times, in particular
bounds for fractional moments and the distribution
function. Interpreting the $GI(n)/GI(n)/\infty$ system or
equivalently the $GI(n)/GI$ system under state-dependent
processor sharing as an infinite-server system with
random states given by the number $n$ of requests
in the system provides results for sojourn times
of virtual requests. In case of $M(n)/GI(n)/\infty$,
the sojourn times of arriving and added requests are
equal in distribution to sojourn times of virtual
requests in modified systems, which yields many results
for the sojourn times of arriving and added requests.
In case of integer moments, the bounds generalize
earlier results for $M/GI(n)/\infty$. In particular,
the mean sojourn times of arriving and added requests
in $M(n)/GI(n)/\infty$ are proportional to the required
service time, generalizing Cohen's famous result
for $M/GI(n)/\infty$.
In this paper, we present a new, optimization-based method to exhibit cyclic behavior in non-reversible stochastic processes. While our method is general, it is strongly motivated by discrete simulations of ordinary differential equations representing non-reversible biological processes, in particular molecular simulations. Here, the discrete time steps of the simulation are often very small compared to the time scale of interest, i.e., of the whole process. In this setting, the detection of a global cyclic behavior of the process becomes difficult because transitions between individual states may appear almost reversible on the small time scale of the simulation. We address this difficulty using a mixed-integer programming model that allows us to compute a cycle of clusters with maximum net flow, i.e., large forward and small backward probability. For a synthetic genetic regulatory network consisting of a ring-oscillator with three genes, we show that this approach can detect the most productive overall cycle, outperforming classical spectral analysis methods. Our method applies to general non-equilibrium steady state systems such as catalytic reactions, for which the objective value computes the effectiveness of the catalyst.
In this paper, we present a new, optimization-based method to exhibit cyclic behavior in non-reversible stochastic processes. While our method is general, it is strongly motivated by discrete simulations of ordinary differential equations representing non-reversible biological processes, in particular molecular simulations. Here, the discrete time steps of the simulation are often very small compared to the time scale of interest, i.e., of the whole process. In this setting, the detection of a global cyclic behavior of the process becomes difficult because transitions between individual states may appear almost reversible on the small time scale of the simulation. We address this difficulty using a mixed-integer programming model that allows us to compute a cycle of clusters with maximum net flow, i.e., large forward and small backward probability. For a synthetic genetic regulatory network consisting of a ring-oscillator with three genes, we show that this approach can detect the most productive overall cycle, outperforming classical spectral analysis methods. Our method applies to general non-equilibrium steady state systems such as catalytic reactions, for which the objective value computes the effectiveness of the catalyst.