Asymptotic equivalence and contiguity of some random graphs

Svante Janson

Introduction

There are many different models of random graphs. Sometimes, the differences are minor, and it can be guessed that the asymptotic behaviour of two models are the same (for all or at least for some interesting properties). This note concerns some cases where it is possible to actually prove such results in a strong form. We begin by defining the two types of asymptotic equality that we will study. All unspecified limits are as n→∞{n\to\infty}.

Let (Xn,An)({\mathcal{X}}_{n},\mathcal{A}_{n}), n≥1n\geq 1, be a sequence of arbitrary measurable spaces and let PnP_{n} and QnQ_{n} be two probability measures on (Xn,An)({\mathcal{X}}_{n},\mathcal{A}_{n}).

The sequence (Pn)n(P_{n})_{n} is asymptotically equivalent to (Qn)n(Q_{n})_{n}, denoted by (Pn)n≅(Qn)n(P_{n})_{n}\cong(Q_{n})_{n}, if for every sequence of measurable sets AnA_{n} (i.e., An∈AnA_{n}\in\mathcal{A}_{n}), we have Pn(An)−Qn(An)→0P_{n}(A_{n})-Q_{n}(A_{n})\to 0.

The sequence (Pn)n(P_{n})_{n} is contiguous with respect to (Qn)n(Q_{n})_{n}, denoted by (Pn)n⊲(Qn)n(P_{n})_{n}\vartriangleleft(Q_{n})_{n}, if for every sequence of measurable sets AnA_{n} such that Qn(An)→0Q_{n}(A_{n})\to 0, we also have Pn(An)→0P_{n}(A_{n})\to 0.

Note that asymptotic equivalence is a symmetric relation while contiguity is not; we say that (Pn)n(P_{n})_{n} and (Qn)n(Q_{n})_{n} are (mutually) contiguous, (Pn)n⊲⊳(Qn)n(P_{n})_{n}\vartriangleleft\vartriangleright(Q_{n})_{n}, if both (Pn)n⊲(Qn)n(P_{n})_{n}\vartriangleleft(Q_{n})_{n} and (Qn)n⊲(Pn)n(Q_{n})_{n}\vartriangleleft(P_{n})_{n}, i.e., if Pn(An)→0  ⟺  Qn(An)→0P_{n}(A_{n})\to 0\iff Q_{n}(A_{n})\to 0 for any sequence of measurable sets An⊆XnA_{n}\subseteq{\mathcal{X}}_{n}. (And similarly for sequences of random variables XnX_{n} and YnY_{n}.)

Asymptotic equivalence implies contiguity, but not conversely (see e.g. Example 1.2 and Remark 1.6), so contiguity is a weaker property.

We illustrate these notions by two simple examples.

In the special case of two constant sequences, Pn=PP_{n}=P and Qn=QQ_{n}=Q where PP and QQ are two probability measures defined on the same space (Xn,An)=(X;A)({\mathcal{X}}_{n},\mathcal{A}_{n})=({\mathcal{X}};\mathcal{A}), (Pn)n≅(Qn)n(P_{n})_{n}\cong(Q_{n})_{n} if and only P=QP=Q, and (Pn)n⊲(Qn)n(P_{n})_{n}\vartriangleleft(Q_{n})_{n} if and only if P≪QP\ll Q, i.e., PP is absolutely continuous with respect to QQ. Hence asymptotic equivalence and contiguity can be thought of as asymptotic versions of equality and absolute continuity, respectively.

and thus (Yn)n⊲(Xn)n(Y_{n})_{n}\vartriangleleft(X_{n})_{n}

In particular, we will study random graphs of the following type. If pijp_{ij}, 1≤i<j≤n1\leq i<j\leq n, are given probabilities in , let G(n,{pij})G(n,\{p_{ij}\}) be the random graph on [n][n] where the edge ijij appears with probability pijp_{ij} and the indicators I_{ij}:=\boldsymbol{1}[\text{edgeijappears}], 1≤i<j≤n1\leq i<j\leq n, are independent. (We will later also consider an extension to random pijp_{ij}, see Section 2.)

Consequently, G(n,{pij})≅G(n,{pij′})G(n,\{p_{ij}\})\cong G(n,\{p_{ij}^{\prime}\}) if ∑i<j∣pij−pij′∣→0\sum_{i<j}|p_{ij}-p_{ij}^{\prime}|\to 0; a simple fact that has been used by many authors. It may be believed that this is essentially best possible, but, somewhat surprisingly, this is not so. In fact, by Corollary 2.12 below, provided pij≤0.9p_{ij}\leq 0.9, say, G(n,{pij})≅G(n,{pij′})G(n,\{p_{ij}\})\cong G(n,\{p_{ij}^{\prime}\}) if ∑i<j(pij−pij′)2/pij→0\sum_{i<j}(p_{ij}-p_{ij}^{\prime})^{2}/p_{ij}\to 0. (Moreover, Theorem 2.2(i) shows that this is best possible if, for example, pij′≤2pijp_{ij}^{\prime}\leq 2p_{ij}.)

For a particular case, suppose that pij′=pij+O(pij2)p_{ij}^{\prime}=p_{ij}+O(p_{ij}^{2}), Then Corollary 2.13 shows that G(n,{pij})≅G(n,{pij′})G(n,\{p_{ij}\})\cong G(n,\{p_{ij}^{\prime}\}) if ∑i<jpij3→0\sum_{i<j}p_{ij}^{3}\to 0, while (1.1) implies this only under the stronger condition ∑i<jpij2→0\sum_{i<j}p_{ij}^{2}\to 0. A typical case where this is an important improvement is when all pij=Θ(1/n)p_{ij}=\Theta(1/n) and ∣pij′−pij∣=Θ(1/n2)|p_{ij}^{\prime}-p_{ij}|=\Theta(1/n^{2}). See further the examples in Section 3.

To see that such an improvement of (1.1) might be possible at all, consider as an example the case when all pijp_{ij} are the same, so we consider the random graph G(n,p)G(n,p):

It follows, using Theorem 4.2 again, that G(n,p)≅G(n,p′)G(n,p)\cong G(n,p^{\prime}) if and only if α=0\alpha=0, i.e. N(p′−p)/Np→0N(p^{\prime}-p)/\sqrt{Np}\to 0, which is equivalent to ∑i<j(p′−p)2/p=N(p′−p)2/p→0\sum_{i<j}(p^{\prime}-p)^{2}/p=N(p^{\prime}-p)^{2}/p\to 0. For example, if p=1/np=1/n and p′=1−e−1/n=p−12n−2+O(n−3)p^{\prime}=1-e^{-1/n}=p-\frac{1}{2}n^{-2}+O(n^{-3}), then N(p′−p)2/p=O(1/n)N(p^{\prime}-p)^{2}/p=O(1/n), so G(n,p)≅G(n,p′)G(n,p)\cong G(n,p^{\prime}), but N∣p−p′∣→1/4N|p-p^{\prime}|\to 1/4, so (1.1) is not enough to show this.

We see that in this example, the trick to improve the simple and ’obvious’ edgewise coupling used in (1.1), is to first ignore the positions of the edges and couple their numbers only; this then is extended to a coupling of the random graphs by randomly reinserting the positions. Corollary 2.12(i) shows that couplings improving the simple edgewise coupling exist also when the edge probabilities are unequal, but in that case we do not know any explicit construction of such couplings.

We give the main results in Section 2 and a number of examples in Section 3; this includes an application to a recent result by van den Esker, van der Hofstad and Hooghiemstra (Example 3.6). Proofs are given in Section 5, after some preliminaries in Section 4.

We use the standard notations opo_{p} and OpO_{p}, see e.g. [14, Section 1.2], and we write whp (with high probability) for events with probability tending to 1 as n→∞{n\to\infty}.

There are also interesting examples of contiguity among random graphs of other types than G(n,{pij})G(n,\{p_{ij}\}). In particular, several different constructions of random regular graphs (or multigraphs) are known to yield distributions that are (mutually) contiguous but not asymptotically equivalent, see e.g. , [14, Section 9.5], . These examples are not covered by the present paper.

Results

We defined above the random graph G(n,p)G(n,{\mathbf{p}}), where p={pij}1≤i<j≤n{\mathbf{p}}=\{p_{ij}\}_{1\leq i<j\leq n} is a vector of probabilities. We extend the definition of G(n,p)G(n,{\mathbf{p}}) to the case when p{\mathbf{p}} is a random vector (with entries in ) by conditioning on p{\mathbf{p}}, i.e., given p={pij}{\mathbf{p}}=\{p_{ij}\}, the edge indicators IijI_{ij} are independent with Iij∼Be⁡(pij)I_{ij}\sim\operatorname{Be}(p_{ij}). Random graphs of this type have been studied in many papers, see for example Bollobás, Janson and Riordan 2007 and the further references given there.

We now state our main results on asymptotically equivalent and contiguity of such random graphs. Actually, the results have nothing to do with the graph structure and the way the indicator variables are indexed by pairs ijij. It therefore seems more natural to consider the more general situation of a (finite or infinite) sequence (Ii)1N(I_{i})_{1}^{N} of indicator variables. (An indicator variable is a random variable with values in {0,1}\{0,1\}, i.e. a random variable with a Bernoulli distribution Be⁡(p)\operatorname{Be}(p) for some p∈p\in.) The results for random graphs then follow by relabelling the indicators.

We define a function ρ:2→[0,∞)\rho:^{2}\to[0,\infty) in Definition 2.1, where we also give some equivalent (within constant factors) alternative formulas that often are more convenient. Since the results below are not affected by changing ρ\rho within constant factors, we could use any of these alternative formulas (and several other similar ones) as our definition. (The motivation for the definition comes in Lemma 4.3.)

We write x≍yx\asymp y (where x,y≥0x,y\geq 0) to denote that cy≤x≤Cycy\leq x\leq Cy for some positive constants c,Cc,C, i.e., that x=Θ(y)x=\Theta(y) (or, equivalently, x=O(y)x=O(y) and y=O(x)y=O(x)). Further, we use x∨yx\vee y and x∧yx\wedge y for the maximum and minimum, respectively, of xx and yy. We interpret 0/00/0 as 0.

Of course, the constant 0.90.9 here and below is arbitrary and could be replaced by any number <1<1.

together with the similar result with 1−p1-p and 1−q1-q. The second follows from x+y≍x∨yx+y\asymp x\vee y for x,y≥0x,y\geq 0 (used thrice). The third is equivalent to

which is easily verified, for example by assuming (by the symmetry p↦1−pp\mapsto 1-p, q↦1−qq\mapsto 1-q) that p≤1/2p\leq 1/2, in which case (2.6) easily reduces to p∨q≍p∨∣p−q∣p\vee q\asymp p\vee|p-q|. ∎

We state our results first for the simpler case of sequences of independent indicator variables with given (non-random) probabilities. The following theorem gives necessary and sufficient conditions for asymptotical equivalence and contiguity. (The asymptotical equivalence criterion follows by a simple and standard type of calculation with Hellinger distances, see the proof in Section 5 and, e.g., [16, p. 158], although we have not seen it stated in this form before. The contiguity criterion is a special case of a result by Oosterhoff and van Zwet 1979 for general sequences of independent variables.) The proofs of the theorems are given in Section 5.

Let 1≤N(n)≤∞1\leq N(n)\leq\infty and let Xn=(Ini)i=1N(n)X_{n}=(I_{ni})_{i=1}^{N(n)} and Xn′=(Ini′)i=1N(n)X^{\prime}_{n}=(I_{ni}^{\prime})_{i=1}^{N(n)} be finite or infinite random vectors consisting of independent indicator variables Ini∼Be⁡(pni)I_{ni}\sim\operatorname{Be}(p_{ni}) and Ini′∼Be⁡(pni′)I_{ni}^{\prime}\sim\operatorname{Be}(p_{ni}^{\prime}).

Xn≅Xn′X_{n}\cong X_{n}^{\prime} if and only if

Xn⊲Xn′X_{n}\vartriangleleft X_{n}^{\prime} if and only if

and, with qni:=1−pniq_{ni}:=1-p_{ni} and qni′:=1−pni′q_{ni}^{\prime}:=1-p_{ni}^{\prime},

By symmetry, Xn⊳Xn′X_{n}\vartriangleright X_{n}^{\prime} is equivalent to (2.8) and

and thus Xn⊲⊳Xn′X_{n}\vartriangleleft\vartriangleright X_{n}^{\prime} is equivalent to (2.8), (2.9) and (2.10).

Often pni≤0.9p_{ni}\leq 0.9 for all nn and ii, and then the second sum in (2.9) vanishes for C>10C>10 and can thus be omitted.

The condition (2.9) is only needed to take care of cases when pnip_{ni} and pni′p_{ni}^{\prime} (or qniq_{ni} and qni′q_{ni}^{\prime}, in case pnip_{ni} and pni′p_{ni}^{\prime} are close to 1) are not of the same order. If no such pnip_{ni} and pni′p_{ni}^{\prime} appear, which is the typical case, then (2.8) is thus enough.

We may rewrite (2.9) in several ways. For example, it is equivalent to (following the formulation in in a more general case): for every sequence Cn→∞C_{n}\to\infty,

It is also equivalent to: For every ε>0\varepsilon>0, there exist CC and n0n_{0} such that if n≥n0n\geq n_{0}, then

As pointed out by Oosterhoff and van Zwet 1979, (2.8) does not imply (2.9) in general. A simple counter example is provided by N(n)=nN(n)=n, pni=n−1p_{ni}=n^{-1}, pni′=n−2p_{ni}^{\prime}=n^{-2}. (On the other hand, it is easy to see, and also follows by the theorem, that (2.7) implies (2.9) and (2.10).)

We have stated Theorem 2.2 in terms of sequences of pairs of random vectors. It is possible (at least partly) to rephrase it in terms of estimates for a single pair (X,X′)(X,X^{\prime}), see Lemmas 5.1 and 5.2 below. Similar reformulations may be made for Theorem 2.9, but we leave these to the reader.

Let 1≤N(n)≤∞1\leq N(n)\leq\infty and suppose that pn={pni}{\mathbf{p}}_{n}=\{p_{ni}\} and pn′={pni′}{\mathbf{p}}_{n}^{\prime}=\{p_{ni}^{\prime}\} are random vectors in N(n)^{N(n)}. Let Xn=(Ini)i=1N(n)X_{n}=(I_{ni})_{i=1}^{N(n)} and Xn′=(Ini′)i=1N(n)X^{\prime}_{n}=(I_{ni}^{\prime})_{i=1}^{N(n)} be random vectors of indicator variables such that the conditioned random vectors (Xn∣pn)(X_{n}\mid{\mathbf{p}}_{n}) and (Xn′∣pn′)(X_{n}^{\prime}\mid{\mathbf{p}}_{n}^{\prime}) are sequences of independent indicator variables with (Ini∣pn)∼Be⁡(pni)(I_{ni}\mid{\mathbf{p}}_{n})\sim\operatorname{Be}(p_{ni}) and (Ini′∣pn′)∼Be⁡(pni′)(I_{ni}^{\prime}\mid{\mathbf{p}}_{n}^{\prime})\sim\operatorname{Be}(p_{ni}^{\prime}).

and, with qni:=1−pniq_{ni}:=1-p_{ni} and qni′:=1−pni′q_{ni}^{\prime}:=1-p_{ni}^{\prime}, for every ε>0\varepsilon>0,

then Xn⊲Xn′X_{n}\vartriangleleft X_{n}^{\prime}.

In analogy to (2.11), the condition (2.15) is equivalent to: For every sequence Cn→∞C_{n}\to\infty,

As said above, Theorems 2.2 and 2.9 apply immediately to random graphs G(n,p)G(n,{\mathbf{p}}). We state a version of Theorem 2.9 for this case, where we have added some simplifying assumptions. Recall that pijp_{ij} and pij′p_{ij}^{\prime} may (and typically do) depend on nn, although we do not show that in our notation.

Let, for each nn, p={pij}{\mathbf{p}}=\{p_{ij}\} and p′={pij′}{\mathbf{p}}^{\prime}=\{p_{ij}^{\prime}\} be random vectors of probabilities and suppose that whp max⁡i,jpij≤0.9\max_{i,j}p_{ij}\leq 0.9.

then G(n,p)≅G(n,p′)G(n,{\mathbf{p}})\cong G(n,{\mathbf{p}}^{\prime}).

then G(n,p)⊳G(n,p′)G(n,{\mathbf{p}})\vartriangleright G(n,{\mathbf{p}}^{\prime}).

If (2.18) holds, and further, for some constant c>0c>0, whp cpij≤pij′≤0.9cp_{ij}\leq p_{ij}^{\prime}\leq 0.9 for all i,ji,j, then G(n,p)⊲⊳G(n,p′)G(n,{\mathbf{p}})\vartriangleleft\vartriangleright G(n,{\mathbf{p}}^{\prime}).

We specialize further to an important case.

Let, for each nn, p={pij}{\mathbf{p}}=\{p_{ij}\} and p′={pij′}{\mathbf{p}}^{\prime}=\{p_{ij}^{\prime}\} be random vectors of probabilities and suppose that pij′=pij+O(pij2)p_{ij}^{\prime}=p_{ij}+O(p_{ij}^{2}).

If ∑i<jpij3=op(1)\sum_{i<j}p_{ij}^{3}=o_{p}(1), then G(n,p)≅G(n,p′)G(n,{\mathbf{p}})\cong G(n,{\mathbf{p}}^{\prime}).

If ∑i<jpij3=Op(1)\sum_{i<j}p_{ij}^{3}=O_{p}(1), and further, for some constant c>0c>0, whp max⁡i,jpij≤0.9\max_{i,j}p_{ij}\leq 0.9, max⁡i,jpij′≤0.9\max_{i,j}p_{ij}^{\prime}\leq 0.9 and pij′≥cpijp_{ij}^{\prime}\geq cp_{ij} for all i,ji,j, then G(n,p)⊲⊳G(n,p′)G(n,{\mathbf{p}})\vartriangleleft\vartriangleright G(n,{\mathbf{p}}^{\prime}).

Examples

Bollobás, Janson and Riordan 2007 study a general class of sparse random graphs G(n,κ)G(n,\kappa) which include many cases studied earlier by various authors. These random graphs are defined as G(n,{pij})G(n,\{p_{ij}\}) with

where κ:S×S→[0,∞)\kappa:{\mathcal{S}}\times{\mathcal{S}}\to[0,\infty) is a symmetric measurable function defined on some measurable space S{\mathcal{S}} and x1,…,xnx_{1},\dots,x_{n} is a random sequence of elements of S{\mathcal{S}}, not necessarily i.i.d. but such that the empirical distribution of x1,…,xnx_{1},\dots,x_{n} converges to a probability measure μ\mu on S{\mathcal{S}}; see for details. (Some further technical conditions are assumed in ; they are not relevant here.) Typically, whp κ(xi,xj)≤n\kappa(x_{i},x_{j})\leq n for all i,ji,j, and then pijp_{ij} equals the simpler p^ij\hat{p}_{ij}. Two natural variations, also treated in and used in various cases by various authors, are obtained by replacing (3.1) by

It was shown in that the same asymptotic results hold for these three versions for the properties studied there. We can now show that, under an extra condition, the three versions are asymptotically equivalent, and thus have the same asymptotic behaviour for any property. Indeed, Corollary 2.13 applies immediately and shows that if

We study some special cases in the following examples.

One common case of the construction in Example 3.1 uses x1,…,xnx_{1},\dots,x_{n} that are i.i.d. on S{\mathcal{S}} with distribution μ\mu. In this case, we show that the condition

implies (3.6) and thus asymptotic equivalence of the three versions. In particular, this holds if ∫S×Sκ(x,y)2 dμ(x) dμ(y)<∞\int_{{\mathcal{S}}\times{\mathcal{S}}}\kappa(x,y)^{2}\,d\mu(x)\,d\mu(y)<\infty.

In fact, if G(t):=μ×μ{(x,y):κ(x,y)>t)}=o(t−2)G(t):=\mu\times\mu\{(x,y):\kappa(x,y)>t)\}=o(t^{-2}), then

Similarly, we can easily shown that (3.7), and thus at least partial contiguity, follows from

In this case, given ε>0\varepsilon>0, there exists C1C_{1} such that

Another case of the construction in Example 3.1 uses S=(0,1]{\mathcal{S}}=(0,1] with μ\mu = Lebesgue measure and the deterministic xi=i/nx_{i}=i/n, i=1,…,ni=1,\dots,n. The homogeneous case κ(x,y)=c/(x∨y)\kappa(x,y)=c/(x\vee y) yielding p^ij=c/(i∨j)\hat{p}_{ij}=c/(i\vee j), where c>0c>0 is a constant, is particularly interesting and related to the CHKNS model, see Bollobás, Janson and Riordan 2007, Durrett 2003; Durrett 2007 and Riordan 2005 and the references given there.

In this case, ∑1≤i<j<∞pij3≤c3∑j≥2j⋅j−3<∞\sum_{1\leq i<j<\infty}p_{ij}^{3}\leq c^{3}\sum_{j\geq 2}j\cdot j^{-3}<\infty, and thus ∑i<jpij3=O(1)\sum_{i<j}p_{ij}^{3}=O(1); if we further for simplicity assume c<2c<2 and thus max⁡ijp^ij<1\max_{ij}\hat{p}_{ij}\allowbreak<1, then Corollary 2.13(iii) implies that G(n,pij(1))⊲⊳G(n,pij(2))⊲⊳G(n,pij(3))G(n,p_{ij}^{(1)})\vartriangleleft\vartriangleright G(n,p_{ij}^{(2)})\vartriangleleft\vartriangleright G(n,p_{ij}^{(3)}).

Note that in this case, p12(1)p_{12}^{(1)}, p12(2)p_{12}^{(2)} and p12(3)p_{12}^{(3)} are constant and different, which shows that the three random graphs are not asymptotically equivalent (for a trivial reason).

We have pij(3)=c/(i∨j+c)p_{ij}^{(3)}=c/(i\vee j+c); the same results hold for the further variation pij=c/(i∨j+d)p_{ij}=c/(i\vee j+d) for any d>c−2d>c-2.

A related case uses the same S=(0,1]{\mathcal{S}}=(0,1], μ\mu = Lebesgue measure and xi=i/nx_{i}=i/n, i=1,…,ni=1,\dots,n, as Example 3.3, now with the homogeneous κ(x,y)=c/xy\kappa(x,y)=c/\sqrt{xy} yielding p^ij=c/ij\hat{p}_{ij}=c/\sqrt{ij}; this case is a mean-field version of the preferential attachment model by Barabási and Albert , see Bollobás, Janson and Riordan 2007 and Riordan 2005 and the references given there.

Also in this case, ∑1≤i<j<∞pij3<∞\sum_{1\leq i<j<\infty}p_{ij}^{3}<\infty, and thus ∑i<jpij3=O(1)\sum_{i<j}p_{ij}^{3}=O(1) (in spite of the fact that (3.10) does not hold); if we further for simplicity assume c<2c<\sqrt{2}, and thus max⁡ijp^ij<1\max_{ij}\hat{p}_{ij}\allowbreak<1, we obtain the same results as in Example 3.3.

A common case of Example 3.1 is when κ(x,y)=ψ(x)ψ(y)\kappa(x,y)=\psi(x)\psi(y) for some function ψ:S→[0,∞)\psi:{\mathcal{S}}\to[0,\infty), see [4, Section 16.4] for discussion and references to previous papers.

In this case, ∑i<jκ(xi,xj)3≤(∑iψ(xi)3)2\sum_{i<j}\kappa(x_{i},x_{j})^{3}\leq\bigl(\sum_{i}\psi(x_{i})^{3}\bigr)^{2}, so (3.6) and (3.7) may be replaced by

If we combine this choice of κ\kappa with the i.i.d. choice of xix_{i} in Example 3.2, it is easily seen, arguing as in (3.9) but now with ∑i(ψ(xi)∧n1/2)\sum_{i}\bigl(\psi(x_{i})\wedge n^{1/2}\bigr), that

implies (3.6) and thus asymptotic equivalence of the three versions; in particular this holds if ∫ψ(x)2 dμ(x)<∞\int\psi(x)^{2}\,d\mu(x)<\infty. Similarly,

implies (3.7) and thus at least partial contiguity.

van den Esker, van der Hofstad and Hooghiemstra 2008+ study a minor variation of the construction in Example 3.5; they let Λ1,…,Λn\Lambda_{1},\dots,\Lambda_{n} be positive i.i.d. random variables with some fixed distribution and define pijp_{ij} by (in our notation) (3.1), (3.3), (3.4) or more generally (3.5) with

(This too can be seen as an instance of the general construction in Example 3.1, see [4, Section 16.4].)

Our results are stated for graphs with a deterministic number of vertices, but can be extended to graphs with random vertex set too by conditioning on the vertex set. One interesting such case is obtained from Example 3.1 by letting x1,…,xnx_{1},\dots,x_{n} be the points of a Poisson process on S{\mathcal{S}} with intensity λμ\lambda\mu, where λ>0\lambda>0 is our parameter and we consider asymptotics as λ→∞\lambda\to\infty; thus nn is random with the distribution Po⁡(λ)\operatorname{Po}(\lambda).

Conditioned on nn, we have the situation in Example 3.2. It follows, for example, that if (3.8) holds, then the random graphs defined in this way using (3.1), (3.3) and (3.4) are asymptotically equivalent; we omit the details.

Bollobás, Janson and Riordan 2007+ study a generalization of the model in Example 3.1 where small sets of edges are added at once, thus allowing a certain degree of clustering; more precisely, for every subgraph FF of the complete graph KnK_{n}, we have a certain probability of adding (the edges of) FF, and these events are independent for different FF. While this introduces dependencies between the edge indicators, the results of the present paper are still applicable to the sequence of indicators IFI_{F} describing the added sets of edges, and asymptotic equivalence or contiguity for two versions of this sequence obviously implies asymptotic equivalence or contiguity for the resulting random graphs too.

We leave the explicit statement of results in this case to the reader.

In this final example, let us return to the case of deterministic p={pij}{\mathbf{p}}=\{p_{ij}\} and let us change all pijp_{ij} proportionately to pij′:=(1+δn)pijp_{ij}^{\prime}:=(1+\delta_{n})p_{ij} for some δn\delta_{n}. Assume for simplicity that all pij≤0.9p_{ij}\leq 0.9 and that ∣δn∣≤1|\delta_{n}|\leq 1.

By Corollary 2.12(i), if δn2∑i<jpij→0\delta_{n}^{2}\sum_{i<j}p_{ij}\to 0, then G(n,p)≅G(n,p′)G(n,{\mathbf{p}})\cong G(n,{\mathbf{p}}^{\prime}). Further, by Corollary 2.12(iii), if δn2∑i<jpij=O(1)\delta_{n}^{2}\sum_{i<j}p_{ij}=O(1) and, for simplicity, δn→0\delta_{n}\to 0, then G(n,p)⊲⊳G(n,p′)G(n,{\mathbf{p}})\vartriangleleft\vartriangleright G(n,{\mathbf{p}}^{\prime}).

In fact, by (2.5), ρ(pij,pij′)≍δn2pij∧∣δn∣pij=δn2pij\rho(p_{ij},p_{ij}^{\prime})\asymp\delta_{n}^{2}p_{ij}\wedge|\delta_{n}|p_{ij}=\delta_{n}^{2}p_{ij}, and thus by Theorem 2.2 the conditions δn2∑i<jpij→0\delta_{n}^{2}\sum_{i<j}p_{ij}\to 0 and δn2∑i<jpij=O(1)\delta_{n}^{2}\sum_{i<j}p_{ij}=O(1) are necessary too for asymptotic equivalence and contiguity, respectively. (The necessity can also be checked by considering the total number of edges, as in the special case in Example 1.5.)

As in Example 3.9, necessity in Theorem 2.2 can in many cases where pij′≤pijp_{ij}^{\prime}\leq p_{ij} for all ii and jj (or conversely) be proved by considering the total numbers ∑iIni\sum_{i}I_{ni} and ∑iIni′\sum_{i}I_{ni}^{\prime}, but this method does not suffice in all cases. A simple counter example is given by N(n)=n4+n8N(n)=n^{4}+n^{8}, pni=n−1p_{ni}=n^{-1} for 1≤i≤n41\leq i\leq n^{4} and pni=n−3p_{ni}=n^{-3} for i>n4i>n^{4}, and pni′=pni−pni2p_{ni}^{\prime}=p_{ni}-p_{ni}^{2}; it is easily checked that then (2.7) and (2.8) do not hold, and thus we do not have asymptotic equivalence or even contiguity, but, using [2, Theorems 2.M and 1.C],

More on asymptotic equivalence and contiguity

We will use two metrics to measure the distance between probability distributions. We state some well-known definitions and facts, see e.g. [2, Appendix A.1] and [11, Chapter IV.1 and V.4a].

If PP and QQ are two probability measures on the same measurable space (X,A)({\mathcal{X}},\mathcal{A}), and RR is any σ\sigma-finite measure on (X,A)({\mathcal{X}},\mathcal{A}) such that P≪RP\ll R and Q≪RQ\ll R, define the total variation distance

(We can, at least symbolically, write (4.3) as H(P,Q):=∫XdP dQH(P,Q):=\int_{{\mathcal{X}}}\sqrt{dP\,dQ}.) Note that these quantities do not depend on the choice of RR. (We may thus take, e.g., R=P+QR=P+Q.)

taking the minimum over all couplings (X′,Y′)(X^{\prime},Y^{\prime}) of XX and YY.

Let XnX_{n} and YnY_{n} be random variables with values in Xn{\mathcal{X}}_{n}. Then the following are equivalent.

This too is well-known and easy: (i)  ⟺  \iff(ii) by (4.4) and Definition 1.1; (ii)  ⟺  \iff(iii) by (4.5); (iii)  ⟺  \iff(iv) by (4.2); (ii)  ⟺  \iff(v) by (4.6). ∎

We calculate the Hellinger distance and integral for two Bernoulli distributions. (This is the origin of our function ρ\rho in Definition 2.1.)

Use (4.2) with P=Be⁡(p)=pδ0+(1−p)δ1P=\operatorname{Be}(p)=p\delta_{0}+(1-p)\delta_{1}, Q=Be⁡(q)=qδ0+(1−q)δ1Q=\operatorname{Be}(q)=q\delta_{0}+(1-q)\delta_{1} and R=δ0+δ1R=\delta_{0}+\delta_{1}, together with the definition (2.1). ∎

Let 1≤N≤∞1\leq N\leq\infty and let, for i∈[N]i\in[N], PiP_{i} and QiQ_{i} be probability measures on the same measurable space (Xi,Ai)({\mathcal{X}}_{i},\mathcal{A}_{i}). If P=∏i=1NPiP=\prod_{i=1}^{N}P_{i} and Q=∏i=1NQiQ=\prod_{i=1}^{N}Q_{i}, then H(P,Q)=∏i=1NH(Pi,Qi)H(P,Q)=\prod_{i=1}^{N}H(P_{i},Q_{i}).

This is stated in, e.g., [11, Proposition IV.1.73], but for completeness we give the simple proof.

If N<∞N<\infty, the result is an immediate consequence of (4.3) and Fubini’s theorem, choosing e.g. Ri=Pi+QiR_{i}=P_{i}+Q_{i} and R=∏i=1NRiR=\prod_{i=1}^{N}R_{i}.

If N=∞N=\infty, let Fn\mathcal{F}_{n} be the σ\sigma-field on ∏i=1∞Xi\prod_{i=1}^{\infty}{\mathcal{X}}_{i} given by {A×∏n+1∞Xi:A∈∏1nAi}\{A\times\prod_{n+1}^{\infty}{\mathcal{X}}_{i}:A\in\prod_{1}^{n}\mathcal{A}_{i}\}, and let P‾n:=P∣Fn\overline{P}_{n}:=P|_{\mathcal{F}_{n}} and Q‾n:=Q∣Fn\overline{Q}_{n}:=Q|_{\mathcal{F}_{n}}. Then, using the finite case,

Proofs

and thus H(Xn,Xn′)→1  ⟺  ∑1N(n)ρ(pni,pni′)→0H(X_{n},X_{n}^{\prime})\to 1\iff\sum_{1}^{N(n)}\rho(p_{ni},p_{ni}^{\prime})\to 0, which yields the result by Theorem 4.2.

(ii): This is, in view of Lemma 4.3 and the equivalence of (2.9) and (2.11), a special case of [16, Theorem 1], to which we refer for a complete proof. Nevertheless, for completeness, we sketch a proof of the more important “if” direction.

First, we can by a simpler version of the argument in the proof of Theorem 2.9 below assume that pni≤C2pni′p_{ni}\leq C_{2}p_{ni}^{\prime} and qni≤C2qni′q_{ni}\leq C_{2}q_{ni}^{\prime} for some constant C2C_{2}. (We define pni′′′p_{ni}^{\prime\prime\prime} by (5.3) with pni′′:=pnip_{ni}^{\prime\prime}:=p_{ni} and use (5.4)–(5.5).) Under this assumption, if we let Pni:=L(Ini)=Be⁡(pni)P_{ni}:={\mathcal{L}}(I_{ni})=\operatorname{Be}(p_{ni}), Pn:=∏iPniP_{n}:=\prod_{i}P_{ni}, Pni′:=L(Ini′)=Be⁡(pni′)P_{ni}^{\prime}:={\mathcal{L}}(I_{ni}^{\prime})=\operatorname{Be}(p_{ni}^{\prime}), Pn′:=∏iPni′P_{n}^{\prime}:=\prod_{i}P_{ni}^{\prime}, we have by Fubini, using ∫(dPni/dPni′) dPni′=1\int(dP_{ni}/dP_{ni}^{\prime})\,dP_{ni}^{\prime}=1 and (2.2),

and thus for any sets AnA_{n}, by the Cauchy–Schwarz inequality,

and thus Pn⊲Pn′P_{n}\vartriangleleft P_{n}^{\prime}, which is the same as Xn⊲Xn′X_{n}\vartriangleleft X_{n}^{\prime}. ∎

We say that a finite or infinite random vectors of indicator variables X=(Ii)i=1NX=(I_{i})_{i=1}^{N} has distribution Be⁡(p)\operatorname{Be}({\mathbf{p}}), where p={pi}i=1N{\mathbf{p}}=\{p_{i}\}_{i=1}^{N} is a deterministic vector with elements in $,iftherandomvariables, if the random variablesI_{i}areindependentindicatorvariableswithare independent indicator variables withI_{i}\sim\operatorname{Be}(p_{i})$.

More generally, if p={pi}i=1N{\mathbf{p}}=\{p_{i}\}_{i=1}^{N} is a random vector with elements in $,with, withN\leq\infty,wesaythatrandomvectorsofindicatorvariables, we say that random vectors of indicator variablesX=(I_{i})_{i=1}^{N}hasdistributionhas distribution\operatorname{Be}({\mathbf{p}})iftheconditionedrandomvectorif the conditioned random vector(X\mid{\mathbf{p}})isasequenceofindependentindicatorvariableswithis a sequence of independent indicator variables with(I_{i}\mid{\mathbf{p}})\sim\operatorname{Be}(p_{i})$.

We next give two results comparing two random vectors with distributions Be⁡(p)\operatorname{Be}({\mathbf{p}}) and Be⁡(p′)\operatorname{Be}({\mathbf{p}}^{\prime}) with deterministic p{\mathbf{p}} and p′{\mathbf{p}}^{\prime}. The first result is easily seen to be equivalent to the “if” direction of Theorem 2.2(i), while the second is equivalent to a special case of the “if” direction of Theorem 2.2(ii).

On the other hand, (2.8) and (2.9) hold for these random vectors (since the sums in (2.9) vanish for any C≥C2C\geq C_{2}), and thus Theorem 2.2(ii) yields Xn⊲Xn′X_{n}\vartriangleleft X_{n}^{\prime}, which is a contradiction. ∎

for every measurable A⊆Xn={0,1}N(n)A\subseteq{\mathcal{X}}_{n}=\{0,1\}^{N(n)}, it follows that

(ii): Let (An)n(A_{n})_{n} be an arbitrary sequence measurable sets with An⊆Xn={0,1}N(n)A_{n}\subseteq{\mathcal{X}}_{n}=\{0,1\}^{N(n)} and let ε>0\varepsilon>0.

in the sequel we consider only n≥n0n\geq n_{0}.

Define pn′′={pni′′}i=1N(n){\mathbf{p}}_{n}^{\prime\prime}=\{p_{ni}^{\prime\prime}\}_{i=1}^{N(n)} by

Moreover, by the construction, with qni′′:=1−pni′′q_{ni}^{\prime\prime}:=1-p_{ni}^{\prime\prime},

Next, define pni′′′={pni′′′}i=1N(n)p_{ni}^{\prime\prime\prime}=\{p_{ni}^{\prime\prime\prime}\}_{i=1}^{N(n)} by

We can construct Xn′′′∼Be⁡(pn′′′)X_{n}^{\prime\prime\prime}\sim\operatorname{Be}({\mathbf{p}}_{n}^{\prime\prime\prime}) using maximal couplings of (Ini′′′∣pn′′′)(I_{ni}^{\prime\prime\prime}\mid{\mathbf{p}}_{n}^{\prime\prime\prime}) and (Ini′′∣pn′′)(I_{ni}^{\prime\prime}\mid{\mathbf{p}}_{n}^{\prime\prime}) so that, using (5.3) and (5.2),

Furthermore, by (5.3), pni′′′≤C2pni′p_{ni}^{\prime\prime\prime}\leq C_{2}p_{ni}^{\prime} and qni′′′:=1−pni′′′≤C2qni′q_{ni}^{\prime\prime\prime}:=1-p_{ni}^{\prime\prime\prime}\leq C_{2}q_{ni}^{\prime} and by (5.3) and (5.2),

In order to apply Theorem 2.9, we reorder {pij}i<j\{p_{ij}\}_{i<j} to {pni}i=1N(n)\{p_{ni}\}_{i=1}^{N(n)}; we do this without further comment. We also let qij:=1−pijq_{ij}:=1-p_{ij} and qij′:=1−pij′q_{ij}^{\prime}:=1-p_{ij}^{\prime}.

(i): By (2.5), whp ρ(pij,pij′)≤C0(pij−pij′)2/pij\rho(p_{ij},p_{ij}^{\prime})\leq C_{0}(p_{ij}-p_{ij}^{\prime})^{2}/p_{ij} for some C0C_{0}, and thus (2.17) implies (2.13), and the conclusion follows by Theorem 2.9(i).

(ii): Similarly, by (2.5) again, (2.18) implies (2.14). Moreover, for any C≥2C\geq 2,

(iii): The extra assumptions allow us to interchange p{\mathbf{p}} and p′{\mathbf{p}}^{\prime} in the assumptions. Hence (ii) yields both G(n,p)⊳G(n,p′)G(n,{\mathbf{p}})\vartriangleright G(n,{\mathbf{p}}^{\prime}) and G(n,p′)⊳G(n,p)G(n,{\mathbf{p}}^{\prime})\vartriangleright G(n,{\mathbf{p}}). ∎

An immediate consequence of Corollary 2.12, since now (pij−pij′)2/pij=O(pij3)(p_{ij}-p_{ij}^{\prime})^{2}/p_{ij}=O(p_{ij}^{3}); note also that the assumption in (i) implies max⁡i,jpij=op(1)\max_{i,j}p_{ij}=o_{p}(1) and thus max⁡i,jpij<0.9\max_{i,j}p_{ij}<0.9 whp. ∎

References