It was recently proven that the correlation function of the stationary version of a reflected Lévy process is nonnegative, nonincreasing, and convex. In another branch of the literature it was established that the mean value of the reflected process starting from zero is nondecreasing and concave. In the present paper it is shown, by putting them in a common framework, that these results extend to substantially more general settings. Indeed, instead of reflected Lévy processes, we consider a class of more general stochastically monotone Markov processes. In this setup we show monotonicity results associated with a supermodular function of two coordinates of our Markov process, from which the above-mentioned monotonicity and convexity/concavity results directly follow, but now for the class of Markov processes considered rather than just reflected Lévy processes. In addition, various results for the transient case (when the Markov process is not in stationarity) are provided. The conditions imposed are natural, in that they are satisfied by various frequently used Markovian models, as illustrated by a series of examples.

This paper is concerned with the optimal number of redundant allocation to n-component coherent systems consisting of heterogeneous dependent components. We assume that the system is built up of L groups of different components, $L\geq 1$, where there are $n_i$ components in group i, and $\sum_{i=1}^{L}n_i=n$. The problem of interest is to allocate $v_i$ active redundant components to each component of type i, $i=1,\dots, L$. To get the optimal values of $v_i$ we propose two cost-based criteria. One of them is introduced based on the costs of renewing the failed components and the costs of refreshing the alive ones at the system failure time. The other criterion is proposed based on the costs of replacing the system at its failure time or at a predetermined time $\tau$, whichever occurs first. The expressions for the proposed functions are derived using the mixture representation of the system reliability function based on the notion of survival signature. We assume that a given copula function models the dependency structure between the components. In the particular case that the system is a series-parallel structure, we provide the formulas for the proposed cost-based functions. The results are discussed numerically for some specific coherent systems.

We explain some key challenges when dealing with a single- or multi-objective optimization problem in practice. To overcome these challenges, we present a mathematical program that optimizes the Nash social welfare function. We refer to this mathematical program as the Nash social welfare program (NSWP). An interesting property of the NSWP is that it can be constructed for any single- or multi-objective optimization problem. We show that solving the NSWP could result in more desirable solutions in practice than its single- or multi-objective counterpart. We also discuss several promising approaches that could be employed to solve the NSWP in practice.

Motivated by applications to wireless networks, cloud computing, data centers, etc., stochastic processing networks have been studied in the literature under various asymptotic regimes. In the heavy traffic regime, the steady-state mean queue length is proved to be $\Theta({1}/{\epsilon})$, where $\epsilon$ is the heavy traffic parameter (which goes to zero in the limit). The focus of this paper is on obtaining queue length bounds on pre-limit systems, thus establishing the rate of convergence to heavy traffic. For the generalized switch, operating under the MaxWeight algorithm, we show that the mean queue length is within $\textrm{O}({\log}({1}/{\epsilon}))$ of its heavy traffic limit. This result holds regardless of the complete resource pooling (CRP) condition being satisfied. Furthermore, when the CRP condition is satisfied, we show that the mean queue length under the MaxWeight algorithm is within $\textrm{O}({\log}({1}/{\epsilon}))$ of the optimal scheduling policy. Finally, we obtain similar results for the rate of convergence to heavy traffic of the total queue length in load balancing systems operating under the ‘join the shortest queue’ routeing algorithm.

We introduce an approach and a software tool for solving coupled energy networks composed of gas and electric power networks. Those networks are coupled to stochastic fluctuations to address possibly fluctuating demand due to fluctuating demands and supplies. Through computational results, the presented approach is tested on networks of realistic size.

A continuum of stochastic dominance rules, also referred to as fractional stochastic dominance (SD), was introduced by Müller, Scarsini, Tsetlin, and Winkler (2017) to cover preferences from first- to second-order SD. Fractional SD can be used to explain many individual behaviors in economics. In this paper we introduce the concept of fractional pure SD, a special case of fractional SD. We investigate further properties of fractional SD, for example the generating processes of fractional pure SD via $\gamma$-transfers of probability, Yaari’s dual characterization by utilizing the special class of distortion functions, the separation theorem in terms of first-order SD and fractional pure SD, Strassen’s representation, and bivariate characterization. We also establish several closure properties of fractional SD under quantile truncation, under comonotonic sums, and under distortion, as well as its equivalence characterization. Examples of distributions ordered in the sense of fractional SD are provided.

Signal-to-interference-plus-noise ratio (SINR) percolation is an infinite-range dependent variant of continuum percolation modeling connections in a telecommunication network. Unlike in earlier works, in the present paper the transmitted signal powers of the devices of the network are assumed random, independent and identically distributed, and possibly unbounded. Additionally, we assume that the devices form a stationary Cox point process, i.e., a Poisson point process with stationary random intensity measure, in two or more dimensions. We present the following main results. First, under suitable moment conditions on the signal powers and the intensity measure, there is percolation in the SINR graph given that the device density is high and interferences are sufficiently reduced, but not vanishing. Second, if the interference cancellation factor $\gamma$ and the SINR threshold $\tau$ satisfy $\gamma \geq 1/(2\tau)$, then there is no percolation for any intensity parameter. Third, in the case of a Poisson point process with constant powers, for any intensity parameter that is supercritical for the underlying Gilbert graph, the SINR graph also percolates with some small but positive interference cancellation factor.

Log-concavity of a joint survival function is proposed as a model for bivariate increasing failure rate (BIFR) distributions. Its connections with or distinctness from other notions of BIFR are discussed. A necessary and sufficient condition for a bivariate survival function to be log-concave (BIFR-LCC) is given that elucidates the impact of dependence between lifetimes on ageing. Illustrative examples are provided to explain BIFR-LCC for both positive and negative dependence.

We investigate vertex levels of containment in a random hypergraph grown in the spirit of a recursive tree. We consider a local profile tracking the evolution of the containment of a particular vertex over time, and a global profile concerned with counts of the number of vertices of a particular containment level.

For the local containment profile, we obtain the exact mean, variance, and probability distribution in terms of standard combinatorial quantities such as generalized harmonic numbers and Stirling numbers of the first kind. Asymptotically, we observe phases: the early vertices have an asymptotically normal distribution, intermediate vertices have a Poisson distribution, and late vertices have a degenerate distribution.

As for the global containment profile, we establish an asymptotically normal distribution for the number of vertices at the smallest containment level as well as their covariances with the number of vertices at the second smallest containment level and the variances of these numbers. The orders in the variance–covariance matrix establish concentration laws.

In this paper we study the allocation problem of relevations in coherent systems. The optimal allocation strategies are obtained by implementing stochastic comparisons of different policies according to the usual stochastic order and the hazard rate order. As special cases of relevations, the load-sharing and minimal repair policies are further investigated. Sufficient (and necessary) conditions are established for various stochastic orderings. Numerical examples are also presented as illustrations.

Random intersection graphs model networks with communities, assuming an underlying bipartite structure of communities and individuals, where these communities may overlap. We generalize the model, allowing for arbitrary community structures within the communities. In our new model, communities may overlap, and they have their own internal structure described by arbitrary finite community graphs. Our model turns out to be tractable. We analyze the overlapping structure of the communities, show local weak convergence (including convergence of subgraph counts), and derive the asymptotic degree distribution and the local clustering coefficient.

In a multitype branching process, it is assumed that immigrants arrive according to a non-homogeneous Poisson or a contagious Poisson process (both processes are formulated as a non-homogeneous birth process with an appropriate choice of transition intensities). We show that the normalized numbers of objects of the various types alive at time t for supercritical, critical, and subcritical cases jointly converge in distribution under those two different arrival processes. Furthermore, we provide some transient expectation results when there are only two types of particles.

We study shot noise processes with cluster arrivals, in which entities in each cluster may experience random delays (possibly correlated), and noises within each cluster may be correlated. We prove functional limit theorems for the process in the large-intensity asymptotic regime, where the arrival rate gets large while the shot shape function, cluster sizes, delays, and noises are unscaled. In the functional central limit theorem, the limit process is a continuous Gaussian process (assuming the arrival process satisfies a functional central limit theorem with a Brownian motion limit). We discuss the impact of the dependence among the random delays and among the noises within each cluster using several examples of dependent structures. We also study infinite-server queues with cluster/batch arrivals where customers in each batch may experience random delays before receiving service, with similar dependence structures.

Copula is one of the widely used techniques to describe the dependency structure between components of a system. Among all existing copulas, the family of Archimedean copulas is the popular one due to its wide range of capturing the dependency structures. In this paper, we consider the systems that are formed by dependent and identically distributed components, where the dependency structures are described by Archimedean copulas. We study some stochastic comparisons results for series, parallel, and general $r$-out-of-$n$ systems. Furthermore, we investigate whether a system of used components performs better than a used system with respect to different stochastic orders. Furthermore, some aging properties of these systems have been studied. Finally, some numerical examples are given to illustrate the proposed results.

We study propagation of avalanches in a certain excitable network. The model is a particular case of the one introduced by Larremore et al. (Phys. Rev. E, 2012) and is mathematically equivalent to an endemic variation of the Reed–Frost epidemic model introduced by Longini (Math. Biosci., 1980). Two types of heuristic approximation are frequently used for models of this type in applications: a branching process for avalanches of a small size at the beginning of the process and a deterministic dynamical system once the avalanche spreads to a significant fraction of a large network. In this paper we prove several results concerning the exact relation between the avalanche model and these limits, including rates of convergence and rigorous bounds for common characteristics of the model.

In this paper we consider a new generalized finite mixture model formed by dependent and identically distributed (d.i.d.) components. We then establish results for the comparisons of lifetimes of two such generalized finite mixture models in two different cases: (i) when the two mixture models are formed from two random vectors $\textbf{X}$ and $\textbf{Y}$ but with the same weights, and (ii) when the two mixture models are formed with the same random vectors but with different weights. Because the lifetimes of k-out-of-n systems and coherent systems are special cases of the mixture model considered, we used the established results to compare the lifetimes of k-out-of-n systems and coherent systems with respect to the reversed hazard rate and hazard rate orderings.

Let $\{Y_{1},\ldots ,Y_{n}\}$ be a collection of interdependent nonnegative random variables, with $Y_{i}$ having an exponentiated location-scale model with location parameter $\mu _i$, scale parameter $\delta _i$ and shape (skewness) parameter $\beta _i$, for $i\in \mathbb {I}_{n}=\{1,\ldots ,n\}$. Furthermore, let $\{L_1^{*},\ldots ,L_n^{*}\}$ be a set of independent Bernoulli random variables, independently of $Y_{i}$'s, with $E(L_{i}^{*})=p_{i}^{*}$, for $i\in \mathbb {I}_{n}.$ Under this setup, the portfolio of risks is the collection $\{T_{1}^{*}=L_{1}^{*}Y_{1},\ldots ,T_{n}^{*}=L_{n}^{*}Y_{n}\}$, wherein $T_{i}^{*}=L_{i}^{*}Y_{i}$ represents the $i$th claim amount. This article then presents several sufficient conditions, under which the smallest claim amounts are compared in terms of the usual stochastic and hazard rate orders. The comparison results are obtained when the dependence structure among the claim severities are modeled by (i) an Archimedean survival copula and (ii) a general survival copula. Several examples are also presented to illustrate the established results.

In this work, we study a new model for continuum line-of-sight percolation in a random environment driven by the Poisson–Voronoi tessellation in the d-dimensional Euclidean space. The edges (one-dimensional facets, or simply 1-facets) of this tessellation are the support of a Cox point process, while the vertices (zero-dimensional facets or simply 0-facets) are the support of a Bernoulli point process. Taking the superposition Z of these two processes, two points of Z are linked by an edge if and only if they are sufficiently close and located on the same edge (1-facet) of the supporting tessellation. We study the percolation of the random graph arising from this construction and prove that a 0–1 law, a subcritical phase, and a supercritical phase exist under general assumptions. Our proofs are based on a coarse-graining argument with some notion of stabilization and asymptotic essential connectedness to investigate continuum percolation for Cox point processes. We also give numerical estimates of the critical parameters of the model in the planar case, where our model is intended to represent telecommunications networks in a random environment with obstructive conditions for signal propagation.

We study the transient and limiting behavior of a queue with a Pólya arrival process. The Pólya process is interesting because it exhibits path-dependent behavior, e.g. it satisfies a non-ergodic law of large numbers: the average number of arrivals over time [0, t] converges almost surely to a nondegenerate limit as $t \rightarrow \infty$. We establish a heavy-traffic diffusion limit for the $\sum_{i=1}^{n} P_i/GI/1$ queue, with arrivals occurring exogenously according to the superposition of n independent and identically distributed Pólya point processes. That limit yields a tractable approximation for the transient queue-length distribution, because the limiting net input process is a Gaussian Markov process with stationary increments. We also provide insight into the long-run performance of queues with path-dependent arrival processes. We show how Little’s law can be stated in this context, and we provide conditions under which there is stability for a queue with a Pólya arrival process.

We study ergodic properties of a class of Markov-modulated general birth–death processes under fast regime switching. The first set of results concerns the ergodic properties of the properly scaled joint Markov process with a parameter that is taken to be large. Under very weak hypotheses, we show that if the averaged process is exponentially ergodic for large values of the parameter, then the same applies to the original joint Markov process. The second set of results concerns steady-state diffusion approximations, under the assumption that the ‘averaged’ fluid limit exists. Here, we establish convergence rates for the moments of the approximating diffusion process to those of the Markov-modulated birth–death process. This is accomplished by comparing the generator of the approximating diffusion and that of the joint Markov process. We also provide several examples which demonstrate how the theory can be applied.