This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

A Partial Characterization of Robinsonian LpL^{p} Graphons

Teddy Mishura [email protected] Department of Mathematics, Kerr Hall South, Toronto Metropolitan University in Toronto, ON
Abstract

We present a characterization of Robinsonian LpL^{p} graphons for p>5p>5. Each LpL^{p} graphon ww is the limit object of a sequence of edge density-normalized simple graphs {Gn/Gn1}\{G_{n}/\|G_{n}\|_{1}\} under the cut distance δ\delta_{\Box}. A graphon ww is Robinson if it satisfies the Robinson property: if xyzx\leq y\leq z, then w(x,z)min{w(x,y),w(y,z)}w(x,z)\leq\min\{w(x,y),w(y,z)\}, and it is Robinsonian if δ(w,u)=0\delta_{\Box}(w,u)=0 for some Robinson uu. In previous work, the author and collaborators introduced a graphon parameter Λ\Lambda that recognizes the Robinson property, where Λ(w)=0\Lambda(w)=0 precisely when ww is Robinson. Using functional analytic arguments, we show here that for p>5p>5, the Robinsonian LpL^{p} graphons ww are precisely those that are the cut distance limit object of graphs GnG_{n} such that Λ(Gn/Gn1)0\Lambda(G_{n}/\|G_{n}\|_{1})\to 0.

keywords:
graphons , Robinson property , weak convergence , cut norm
MSC:
[2020] 47B47 , 05C50 , 54C08

[intoc]

1 Introduction

The analysis of structure in large networks has been an important question since the beginnings of computer science—many real world data structures are modeled as vast, complex graphs, and determining both local and global behavior is more desirable now than ever. In particular, many physical networks are dynamic in nature, growing in size and structure as new data is collected; equally as important in such changing systems is recognizing and characterizing emergent properties such as clustering or connectedness. The language of graph limit theory offers a powerful method to attack these problems in the form of LpL^{p} graphons, symmetric measurable functions from [0,1]2[0,1]^{2}\to{\mathbb{R}} with finite pp-norm p\|\cdot\|_{p}.

First introduced by Lovász and Szegedy [15], graphons are symmetric measurable functions from [0,1]2[0,1][0,1]^{2}\to[0,1], and can be viewed as the limit objects of sequences of graphs under the cut norm \|\cdot\|_{\Box}, where

w=supA,B[0,1]|A×Bw(x,y)𝑑x𝑑y|.\|w\|_{\square}=\sup_{A,B\subseteq[0,1]}\Bigg{|}\iint_{A\times B}w(x,y)dxdy\Bigg{|}.

For any labeled graph GG, the adjacency matrix AGA_{G} can be inserted into the unit square as a {0,1}\{0,1\}–valued step function wGw_{G}; graphons form the completion of these functions under ,\|\cdot\|_{\Box}, and we denote the space of graphons by 𝒲0{\mathcal{W}}_{0}. Given an unlabeled graph HH, where the adjacency matrix is ill defined, we compare its distance to a graphon using the cut distance δ\delta_{\Box}—the minimum difference in cut norm over all labelings of HH. This embedding of the adjacency matrix forces sequences of graphs {Gn}\{G_{n}\} with vanishing edge density to converge to the zero graphon, yielding a trivial result for many classes of networks, such as any graph family whose edges are governed by a power law. Borgs et al. resolved this issue in [1, 2] by instead considering the convergence of the graph sequence {Gn}\{G_{n}\} normalized by its edge density {Gn/Gn1}\{G_{n}/\|G_{n}\|_{1}\}. The limit functions of these sequences need only be bounded in pp-norm and are thus commonly known as LpL^{p} graphons, whose space is denoted 𝒲p.{\mathcal{W}}^{p}.

The main attraction of graph limit theory is the ability to answer combinatorial problems using analytic tools. Thus, when considering the problem of recognizing emergent structure and behaviors in a collection of large networks 𝒢\mathcal{G}, it is tempting to find an LpL^{p} graphon ww that is extremely close to 𝒢\mathcal{G} in cut norm and then test ww for the desired structure or behavior. While this approach is often quite powerful, it is not always the case that a combinatorial property of graphs translates to graphons (see for example [11]). Thus, recognizing and characterizing the properties that do extend from graphs to graphons is very useful for the analysis of growing networks, as determining the underlying graphon of a collection of networks can give insight into the collection as a whole.

We study here the Robinson property: an LpL^{p} graphon ww is Robinson if for all xyzx\leq y\leq z, then

w(x,z)min(w(x,y),w(y,z)).w(x,z)\leq\min(w(x,y),w(y,z)).

These functions were defined as diagonally increasing in [4] because they increase toward the main diagonal along a horizontal or vertical line, and are generalizations of Robinson matrices. Sometimes called R-matrices, Robinson matrices appear in the classic problem of seriation (see [14] for a comprehensive review of seriation) and their study is a field of much interest [3, 5, 6, 13]. Given a sequence of graphs {Gn}\{G_{n}\}, can one recognize if they are sampled from a Robinson graphon? It is simpler to tackle the labeled case first and then extend to the unlabeled case, so the problem becomes the following: Construct a function Φ:𝒲p[0,)\Phi:{\mathcal{W}}^{p}\to[0,\infty) that measures the Robinson property of a graphon. Such functions are called graphon parameters, and can be chosen to recognize desired structure in ww. A graphon parameter Φ\Phi is suitable for Robinson measurement if it is subadditive and satisfies the following three properties:

  • 1.

    (Recognition) Φ(w)=0\Phi(w)=0 if and only if ww is Robinson a.e.

  • 2.

    (Continuity) Φ\Phi is continuous with respect to the cut norm.

  • 3.

    (Stability) Given a graphon ww, there exists a Robinson graphon uu such that wuf(Φ(w))\|w-u\|_{\Box}\leq f(\Phi(w)), where f:[0,)[0,)f:[0,\infty)\to[0,\infty) is a nondecreasing function with limx0+f(x)=0.\lim_{x\to 0^{+}}f(x)=0.

Provided a sequence of associated graphons {wGn}\{w_{G_{n}}\} and a graphon w𝒲pw\in{\mathcal{W}}^{p} such that wGnw0\|w_{G_{n}}-w\|_{\Box}\to 0, recognition and stability combined with subadditivity guarantee that the Robinson graphon uu is close in cut norm to each wGnw_{G_{n}} by some measure of how Robinson ww is, while continuity ensures that the if the {wGn}\{w_{G_{n}}\} are Robinson, then so is ww. Thus, any such Φ\Phi will be able to recognize when a sequence of labeled graphs is sampled from a Robinson graphon or a graphon that is almost Robinson (see [10]). However, extending these results to unlabeled graphs is both difficult and dependent on the density of the graph sequence in question.

In this paper, we resolve this problem for unlabeled graphs that converge to LpL^{p} graphons for p>5p>5 in terms of the graphon parameter Λ\Lambda defined in [10]. All graphons in 𝒲p{\mathcal{W}}^{p} for p>5p>5 are “Λ\Lambda-close” to certain Robinson graphons in the cut norm—thus, if a graph sequence {Gn}\{G_{n}\} whose associated graphons {wGn}\{w_{G_{n}}\} grow increasingly Robinson converges to a graphon ww in cut distance, there must exist a sequence of Robinson graphons that converge to ww in cut distance. Through a combination of the specific properties of these Robinson approximations and functional analytic arguments, we showed that in this case there must exist some Robinson graphon uu such that δ(u/u1,w/w1)=0\delta_{\Box}(u/\|u\|_{1},w/\|w\|_{1})=0. Therefore, ww is Robinsonian precisely when it is the limit object of a sequence of simple graphs {Gn}\{G_{n}\} such that Λ(Gn/Gn1)0\Lambda(G_{n}/\|G_{n}\|_{1})\to 0.

2 Notation and background

In this section, we define basic notation used throughout the paper; afterwards, we present graph limit theory, stating required definitions alongside notable results, then transition into explanation of the Robinson property for both matrices and graphons.

All functions and sets discussed during this paper are assumed to be measurable. The Lebesgue measure of a set AA\subseteq{\mathbb{R}} is denoted by |A||A|. The notation p\|\cdot\|_{p} is used to represent the standard pp-norm; that is, for a function f:Xf:X\to{\mathbb{R}} and 1p<1\leq p<\infty,

fp:=(X|f|p)1p and f:=infa{f(x)a a.e.}.\|f\|_{p}:=\left(\int_{X}|f|^{p}\right)^{\frac{1}{p}}\mbox{ and }\|f\|_{\infty}:=\inf_{a\in{\mathbb{R}}}\{f(x)\leq a\mbox{ a.e.}\}.

Let XX be a normed vector space with inner product ,,\langle\cdot,\cdot\rangle, and let X={ϕ:X|ϕ continuous}X^{*}=\{\phi:X\to{\mathbb{R}}~{}|~{}\phi\text{ continuous}\} be the dual of X.X. If {μn}n1X\{\mu_{n}\}_{n\geq 1}\subset X^{*} and μX\mu\in X^{*}, we say that μn\mu_{n} converges weak to μ\mu, writing μnμ\mu_{n}\to\mu, if for all xXx\in X, we have that

x,μnx,μ.\langle x,\mu_{n}\rangle\to\langle x,\mu\rangle.

In particular, the function space Lp[0,1]2L^{p}[0,1]^{2} has Lq[0,1]2L^{q}[0,1]^{2} as its dual, where q=p/(p1)q=p/(p-1). Thus, for fnf_{n} to converge weak to ff in Lq[0,1]2L^{q}[0,1]^{2}, for all gLp[0,1]2g\in L^{p}[0,1]^{2}, it must be true that

([0,1]2|fn|p|g|p𝑑x𝑑y)1p([0,1]2|f|p|g|p𝑑x𝑑y)1p.\left(\iint_{[0,1]^{2}}|f_{n}|^{p}|g|^{p}dxdy\right)^{\frac{1}{p}}\to\left(\iint_{[0,1]^{2}}|f|^{p}|g|^{p}dxdy\right)^{\frac{1}{p}}.

2.1 Graph limit theory

A graph G=(V,E)G=(V,E) is a collection of vertices V(G)V(G) and edges E(G)E(G). In this paper, we will only consider simple graphs: those with neither loops nor multiple edges between vertices. We also draw a distinction between labeled and unlabeled graphs, where labeled graphs GG have the vertex set V(G)={1,,|V(G)|}V(G)=\{1,\ldots,|V(G)|\}. Any unlabeled graph HH can be labeled by uniquely assigning elements from {1,,|V(H)|}\{1,\ldots,|V(H)|\} to V(H)V(H).

Every labeled graph GG has an adjacency matrix AGA_{G} with (AG)ij=1(A_{G})_{ij}=1 if iji\sim j and 0 otherwise, where iji\sim j denotes that vertices ii and jj are connected by an edge. Furthermore, AGA_{G} can be embedded into [0,1]2[0,1]^{2} as a step function wGw_{G} as follows: Divide [0,1][0,1] into nn intervals I1,,InI_{1},...,I_{n} of equal measure, and for xIix\in I_{i} and yIjy\in I_{j}, let wG(x,y)=(AG)ijw_{G}(x,y)=(A_{G})_{ij}. Importantly, we note that both AGA_{G} and wGw_{G} are dependent on the labeling of V(G)V(G): see Figure 1 for an example. Unless otherwise noted, we always consider a labeled graph GG the same as its associated step function wGw_{G}. These functions wGw_{G} are members of a broader class of functions known as graphons, symmetric functions w:[0,1]2[0,1]w:[0,1]^{2}\to[0,1]. The space of graphons is denoted 𝒲0{\mathcal{W}}_{0}.

Refer to caption
Figure 1: Two isomorphic labeled graphs HH and HH^{\prime} with different associated step functions. Note that wHw_{H} and wHw_{H^{\prime}} are symmetric as functions on [0,1]2,[0,1]^{2}, not as matrices.

The foundational idea of graph limit theory is that a sequence of graphs {Gn}\{G_{n}\} can converge to a graphon ww. We shall view this entirely through an analytic framework, where we consider the convergence of the functions {wGn}\{w_{G_{n}}\} to the function ww in the famous cut norm \|\cdot\|_{\Box}, first defined by Frieze and Kannan [7]. The cut norm of w:[0,1]2w:[0,1]^{2}\to{\mathbb{R}} is given by

w=supS,T[0,1]|S×Tw(x,y)𝑑x𝑑y|.\|w\|_{\Box}=\sup_{S,T\subseteq[0,1]}\Bigg{|}\iint_{S\times T}w(x,y)dxdy\Bigg{|}.

For a labeled graph GG, we define G:=wG\|G\|_{\Box}:=\|w_{G}\|_{\Box}. A sequence of labeled graphs {Gn}\{G_{n}\} converges to a graphon ww if wGnw0\|w_{G_{n}}-w\|_{\Box}\to 0.

For a sequence of unlabeled graphs {Hn}\{H_{n}\} to converge to a graphon ww, something more is needed: Hnw\|H_{n}-w\|_{\Box} is only defined if HnH_{n} is labeled. Thus, the natural extension is to find the minimum cut norm difference over all possible labelings of HnH_{n}. The cut distance generalizes this idea to functions on [0,1]2[0,1]^{2}, where it is now the set [0,1][0,1] that must be relabeled. Let Φ\Phi be the set of all bijections φ:[0,1][0,1]\varphi:[0,1]\to[0,1] such that if A[0,1]A\subset[0,1], then |φ(A)|=|A||\varphi(A)|=|A|. We call such bijections measure preserving. The cut distance between two functions u,w:[0,1]2u,w:[0,1]^{2}\to{\mathbb{R}} is given by

δ(u,w)=infφΦuwφ,\delta_{\square}(u,w)=\inf_{\varphi\in\Phi}\|u-w^{\varphi}\|_{\square},

where wφ(x,y)=w(φ(x),φ(y)).w^{\varphi}(x,y)=w(\varphi(x),\varphi(y)). However, the distance δ\delta_{\square} is only a pseudometric—two nonidentical graphons can have cut distance 0. We fix this by identifying graphons of cut distance 0 from each other to get the set 𝒲~0\widetilde{{\mathcal{W}}}_{0} of unlabeled graphons. It is a classic result of graph limit theory that the metric space (𝒲~0,δ)(\widetilde{{\mathcal{W}}}_{0},\delta_{\Box}) is compact. Unless otherwise stated, the convergence of a graph sequence {Gn}\{G_{n}\} refers to the convergence of the unlabeled graphons {wGn}\{w_{G_{n}}\} in δ\delta_{\Box}. Furthermore, it is clear that for any graphon w𝒲0w\in{\mathcal{W}}_{0},

ww1wpw.\|w\|_{\square}\leq\|w\|_{1}\leq\|w\|_{p}\leq\|w\|_{\infty}.

Thus, if wGn10\|w_{G_{n}}\|_{1}\to 0 for some sequence of graphs {Gn}\{G_{n}\}, that sequence must trivially converge to 0. These sequences are called sparse; graph sequences that are not sparse are called dense. We note that wG1=|E(G)|/|V(G)|2\|w_{G}\|_{1}=|E(G)|/|V(G)|^{2} is the edge density of a graph, and so sparse graph sequences are exactly those with a subquadratic number of edges. Such graph sequences were shown in [1] to have nontrivial limits when normalized by their edge density—or, for the associated graphons, normalized by their L1L^{1} norm. Thus, for sparse graph sequences {Gn}\{G_{n}\} to converge, there must be some limit graphon ww such that

wGnwGn1w0 or δ(wGnwGn1,w)0\left\|\frac{w_{G_{n}}}{\|w_{G_{n}}\|_{1}}-w\right\|_{\Box}\to 0~{}\text{ or }~{}\delta_{\Box}\left(\frac{w_{G_{n}}}{\|w_{G_{n}}\|_{1}},w\right)\to 0

as opposed to {wGn}\{w_{G_{n}}\}. However, the limit objects of such sequences must be unbounded, and so we introduce LpL^{p} graphons: symmetric, measurable functions wLp([0,1]2)w\in L^{p}([0,1]^{2}). We refer to the space of LpL^{p} graphons as 𝒲p{\mathcal{W}}^{p}, noting that 𝒲0𝒲𝒲p{\mathcal{W}}_{0}\subset{\mathcal{W}}^{\infty}\subseteq{\mathcal{W}}^{p} for p1p\geq 1. This generalization of the notion of a graphon still brings with it powerful results from the dense theory: The metric space (𝒲~p,δ)(\widetilde{{\mathcal{W}}}^{p},\delta_{\Box}) has a compact unit ball.

As every graph can construct a graphon, so can every graphon construct a graph; sampling from a graphon refers to creating one of these random graphs. Given a graphon ww and a set S={x1,,xn}S=\{x_{1},...,x_{n}\} where each xi[0,1]x_{i}\in[0,1], create a random simple graph 𝒢(n,w){\mathcal{G}}(n,w) on V(G)={1,,n}V(G)=\{1,\ldots,n\} as follows: Connect vertices ii and jj with probability w(xi,xj)w(x_{i},x_{j}), making independent decisions for distinct pairs (i,j)(i,j), where i,j[n]i,j\in[n] and iji\neq j. We refer to 𝒢(n,w){\mathcal{G}}(n,w) as a w-random graph. Importantly, for large enough samples, the cut distance between a graphon and a sample is small with high probability.

It is often important to consider the average value of graphons over sets of the form A×BA\times B. Let w𝒲1w\in{\mathcal{W}}^{1} and let A,B[0,1]A,B\subseteq[0,1]. The cell average of ww over A×BA\times B is given by

w¯(A×B)=1|A×B|A×Bw𝑑x𝑑y.\overline{w}(A\times B)=\frac{1}{|A\times B|}\iint_{A\times B}w~{}dxdy.

Given a fixed LpL^{p} graphon w𝒲pw\in{\mathcal{W}}^{p} and a partition 𝒫={S1,,Sk}\mathcal{P}=\{S_{1},\ldots,S_{k}\} of [0,1][0,1] into nonempty measurable sets each of measure greater than 0, we define the step function w𝒫w_{\mathcal{P}} by

w𝒫(x,y)=1|Si×Sj|Si×Sjw(x,y)𝑑x𝑑y=w¯(Si×Sj)(xSi,ySj),w_{\mathcal{P}}(x,y)=\frac{1}{|S_{i}\times S_{j}|}\iint_{S_{i}\times S_{j}}w(x,y)dxdy=\overline{w}(S_{i}\times S_{j})\quad(x\in S_{i},y\in S_{j}),

where the operator ww𝒫w\mapsto w_{\mathcal{P}} is called the stepping operator. The stepping operator does not increase either the pp-norm or cut norm [2] of any graphon ww it is applied to; indeed, if w𝒲1w\in{\mathcal{W}}^{1}, 𝒫\mathcal{P} is any partition of [0,1][0,1] into a finite number of nonempty measurable sets, and p1p\geq 1, then

w𝒫pwp and w𝒫w.\|w_{\mathcal{P}}\|_{p}\leq\|w\|_{p}~{}\text{ and }~{}\|w_{\mathcal{P}}\|_{\square}\leq\|w\|_{\square}.

We note that this property is sometimes referred to as being contractive.

2.2 The Robinson property

The Robinson property for matrices was first introduced in the study of the classical seriation problem [16], whose objective is to order a set of items so that similar items are placed close to one another. A symmetric matrix A=[aij]A=[a_{ij}] is Robinson if

ijkaikmin{aij,ajk},i\leq j\leq k\implies a_{ik}\leq\min\{a_{ij},a_{jk}\},

and is Robinsonian if it becomes a Robinson matrix after simultaneous application of a permutation π\pi to its rows and columns. In that case, the permutation π\pi is called a Robinson ordering of AA. Graphs GG with Robinsonian adjacency matrices AGA_{G} are the well known unit interval graphs, sometimes also referred to as 11dimensional geometric graphs. A graph GG is a unit interval graph if and only if there exists an ordering \prec on V(G)V(G) such that

v,z,wV(G),vzw and vwzv and zw.\forall v,z,w\in V(G),~{}v\prec z\prec w\text{ and }v\sim w\implies z\sim v\text{ and }z\sim w.

This ordering \prec is precisely the Robinson ordering of AGA_{G}. Such Robinsonian matrices that are {0,1}\{0,1\}-valued are also studied under the name of the “consecutive 1s property” (see [8]).

When measuring large quantities of matrices or graphs for the Robinson property, such as in quality analysis of data, it is often useful to determine if these objects can be viewed as samples of some process ww. Were this the case, only the Robinson property of the process ww would need to be analyzed for quantitative results concerning the data as a whole. A graphon w𝒲pw\in{\mathcal{W}}^{p} is said to be Robinson if

xyzw(x,z)min{w(x,y),w(y,z)}.x\leq y\leq z\implies w(x,z)\leq\min\{w(x,y),w(y,z)\}.

Robinson graphons were first considered in [4] under the name diagonally increasing graphons. We call a graphon Robinson almost everywhere, or Robinson a.e. for short, if it is equal a.e. to a Robinson graphon. A graphon w𝒲pw\in{\mathcal{W}}^{p} is called Robinsonian if there exists a Robinson LpL^{p} graphon u𝒲pu\in{\mathcal{W}}^{p} such that δ(w,u)=0\delta_{\Box}(w,u)=0. This natural extension of the Robinson property to the world of graphons results in the addition of significant complexity; thus, much effort has been focused on measuring and approximating the Robinson property—if a graphon is close enough to being Robinson, its samples should share the same behavior.

2.3 Sampling from Robinson graphons

In [4], the authors introduced a (labeled) graphon parameter Γ\Gamma to estimate the Robinsonicity of a given graphon. We recall its definition below.

Definition 2.1 (The parameter Γ\Gamma).

For w𝒲0w\in{\mathcal{W}}_{0} and a measurable subset AA of [0,1][0,1], let

Γ(w,A):=y<z[xA[0,y](w(x,z)w(x,y))𝑑x]+𝑑y𝑑z+y<z[xA[z,1](w(x,y)w(x,z))𝑑x]+𝑑y𝑑z,\Gamma(w,A):=\iint_{y<z}\left[\int_{x\in A\cap[0,y]}(w(x,z)-w(x,y))dx\right]_{+}dydz+\iint_{y<z}\left[\int_{x\in A\cap[z,1]}(w(x,y)-w(x,z))dx\right]_{+}dydz,

with [x]+=max(x,0)[x]_{+}=\max(x,0). Moreover, let Γ(w):=supA𝒜Γ(w,A),\Gamma(w):=\sup_{A\in\mathcal{A}}\Gamma(w,A), where the supremum is taken over all measurable subsets of [0,1][0,1]. If GG is a simple, labeled graph, then we define Γ(G):=Γ(wG)\Gamma(G):=\Gamma(w_{G}).

It was shown in [4, 9] that Γ\Gamma is suitable for Robinson measurement of graphons; that is, it possesses recognition, continuity, and stability. Thus, as Γ\Gamma is continuous with respect to graphons in the cut norm, we hope that we can pass that continuity down to graph sequences. Namely, given a convergent sequence of dense graphs {Gn}\{G_{n}\} and a graphon w𝒲0w\in{\mathcal{W}}_{0} where δ(Gn,w)0\delta_{\Box}(G_{n},w)\to 0, we desire a result of the form Γ(Gn)Γ(w)\Gamma(G_{n})\to\Gamma(w). The authors showed a form of this independent of vertex labeling, which we recall below.

Proposition 2.2 (Corollary 6.5, [4]).

Let {Gn}\{G_{n}\} be a sequence of simple graphs with |V(Gn)||V(G_{n})|\to\infty and let w𝒲0w\in{\mathcal{W}}_{0} such that δ(Gn,w)0\delta_{\Box}(G_{n},w)\to 0. If Φn\Phi_{n} is the set of all vertex labelings of GnG_{n} and Φ\Phi is the set of all measure preserving bijections of [0,1][0,1], then

minϕnΦnΓ(Gnϕn)infψΦΓ(wψ).\min_{\phi_{n}\in\Phi_{n}}\Gamma(G_{n}^{\phi_{n}})\to\inf_{\psi\in\Phi}\Gamma(w^{\psi}).

In the case where the lefthand side of this limit tends to 0, it is tempting to conclude that the graphs {Gn}\{G_{n}\} are similar to random models sampled from a Robinsonian graphon; however, this does not follow from any of the previous results. Indeed, while infψΦΓ(wψ)=0\inf_{\psi\in\Phi}\Gamma(w^{\psi})=0, this would only imply that the {Gn}\{G_{n}\} are arbitrarily close in δ\delta_{\Box} to Robinson graphons {un}\{u_{n}\}—and a δ\delta_{\Box}–equivalence class may contain more than one Robinson graphon. The proof of stability of Γ\Gamma in [9] combined with a functional analytic argument was necessary to show that graphs {Gn}\{G_{n}\} as mentioned above are indeed sampled from a unique Robinsonian graphon.

Theorem 2.3 (Theorem 5.3, [9]).

Let {Gn}\{G_{n}\} be a sequence of simple graphs converging to a graphon w𝒲0w\in{\mathcal{W}}_{0} in δ\delta_{\Box}. Then ww is Robinsonian if and only if there exist vertex labelings {ϕn}n1\{\phi_{n}\}_{n\geq 1} such that Γ(Gnϕn)0\Gamma(G_{n}^{\phi_{n}})\to 0.

3 Sparse Robinsonian graph sequences

In this section we state and prove our main result regarding the characterization of Robinsonian LpL^{p} graphons. Specifically, we shall prove a generalization of Theorem 2.3 where w𝒲pw\in\mathcal{W}^{p} and where the simple graphs {Gn}n1\{G_{n}\}_{n\geq 1} that converge to ww do so by satisfying δ(Gn/Gn1,w/w1)0.\delta_{\Box}(G_{n}/\|G_{n}\|_{1},w/\|w\|_{1})\to 0. However, while Γ\Gamma is a suitable measurement of the Robinson property for graphons in the space 𝒲0{\mathcal{W}}_{0}, the techniques used to show continuity and stability of Γ\Gamma on 𝒲0{\mathcal{W}}_{0} are unable to be extended to the more general space 𝒲p{\mathcal{W}}^{p}. Thus, we shall instead prove our characterization result in terms of the graphon parameter Λ\Lambda, which was introduced in [10] to address this discrepancy between 𝒲0{\mathcal{W}}_{0} and 𝒲p{\mathcal{W}}^{p}. We recall its definition and relevant properties here.

Definition 3.1 (The parameter Λ\Lambda).

Let w𝒲1w\in\mathcal{W}^{1}. Define

Λ(w)=12supABC,|A|=|B|=|C|[A×Cw𝑑x𝑑yB×Cw𝑑x𝑑y]+12supXYZ,|X|=|Y|=|Z|[X×Zw𝑑x𝑑yX×Yw𝑑x𝑑y]\Lambda(w)=\frac{1}{2}\sup_{\begin{subarray}{c}A\leq B\leq C,\\ |A|=|B|=|C|\end{subarray}}\bigg{[}\iint_{A\times C}w~{}dxdy-\iint_{B\times C}w~{}dxdy\bigg{]}\\ +\frac{1}{2}\sup_{\begin{subarray}{c}X\leq Y\leq Z,\\ |X|=|Y|=|Z|\end{subarray}}\bigg{[}\iint_{X\times Z}w~{}dxdy-\iint_{X\times Y}w~{}dxdy\bigg{]}

where A,B,CA,B,C and X,Y,ZX,Y,Z are measurable subsets of [0,1][0,1] and ABA\leq B if aba\leq b for all aAa\in A and bBb\in B. If GG is a labeled graph, we define Λ(G):=Λ(wG)\Lambda(G):=\Lambda(w_{G}).

Proposition 3.2.

Let 1p1\leq p\leq\infty. Suppose w,u𝒲pw,u\in{\mathcal{W}}^{p}. Then we have

  • (i)

    (Recognition) Λ\Lambda characterizes Robinson LpL^{p}-graphons, i.e., ww is Robinson a.e. if and only if Λ(w)=0\Lambda(w)=0.

  • (ii)

    (Continuity) Λ\Lambda is continuous with respect to cut norm, i.e., |Λ(w)Λ(u)|2wu|\Lambda(w)-\Lambda(u)|\leq 2\|w-u\|_{\Box}.

  • (iii)

    (Stability) For every w𝒲pw\in{\mathcal{W}}^{p} where p>5p>5, there exists a Robinson graphon u𝒲pu\in{\mathcal{W}}^{p} satisfying

    uw78Λ(w)p55p5.\|u-w\|_{\Box}\leq 78\Lambda(w)^{\frac{p-5}{5p-5}}.

Informally, Λ\Lambda can be thought of as extending the Robinson property from pointwise comparison to comparison of blocks. We make special note of item (iii)(iii), as this result is constructive—the Robinson graphon uu is the α\alpha-Robinson approximation of ww for some chosen α\alpha, which is a generalized version of the Robinson construction introduced in [9].

Definition 3.3 (α\alpha-Robinson approximation for graphons).

Let p1p\geq 1, and fix a parameter 0<α<10<\alpha<1. Given a graphon w𝒲pw\in{\mathcal{W}}^{p}, the α\alpha-Robinson approximation RwαR_{w}^{\alpha} of ww is the symmetric function defined as follows: For all 0xy10\leq x\leq y\leq 1, define

Rwα(x,y)=sup{1|S×T|S×Tw:S×T[0,x]×[y,1],|S|=|T|=α},R_{w}^{\alpha}(x,y)=\sup\left\{\frac{1}{|S\times T|}\iint_{S\times T}w\,:\,S\times T\subseteq[0,x]\times[y,1],\,|S|=|T|=\alpha\right\},

taking the convention that sup=0\sup\emptyset=0. Moreover, we set Rwα=wR_{w}^{\alpha}=w if α=0\alpha=0 and ww is Robinson.

The goal is thus to show for ww and {Gn}\{G_{n}\} where δ(Gn/Gn1,w/w1)0\delta_{\Box}(G_{n}/\|G_{n}\|_{1},w/\|w\|_{1})\to 0 that ww is Robinsonian if and only if there exists a sequence of vertex labelings {ϕn}n1\{\phi_{n}\}_{n\geq 1} such that Λ(Gnϕn/Gn1)0.\Lambda(G_{n}^{\phi_{n}}/\|G_{n}\|_{1})\to 0. While the forward direction is rather direct, the reverse direction requires some heavier machinery. There are two ways to show that ww is Robinsonian: Firstly, find a measure preserving bijection ψ\psi such that Λ(wψ)=0,\Lambda(w^{\psi})=0, and secondly, find a Robinson graphon u𝒲pu\in{\mathcal{W}}^{p} such that δ(u/u1,w/w1)=0\delta_{\Box}(u/\|u\|_{1},w/\|w\|_{1})=0.

Finding such a ψ\psi is not obvious—indeed, while it is tempting to think of ψ\psi as the limit object of the {ϕn}\{\phi_{n}\}, it is not clear that such an object need exist. Proposition 3.2 (iii) provides an avenue of attack for the latter option of showing that ww is Robinsonian, as we can use these Robinson approximations of the associated graphons wGnw_{G_{n}} to find a Robinson graphon uu of cut distance 0 from ww. We proceed with this plan below, beginning with a technical lemma showing that for certain partitions 𝒫\mathcal{P} of [0,1][0,1], the step graphon w𝒫w_{\mathcal{P}} is Robinsonian.

Lemma 3.4.

Let η>0\eta>0 and let 𝒫={P1,,Pm}\mathcal{P}=\{P_{1},\ldots,P_{m}\} be a partition of [0,1][0,1] into measurable sets each of measure at least η.\eta. If w𝒲pw\in{\mathcal{W}}^{p} with p>5p>5 is a graphon such that there exists a sequence of graphs {Gn}n1\{G_{n}\}_{n\geq 1} and vertex labelings {ϕn}n1\{\phi_{n}\}_{n\geq 1} where

δ(GnGn1,ww1)0andΛ(GnϕnGn1)0,\delta_{\Box}\left(\frac{G_{n}}{\|G_{n}\|_{1}},\frac{w}{\|w\|_{1}}\right)\to 0\quad\quad\quad\text{and}\quad\quad\quad\Lambda\left(\frac{G_{n}^{\phi_{n}}}{\|G_{n}\|_{1}}\right)\to 0,

then the step graphon w𝒫w_{\mathcal{P}} is Robinsonian.

Proof.

Fix η>0\eta>0, let 𝒫={P1,,Pm}\mathcal{P}=\{P_{1},\ldots,P_{m}\} be a partition of [0,1][0,1] into measurable sets each of measure at least η\eta, and let wn:=wGnϕn/wGn1w_{n}:=w_{G_{n}^{\phi_{n}}}/\|w_{G_{n}}\|_{1} and w:=w/w1w^{\prime}:=w/\|w\|_{1} for ease of reading. By the assumption of convergence of wnww_{n}\to w^{\prime} in cut distance, we have that δ((wn)𝒫,(w)𝒫)0\delta_{\Box}((w_{n})_{\mathcal{P}},(w^{\prime})_{\mathcal{P}})\to 0. Furthermore, by Proposition 3.2 (iii), for every nn\in{\mathbb{N}}, there must exist a Robinson graphon un𝒲pu_{n}\in{\mathcal{W}}^{p} such that

un(wn)𝒫78Λ((wn)𝒫)p55p5.\|u_{n}-(w_{n})_{\mathcal{P}}\|_{\Box}\leq 78\Lambda((w_{n})_{\mathcal{P}})^{\frac{p-5}{5p-5}}.

We note that as Λ(wn)0\Lambda(w_{n})\to 0, it must be the case that Λ((wn)𝒫)0\Lambda((w_{n})_{\mathcal{P}})\to 0 as well. This further implies that δ(un,(w)𝒫)0\delta_{\Box}(u_{n},(w^{\prime})_{\mathcal{P}})\to 0 as nn\to\infty, as

δ(un,(w)𝒫)δ(un,(wn)𝒫)+δ((wn)𝒫,(w)𝒫)un(wn)𝒫+δ((wn)𝒫,(w)𝒫).\delta_{\Box}(u_{n},(w^{\prime})_{\mathcal{P}})\leq\delta_{\Box}(u_{n},(w_{n})_{\mathcal{P}})+\delta_{\Box}((w_{n})_{\mathcal{P}},(w^{\prime})_{\mathcal{P}})\leq\|u_{n}-(w_{n})_{\mathcal{P}}\|_{\Box}+\delta_{\Box}((w_{n})_{\mathcal{P}},(w^{\prime})_{\mathcal{P}}).

Before we continue, we shall show that the functions unu_{n} are uniformly bounded above; indeed, note that

|un(x,y)|=|R(wn)𝒫α(x,y)|\displaystyle|u_{n}(x,y)|=|R^{\alpha}_{(w_{n})_{\mathcal{P}}}(x,y)| =supA,B[0,x]×[y,1]|A|=|B|=α|1α2A×B(wn)𝒫|\displaystyle=\sup_{\begin{subarray}{c}A,B\subseteq[0,x]\times[y,1]\\ |A|=|B|=\alpha\end{subarray}}\left|\frac{1}{\alpha^{2}}\iint_{A\times B}(w_{n})_{\mathcal{P}}\right|
(wn)𝒫=sup1i,jm|1|Pi×Pj|Pi×Pjwn|\displaystyle\leq\|(w_{n})_{\mathcal{P}}\|_{\infty}=\sup_{1\leq i,j\leq m}\left|\frac{1}{|P_{i}\times P_{j}|}\iint_{P_{i}\times P_{j}}w_{n}\right|
sup1i,jm1|Pi×Pj|wn1𝟙Pi×Pj\displaystyle\leq\sup_{1\leq i,j\leq m}\frac{1}{|P_{i}\times P_{j}|}\|w_{n}\|_{1}\|\mathbbm{1}_{P_{i}\times P_{j}}\|_{\infty}
=sup1i,jm1|Pi×Pj|1η2.\displaystyle=\sup_{1\leq i,j\leq m}\frac{1}{|P_{i}\times P_{j}|}\leq\frac{1}{\eta^{2}}.

Now, for all nn\in{\mathbb{N}}, the functions unu_{n} are elements of the space B1/η2(L[0,1]2)={fL[0,1]2:f1η2}B_{1/\eta^{2}}(L^{\infty}[0,1]^{2})=\{f\in L^{\infty}[0,1]^{2}:\|f\|_{\infty}\leq\frac{1}{\eta^{2}}\}, and because the Banach space L[0,1]2L^{\infty}[0,1]^{2} is isometrically isomorphic to the Banach space dual of L1[0,1]2L^{1}[0,1]^{2}, it is possible to equip B1/η2(L[0,1]2)B_{1/\eta^{2}}(L^{\infty}[0,1]^{2}) with the weak topology induced by this duality. By the Banach-Alaglou theorem, B1/η2(L[0,1]2)B_{1/\eta^{2}}(L^{\infty}[0,1]^{2}) is compact in this topology; additionally, as L1[0,1]2L^{1}[0,1]^{2} is separable, B1/η2(L[0,1]2)B_{1/\eta^{2}}(L^{\infty}[0,1]^{2}) is metrizable in the weak topology and thus is sequentially compact as well. Therefore, {un}n\{u_{n}\}_{n\in{\mathbb{N}}} has a weak convergent subsequence.

Without loss of generality, we can thus assume that {un}\{u_{n}\} converges to some zB1/η2(L[0,1]2)z\in B_{1/\eta^{2}}(L^{\infty}[0,1]^{2}) in the weak topology; i.e. for every hL1[0,1]2h\in L^{1}[0,1]^{2} we have that [0,1]2unh[0,1]2zh\iint_{[0,1]^{2}}u_{n}h\to\iint_{[0,1]^{2}}zh as n.n\to\infty. In particular, for all measurable sets S,T[0,1]S,T\subseteq[0,1], we have that

S×TunS×Tz as n.\iint_{S\times T}u_{n}\to\iint_{S\times T}z~{}\text{ as }n\to\infty. (1)

By (1), we get that for any partition 𝒫\mathcal{P} of [0,1][0,1] whose elements are all bounded below in measure, (un)𝒫(u_{n})_{\mathcal{P}} converges pointwise to z𝒫z_{\mathcal{P}}. Furthermore, as every (un)𝒫(u_{n})_{\mathcal{P}} is Robinson, it must be the case that z𝒫z_{\mathcal{P}} is Robinson as well. We can also show that zz is Robinson a.e. by letting 𝒫n\mathcal{P}_{n} be any partition of [0,1][0,1] into sets of exactly measure 1/n1/n and noting that zz𝒫n10\|z-z_{\mathcal{P}_{n}}\|_{1}\to 0. Additionally, again using (1), it is true that δ(un,z)0\delta_{\Box}(u_{n},z)\to 0 as nn\to\infty. Therefore,

δ((w)𝒫,z)δ((w)𝒫,un)+δ(un,z)0\delta_{\Box}((w^{\prime})_{\mathcal{P}},z)\leq\delta_{\Box}((w^{\prime})_{\mathcal{P}},u_{n})+\delta_{\Box}(u_{n},z)\to 0

as nn\to\infty, showing that (w)𝒫(w^{\prime})_{\mathcal{P}} is Robinsonian and thus that w𝒫w_{\mathcal{P}} is Robinson as well, proving the claim. ∎

With this lemma at our disposal, we now state and prove our main result about the partial characterization of Robinsonian LpL^{p} graphons.

Theorem 3.5.

Let {Gn}\{G_{n}\} be a sequence of simple labeled graphs and let w𝒲pw\in{\mathcal{W}}^{p} with p>5p>5 such that δ(Gn/Gn1,w/w1)0\delta_{\Box}(G_{n}/\|G_{n}\|_{1},w/\|w\|_{1})\to 0. Then, ww is Robinsonian if and only if Λ(Gnϕn/Gn1)0\Lambda(G^{\phi_{n}}_{n}/\|G_{n}\|_{1})\to 0 for some vertex labelings {ϕn}n1\{\phi_{n}\}_{n\geq 1}.

Proof.

We start with the forward direction; assume that ww is Robinsonian. Then there exists a Robinson graphon u𝒲pu\in{\mathcal{W}}^{p} such that δ(w/w1,u/u1)=0,\delta_{\Box}(w/\|w\|_{1},u/\|u\|_{1})=0, from which it follows by using the triangle inequality that

δ(GnGn1,uu1)δ(GnGn1,ww1)+δ(ww1,uu1)0.\delta_{\Box}\left(\frac{G_{n}}{\|G_{n}\|_{1}},\frac{u}{\|u\|_{1}}\right)\leq\delta_{\Box}\left(\frac{G_{n}}{\|G_{n}\|_{1}},\frac{w}{\|w\|_{1}}\right)+\delta_{\Box}\left(\frac{w}{\|w\|_{1}},\frac{u}{\|u\|_{1}}\right)\to 0.

So there must exist a sequence of vertex labelings {ϕn}n1\{\phi_{n}\}_{n\geq 1} such that Gnϕn/Gn1u/u10\|G_{n}^{\phi_{n}}/\|G_{n}\|_{1}-u/\|u\|_{1}\|_{\Box}\to 0 as nn tends to infinity. Noting that Λ(u/u1)=0\Lambda(u/\|u\|_{1})=0 because uu is Robinson, this implies that

Λ(GnϕnGn1)=|Λ(GnϕnGn1)Λ(uu1)|GnϕnGn1uu10\Lambda\left(\frac{G_{n}^{\phi_{n}}}{\|G_{n}\|_{1}}\right)=\left|\Lambda\left(\frac{G_{n}^{\phi_{n}}}{\|G_{n}\|_{1}}\right)-\Lambda\left(\frac{u}{\|u\|_{1}}\right)\right|\leq\left\|\frac{G_{n}^{\phi_{n}}}{\|G_{n}\|_{1}}-\frac{u}{\|u\|_{1}}\right\|_{\Box}\to 0

as required, proving the forward direction.

For the reverse direction, assume that each GnG_{n} has a vertex labeling ϕn\phi_{n} such that Λ(Gnϕn)0\Lambda(G_{n}^{\phi_{n}})\to 0; note also that δ(Gn/Gn1,w/w1)0\delta_{\Box}(G_{n}/\|G_{n}\|_{1},w/\|w\|_{1})\to 0 by assumption. For all nn\in{\mathbb{N}}, let 𝒫n\mathcal{P}_{n} be any partition of [0,1][0,1] into measurable sets each of measure exactly 1/n1/n, and consider the graphons {w𝒫n}n1.\{w_{\mathcal{P}_{n}}\}_{n\geq 1}. By Lemma 3.4, for all nn\in{\mathbb{N}}, there exists a sequence {zn}n\{z_{n}\}_{n\in{\mathbb{N}}} of Robinson graphons and a sequence {ψn}\{\psi_{n}\} of measure preserving bijections of [0,1][0,1] such that

w𝒫nψnzn=0.\|w_{\mathcal{P}_{n}}^{\psi_{n}}-z_{n}\|_{\Box}=0.

We consider now the sequence of LpL^{p} graphons {w𝒫nψn}n\{w^{\psi_{n}}_{\mathcal{P}_{n}}\}_{n\in{\mathbb{N}}}. These are elements of a bounded ball of radius wp\|w\|_{p} in the space Lp[0,1]2L^{p}[0,1]^{2}, which is compact and sequentially compact in the weak topology induced by Lq[0,1]2L^{q}[0,1]^{2}. Thus, {w𝒫nψn}n\{w^{\psi_{n}}_{\mathcal{P}_{n}}\}_{n\in{\mathbb{N}}} has a weak convergent subsequence in this topology and we can say without loss of generality that there exists some yLp[0,1]2y\in L^{p}[0,1]^{2} such that ypwp\|y\|_{p}\leq\|w\|_{p} and that for all measurable S,T[0,1]S,T\subseteq[0,1] we have

S×T(w𝒫nψn)pS×Typ as n.\iint_{S\times T}\left(w^{\psi_{n}}_{\mathcal{P}_{n}}\right)^{p}\to\iint_{S\times T}y^{p}~{}\text{ as }n\to\infty.

As every w𝒫nψnw^{\psi_{n}}_{\mathcal{P}_{n}} is Robinson a.e., it must be that yy is Robinson a.e. as well. We can therefore show that

δ(w,y)δ(w,w𝒫n)+δ(y,w𝒫n)0\delta_{\Box}(w,y)\leq\delta_{\Box}(w,w_{\mathcal{P}_{n}})+\delta_{\Box}(y,w_{\mathcal{P}_{n}})\to 0

as nn\to\infty, demonstrating that ww is Robinsonian and finishing the proof. ∎

4 Further Directions

Several natural questions present themselves as extensions to this work, the most immediate of which is whether the statement of Theorem 3.5 holds for 1p51\leq p\leq 5. As the requirement of p>5p>5 is an artifact of the proof method for Proposition 3.2 (iii), an obvious method of attack for this question would be to improve the techniques used to show this result. It is additionally interesting to consider whether this characterization result could be shown for a more general concept of sparse graphons, such as the kk-shapes introduced in [12], especially as these objects are defined independently of the vertex labeling of the graphs they are associated with. It is likely that a new parameter would need to be defined to quantify these approximation results in this extended setting, and it is of significant interest to the author to see whether similar claims could hold true for these more general objects. Finally, it is natural to wonder whether there are other graph properties for which such characterization results hold—monotonicity is a likely candidate due to its similarity to Robinsonicity, and categorizing other such properties seems a fruitful direction for future work.

5 Acknowledgements

The author would like to acknowledge Dr. Mahya Ghandehari, Charli Klein, and Sarah Tillman for their help in editing this paper and the Department of Mathematics at Toronto Metropolitan University for its continued support throughout the process of this research. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

References

  • Borgs et al. [2018] C. Borgs, J. Chayes, H. Cohn, and Y. Zhao. An LpL^{p} theory of sparse graph convergence II: LD convergence, quotients and right convergence. The Annals of Probability, 46(1), Jan 2018. ISSN 0091-1798.
  • Borgs et al. [2019] C. Borgs, J. Chayes, H. Cohn, and Y. Zhao. An LpL^{p} theory of sparse graph convergence I: Limits, sparse random graph models, and power law distributions. Transactions of the American Mathematical Society, 372(5):3019–3062, May 2019. ISSN 1088-6850.
  • Chepoi and Seston [2009] V. Chepoi and M. Seston. Seriation in the Presence of Errors: A Factor 16 Approximation Algorithm for ll_{\infty}-Fitting Robinson Structures to Distances. Algorithmica, 59:521–568, 02 2009.
  • Chuangpishit et al. [2015] H. Chuangpishit, M. Ghandehari, M. Hurshman, J. Janssen, and N. Kalyaniwalla. Linear embeddings of graphs and graph limits. Journal of Combinatorial Theory, Series B, 113:162–184, Jul 2015. ISSN 0095-8956.
  • Flammarion et al. [2019] N. Flammarion, C. Mao, and P. Rigollet. Optimal rates of statistical seriation. Bernoulli, 25(1):623–653, 2019.
  • Fogel et al. [2016] F. Fogel, A. d’Aspremont, and M. Vojnovic. Spectral ranking using seriation. J. Mach. Learn. Res., 17(1):3013–3057, January 2016. ISSN 1532-4435.
  • Frieze and Kannan [1999] A. Frieze and R. Kannan. Quick approximation to matrices and applications. Combinatorica, 19:175–220, 02 1999.
  • Fulkerson and Gross [1965] D. R. Fulkerson and O. A. Gross. Incidence matrices and interval graphs. Pacific J. Math., 15:835–855, 1965. ISSN 0030-8730,1945-5844.
  • Ghandehari and J. [2024] M. Ghandehari and Jeannette J. Graph sequences sampled from Robinson graphons. European J. Combin., 116:Paper No. 103859, 31, 2024. ISSN 0195-6698,1095-9971.
  • Ghandehari and Mishura [2024] M. Ghandehari and T. Mishura. Robust Recovery of Robinson Property in LpL^{p}-Graphons: A Cut-Norm Approach. Electron. J. Comb., 31, 2024.
  • Hladký and Rocha [2020] J. Hladký and I. Rocha. Independent sets, cliques, and colorings in graphons. European Journal of Combinatorics, 88:103108, 2020. ISSN 0195-6698. Selected papers of EuroComb17.
  • Kunszenti-Kovács et al. [2019] D. Kunszenti-Kovács, L. Lovász, and B. Szegedy. Measures on the square as sparse graph limits. Journal of Combinatorial Theory, Series B, 138:1–40, 2019. ISSN 0095-8956.
  • Laurent et al. [2017] M. Laurent, M. Seminaroti, and S. Tanigawa. A structural characterization for certifying robinsonian matrices. Electronic Journal of Combinatorics, 24, 01 2017.
  • Liiv [2010] I. Liiv. Seriation and matrix reordering methods: An historical overview. Statistical Analysis and Data Mining: The ASA Data Science Journal, 3(2):70–91, 2010.
  • Lovász and Szegedy [2006] L. Lovász and B. Szegedy. Limits of dense graph sequences. Journal of Combinatorial Theory, Series B, 96(6):933–957, 2006. ISSN 0095-8956.
  • Robinson [1951] W. S. Robinson. A method for chronologically ordering archaeological deposits. American Antiquity, 16(4):293–301, 1951.