Bifurcating Markov chains (BMCs) are Markov chains indexed by a full binary tree representing the evolution of a trait along a population where each individual has two children. We provide a central limit theorem for additive functionals of BMCs under $L^2$ -ergodic conditions with three different regimes. This completes the pointwise approach developed in a previous work. As an application, we study the elementary case of a symmetric bifurcating autoregressive process, which justifies the nontrivial hypothesis considered on the kernel transition of the BMCs. We illustrate in this example the phase transition observed in the fluctuations.

Empirical studies (e.g. Jiang et al. (2015) and Mislove et al. (2007)) show that online social networks have not only in- and out-degree distributions with Pareto-like tails, but also a high proportion of reciprocal edges. A classical directed preferential attachment (PA) model generates in- and out-degree distributions with power-law tails, but the theoretical properties of the reciprocity feature in this model have not yet been studied. We derive asymptotic results on the number of reciprocal edges between two fixed nodes, as well as the proportion of reciprocal edges in the entire PA network. We see that with certain choices of parameters, the proportion of reciprocal edges in a directed PA network is close to 0, which differs from the empirical observation. This points out one potential problem of fitting a classical PA model to a given network dataset with high reciprocity, and indicates that alternative models need to be considered.

In this paper we propose a general framework for modeling an insurance liability cash flow in continuous time, by generalizing the reduced-form framework for credit risk and life insurance. In particular, we assume a nontrivial dependence structure between the reference filtration and the insurance internal filtration. We apply these results for pricing and hedging non-life insurance liabilities in hybrid financial and insurance markets, while taking into account the role of inflation under the benchmarked risk-minimization approach. This framework offers at the same time a general and flexible structure, and an explicit and treatable pricing-hedging formula.

We study the integrated telegraph process $X_t$ under the assumption of general distribution for the random times between consecutive reversals of direction. Specifically, $X_t$ represents the position, at time t, of a particle moving U time units upwards with velocity c and D time units downwards with velocity $-c$ . The latter motions are repeated cyclically, according to independent alternating renewals. Explicit expressions for the probability law of $X_t$ are given in the following cases: (i) (U, D) gamma-distributed; (ii) U exponentially distributed and D gamma-distributed. For certain values of the parameters involved, the probability law of $X_t$ is provided in a closed form. Some expressions for the moment generating function of $X_t$ and its Laplace transform are also obtained. The latter allows us to prove the existence of a Kac-type condition under which the probability density function of the integrated telegraph process, with identically distributed gamma intertimes, converges to that of the standard Brownian motion.

Finally, we consider the square of $X_t$ and disclose its distribution function, specifying the expression for some choices of the distribution of (U, D).

We apply general moment identities for Poisson stochastic integrals with random integrands to the computation of the moments of Markovian growth–collapse processes. This extends existing formulas for mean and variance available in the literature to closed-form moment expressions of all orders. In comparison with other methods based on differential equations, our approach yields explicit summations in terms of the time parameter. We also treat the case of the associated embedded chain, and provide recursive codes in Maple and Mathematica for the computation of moments and cumulants of any order with arbitrary cut-off moment sequences and jump size functions.

We consider residue expansions for survival and density/mass functions of first-passage distributions in finite-state semi-Markov processes (SMPs) in continuous and integer time. Conditions are given which guarantee that the residue expansions for these functions have a dominant exponential/geometric term. The key condition assumes that the relevant states for first passage contain an irreducible class, thus ensuring the same sort of dominant exponential/geometric terms as one gets for phase-type distributions in Markov processes. Essentially, the presence of an irreducible class along with some other conditions ensures that the boundary singularity b for the moment generating function (MGF) of the first-passage-time distribution is a simple pole. In the continuous-time setting we prove that b is a dominant pole, in that the MGF has no other pole on the vertical line $\{\text{Re}(s)=b\}.$ In integer time we prove that b is dominant if all holding-time mass functions for the SMP are aperiodic and non-degenerate. The expansions and pole characterisations address first passage to a single new state or a subset of new states, and first return to the starting state. Numerical examples demonstrate that the residue expansions are considerably more accurate than saddlepoint approximations and can provide a substitute for exact computation above the 75th percentile.

We study an intertemporal consumption and portfolio choice problem under Knightian uncertainty in which agent’s preferences exhibit local intertemporal substitution. We also allow for market frictions in the sense that the pricing functional is nonlinear. We prove existence and uniqueness of the optimal consumption plan, and we derive a set of sufficient first-order conditions for optimality. With the help of a backward equation, we are able to determine the structure of optimal consumption plans. We obtain explicit solutions in a stationary setting in which the financial market has different risk premia for short and long positions.

We consider two-dimensional Lévy processes reflected to stay in the positive quadrant. Our focus is on the non-standard regime when the mean of the free process is negative but the reflection vectors point away from the origin, so that the reflected process escapes to infinity along one of the axes. Under rather general conditions, it is shown that such behaviour is certain and each component can dominate the other with positive probability for any given starting position. Additionally, we establish the corresponding invariance principle providing justification for the use of the reflected Brownian motion as an approximate model. Focusing on the probability that the first component dominates, we derive a kernel equation for the respective Laplace transform in the starting position. This is done for the compound Poisson model with negative exponential jumps and, by means of approximation, for the Brownian model. Both equations are solved via boundary value problem analysis, which also yields the domination probability when starting at the origin. Finally, certain asymptotic analysis and numerical results are presented.

Restless bandits are a class of sequential resource allocation problems concerned with allocating one or more resources among several alternative processes where the evolution of the processes depends on the resources allocated to them. Such models capture the fundamental trade-offs between exploration and exploitation. In 1988, Whittle developed an index heuristic for restless bandit problems which has emerged as a popular solution approach because of its simplicity and strong empirical performance. The Whittle index heuristic is applicable if the model satisfies a technical condition known as indexability. In this paper, we present two general sufficient conditions for indexability and identify simpler-to-verify refinements of these conditions. We then revisit a previously proposed algorithm called the adaptive greedy algorithm which is known to compute the Whittle index for a sub-class of restless bandits. We show that a generalization of the adaptive greedy algorithm computes the Whittle index for all indexable restless bandits. We present an efficient implementation of this algorithm which can compute the Whittle index of a restless bandit with K states in $\mathcal{O}\!\left(K^3\right)$ computations. Finally, we present a detailed numerical study which affirms the strong performance of the Whittle index heuristic.

The hyperbolic random geometric graph was introduced by Krioukov et al. (Phys. Rev. E82, 2010). Among many equivalent models for the hyperbolic space, we study the d-dimensional Poincaré ball ( $d\ge 2$ ), with a general connectivity radius. While many phase transitions are known for the expectation asymptotics of certain subgraph counts, very little is known about the second-order results. Two of the distinguishing characteristics of geometric graphs on the hyperbolic space are the presence of tree-like hierarchical structures and the power-law behaviour of the degree distribution. We aim to reveal such characteristics in detail by investigating the behaviour of sub-tree counts. We show multiple phase transitions for expectation and variance in the resulting hyperbolic geometric graph. In particular, the expectation and variance of the sub-tree counts exhibit an intricate dependence on the degree sequence of the tree under consideration. Additionally, unlike the thermodynamic regime of the Euclidean random geometric graph, the expectation and variance may exhibit different growth rates, which is indicative of power-law behaviour. Finally, we also prove a normal approximation for sub-tree counts using the Malliavin–Stein method of Last et al. (Prob. Theory Relat. Fields165, 2016), along with the Palm calculus for Poisson point processes.

The study of multivariate extremes is dominated by multivariate regular variation, although it is well known that this approach does not provide adequate distinction between random vectors whose components are not always simultaneously large. Various alternative dependence measures and representations have been proposed, with the most well-known being hidden regular variation and the conditional extreme value model. These varying depictions of extremal dependence arise through consideration of different parts of the multivariate domain, and particularly through exploring what happens when extremes of one variable may grow at different rates from other variables. Thus far, these alternative representations have come from distinct sources, and links between them are limited. In this work we elucidate many of the relevant connections through a geometrical approach. In particular, the shape of the limit set of scaled sample clouds in light-tailed margins is shown to provide a description of several different extremal dependence representations.

Let $X_t^\sharp$ be a multivariate process of the form $X_t =Y_t - Z_t$ , $X_0=x$ , killed at some terminal time T, where $Y_t$ is a Markov process having only jumps of length smaller than $\delta$ , and $Z_t$ is a compound Poisson process with jumps of length bigger than $\delta$ , for some fixed $\delta>0$ . Under the assumptions that the summands in $Z_t$ are subexponential, we investigate the asymptotic behaviour of the potential function $u(x)= \mathbb{E}^x \int_0^\infty \ell\big(X_s^\sharp\big)ds$ . The case of heavy-tailed entries in $Z_t$ corresponds to the case of ‘big claims’ in insurance models and is of practical interest. The main approach is based on the fact that u(x) satisfies a certain renewal equation.

The practically important classes of equal-input and of monotone Markov matrices are revisited, with special focus on embeddability, infinite divisibility, and mutual relations. Several uniqueness results for the classic Markov embedding problem are obtained in the process. To achieve our results, we need to employ various algebraic and geometric tools, including commutativity, permutation invariance, and convexity. Of particular relevance in several demarcation results are Markov matrices that are idempotents.

A continuous-state branching process with immigration having branching mechanism $\Psi$ and immigration mechanism $\Phi$ , a CBI $(\Psi,\Phi)$ process for short, may have either of two different asymptotic regimes, depending on whether $\int_{0}\frac{\Phi(u)}{|\Psi(u)|}\textrm{d} u<\infty$ or $\int_{0}\frac{\Phi(u)}{|\Psi(u)|}\textrm{d} u=\infty$ . When $\int_{0}\frac{\Phi(u)}{|\Psi(u)|}\textrm{d} u<\infty$ , the CBI process has either a limit distribution or a growth rate dictated by the branching dynamics. When $\scriptstyle\int_{0}\tfrac{\Phi(u)}{|\Psi(u)|}\textrm{d} u=\infty$ , immigration overwhelms branching dynamics. Asymptotics in the latter case are studied via a nonlinear time-dependent renormalization in law. Three regimes of weak convergence are exhibited. Processes with critical branching mechanisms subject to a regular variation assumption are studied. This article proves and extends results stated by M. Pinsky in ‘Limit theorems for continuous state branching processes with immigration’ (Bull. Amer. Math. Soc.78, 1972).

We correct typographical errors in our original paper.

The generalized perturbative approach is an all-purpose variant of Stein’s method used to obtain rates of normal approximation. Originally developed for functions of independent random variables, this method is here extended to functions of the realization of a hidden Markov model. In this dependent setting, rates of convergence are provided in some applications, leading, in each instance, to an extra log-factor vis-à-vis the rate in the independent case.

We solve non-Markovian optimal switching problems in discrete time on an infinite horizon, when the decision-maker is risk-aware and the filtration is general, and establish existence and uniqueness of solutions for the associated reflected backward stochastic difference equations. An example application to hydropower planning is provided.

In this paper, we construct operator fractional Lévy motion (ofLm), a broad class of infinitely divisible stochastic processes that are covariance operator self-similar and have wide-sense stationary increments. The ofLm class generalizes the univariate fractional Lévy motion as well as the multivariate operator fractional Brownian motion (ofBm). OfLm can be divided into two types, namely, moving average (maofLm) and real harmonizable (rhofLm), both of which share the covariance structure of ofBm under assumptions. We show that maofLm and rhofLm admit stochastic integral representations in the time and Fourier domains, and establish their distinct small- and large-scale limiting behavior. We also characterize time-reversibility for ofLm through parametric conditions related to its Lévy measure. In particular, we show that, under non-Gaussianity, the parametric conditions for time-reversibility are generally more restrictive than those for the Gaussian case (ofBm).