We analyze the optimal policy for the sequential selection of an alternating subsequence from a sequence of n independent observations from a continuous distribution F, and we prove a central limit theorem for the number of selections made by that policy. The proof exploits the backward recursion of dynamic programming and assembles a detailed understanding of the associated value functions and selection rules.

We show that the total number of collisions in the exchangeable coalescent process driven by the beta (1, b) measure converges in distribution to a 1-stable law, as the initial number of particles goes to ∞. The stable limit law is also shown for the total branch length of the coalescent tree. These results were known previously for the instance b = 1, which corresponds to the Bolthausen-Sznitman coalescent. The approach we take is based on estimating the quality of a renewal approximation to the coalescent in terms of a suitable Wasserstein distance. Application of the method to beta (a, b)-coalescents with 0 < a < 1 leads to a simplified derivation of the known (2 - a)-stable limit. We furthermore derive asymptotic expansions for the moments of the number of collisions and of the total branch length for the beta (1, b)-coalescent by exploiting the method of sequential approximations.

In this paper we study a reinsurance game between two insurers whose surplus processes are modeled by arithmetic Brownian motions. We assume a minimax criterion in the game. One insurer tries to maximize the probability of absolute dominance while the other tries to minimize it through reinsurance control. Here absolute dominance is defined as the event that liminf of the difference of the surplus levels tends to -∞. Under suitable parameter conditions, the game is solved with the value function and the Nash equilibrium strategy given in explicit form.

Let T be a stopping time associated with a sequence of independent and identically distributed or exchangeable random variables taking values in {0, 1, 2, …, m}, and let ST,i be the stopped sum denoting the number of appearances of outcome 'i' in X1, …, XT, 0 ≤ i ≤ m. In this paper we present results revealing that, if the distribution of T is known, then we can also derive the joint distribution of (T, ST,0, ST,1, …, ST,m). Two applications, which have independent interest, are offered to illustrate the applicability and the usefulness of the main results.

In this paper a method based on a Markov chain Monte Carlo (MCMC) algorithm is proposed to compute the probability of a rare event. The conditional distribution of the underlying process given that the rare event occurs has the probability of the rare event as its normalizing constant. Using the MCMC methodology, a Markov chain is simulated, with the aforementioned conditional distribution as its invariant distribution, and information about the normalizing constant is extracted from its trajectory. The algorithm is described in full generality and applied to the problem of computing the probability that a heavy-tailed random walk exceeds a high threshold. An unbiased estimator of the reciprocal probability is constructed whose normalized variance vanishes asymptotically. The algorithm is extended to random sums and its performance is illustrated numerically and compared to existing importance sampling algorithms.

We consider a two-dimensional reflecting random walk on the nonnegative integer quadrant. This random walk is assumed to be skip free in the direction to the boundary of the quadrant, but may have unbounded jumps in the opposite direction, which are referred to as upward jumps. We are interested in the tail asymptotic behavior of its stationary distribution, provided it exists. Assuming that the upward jump size distributions have light tails, we find the rough tail asymptotics of the marginal stationary distributions in all directions. This generalizes the corresponding results for the skip-free reflecting random walk in Miyazawa (2009). We exemplify these results for a two-node queueing network with exogenous batch arrivals.

The yeast Saccharomyces cerevisiae has emerged as an ideal model system to study the dynamics of prion proteins which are responsible for a number of fatal neurodegenerative diseases in humans. Within an infected cell, prion proteins aggregate in complexes which may increase in size or be fragmented and are transmitted upon cell division. Recent work in yeast suggests that only aggregates below a critical size are transmitted efficiently. We formulate a continuous-time branching process model of a yeast colony under conditions of prion curing. We generalize previous approaches by providing an explicit formula approximating prion loss as influenced by both aggregate growth and size-dependent transmission.

We construct a flow of continuous-time and discrete-state branching processes. Some scaling limit theorems for the flow are proved, which lead to the path-valued branching processes and nonlocal branching superprocesses, over the positive half line, studied in Li (2014).

In this paper we provide the basis for new methods of inference for max-stable processes ξ on general spaces that admit a certain incremental representation, which, in important cases, has a much simpler structure than the max-stable process itself. A corresponding peaks-over-threshold approach will incorporate all single events that are extreme in some sense and will therefore rely on a substantially larger amount of data in comparison to estimation procedures based on block maxima. Conditioning a process η in the max-domain of attraction of ξ on being extremal, several convergence results for the increments of η are proved. In a similar way, the shape functions of mixed moving maxima (M3) processes can be extracted from suitably conditioned single events η. Connecting the two approaches, transformation formulae for processes that admit both an incremental and an M3 representation are identified.

In Chaudhuri and Dasgupta's 2006 paper a certain stochastic model for ‘replicating character strings’ (such as in DNA sequences) was studied. In their model, a random ‘input’ sequence was subjected to random mutations, insertions, and deletions, resulting in a random ‘output’ sequence. In this paper their model will be set up in a slightly different way, in an effort to facilitate further development of the theory for their model. In their 2006 paper, Chaudhuri and Dasgupta showed that, under certain conditions, strict stationarity of the ‘input’ sequence would be preserved by the ‘output’ sequence, and they proved a similar ‘preservation’ result for the property of strong mixing with exponential mixing rate. In our setup, we will in spirit slightly extend their ‘preservation of stationarity’ result, and also prove a ‘preservation’ result for the property of absolute regularity with summable mixing rate.

As the name suggests, the family of general error distributions has been used to model nonnormal errors in a variety of situations. In this article we show that the asymptotic distribution of linearly normalized partial maxima of random observations from the general error distributions is Gumbel when the parameter of these distributions lies in the interval (0, 1). Our result fills a gap in the literature. We also establish the corresponding density convergence, obtain an asymptotic distribution of the partial maxima under power normalization, and state and prove a strong law. We also study the asymptotic behaviour of observations near the partial maxima and the sum of such observations.

The paper considers a statistical concept of causality in continuous time between filtered probability spaces, based on Granger’s definition of causality. This causality concept is connected with the preservation of the martingale representation property when the filtration is getting smaller. We also give conditions, in terms of causality, for every martingale to be a continuous semimartingale, and we consider the equivalence between the concept of causality and the preservation of the martingale representation property under change of measure. In addition, we apply these results to weak solutions of stochastic differential equations. The results can be applied to the economics of securities trading.

There is much interest within the mathematical biology and statistical physics community in converting stochastic agent-based models for random walkers into a partial differential equation description for the average agent density. Here a collection of noninteracting biased random walkers on a one-dimensional lattice is considered. The usual master equation approach requires that two continuum limits, involving three parameters, namely step length, time step and the random walk bias, approach zero in a specific way. We are interested in the case where the two limits are not consistent. New results are obtained using a Fokker–Planck equation and the results are highly dependent on the simulation update schemes. The theoretical results are confirmed with examples. These findings provide insight into the importance of updating schemes to an accurate macroscopic description of stochastic local movement rules in agent-based models when the lattice spacing represents a physical object such as cell diameter.

Consider a one-dimensional diffusion process on the diffusion interval I originated in x0 ∈ I. Let a(t) and b(t) be two continuous functions of t, t > t0, with bounded derivatives, a(t) < b(t), and a(t), b(t) ∈ I, for all t > t0. We study the joint distribution of the two random variables Ta and Tb, the first hitting times of the diffusion process through the two boundaries a(t) and b(t), respectively. We express the joint distribution of Ta and Tb in terms of ℙ(Ta < t, Ta < Tb) and ℙ(Tb < t, Ta > Tb), and we determine a system of integral equations verified by these last probabilities. We propose a numerical algorithm to solve this system and we prove its convergence properties. Examples and modeling motivation for this study are also discussed.

The recently introduced two-parameter infinitely-many-neutral-alleles model extends the celebrated one-parameter version (which is related to Kingman's distribution) to diffusive two-parameter Poisson-Dirichlet frequencies. In this paper we investigate the dynamics driving the species heterogeneity underlying the two-parameter model. First we show that a suitable normalization of the number of species is driven by a critical continuous-state branching process with immigration. Secondly, we provide a finite-dimensional construction of the two-parameter model, obtained by means of a sequence of Feller diffusions of Wright-Fisher flavor which feature finitely many types and inhomogeneous mutation rates. Both results provide insight into the mathematical properties and biological interpretation of the two-parameter model, showing that it is structurally different from the one-parameter case in that the frequency dynamics are driven by state-dependent rather than constant quantities.

We consider an investment problem where observing and trading are only possible at random times. In addition, we introduce drawdown constraints which require that the investor's wealth does not fall under a prior fixed percentage of its running maximum. The financial market consists of a riskless bond and a stock which is driven by a Lévy process. Moreover, a general utility function is assumed. In this setting we solve the investment problem using a related limsup Markov decision process. We show that the value function can be characterized as the unique fixed point of the Bellman equation and verify the existence of an optimal stationary policy. Under some mild assumptions the value function can be approximated by the value function of a contracting Markov decision process. We are able to use Howard's policy improvement algorithm for computing the value function as well as an optimal policy. These results are illustrated in a numerical example.

In this paper, bounded variation control of one-dimensional diffusion processes is considered. We assume that the agent is allowed to control the diffusion only at the jump times of an observable, independent Poisson process. The agent's objective is to maximize the expected present value of the cumulative payoff generated by the controlled diffusion over its lifetime. We propose a relatively weak set of assumptions on the underlying diffusion and the instantaneous payoff structure, under which we solve the problem in closed form. Moreover, we illustrate the main results with an explicit example.

We derive explicit lower and upper bounds for the probability generating functional of a stationary locally stable Gibbs point process, which can be applied to summary statistics such as the F function. For pairwise interaction processes we obtain further estimates for the G and K functions, the intensity, and higher-order correlation functions. The proof of the main result is based on Stein's method for Poisson point process approximation.

We construct a wavelet-based almost-sure uniform approximation of fractional Brownian motion (FBM) (Bt(H))_t∈[0,1] of Hurst index H ∈ (0, 1). Our results show that, by Haar wavelets which merely have one vanishing moment, an almost-sure uniform expansion of FBM for H ∈ (0, 1) can be established. The convergence rate of our approximation is derived. We also describe a parallel algorithm that generates sample paths of an FBM efficiently.

Pairs trading is a trading strategy which is used very frequently in the financial industry. An investment opportunity arises when the spread between two assets, which historically have exhibited autoregressive behavior, deviates from its recent history. In this case, the investor takes a long position in the asset which is expected to outperform going forward and finances this by taking a short position in the other one. If the spread converges, the investor can close both positions to generate a profit. We model the spread between two assets as an Ornstein-Uhlenbeck process and assume a constant opportunity cost. We then study the optimal liquidation strategy for an investor who wants to optimize profit in excess of the opportunity cost. Including this cost is important from an applied perspective, as the performance of any investment is always evaluated relative to the performance of the opportunity set.