Showing 1–29 of 29 results for author: Olvera-Cravioto, M

Opinion dynamics on non-sparse networks with community structure

Authors: Panagiotis Andreou, Mariana Olvera-Cravioto

Abstract: We study the evolution of opinions on a directed network with community structure. Individuals update their opinions synchronously based on a weighted average of their neighbors' opinions, their own previous opinions, and external media signals. Our model is akin to the popular Friedkin-Johnsen model, and is able to incorporate factors such as stubbornness, confirmation bias, selective exposure, a… ▽ More We study the evolution of opinions on a directed network with community structure. Individuals update their opinions synchronously based on a weighted average of their neighbors' opinions, their own previous opinions, and external media signals. Our model is akin to the popular Friedkin-Johnsen model, and is able to incorporate factors such as stubbornness, confirmation bias, selective exposure, and multiple topics, which are believed to play an important role in the formation of opinions. Our main result shows that, in the large graph limit, the opinion process concentrates around its mean-field approximation for any level of edge density, provided the average degree grows to infinity. Moreover, we show that the opinion process exhibits propagation of chaos. We also give results for the trajectories of individual vertices and the stationary version of the opinion process, and prove that the limits in time and in the size of the network commute. The mean-field approximation is explicit and can be used to quantify consensus and polarization. △ Less

Submitted 16 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: Added references

PageRank Nibble on the sparse directed stochastic block model

Authors: Sayan Banerjee, Prabhanka Deka, Mariana Olvera-Cravioto

Abstract: We present new results on community recovery based on the PageRank Nibble algorithm on a sparse directed stochastic block model (dSBM). Our results are based on a characterization of the local weak limit of the dSBM and the limiting PageRank distribution. This characterization allows us to estimate the probability of misclassification for any given connection kernel and any given number of seeds (… ▽ More We present new results on community recovery based on the PageRank Nibble algorithm on a sparse directed stochastic block model (dSBM). Our results are based on a characterization of the local weak limit of the dSBM and the limiting PageRank distribution. This characterization allows us to estimate the probability of misclassification for any given connection kernel and any given number of seeds (vertices whose community label is known). The fact that PageRank is a local algorithm that can be efficiently computed in both a distributed and asynchronous fashion, makes it an appealing method for identifying members of a given community in very large networks where the identity of some vertices is known. △ Less

Submitted 12 March, 2023; originally announced March 2023.

Comments: PageRank Nibble, directed stochastic block model, local weak convergence, community detection

MSC Class: 0580; 05C81; 05C85; 60J80

Local weak limits for collapsed branching processes with random out-degrees

Authors: Sayan Banerjee, Prabhanka Deka, Mariana Olvera-Cravioto

Abstract: We obtain local weak limits in probability for Collapsed Branching Processes (CBP), which are directed random networks obtained by collapsing random-sized families of individuals in a general continuous-time branching process. The local weak limit of a given CBP, as the network grows, is shown to be a related continuous-time branching process stopped at an independent exponential time. This is don… ▽ More We obtain local weak limits in probability for Collapsed Branching Processes (CBP), which are directed random networks obtained by collapsing random-sized families of individuals in a general continuous-time branching process. The local weak limit of a given CBP, as the network grows, is shown to be a related continuous-time branching process stopped at an independent exponential time. This is done through an explicit coupling of the in-components of vertices with the limiting object. We also show that the in-components of a finite collection of uniformly chosen vertices locally weakly converge (in probability) to i.i.d. copies of the above limit, reminiscent of propagation of chaos in interacting particle systems. We obtain as special cases novel descriptions of the local weak limits of directed preferential and uniform attachment models. We also outline some applications of our results for analyzing the limiting in-degree and PageRank distributions. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 23 pages, 3 figures

Opinion dynamics on directed complex networks

Authors: Nicolas Fraiman, Tzu-Chi Lin, Mariana Olvera-Cravioto

Abstract: We propose and analyze a mathematical model for the evolution of opinions on directed complex networks. Our model generalizes the popular DeGroot and Friedkin-Johnsen models by allowing vertices to have attributes that may influence the opinion dynamics. We start by establishing sufficient conditions for the existence of a stationary opinion distribution on any fixed graph, and then provide an inc… ▽ More We propose and analyze a mathematical model for the evolution of opinions on directed complex networks. Our model generalizes the popular DeGroot and Friedkin-Johnsen models by allowing vertices to have attributes that may influence the opinion dynamics. We start by establishing sufficient conditions for the existence of a stationary opinion distribution on any fixed graph, and then provide an increasingly detailed characterization of its behavior by considering a sequence of directed random graphs having a local weak limit. Our most explicit results are obtained for graph sequences whose local weak limit is a marked Galton-Watson tree, in which case our model can be used to explain a variety of phenomena, e.g., conditions under which consensus can be achieved, mechanisms in which opinions can become polarized, and the effect of disruptive stubborn agents on the formation of opinions. △ Less

Submitted 10 January, 2024; v1 submitted 2 September, 2022; originally announced September 2022.

MSC Class: 60G10; 05C80; 60J05; 60J85

Strong couplings for static locally tree-like random graphs

Authors: Mariana Olvera-Cravioto

Abstract: The goal of this paper is to provide a general purpose result for the coupling of exploration processes of random graphs, both undirected and directed, with their local weak limits when this limit is a marked Galton-Watson process. This class includes in particular the configuration model and the family of inhomogeneous random graphs with rank-1 kernel. Vertices in the graph are allowed to have at… ▽ More The goal of this paper is to provide a general purpose result for the coupling of exploration processes of random graphs, both undirected and directed, with their local weak limits when this limit is a marked Galton-Watson process. This class includes in particular the configuration model and the family of inhomogeneous random graphs with rank-1 kernel. Vertices in the graph are allowed to have attributes on a general separable metric space and can potentially influence the construction of the graph itself. The coupling holds for any fixed depth of a breadth-first exploration process. △ Less

Submitted 31 January, 2022; v1 submitted 21 February, 2021; originally announced February 2021.

MSC Class: 05C80; 60J80

PageRank Asymptotics on Directed Preferential Attachment Networks

Authors: Sayan Banerjee, Mariana Olvera-Cravioto

Abstract: We characterize the tail behavior of the distribution of the PageRank of a uniformly chosen vertex in a directed preferential attachment graph and show that it decays as a power law with an explicit exponent that is described in terms of the model parameters. Interestingly, this power law is heavier than the tail of the limiting in-degree distribution, which goes against the commonly accepted {\em… ▽ More We characterize the tail behavior of the distribution of the PageRank of a uniformly chosen vertex in a directed preferential attachment graph and show that it decays as a power law with an explicit exponent that is described in terms of the model parameters. Interestingly, this power law is heavier than the tail of the limiting in-degree distribution, which goes against the commonly accepted {\em power law hypothesis}. This deviation from the power law hypothesis points at the structural differences between the inbound neighborhoods of typical vertices in a preferential attachment graph versus those in static random graph models where the power law hypothesis has been proven to hold (e.g., directed configuration models and inhomogeneous random digraphs). In addition to characterizing the PageRank distribution of a typical vertex, we also characterize the explicit growth rate of the PageRank of the oldest vertex as the network size grows. △ Less

Submitted 18 October, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: 26 pages

MSC Class: Primary:05C80; Secondary: 60J80; 68P10; 41A60; 60B10

Stochastic recursions on directed random graphs

Authors: Nicolas Fraiman, Tzu-Chi Lin, Mariana Olvera-Cravioto

Abstract: For a directed graph $G(V_n, E_n)$ on the vertices $V_n = \{1,2, \dots, n\}$, we study the distribution of a Markov chain $\{ {\bf R}^{(k)}: k \geq 0\}$ on $\mathbb{R}^n$ such that the $i$th component of ${\bf R}^{(k)}$, denoted $R_i^{(k)}$, corresponds to the value of the process on vertex $i$ at time $k$. We focus on processes $\{ {\bf R}^{(k)}: k \geq 0\}$ where the value of $R_i^{(k+1)}$ depen… ▽ More For a directed graph $G(V_n, E_n)$ on the vertices $V_n = \{1,2, \dots, n\}$, we study the distribution of a Markov chain $\{ {\bf R}^{(k)}: k \geq 0\}$ on $\mathbb{R}^n$ such that the $i$th component of ${\bf R}^{(k)}$, denoted $R_i^{(k)}$, corresponds to the value of the process on vertex $i$ at time $k$. We focus on processes $\{ {\bf R}^{(k)}: k \geq 0\}$ where the value of $R_i^{(k+1)}$ depends only on the values $\{ R_j^{(k)}: j \to i\}$ of its inbound neighbors, and possibly on vertex attributes. We then show that, provided $G(V_n, E_n)$ converges in the local weak sense to a marked Galton-Watson process, the dynamics of the process for a uniformly chosen vertex in $V_n$ can be coupled, for any fixed $k$, to a process $\{ \mathcal{R}_\emptyset^{(r)}: 0 \leq r \leq k\}$ constructed on the limiting marked Galton-Watson tree. Moreover, we derive sufficient conditions under which $\mathcal{R}^{(k)}_\emptyset$ converges, as $k \to \infty$, to a random variable $\mathcal{R}^*$ that can be characterized in terms of the attracting endogenous solution to a branching distributional fixed-point equation. Our framework can also be applied to processes $\{ {\bf R}^{(k)}: k \geq 0\}$ whose only source of randomness comes from the realization of the graph $G(V_n, E_n)$. △ Less

Submitted 26 October, 2022; v1 submitted 19 October, 2020; originally announced October 2020.

MSC Class: 60G10; 05C80; 60J05 (Primary); 60J85 (Secondary)

Importance sampling for maxima on trees

Authors: Bojan Basrak, Michael Conroy, Mariana Olvera-Cravioto, Zbigniew Palmowski

Abstract: We consider the distributional fixed-point equation: $$R \stackrel{\mathcal{D}}{=} Q \vee \left( \bigvee_{i=1}^N C_i R_i \right),$$ where the $\{R_i\}$ are i.i.d.~copies of $R$, independent of the vector $(Q, N, \{C_i\})$, where $N \in \mathbb{N}$, $Q, \{C_i\} \geq 0$ and $P(Q > 0) > 0$. By setting $W = \log R$, $X_i = \log C_i$, $Y = \log Q$ it is equivalent to the high-order Lindley equation… ▽ More We consider the distributional fixed-point equation: $$R \stackrel{\mathcal{D}}{=} Q \vee \left( \bigvee_{i=1}^N C_i R_i \right),$$ where the $\{R_i\}$ are i.i.d.~copies of $R$, independent of the vector $(Q, N, \{C_i\})$, where $N \in \mathbb{N}$, $Q, \{C_i\} \geq 0$ and $P(Q > 0) > 0$. By setting $W = \log R$, $X_i = \log C_i$, $Y = \log Q$ it is equivalent to the high-order Lindley equation $$W \stackrel{\mathcal{D}}{=} \max\left\{ Y, \, \max_{1 \leq i \leq N} (X_i + W_i) \right\}.$$ It is known that under Kesten assumptions, $$P(W > t) \sim H e^{-αt}, \qquad t \to \infty,$$ where $α>0$ solves the Cramér-Lundberg equation $E \left[ \sum_{j=1}^N C_i ^α\right] = E\left[ \sum_{i=1}^N e^{αX_i} \right] = 1$. The main goal of this paper is to provide an explicit representation for $P(W > t)$, which can be directly connected to the underlying weighted branching process where $W$ is constructed and that can be used to construct unbiased and strongly efficient estimators for all $t$. Furthermore, we show how this new representation can be directly analyzed using Alsmeyer's Markov renewal theorem, yielding an alternative representation for the constant $H$. We provide numerical examples illustrating the use of this new algorithm. △ Less

Submitted 11 September, 2020; v1 submitted 19 April, 2020; originally announced April 2020.

PageRank's behavior under degree-degree correlations

Authors: Mariana Olvera-Cravioto

Abstract: The focus of this work is the asymptotic analysis of the tail distribution of Google's PageRank algorithm on large scale-free directed networks. In particular, the main theorem provides the convergence, in the Kantorovich-Rubinstein metric, of the rank of a randomly chosen vertex in graphs generated via either a directed configuration model or an inhomogeneous random digraph. The theorem fully cha… ▽ More The focus of this work is the asymptotic analysis of the tail distribution of Google's PageRank algorithm on large scale-free directed networks. In particular, the main theorem provides the convergence, in the Kantorovich-Rubinstein metric, of the rank of a randomly chosen vertex in graphs generated via either a directed configuration model or an inhomogeneous random digraph. The theorem fully characterizes the limiting distribution by expressing it as a random sum of i.i.d.~copies of the attracting endogenous solution to a branching distributional fixed-point equation. In addition, we provide the asymptotic tail behavior of the limit and use it to explain the effect that in-degree/out-degree correlations in the underlying graph can have on the qualitative performance of PageRank. △ Less

Submitted 20 September, 2019; originally announced September 2019.

MSC Class: 05C80; 60J80; 68P20; 41A60; 37A30; 60B10

Last-mile shared delivery: A discrete sequential packing approach

Authors: Junyu Cao, Mariana Olvera-Cravioto, Zuo-Jun Shen

Abstract: We propose a model for optimizing the last-mile delivery of n packages, from a distribution center to their final recipients, using a strategy that combines the use of ride-sharing platforms (e.g., Uber or Lyft) with traditional in-house van delivery systems. The main objective is to compute the optimal reward offered to private drivers for each of the n packages, such that the total expected cost… ▽ More We propose a model for optimizing the last-mile delivery of n packages, from a distribution center to their final recipients, using a strategy that combines the use of ride-sharing platforms (e.g., Uber or Lyft) with traditional in-house van delivery systems. The main objective is to compute the optimal reward offered to private drivers for each of the n packages, such that the total expected cost of delivering all packages is minimized. Our technical approach is based on the formulation of a discrete sequential packing problem, where bundles of packages are picked up from the warehouse at random times during the interval [0, T]. Our theoretical results include both exact and asymptotic (as $n \to \infty$) expressions for the expected number of packages that will be picked up by time T, and are closely related to the classical Renyi's parking/packing problem. △ Less

Submitted 14 May, 2018; originally announced May 2018.

MSC Class: Primary: 60F10; secondary: 60C05; 60G55

On a general class of inhomogeneous random digraphs

Authors: Junyu Cao, Mariana Olvera-Cravioto

Abstract: We study a family of directed random graphs whose arcs are sampled independently of each other, and are present in the graph with a probability that depends on the attributes of the vertices involved. In particular, this family of models includes as special cases the directed versions of the Erdos-Renyi model, graphs with given expected degrees, the generalized random graph, and the Poissonian ran… ▽ More We study a family of directed random graphs whose arcs are sampled independently of each other, and are present in the graph with a probability that depends on the attributes of the vertices involved. In particular, this family of models includes as special cases the directed versions of the Erdos-Renyi model, graphs with given expected degrees, the generalized random graph, and the Poissonian random graph. We establish the phase transition for the existence of a giant strongly connected component and provide some other basic properties, including the limiting joint distribution of the degrees and the mean number of arcs. In particular, we show that by choosing the joint distribution of the vertex attributes according to a multivariate regularly varying distribution, one can obtain scale-free graphs with arbitrary in-degree/out-degree dependence. △ Less

Submitted 8 December, 2017; originally announced December 2017.

MSC Class: Primary 05C80; Secondary 90B15; 60C05

Likelihood Ratio Gradient Estimation for Steady-State Parameters

Authors: Peter W. Glynn, Mariana Olvera-Cravioto

Abstract: We consider a discrete-time Markov chain $\boldsymbolΦ$ on a general state-space ${\sf X}$, whose transition probabilities are parameterized by a real-valued vector $\boldsymbolθ$. Under the assumption that $\boldsymbolΦ$ is geometrically ergodic with corresponding stationary distribution $π(\boldsymbolθ)$, we are interested in estimating the gradient $\nabla α(\boldsymbolθ)$ of the steady-state e… ▽ More We consider a discrete-time Markov chain $\boldsymbolΦ$ on a general state-space ${\sf X}$, whose transition probabilities are parameterized by a real-valued vector $\boldsymbolθ$. Under the assumption that $\boldsymbolΦ$ is geometrically ergodic with corresponding stationary distribution $π(\boldsymbolθ)$, we are interested in estimating the gradient $\nabla α(\boldsymbolθ)$ of the steady-state expectation $$α(\boldsymbolθ) = π( \boldsymbolθ) f.$$ To this end, we first give sufficient conditions for the differentiability of $α(\boldsymbolθ)$ and for the calculation of its gradient via a sequence of finite horizon expectations. We then propose two different likelihood ratio estimators and analyze their limiting behavior. △ Less

Submitted 9 March, 2018; v1 submitted 9 July, 2017; originally announced July 2017.

MSC Class: Primary: 65C05; 60J22; Secondary: 60J10; 60G42

PageRank on inhomogeneous random digraphs

Authors: Jiung Lee, Mariana Olvera-Cravioto

Abstract: We study the typical behavior of a generalized version of Google's PageRank algorithm on a large family of inhomogeneous random digraphs. This family includes as special cases directed versions of classical models such as the Erdös-Rényi model, the Chung-Lu model, the Poissonian random graph and the generalized random graph, and is suitable for modeling scale-free directed complex networks where t… ▽ More We study the typical behavior of a generalized version of Google's PageRank algorithm on a large family of inhomogeneous random digraphs. This family includes as special cases directed versions of classical models such as the Erdös-Rényi model, the Chung-Lu model, the Poissonian random graph and the generalized random graph, and is suitable for modeling scale-free directed complex networks where the number of neighbors a vertex has is related to its attributes. In particular, we show that the rank of a randomly chosen node in a graph from this family converges weakly to the attracting endogenous solution to the stochastic fixed-point equation $$\mathcal{R} \stackrel{\mathcal{D}}{=} \sum_{i=1}^{\mathcal{N}} \mathcal{C}_i \mathcal{R}_i + \mathcal{Q},$$ where $(\mathcal{N}, \mathcal{Q}, \{\mathcal{C}_i\}_{i \geq 1})$ is a real-valued vector with $\mathcal{N}\in \{0,1,2,...\}$, the $\{ \mathcal{R}_i\}$ are i.i.d.~copies of $\mathcal{R}$, independent of $(\mathcal{N}, \mathcal{Q}, \{\mathcal{C}_i\}_{i \geq 1})$, with $\{ \mathcal{C}_i\}$ i.i.d.~and independent of $(\mathcal{N}, \mathcal{Q})$; $\stackrel{\mathcal{D}}{=}$ denotes equality in distribution. This result can then be used to provide further evidence of the power-law behavior of PageRank on scale-free graphs. △ Less

Submitted 11 November, 2020; v1 submitted 8 July, 2017; originally announced July 2017.

MSC Class: Primary: 05C80; 60J80; 68P20. Secondary: 41A60; 37A30; 60B10

Journal ref: Stochastic Processes and their Applications, Vol.130, No. 4, pp. 2312-2348, 2020

Convergence of the Population Dynamics algorithm in the Wasserstein metric

Authors: Mariana Olvera-Cravioto

Abstract: We study the convergence of the population dynamics algorithm, which produces sample pools of random variables having a distribution that closely approximates that of the {\em special endogenous solution} to a stochastic fixed-point equation of the form: $$R\stackrel{\mathcal D}{=} Φ( Q, N, \{ C_i \}, \{R_i\}),$$ where $(Q, N, \{C_i\})$ is a real-valued random vector with $N \in \mathbb{N}$, and… ▽ More We study the convergence of the population dynamics algorithm, which produces sample pools of random variables having a distribution that closely approximates that of the {\em special endogenous solution} to a stochastic fixed-point equation of the form: $$R\stackrel{\mathcal D}{=} Φ( Q, N, \{ C_i \}, \{R_i\}),$$ where $(Q, N, \{C_i\})$ is a real-valued random vector with $N \in \mathbb{N}$, and $\{R_i\}_{i \in \mathbb{N}}$ is a sequence of i.i.d. copies of $R$, independent of $(Q, N, \{C_i\})$; the symbol $\stackrel{\mathcal{D}}{=}$ denotes equality in distribution. Specifically, we show its convergence in the Wasserstein metric of order $p$ ($p \geq 1$) and prove the consistency of estimators based on the sample pool produced by the algorithm. △ Less

Submitted 8 February, 2018; v1 submitted 26 May, 2017; originally announced May 2017.

Typical distances in the directed configuration model

Authors: Pim van der Hoorn, Mariana Olvera-Cravioto

Abstract: We analyze the distribution of the distance between two nodes, sampled uniformly at random, in digraphs generated via the directed configuration model, in the supercritical regime. Under the assumption that the covariance between the in-degree and out-degree is finite, we show that the distance grows logarithmically in the size of the graph. In contrast with the undirected case, this can happen ev… ▽ More We analyze the distribution of the distance between two nodes, sampled uniformly at random, in digraphs generated via the directed configuration model, in the supercritical regime. Under the assumption that the covariance between the in-degree and out-degree is finite, we show that the distance grows logarithmically in the size of the graph. In contrast with the undirected case, this can happen even when the variance of the degrees is infinite. The main tool in the analysis is a new coupling between a breadth-first graph exploration process and a suitable branching process based on the Kantorovich-Rubinstein metric. This coupling holds uniformly for a much larger number of steps in the exploration process than existing ones, and is therefore of independent interest. △ Less

Submitted 21 April, 2017; v1 submitted 14 November, 2015; originally announced November 2015.

MSC Class: 05C80; 60B10

Efficient Simulation for Branching Linear Recursions

Authors: Ningyuan Chen, Mariana Olvera-Cravioto

Abstract: We consider a linear recursion of the form $$R^{(k+1)}\stackrel{\mathcal D}{=}\sum_{i=1}^{N}C_iR^{(k)}_i+Q,$$ where $(Q,N,C_1,C_2,\dots)$ is a real-valued random vector with $N\in\mathbb{N}=\{0, 1, 2, \dots\}$, $\{R^{(k)}_i\}_{i\in\mathbb{N}}$ is a sequence of i.i.d. copies of $R^{(k)}$, independent of $(Q,N,C_1,C_2,\dots)$, and $\stackrel{\mathcal{D}}{=}$ denotes equality in distribution. For sui… ▽ More We consider a linear recursion of the form $$R^{(k+1)}\stackrel{\mathcal D}{=}\sum_{i=1}^{N}C_iR^{(k)}_i+Q,$$ where $(Q,N,C_1,C_2,\dots)$ is a real-valued random vector with $N\in\mathbb{N}=\{0, 1, 2, \dots\}$, $\{R^{(k)}_i\}_{i\in\mathbb{N}}$ is a sequence of i.i.d. copies of $R^{(k)}$, independent of $(Q,N,C_1,C_2,\dots)$, and $\stackrel{\mathcal{D}}{=}$ denotes equality in distribution. For suitable vectors $(Q,N,C_1,C_2,\dots)$ and provided the initial distribution of $R^{(0)}$ is well-behaved, the process $R^{(k)}$ is known to converge to the endogenous solution of the corresponding stochastic fixed-point equation, which appears in the analysis of information ranking algorithms, e.g., PageRank, and in the complexity analysis of divide and conquer algorithms, e.g. Quicksort. Naive Monte Carlo simulation of $R^{(k)}$ based on the branching recursion has exponential complexity in $k$, and therefore the need for efficient methods. We propose in this paper an iterative bootstrap algorithm that has linear complexity and can be used to approximately sample $R^{(k)}$. We show the consistency of estimators based on our proposed algorithm. △ Less

Submitted 31 March, 2015; originally announced March 2015.

Comments: submitted to WSC 2015

MSC Class: 68U20; 62F40; 60B10

Parallel queues with synchronization

Authors: Mariana Olvera-Cravioto, Octavio Ruiz-Lacedelli

Abstract: Motivated by the growing interest in today's massive parallel computing capabilities we analyze a queueing network with many servers in parallel to which jobs arrive a according to a Poisson process. Each job, upon arrival, is split into several pieces which are randomly routed to specific servers in the network, without centralized information about the status of the servers' individual queues. T… ▽ More Motivated by the growing interest in today's massive parallel computing capabilities we analyze a queueing network with many servers in parallel to which jobs arrive a according to a Poisson process. Each job, upon arrival, is split into several pieces which are randomly routed to specific servers in the network, without centralized information about the status of the servers' individual queues. The main feature of this system is that the different pieces of a job must initiate their service in a synchronized fashion. Moreover, the system operates in a FCFS basis. The synchronization and service discipline create blocking and idleness among the servers, which is compensated by the fast service time attained through the parallelization of the work. We analyze the stationary waiting time distribution of jobs under a many servers limit and provide exact tail asymptotics; these asymptotics generalize the celebrated Cramér-Lundberg approximation for the single-server queue. △ Less

Submitted 16 July, 2015; v1 submitted 31 December, 2014; originally announced January 2015.

MSC Class: 60K25; 60J85; 68M20; 68M14

Coupling on weighted branching trees

Authors: Ningyuan Chen, Mariana Olvera-Cravioto

Abstract: This paper considers linear functions constructed on two different weighted branching processes and provides explicit bounds for their Kantorovich-Rubinstein distance in terms of couplings of their corresponding generic branching vectors. Motivated by applications to the analysis of random graphs, we also consider a variation of the weighted branching process where the generic branching vector has… ▽ More This paper considers linear functions constructed on two different weighted branching processes and provides explicit bounds for their Kantorovich-Rubinstein distance in terms of couplings of their corresponding generic branching vectors. Motivated by applications to the analysis of random graphs, we also consider a variation of the weighted branching process where the generic branching vector has a different dependence structure from the usual one. By applying the bounds to sequences of weighted branching processes, we derive sufficient conditions for the convergence in the Kantorovich-Rubinstein distance of linear functions. We focus on the case where the limits are endogenous fixed points of suitable smoothing transformations. △ Less

Submitted 25 June, 2015; v1 submitted 4 October, 2014; originally announced October 2014.

MSC Class: 60J80; 60B10; 60H25

Ranking algorithms on directed configuration networks

Authors: Ningyuan Chen, Nelly Litvak, Mariana Olvera-Cravioto

Abstract: This paper studies the distribution of a family of rankings, which includes Google's PageRank, on a directed configuration model. In particular, it is shown that the distribution of the rank of a randomly chosen node in the graph converges in distribution to a finite random variable $\mathcal{R}^*$ that can be written as a linear combination of i.i.d. copies of the endogenous solution to a stochas… ▽ More This paper studies the distribution of a family of rankings, which includes Google's PageRank, on a directed configuration model. In particular, it is shown that the distribution of the rank of a randomly chosen node in the graph converges in distribution to a finite random variable $\mathcal{R}^*$ that can be written as a linear combination of i.i.d. copies of the endogenous solution to a stochastic fixed point equation of the form $$\mathcal{R} \stackrel{\mathcal{D}}{=} \sum_{i=1}^{\mathcal{N}} \mathcal{C}_i \mathcal{R}_i + \mathcal{Q},$$ where $(\mathcal{Q}, \mathcal{N}, \{ \mathcal{C}_i\})$ is a real-valued vector with $\mathcal{N} \in \{0,1,2,\dots\}$, $P(|\mathcal{Q}| > 0) > 0$, and the $\{\mathcal{R}_i\}$ are i.i.d. copies of $\mathcal{R}$, independent of $(\mathcal{Q}, \mathcal{N}, \{ \mathcal{C}_i\})$. Moreover, we provide precise asymptotics for the limit $\mathcal{R}^*$, which when the in-degree distribution in the directed configuration model has a power law imply a power law distribution for $\mathcal{R}^*$ with the same exponent. △ Less

Submitted 12 October, 2014; v1 submitted 25 September, 2014; originally announced September 2014.

MSC Class: 05C80; 60J80; 68P20

PageRank in scale-free random graphs

Authors: Ningyuan Chen, Nelly Litvak, Mariana Olvera-Cravioto

Abstract: We analyze the distribution of PageRank on a directed configuration model and show that as the size of the graph grows to infinity it can be closely approximated by the PageRank of the root node of an appropriately constructed tree. This tree approximation is in turn related to the solution of a linear stochastic fixed point equation that has been thoroughly studied in the recent literature. We analyze the distribution of PageRank on a directed configuration model and show that as the size of the graph grows to infinity it can be closely approximated by the PageRank of the root node of an appropriately constructed tree. This tree approximation is in turn related to the solution of a linear stochastic fixed point equation that has been thoroughly studied in the recent literature. △ Less

Submitted 15 August, 2014; originally announced August 2014.

Maximums on Trees

Authors: Predrag R. Jelenkovic, Mariana Olvera-Cravioto

Abstract: We study the minimal/endogenous solution $R$ to the maximum recursion on weighted branching trees given by $$R\stackrel{\mathcal{D}}{=}\left(\bigvee_{i=1}^NC_iR_i \right)\vee Q,$$ where $(Q,N,C_1,C_2,\dots)$ is a random vector with $N\in \mathbb{N}\cup\{\infty\}$, $P(|Q|>0)>0$ and nonnegative weights $\{C_i\}$, and $\{R_i\}_{i\in\mathbb{N}}$ is a sequence of i.i.d. copies of $R$ independent of… ▽ More We study the minimal/endogenous solution $R$ to the maximum recursion on weighted branching trees given by $$R\stackrel{\mathcal{D}}{=}\left(\bigvee_{i=1}^NC_iR_i \right)\vee Q,$$ where $(Q,N,C_1,C_2,\dots)$ is a random vector with $N\in \mathbb{N}\cup\{\infty\}$, $P(|Q|>0)>0$ and nonnegative weights $\{C_i\}$, and $\{R_i\}_{i\in\mathbb{N}}$ is a sequence of i.i.d. copies of $R$ independent of $(Q,N,C_1,C_2,\dots)$; $\stackrel{\mathcal{D}}{=}$ denotes equality in distribution. Furthermore, when $Q>0$ this recursion can be transformed into its additive equivalent, which corresponds to the maximum of a branching random walk and is also known as a high-order Lindley equation. We show that, under natural conditions, the asymptotic behavior of $R$ is power-law, i.e., $P(|R|>x)\sim Hx^{-α}$, for some $α>0$ and $H>0$. This has direct implications for the tail behavior of other well known branching recursions. △ Less

Submitted 23 May, 2014; originally announced May 2014.

Comments: arXiv admin note: text overlap with arXiv:1006.3295

MSC Class: 60H25; 60J80; 60F10; 60K05

Directed random graphs with given degree distributions

Authors: Ningyuan Chen, Mariana Olvera-Cravioto

Abstract: Given two distributions F and G on the nonnegative integers we propose an algorithm to construct in- and out-degree sequences from samples of i.i.d. observations from F and G, respectively, that with high probability will be graphical, that is, from which a simple directed graph can be drawn. We then analyze a directed version of the configuration model and show that, provided that F and G have fi… ▽ More Given two distributions F and G on the nonnegative integers we propose an algorithm to construct in- and out-degree sequences from samples of i.i.d. observations from F and G, respectively, that with high probability will be graphical, that is, from which a simple directed graph can be drawn. We then analyze a directed version of the configuration model and show that, provided that F and G have finite variance, the probability of obtaining a simple graph is bounded away from zero as the number of nodes grows. We show that conditional on the resulting graph being simple, the in- and out-degree distributions are (approximately) F and G for large size graphs. Moreover, when the degree distributions have only finite mean we show that the elimination of self-loops and multiple edges does not significantly change the degree distributions in the resulting simple graph. △ Less

Submitted 10 July, 2012; originally announced July 2012.

MSC Class: 05C80; 60C05

Tail behavior of solutions of linear recursions on trees

Authors: Mariana Olvera-Cravioto

Abstract: Consider the linear nonhomogeneous fixed point equation R =_d sum_{i=1}^N C_i R_i + Q, where (Q,N,C_1,...,C_N) is a random vector with N in{0,1,2,3,...}U{infty}, {C_i}_{i=1}^N >= 0, P(|Q|>0) > 0, and {R_i}_{i=1}^N is a sequence of i.i.d. random variables independent of (Q,N,C_1,...,C_N) having the same distribution as R. It is known that R will have a heavy-tailed distribution under several diff… ▽ More Consider the linear nonhomogeneous fixed point equation R =_d sum_{i=1}^N C_i R_i + Q, where (Q,N,C_1,...,C_N) is a random vector with N in{0,1,2,3,...}U{infty}, {C_i}_{i=1}^N >= 0, P(|Q|>0) > 0, and {R_i}_{i=1}^N is a sequence of i.i.d. random variables independent of (Q,N,C_1,...,C_N) having the same distribution as R. It is known that R will have a heavy-tailed distribution under several different sets of assumptions on the vector (Q,N,C_1,...,C_N). This paper investigates the settings where either Z_N = sum_{i=1}^N C_i or Q are regularly varying with index -alpha < -1 and E[sum_{i=1}^N C_i^alpha] < 1. This work complements previous results showing that P(R>t) Ht^{-alpha} provided there exists a solution alpha > 0 to the equation E[sum_{i=1}^N|C_i|^alpha] = 1, and both Q and Z_N have lighter tails. △ Less

Submitted 18 August, 2011; originally announced August 2011.

MSC Class: 60H25; 60J80; 60F10; 60K05

Asymptotics for Weighted Random Sums

Authors: Mariana Olvera-Cravioto

Abstract: Let $\{X_i\}$ be a sequence of independent identically distributed random variables with an intermediate regularly varying (IR) right tail $\bar{F}$. Let $(N, C_1, ..., C_N)$ be a nonnegative random vector independent of the $\{X_i\}$ with $N \in \mathbb{N} \cup \{\infty\}$. We study the weighted random sum $S_N = \sum_{i=1}^N C_i X_i$, and its maximum,… ▽ More Let $\{X_i\}$ be a sequence of independent identically distributed random variables with an intermediate regularly varying (IR) right tail $\bar{F}$. Let $(N, C_1, ..., C_N)$ be a nonnegative random vector independent of the $\{X_i\}$ with $N \in \mathbb{N} \cup \{\infty\}$. We study the weighted random sum $S_N = \sum_{i=1}^N C_i X_i$, and its maximum, $M_N = \sup_{1 \leq k < N+1} \sum_{i=1}^k C_i X_i$. These type of sums appear in the analysis of stochastic recursions, including weighted branching processes and autoregressive processes. In particular, we derive conditions under which $$P(M_N > x) \sim P(S_N > x) \sim E[\sum_{i=1}^N \bar{F}(x/C_i)],$$ as $x \to \infty$. When $E[X_1] > 0$ and the distribution of $Z_N = \sum_{i=1}^N C_i$ is also IR, we obtain the asymptotics $$P(M_N > x) \sim P(S_N > x) \sim E[\sum_{i=1}^N \bar{F}(x/C_i)] + P(Z_N > x/E[X_1]).$$ For completeness, when the distribution of $Z_N$ is IR and heavier than $\bar{F}$, we also obtain conditions under which the asymptotic relations $$P(M_N > x) \sim P(S_N > x) \sim P(Z_N > x/E[X_1])$$ hold. △ Less

Submitted 17 April, 2012; v1 submitted 1 February, 2011; originally announced February 2011.

MSC Class: 60G50; 60F10; 60J80; 60G70

Implicit Renewal Theorem for Trees with General Weights

Authors: Predrag R. Jelenkovic, Mariana Olvera-Cravioto

Abstract: Consider distributional fixed point equations of the form R =d f(C_i, R_i, 1 <= i <= N), where f(.) is a possibly random real valued function, N in {0, 1, 2, 3,...} U {infty}, {C_i}_{i=1}^N are real valued random weights and {R_i}_{i >= 1} are iid copies of R, independent of (N, C_1,..., C_N); =d represents equality in distribution. Fixed point equations of this type are of utmost importance for s… ▽ More Consider distributional fixed point equations of the form R =d f(C_i, R_i, 1 <= i <= N), where f(.) is a possibly random real valued function, N in {0, 1, 2, 3,...} U {infty}, {C_i}_{i=1}^N are real valued random weights and {R_i}_{i >= 1} are iid copies of R, independent of (N, C_1,..., C_N); =d represents equality in distribution. Fixed point equations of this type are of utmost importance for solving many applied probability problems, ranging from average case analysis of algorithms to statistical physics. We develop an Implicit Renewal Theorem that enables the characterization of the power tail behavior of the solutions R to many equations of multiplicative nature that fall in this category. This result extends the prior work in Jelenkovic and Olvera-Cravioto (2010), which assumed nonnegative weights {C_i}, to general real valued weights. We illustrate the developed theorem by deriving the power tail asymptotics of the solution R to the linear equation R =d sum_{i=1}^N C_i R_i + Q. △ Less

Submitted 19 October, 2011; v1 submitted 9 December, 2010; originally announced December 2010.

Comments: 1 figure, 30 pages

MSC Class: 60H26; 60F10; 60K05; 60J80

On the transition from heavy traffic to heavy tails for the M/G/1 queue: The regularly varying case

Authors: Mariana Olvera-Cravioto, Jose Blanchet, Peter Glynn

Abstract: Two of the most popular approximations for the distribution of the steady-state waiting time, $W_{\infty}$, of the M/G/1 queue are the so-called heavy-traffic approximation and heavy-tailed asymptotic, respectively. If the traffic intensity, $ρ$, is close to 1 and the processing times have finite variance, the heavy-traffic approximation states that the distribution of $W_{\infty}$ is roughly expo… ▽ More Two of the most popular approximations for the distribution of the steady-state waiting time, $W_{\infty}$, of the M/G/1 queue are the so-called heavy-traffic approximation and heavy-tailed asymptotic, respectively. If the traffic intensity, $ρ$, is close to 1 and the processing times have finite variance, the heavy-traffic approximation states that the distribution of $W_{\infty}$ is roughly exponential at scale $O((1-ρ)^{-1})$, while the heavy tailed asymptotic describes power law decay in the tail of the distribution of $W_{\infty}$ for a fixed traffic intensity. In this paper, we assume a regularly varying processing time distribution and obtain a sharp threshold in terms of the tail value, or equivalently in terms of $(1-ρ)$, that describes the point at which the tail behavior transitions from the heavy-traffic regime to the heavy-tailed asymptotic. We also provide new approximations that are either uniform in the traffic intensity, or uniform on the positive axis, that avoid the need to use different expressions on the two regions defined by the threshold. △ Less

Submitted 7 April, 2011; v1 submitted 27 September, 2010; originally announced September 2010.

Comments: Published in at http://dx.doi.org/10.1214/10-AAP707 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AAP-AAP707

Journal ref: Annals of Applied Probability 2011, Vol. 21, No. 2, 645-668

Uniform Approximations for the M/G/1 Queue with Subexponential Processing Times

Authors: Mariana Olvera-Cravioto, Peter W. Glynn

Abstract: This paper studies the asymptotic behavior of the steady-state waiting time, W_infty, of the M/G/1 queue with subexponenential processing times for different combinations of traffic intensities and overflow levels. In particular, we provide insights into the regions of large deviations where the so-called heavy traffic approximation and heavy tail asymptotic hold. For queues whose service time dis… ▽ More This paper studies the asymptotic behavior of the steady-state waiting time, W_infty, of the M/G/1 queue with subexponenential processing times for different combinations of traffic intensities and overflow levels. In particular, we provide insights into the regions of large deviations where the so-called heavy traffic approximation and heavy tail asymptotic hold. For queues whose service time distribution decays slower than e^{-sqrt{t}} we identify a third region of asymptotics where neither the heavy traffic nor the heavy tailed approximations are valid. These results are obtained by deriving approximations for P(W_infty > x) that are either uniform in the traffic intensity as the tail value goes to infinity or uniform on the positive axis as the traffic intensity converges to one. Our approach makes clear the connection between the asymptotic behavior of the steady-state waiting time distribution and that of an associated random walk. △ Less

Submitted 19 March, 2011; v1 submitted 27 September, 2010; originally announced September 2010.

MSC Class: 60K25; 68M20; 60F10

Implicit Renewal Theory and Power Tails on Trees

Authors: Predrag R. Jelenković, Mariana Olvera-Cravioto

Abstract: We extend Goldie's (1991) Implicit Renewal Theorem to enable the analysis of recursions on weighted branching trees. We illustrate the developed method by deriving the power tail asymptotics of the distributions of the solutions R to: R =_D sum_{i=1}^N C_i R_i + Q, R =_D max(max_{i=1}^N C_i R_i, Q), and similar recursions, where (Q, N, C_1,..., C_N) is a nonnegative random vector with N in {0, 1,… ▽ More We extend Goldie's (1991) Implicit Renewal Theorem to enable the analysis of recursions on weighted branching trees. We illustrate the developed method by deriving the power tail asymptotics of the distributions of the solutions R to: R =_D sum_{i=1}^N C_i R_i + Q, R =_D max(max_{i=1}^N C_i R_i, Q), and similar recursions, where (Q, N, C_1,..., C_N) is a nonnegative random vector with N in {0, 1, 2, 3, ..., infinity}, and {R_i}_{i >= 1} are iid copies of R, independent of (Q, N, C_1,..., C_N); =_D denotes the equality in distribution. △ Less

Submitted 31 May, 2012; v1 submitted 16 June, 2010; originally announced June 2010.

Comments: arXiv admin note: substantial text overlap with arXiv:1012.2165

MSC Class: 60H25 (Primary) 60J80; 60F10; 60K05 (Secondary)

Information Ranking and Power Laws on Trees

Authors: Predrag R. Jelenkovic, Mariana Olvera-Cravioto

Abstract: We study the situations when the solution to a weighted stochastic recursion has a power law tail. To this end, we develop two complementary approaches, the first one extends Goldie's (1991) implicit renewal theorem to cover recursions on trees; and the second one is based on a direct sample path large deviations analysis of weighted recursive random sums. We believe that these methods may be of i… ▽ More We study the situations when the solution to a weighted stochastic recursion has a power law tail. To this end, we develop two complementary approaches, the first one extends Goldie's (1991) implicit renewal theorem to cover recursions on trees; and the second one is based on a direct sample path large deviations analysis of weighted recursive random sums. We believe that these methods may be of independent interest in the analysis of more general weighted branching processes as well as in the analysis of algorithms. △ Less

Submitted 28 July, 2010; v1 submitted 11 May, 2009; originally announced May 2009.

MSC Class: 60H25; 60J80