We introduce and study a game-theoretic model to understand the spread of an epidemic in a homogeneous population. A discrete-time stochastic process is considered where, in each epoch, first, a randomly chosen agent updates their action trying to maximize a proposed utility function, and then agents who have viral exposures beyond their immunity get infected. Our main results discuss asymptotic limiting distributions of both the cardinality of the subset of infected agents and the action profile, considered under various values of two parameters (initial action and immunity profile). We also show that the theoretical distributions are almost always achieved in the first few epochs.

Let $\mathbb{P}_\kappa(n)$ be the probability that n points $z_1,\ldots,z_n$ picked uniformly and independently in $\mathfrak{C}_\kappa$, a regular $\kappa$-gon with area 1, are in convex position, that is, form the vertex set of a convex polygon. In this paper, we compute $\mathbb{P}_\kappa(n)$ up to asymptotic equivalence, as $n\to+\infty$, for all $\kappa\geq 3$, which improves on a famous result of Bárány (Ann. Prob. 27, 1999). The second purpose of this paper is to establish a limit theorem which describes the fluctuations around the limit shape of an n-tuple of points in convex position when $n\to+\infty$. Finally, we give an asymptotically exact algorithm for the random generation of $z_1,\ldots,z_n$, conditioned to be in convex position in $\mathfrak{C}_\kappa$.

The study of many population growth models is complicated by only partial observation of the underlying stochastic process driving the model. For example, in an epidemic outbreak we might know when individuals show symptoms to a disease and are removed, but not when individuals are infected. Motivated by the above example and the long-established approximation of epidemic processes by branching processes, we explore the number of individuals alive in a time-inhomogeneous branching process with a general phase-type lifetime distribution given only (partial) information on the times of deaths of individuals. Deaths are detected independently with a detection probability that can vary with time and type. We show that the number of individuals alive immediately after the kth detected death can be expressed as the mixture of random variables each of which consists of the sum of k independent zero-modified geometric distributions. Furthermore, in the case of an Erlang lifetime distribution, we derive an easy-to-compute mixture of negative binomial distributions as an approximation of the number of individuals alive immediately after the kth detected death.

This paper investigates structural changes in the parameters of first-order autoregressive (AR) models by analyzing the edge eigenvalues of the precision matrices. Specifically, edge eigenvalues in the precision matrix are observed if and only if there is a structural change in the AR coefficients. We show that these edge eigenvalues correspond to the zeros of a determinantal equation. Additionally, we propose a consistent estimator for detecting outliers within the panel time series framework, supported by numerical experiments.

The preferential attachment model is a natural and popular random graph model for a growing network that contains very well-connected ‘hubs’. We study the higher-order connectivity of such a network by investigating the topological properties of its clique complex. We concentrate on the Betti numbers, a sequence of topological invariants of the complex related to the numbers of holes (equivalently, repeated connections) of different dimensions. We prove that the expected Betti numbers grow sublinearly fast, with the trivial exceptions of those at dimensions 0 and 1. Our result also shows that preferential attachment graphs undergo infinitely many phase transitions within the parameter regime where the limiting degree distribution has an infinite variance. Regarding higher-order connectivity, our result shows that preferential attachment favors higher-order connectivity. We illustrate our theoretical results with simulations.

This paper characterizes irreducible phase-type representations for exponential distributions. Bean and Green (2000) gave a set of necessary and sufficient conditions for a phase-type distribution with an irreducible generator matrix to be exponential. We extend these conditions to irreducible representations, and we thus give a characterization of all irreducible phase-type representations for exponential distributions. We consider the results in relation to time-reversal of phase-type distributions, PH-simplicity, and the algebraic degree of a phase-type distribution, and we give applications of the results. In particular we give the conditions under which a Coxian distribution becomes exponential, and we construct bivariate exponential distributions. Finally, we translate the main findings to the discrete case of geometric distributions.

Let $(X_t)_{t \geqslant 0}$ be the solution of the stochastic differential equation

\[\mathrm{d} X_t = b(X_t) \,\mathrm{d} t+A\,\mathrm{d} Z_t, \quad X_{0}=x,\]

$b\,:\, \mathbb{R}^d \rightarrow \mathbb{R}^d$

$A \in \mathbb{R}^{d \times d}$

$(Z_t)_{t\geqslant 0}$

$\alpha$

$\alpha \in (1,2)$

$x\in\mathbb{R}^{d}$

$\Gamma = (\gamma_n)_{n\in \mathbb{N}}$

$(X_t)_{t \geqslant 0}$

$\alpha$

$\gamma^{\frac{1}{\alpha}}_n$

$\gamma^{1+\frac {1}{\alpha}-\frac{1}{\kappa}}_n$

$\kappa \in [1,\alpha)$

$\gamma^{\frac{2-\alpha}{\alpha}}_n$

$\alpha$

We study continuous-time Markov chains on the nonnegative integers under mild regularity conditions (in particular, the set of jump vectors is finite and both forward and backward jumps are possible). Based on the so-called flux balance equation, we derive an iterative formula for calculating stationary measures. Specifically, a stationary measure $\pi(x)$ evaluated at $x\in\mathbb{N}_0$ is represented as a linear combination of a few generating terms, similarly to the characterization of a stationary measure of a birth–death process, where there is only one generating term, $\pi(0)$. The coefficients of the linear combination are recursively determined in terms of the transition rates of the Markov chain. For the class of Markov chains we consider, there is always at least one stationary measure (up to a scaling constant). We give various results pertaining to uniqueness and nonuniqueness of stationary measures, and show that the dimension of the linear space of signed invariant measures is at most the number of generating terms. A minimization problem is constructed in order to compute stationary measures numerically. Moreover, a heuristic linear approximation scheme is suggested for the same purpose by first approximating the generating terms. The correctness of the linear approximation scheme is justified in some special cases. Furthermore, a decomposition of the state space into different types of states (open and closed irreducible classes, and trapping, escaping and neutral states) is presented. Applications to stochastic reaction networks are well illustrated.

We use the framework of multivariate regular variation to analyse the extremal behaviour of preferential attachment models. To this end, we follow a directed linear preferential attachment model for a random, heavy-tailed number of steps in time and treat the incoming edge count of all existing nodes as a random vector of random length. By combining martingale properties, moment bounds and a Breiman type theorem we show that the resulting quantity is multivariate regularly varying, both as a vector of fixed length formed by the edge counts of a finite number of oldest nodes, and also as a vector of random length viewed in sequence space. A Pólya urn representation allows us to explicitly describe the extremal dependence between the degrees with the help of Dirichlet distributions. As a by-product of our analysis we establish new results for almost sure convergence of the edge counts in sequence space as the number of nodes goes to infinity.

For a continuous-time phase-type (PH) distribution, starting with its Laplace–Stieltjes transform, we obtain a necessary and sufficient condition for its minimal PH representation to have the same order as its algebraic degree. To facilitate finding this minimal representation, we transform this condition equivalently into a non-convex optimization problem, which can be effectively addressed using an alternating minimization algorithm. The algorithm convergence is also proved. Moreover, the method we develop for the continuous-time PH distributions can be used directly for the discrete-time PH distributions after establishing an equivalence between the minimal representation problems for continuous-time and discrete-time PH distributions.

Motivated by recent developments of quasi-stationary Monte Carlo methods, we investigate the stability of quasi-stationary distributions of killed Markov processes under perturbations of the generator. We first consider a general bounded self-adjoint perturbation operator, and then study a particular unbounded perturbation corresponding to truncation of the killing rate. In both scenarios, we quantify the difference between eigenfunctions of the smallest eigenvalue of the perturbed and unperturbed generators in a Hilbert space norm. As a consequence, $\mathcal{L}^1$-norm estimates of the difference of the resulting quasi-stationary distributions in terms of the perturbation are provided.