Abstract We solve the non-discounted, finite-horizon optimal stopping problem of a Gauss–Markov bridge by using a time-space transformation approach. The associated optimal stopping boundary is proved to be Lipschitz continuous on any closed interval that excludes the horizon, and it is characterized by the unique solution of an integral equation. A Picard iteration algorithm is discussed and implemented to exemplify the numerical computation and geometry of the optimal stopping boundary for some illustrative cases.

Abstract We present a closed-form solution to a discounted optimal stopping zero-sum game in a model based on a generalised geometric Brownian motion with coefficients depending on its running maximum and minimum processes. The optimal stopping times forming a Nash equilibrium are shown to be the first times at which the original process hits certain boundaries depending on the running values of the associated maximum and minimum processes. The proof is based on the reduction of the original game to the equivalent coupled free-boundary problem and the solution of the latter problem by means of the smooth-fit and normal-reflection conditions. We show that the optimal stopping boundaries are partially determined as either unique solutions to the appropriate system of arithmetic equations or unique solutions to the appropriate first-order nonlinear ordinary differential equations. The results obtained are related to the valuation of the perpetual lookback game options with floating strikes in the appropriate diffusion-type extension of the Black–Merton–Scholes model.

Abstract We derive some key extremal features for stationary kth-order Markov chains that can be used to understand how the process moves between an extreme state and the body of the process. The chains are studied given that there is an exceedance of a threshold, as the threshold tends to the upper endpoint of the distribution. Unlike previous studies with $k>1$ , we consider processes where standard limit theory describes each extreme event as a single observation without any information about the transition to and from the body of the distribution. Our work uses different asymptotic theory which results in non-degenerate limit laws for such processes. We study the extremal properties of the initial distribution and the transition probability kernel of the Markov chain under weak assumptions for broad classes of extremal dependence structures that cover both asymptotically dependent and asymptotically independent Markov chains. For chains with $k>1$ , the transition of the chain away from the exceedance involves novel functions of the k previous states, in comparison to just the single value, when $k=1$ . This leads to an increase in the complexity of determining the form of this class of functions, their properties, and the method of their derivation in applications. We find that it is possible to derive an affine normalization, dependent on the threshold excess, such that non-degenerate limiting behaviour of the process, in the neighbourhood of the threshold excess, is assured for all lags. We find that these normalization functions have an attractive structure that has parallels to the Yule–Walker equations. Furthermore, the limiting process is always linear in the innovations. We illustrate the results with the study of kth-order stationary Markov chains with exponential margins based on widely studied families of copula dependence structures.

Abstract We analyze the process M(t) representing the maximum of the one-dimensional telegraph process X(t) with exponentially distributed upward random times and generally distributed downward random times. The evolution of M(t) is governed by an alternating renewal of two phases: a rising phase R and a constant phase C. During a rising phase, X(t) moves upward, whereas, during a constant phase, it moves upward and downward, continuing to move until it attains the maximal level previously reached. Under some choices of the distribution of the downward times, we are able to determine the distribution of C, which allows us to obtain some bounds for the survival function of M(t). In the particular case of exponential downward random times, we derive an explicit expression for the survival function of M(t). Finally, the moments of the first passage time $\Theta_w$ of the process X(t) through a fixed boundary $w>0$ are analyzed.

Abstract In this paper, we explore the applications of Tail Variance (TV) as a measure of tail riskiness and the confidence level of using Tail Conditional Expectation (TCE)-based risk capital. While TCE measures the expected loss of a risk that exceeds a certain threshold, TV measures the variability of risk along its tails. We first derive analytical formulas of TV and TCE for a large variety of probability distributions. These formulas are useful instruments for relevant research works on tail risk measures. We then propose a distribution-free approach utilizing TV to estimate the lower bounds of the confidence level of using TCE-based risk capital. In doing so, we introduce sharpened conditional probability inequalities, which halve the bounds of conventional Markov and Cantelli inequalities. Such an approach is easy to implement. We further investigate the characterization of tail risks by TV through an exploration of TV’s asymptotics. A distribution-free limit formula is derived for the asymptotics of TV. To further investigate the asymptotic properties, we consider two broad distribution families defined on tails, namely, the polynomial-tailed distributions and the exponential-tailed distributions. The two distribution families are found to exhibit an asymptotic equivalence between TV and the reciprocal square of the hazard rate. We also establish asymptotic relationships between TCE and VaR for the two families. Our asymptotic analysis contributes to the existing research by unifying the asymptotic expressions and the convergence rate of TV for Student-t distributions, exponential distributions, and normal distributions, which complements the discussion on the convergence rate of univariate cases in [28]. To show the usefulness of our results, we present two case studies based on real data from the industry. We first show how to use conditional inequalities to assess the confidence of using TCE-based risk capital for different types of insurance businesses. Then, for financial data, we provide alternative evidence for the relationship between the data frequency and the tail categorization by the asymptotics of TV.

Abstract Continuous-time Markov chains are frequently used to model the stochastic dynamics of (bio)chemical reaction networks. However, except in very special cases, they cannot be analyzed exactly. Additionally, simulation can be computationally intensive. An approach to address these challenges is to consider a more tractable diffusion approximation. Leite and Williams (Ann. Appl. Prob.29, 2019) proposed a reflected diffusion as an approximation for (bio)chemical reaction networks, which they called the constrained Langevin approximation (CLA) as it extends the usual Langevin approximation beyond the first time some chemical species becomes zero in number. Further explanation and examples of the CLA can be found in Anderson et al. (SIAM Multiscale Modeling Simul.17, 2019).In this paper, we extend the approximation of Leite and Williams to (nearly) density-dependent Markov chains, as a first step to obtaining error estimates for the CLA when the diffusion state space is one-dimensional, and we provide a bound for the error in a strong approximation. We discuss some applications for chemical reaction networks and epidemic models, and illustrate these with examples. Our method of proof is designed to generalize to higher dimensions, provided there is a Lipschitz Skorokhod map defining the reflected diffusion process. The existence of such a Lipschitz map is an open problem in dimensions more than one.

Abstract We show joint convergence of the Łukasiewicz path and height process for slightly supercritical Galton–Watson forests. This shows that the height processes for supercritical continuous-state branching processes as constructed by Lambert (2002) are the limit under rescaling of their discrete counterparts. Unlike for (sub-)critical Galton–Watson forests, the height process does not encode the entire metric structure of a supercritical Galton–Watson forest. We demonstrate that this result is nonetheless useful, by applying it to the configuration model with an independent and identically distributed power-law degree sequence in the critical window, of which we obtain the metric space scaling limit in the product Gromov–Hausdorff–Prokhorov topology, which is of independent interest.

Abstract Inaccuracy and information measures based on cumulative residual entropy are quite useful and have received considerable attention in many fields, such as statistics, probability, and reliability theory. In particular, many authors have studied cumulative residual inaccuracy between coherent systems based on system lifetimes. In a previous paper (Bueno and Balakrishnan,Prob. Eng. Inf. Sci.36, 2022), we discussed a cumulative residual inaccuracy measure for coherent systems at component level, that is, based on the common, stochastically dependent component lifetimes observed under a non-homogeneous Poisson process. In this paper, using a point process martingale approach, we extend this concept to a cumulative residual inaccuracy measure between non-explosive point processes and then specialize the results to Markov occurrence times. If the processes satisfy the proportional risk hazard process property, then the measure determines the Markov chain uniquely. Several examples are presented, including birth-and-death processes and pure birth process, and then the results are applied to coherent systems at component level subject to Markov failure and repair processes.

Abstract We consider Poisson hail models and characterize up to boundaries the collection of critical moments which guarantee stability. In particular, we treat the case of infinite speed of propagation.

Abstract The Minkowski functionals, including the Euler characteristic statistics, are standard tools for morphological analysis in cosmology. Motivated by cosmic research, we examine the Minkowski functional of the excursion set for an isotropic central limit random field, whose k-point correlation functions (kth-order cumulants) have the same structure as that assumed in cosmic research. Using 3- and 4-point correlation functions, we derive the asymptotic expansions of the Euler characteristic density, which is the building block of the Minkowski functional. The resulting formula reveals the types of non-Gaussianity that cannot be captured by the Minkowski functionals. As an example, we consider an isotropic chi-squared random field and confirm that the asymptotic expansion accurately approximates the true Euler characteristic density.