Title: PHASE TRANSITIONS FOR THE INFORMATION BOTTLE-NECK IN REPRESENTATION LEARNING

Abstract: In the Information Bottleneck (IB), when tuning the relative strength between compression and prediction terms, how do the two terms behave, and what's their relationship with the dataset and the learned representation? In this paper, we set out to answer these questions by studying multiple phase transitions in the IB objective: IB β [p(z|x)] = I(X; Z)-βI(Y ; Z) defined on the encoding distribution p(z|x) for input X, target Y and representation Z, where sudden jumps of dI(Y ;Z) dβ and prediction accuracy are observed with increasing β. We introduce a definition for IB phase transitions as a qualitative change of the IB loss landscape, and show that the transitions correspond to the onset of learning new classes. Using secondorder calculus of variations, we derive a formula that provides a practical condition for IB phase transitions, and draw its connection with the Fisher information matrix for parameterized models. We provide two perspectives to understand the formula, revealing that each IB phase transition is finding a component of maximum (nonlinear) correlation between X and Y orthogonal to the learned representation, in close analogy with canonical-correlation analysis (CCA) in linear settings. Based on the theory, we present an algorithm for discovering phase transition points. Finally, we verify that our theory and algorithm accurately predict phase transitions in categorical datasets, predict the onset of learning new classes and class difficulty in MNIST, and predict prominent phase transitions in CIFAR10.

1. INTRODUCTION:

The Information Bottleneck (IB) objective (Tishby et al., 2000) :

EQUATION

explicitly trades off model compression (I(X; Z), I(•; •) denoting mutual information) with predictive performance (I(Y ; Z)) using the Lagrange multiplier β, where X, Y are observed random variables, and Z is a learned representation of X. The IB method has proved effective in a variety of scenarios, including improving the robustness against adversarial attacks (Alemi et al., 2016; Fischer, 2018) , learning invariant and disentangled representations (Achille & Soatto, 2018a; b) , underlying information-based geometric clustering (Strouse & Schwab, 2017b) , improving the training and performance in adversarial learning (Peng et al., 2018) , and facilitating skill discovery (Sharma et al., 2019) and learning goal-conditioned policy (Goyal et al., 2019) in reinforcement learning.

From Eq. (1) we see that when β → 0 it will encourage I(X; Z) = 0 which leads to a trivial representation Z that is independent of X, while when β → +∞, it reduces to a maximum likelihood objective1 that does not constrain the information flow. Between these two extremes, how will the IB objective behave? Will prediction and compression performance change smoothly, or do there exist interesting transitions in between? In Wu et al. (2019) , the authors observe and study the learnability transition, i.e. the β value such that the IB objective transitions from a trivial global minimum to learning a nontrivial representation. They also show how this first phase transition relates to the structure of the dataset. However, to answer the full question, we need to consider the full range of β. Motivation. To get a sense of how I(Y ; Z) and I(X; Z) vary with β, we train Variational Information Bottleneck (VIB) models (Alemi et al., 2016) on the CIFAR10 dataset (Krizhevsky & Hinton, 2009) , where each experiment is at a different β and random initialization of the model. Fig. 1 shows the I(X; Z), I(Y ; Z) and accuracy vs. β, as well as I(Y ; Z) vs. I(X; Z) for CIFAR10 with 20% label noise (see Appendix I for details).

From Fig. 1 (b)(c), we see that as we increase β, instead of going up smoothly, both I(X; Z) and I(Y ; Z) show multiple phase transitions, where the slopes dI(X;Z) dβ and dI(Y ;Z) dβ are discontinuous and the accuracy has discrete jumps. The observation lets us refine our question: When do the phase transitions occur, and how do they depend on the structure of the dataset? These questions are important, since answering them will help us gain a better understanding of the IB objective and its close interplay with the dataset and the learned representation.

Moreover, the IB objective belongs to a general form of two-term trade-offs in many machine learning objectives: L = Prediction-loss + β • Complexity, where the complexity term generally takes the form of regularization. Usually, learning is set at a specific β. Many more insights can be gained if we understand the behavior of the prediction loss and model complexity with varying β, and how they depend on the dataset. The techniques developed to address the question in the IB setting may also help us understand the two-term tradeoff in other learning objectives.

Contributions. In this work, we begin to address the above question in IB settings. Specifically:

• We identify a qualitative change of the IB loss landscape w.r.t. p(z|x) for varying β as IB phase transitions (Section 3). • Based on the definition, we introduce a quantity G[p(z|x)] and use it to prove a theorem giving a practical condition for IB phase transitions. We further reveal the connection between G[p(z|x)] and the Fisher information matrix when p(z|x) is parameterized by θ (Section 3). • We reveal the close interplay between the IB objective, the dataset and the learned representation, by showing that in IB, each phase transition corresponds to learning a new nonlinear component of maximum correlation between X and Y , orthogonal to the previously-learned Z, and each with decreasing strength (Section 4).

To the best of our knowledge, our work provides the first theoretical formula to address IB phase transitions in the most general setting. In addition, we present an algorithm for iteratively finding the IB phase transition points (Section 5). We show that our theory and algorithm give tight matches with the observed phase transitions in categorical datasets, predict the onset of learning new classes and class difficulty in MNIST, and predict prominent transitions in CIFAR10 experiments (Section 6).

2. RELATED WORK:

The Information Bottleneck Method (Tishby et al., 2000) provides a tabular method based on the Blahut-Arimoto (BA) Algorithm (Blahut, 1972) to numerically solve the IB functional for the optimal encoder distribution P (Z|X), given the trade-off parameter β and the cardinality of the representation variable Z. This work has been extended in a variety of directions, including to the case where all three variables X, Y, Z are multivariate Gaussians (Chechik et al., 2005) , cases of variational bounds on the IB and related functionals for amortized learning (Alemi et al., 2016; Achille & Soatto, 2018a; Fischer, 2018) , and a more generalized interpretation of the constraint on model complexity as a Kolmogorov Structure Function (Achille et al., 2018) . Previous theoretical analyses of IB include Rey & Roth (2012) , which looks at IB through the lens of copula functions, and Shamir et al. (2010) , which starts to tackle the question of how to bound generalization with IB. We will make practical use of the original IB algorithm, as well as the amortized bounds of the Variational Informormation Bottleneck (Alemi et al., 2016) and the Conditional Entropy Bottleneck (Fischer, 2018) .

Phase transitions, where key quantities change discontinuously with varying relative strength in the two-term trade- 

β c i = 1 1-λi

, where λ i are eigenvalues of the matrix Σ x|y Σ -1

x , and Σ x is the covariance matrix. This work provides valuable insights for IB, but is limited to the special case that X, Y and Z are jointly Gaussian. Phase transitions in the general IB setting have also been observed, which Tishby (2018) describes as "information bifurcation". In Wu et al. (2019) , the authors study the first phase transition, i.e. the learnability phase transition, and provide insights on how the learnability depends on the dataset. Our work is the first work that addresses all the IB phase transitions in the most general setting, and provides theoretical insights on the interplay between the IB objective, its phase transitions, the dataset, and the learned representation.

3. FORMULA FOR IB PHASE TRANSITIONS:



3.1. DEFINITIONS:

Let X ∈ X , Y ∈ Y, Z ∈ Z be random variables denoting the input, target and representation, respectively, having a joint probability distribution p(X, Y, Z), with X × Y × Z its support. X, Y and Z satisfy the Markov chain Z -X -Y , i.e. Y and Z are conditionally independent given X. We assume that the integral (or summing if X, Y or Z are discrete random variables) is on X × Y × Z. We use x, y and z to denote the instances of the respective random variables. The above settings are used throughout the paper. We can view the IB objective IB β [p(z|x)] (Eq. 1) as a functional of the encoding distribution p(z|x). To prepare for the introduction of IB phase transitions, we first define relative perturbation function and second variation, as follows. Definition 1. Relative perturbation function: For p(z|x), its relative perturbation function r(z|x) is a bounded function that maps X × Z to R and satisfies E z∼p(z|x) [r(z|x)] = 0. Formally, define

Q Z|X := {r(z|x) : X × Z → R E z∼p(z|x) [r(z|x)] = 0, and ∃M > 0 s.t. ∀X ∈ X , Z ∈ Z, |r(z|x)| ≤ M }. We have that r(z|x) ∈ Q Z|X iff r(z|x) is a relative perturbation function of p(z|x). The perturbed probability (density) is p (z|x) = p(z|x) (1 + • r(z|x)) for some > 0.

Definition 2. Second variation: Let functional F [f (x)] be defined on some normed linear space R. Let us add a perturbative function • h(x) to f (x), and now the functional

F [f (x) + • h(x)] can be expanded as ∆F [f (x)] = F [f (x) + • h(x)] -F [f (x)] = ϕ 1 [ • h(x)] + ϕ 2 [ • h(x)] + ϕ r [ • h(x)] • h(x) 2 such that lim →0 ϕ r [ • h(x)] = 0, where • denotes the norm, ϕ 1 [ • h(x)] = dF [f (x)] d is a linear functional of •h(x)

, and is called the first variation, denoted as δF

[f (x)]. ϕ 2 [ •h(x)] = 1

We can think of the perturbation function • h(x) as an infinite-dimensional "vector" (x being the indices), with being its amplitude and h(x) its direction. With the above preparations, we define the IB phase transition as a change in the local curvature on the global minimum of IB β [p(z|x)]. Definition 3. IB phase transitions: Let r(z|x) ∈ Q Z|X be a perturbation function of p(z|x), p * β (z|x) denote the optimal solution of IB β [p(z|x)] at β, where the IB functional IB[•] is defined in Eq. (1). The IB phase transitions β c i are the β values satisfying the following two conditions:

(1) ∀r(z|x) ∈ Q Z|X , δ 2 IB β [p(z|x)] p * β (z|x) ≥ 0; (2) lim β →β + inf r(z|x)∈Q Z|X δ 2 IB β [p(z|x)] p * β (z|x) = 0 -.

Here β + and 0 -denote one-sided limits.

We can understand the δ 2 IB β [p(z|x)] as a local "curvature" of the IB objective IB β (Eq. 1) w.r.t. p(z|x), along some relative perturbation r(z|x). A phase transition occurs when the convexity of IB β [p(z|x)] w.r.t. p(z|x) changes from a minimum to a saddle point in the neighborhood of its optimal solution p * β (z|x) as β increases from β c to β c + 0 + . This means that there exists a perturbation to go downhill and find a better minimum. We validate this definition empirically below.

3.2. CONDITION FOR IB PHASE TRANSITIONS:

The definition for IB phase transition (Definition 3) indicates the important role δ 2 IB β [p(z|x)] plays on the optimal solution in providing the condition for phase transitions. To concretize it and prepare for a more practical condition for IB phase transitions, we expand IB β [p(z|x)(1 + • r(z|x))] to the second order of , giving: Lemma 0.1.

For IB β [p(z|x)], the condition of ∀r(z|x) ∈ Q Z|X , δ 2 IB β [p(z|x)] ≥ 0 is equivalent to β ≤ G[p(z|x)]. The threshold function G[p(z|x)] is given by: G[p(z|x)] := inf r(z|x)∈Q Z|X G[r(z|x); p(z|x)] G[r(z|x); p(z|x)] := E x,z∼p(x,z) [r 2 (z|x)] -E z∼p(z) E x∼p(x|z) [r(z|x)] 2 E y,z∼p(y,z) E x∼p(x|y,z) [r(z|x)] 2 -E z∼p(z) E x∼p(x|z) [r(z|x)] 2 (2)

The proof is given in Appendix B, in which we also give Eq. ( 20) for empirical estimation. Note that Lemma 0.1 is very general and can be applied to any p(z|x), not only at the optimal solution p * β (z|x).

The Fisher Information matrix. In practice, the encoder p θ (z|x) is usually parameterized by some parameter vector θ = (θ 1 , θ 2 , ...θ k ) T ∈ Θ, e.g. weights and biases in a neural net, where Θ is the parameter field. An infinitesimal change of θ ← θ + ∆θ induces a relative perturbation • r(z|x) ∆θ T ∂logp θ (z|x) ∂θ on p θ (z|x), from which we can compute the threshold function G Θ [p θ (z|x)]:

Lemma 0.2. For IB β [p θ (z|x)] objective, the condition of ∀∆θ ∈ Θ, δ 2 IB β [p θ (z|x)] ≥ 0 is equiva- lent to β ≤ G Θ [p θ (z|x)],

where

EQUATION

)

where are the conditional Fisher information matrix (Zegers, 2015) of θ for Z conditioned on X and Y , respectively. λ max is the largest eigenvalue of C -1 I Z|Y (θ) -I Z (θ) (C T ) -1 with v max the corresponding eigenvector, where CC T is the Cholesky decomposition of the matrix I Z|X (θ) -I Z (θ), and v max is the eigenvector for λ max . The infimum is attained at ∆θ = (C T ) -1 v max .

I Z (θ) := dzp θ (z) ∂logp θ (z)

The proof is in appendix C. We see that for parameterized encoders p θ (z|x), each term of G[p(z|x)] in Eq. ( 2) can be replaced by a bilinear form with the Fisher information matrix of the respective variables. Although this lemma is not required to understand the more general setting of Lemma 0.1, where the model is described in a functional space, Lemma 0.2 helps understand G[p(z|x)] for parameterized models, which permits directly linking the phase transitions to the model's parameters.

Phase Transitions. Now we introduce Theorem 1 that gives a concrete and practical condition for IB phase transitions, which is the core result of the paper: Theorem 1. The IB phase transition points {β c i } as defined in Definition 3 are given by the roots of the following equation:

G[p * β (z|x)] = β (4)

where G[p(z|x)] is given by Eq. ( 2) and p * β (z|x) is the optimal solution of IB β [p(z|x)] at β.

We can understand Eq. ( 4) as the condition when δ 2 IB β [p(z|x)] is about to be able to be negative at the optimal solution p * β (z|x) for a given β. The proof for Theorem 1 is given in Appendix D. In Section 4, we will analyze Theorem 1 in detail.

4. UNDERSTANDING THE FORMULA FOR IB PHASE TRANSITIONS:

In this section we set out to understand G[p(z|x)] as given by Eq. ( 2) and the phase transition condition as given by Theorem 1, from the perspectives of Jensen's inequality and representational maximum correlation.

4.1. JENSEN'S INEQUALITY:

The condition for IB phase transitions given by Theorem 1 involves

G[p(z|x)] = inf r(z|x)∈Q Z|X G[r(z|x); p(z|x)] which is in itself an optimization problem. We can understand G[p(z|x)] = inf r(z|x)∈Q Z|X A-C B-C in Eq.

(2) using Jensen's inequality:

E x,z∼p(x,z) [r 2 (z|x)] A ≥ E y,z∼p(y,z) E x∼p(x|y,z) [r(z|x)] 2 B ≥ E z∼p(z) E x∼p(x|z) [r(z|x)] 2 C (5)

The equality between A and B holds when the perturbation r(z|x) is constant w.r.t. x for any z; the equality between B and C holds when E x∼p(x|y,z) [r(z|x)] is constant w.r.t. y for any z. Therefore, the minimization of A-C B-C encourages the relative perturbation function r(z|x) to be as constant w.r.t. x as possible (minimizing intra-class difference), but as different w.r.t. different y as possible (maximizing inter-class difference), resulting in a clustering of the values of r(z|x) for different examples x according to their class y. Because of this clustering property in classification problems, we conjecture that there are at most |Y| -1 phase transitions, where |Y| is the number of classes, with each phase transition differentiating one or more classes.

4.2. REPRESENTATIONAL MAXIMUM CORRELATION:

Under certain conditions we can further simplify G[p(z|x)] and gain a deeper understanding of it. Firstly, inspired by maximum correlation (Anantharam et al., 2013) , we introduce two new concepts, representational maximum correlation and conditional maximum correlation, as follows. Definition 4. Given a joint distribution p(X, Y ), and a representation Z satisfying the Markov chain Z -X -Y , the representational maximum correlation ρ r (X, Y ; Z) is defined as

EQUATION

where

S 1 = {(f : X × Z → R, g : Y × Z → R) f, g bounded, and E x∼p(x|z) [f (x, z)] = E y∼p(y|z) [g(y, z)] = 0, E x,z∼p(x,z) [f 2 (x, z)] = E y,z∼p(y,z) [g 2 (y, z)] = 1}.

The conditional maximum correlation ρ m (X, Y |Z) is defined as:

EQUATION

where

S 2 = {(f : X → R, g : Y → R) f, g bounded, and ∀z ∈ Z : E x∼p(x|z) [f (x)] = E y∼p(y|z) [g(y)] = 0, E x∼p(x|z) [f 2 (x)] = E y∼p(y|z) [g 2 (y)] = 1}.

We prove the following Theorem 2, which expresses G[p(z|x)] in terms of representational maximum correlation and related quantities, with proof given in Appendix F.

Theorem 2. Define Q (0) Z|X := {r(z|x) : X × Z → R r bounded}. If Q (0)

Z|X and Q Z|X satisfy:

EQUATION

(ii) The representational maximum correlation and conditional maximum correlation:

EQUATION

(iii) When Z is continuous, an optimal relative perturbation function r(z|x) for G[p(z|x)] is given by

EQUATION

where z * = arg max z∈Z ρ m (X, Y |Z = z), and h * (x) is the optimal solution for the learn- Wu et al. (2019) ).

ability threshold function h * (x) = arg min h(x)∈{h:X →R | h bounded} β 0 [h(x)] with p(X, Y |Z = z * ) (β 0 [h(x)] is given in Theorem 4 of

(iv) For discrete X, Y and Z, we have

EQUATION

where σ 2 (Z) is the second largest singular value of the matrix

Q X,Y |Z := p(x,y|z) √ p(x|z)p(y|z) x,y = p(x,y) √ p(x)p(y) p(z|x) p(z|y) x,y

.

Theorem 2 furthers our understanding of G[p(z|x)] and the phase transition condition (Theorem 1), which we elaborate as follows.

Discovering maximum correlation in the orthogonal space of a learned representation: Intuitively, the representational maximum correlation measures the maximum linear correlation between f (X, Z) and g(Y, Z) among all real-valued functions f, g, under the constraint that f (X, Z) is "orthogonal" to p(X|Z) and g

(Y, Z) is "orthogonal" to p(Y |Z). Theorem 2 (i) reveals that G[p(z|x)]

is the inverse square of this representational maximum correlation. Theorem 2 (ii) further shows that G[p(z|x)] is finding a specific z * on which maximum (nonlinear) correlation between X and Y conditioned on Z can be found. Combined with Theorem 1, we have that when we continuously increase β, for the optimal representation Z * β given by p * β (z|x) at β, ρ r (X, Y ; Z * β ) shall monotonically decrease due to that X and Y has to find their maximum correlation on the orthogonal space of an increasingly better representation Z * β that captures more information about X. A phase transition occurs when ρ r (X, Y ; Z * β ) reduces to 1 √ β , after which as β continues to increase, ρ r (X, Y ; Z * β ) will try to find maximum correlation between X and Y orthogonal to the full previously learned representation. This is reminiscent of canonical-correlation analysis (CCA) (Hotelling, 1992) in linear settings, where components with decreasing linear maximum correlation that are orthogonal to previous components are found one by one. In comparison, we show that in IB, each phase transition corresponds to learning a new nonlinear component of maximum correlation between X and Y in Z, orthogonal to the previously-learned Z. In the case of classification where different classes may have different difficulty (e.g. due to label noise or support overlap), we should expect that classes that are less difficult as measured by a larger maximum correlation between X and Y are learned earlier.

Conspicuous subset conditioned on a single z: Furthermore, we show in (iii) that an optimal relative perturbation function r(z|x) can be decomposed into a product of two factors, a δ(z-z * )

p(z)

factor that only focus on perturbing a specific point z * in the representation space, and an h * (x) factor that is finding the "conspicuous subset" (Wu et al., 2019) , i.e. the most confident, large, typical, and imbalanced subset in the X space for the distribution (X, Y ) ∼ p(X, Y |z * ).

Singular values:

In categorical settings, (iv) reveals a connection between G[p(z|x)] and the singular value of the Q X,Y |Z matrix. Due to the property of SVD, we know that the square of the singular values of Q X,Y |Z equals the non-negative eigenvalue of the matrix

Q T X,Y |Z Q X,Y |Z .

Then the phase transition condition in Theorem 1 is equivalent to a (nonlinear) eigenvalue problem. This is resonant with previous analogy with CCA in linear settings, and is also reminiscent of the linear eigenvalue problem in Gaussian IB (Chechik et al., 2005) .

5. ALGORITHM FOR PHASE TRANSITIONS DISCOVERY IN CLASSIFICATION:

As a consequence of the theoretical analysis above, we are able to derive an algorithm to efficiently estimate the phase transitions for a given model architecture and dataset. This algorithm also permits us to empirically confirm some of our theoretical results in Section 6.

Typically, classification involves high-dimensional inputs X. Without sweeping the full range of β where at each β it is a full learning problem, it is in general a difficult task to estimate the phase transitions. In Algorithm 1, we present a two-stage approach.

In the first stage, we train a single maximum likelihood neural network f θ with the same encoder architecture as in the (variational) IB to estimate p(y|x), and obtain an N × C matrix p(y|x), where N is the number of examples in the dataset and C is the number of classes. In the second stage, we perform an iterative algorithm w.r.t. G and β, alternatively, to converge to a phase transition point.

Specifically, for a given β, we use a Blahut-Arimoto type IB algorithm (Tishby et al., 2000) to efficiently reach IB optimal p * β (z|x) at β, then use SVD (with the formula given in Theorem 2 (iv)) to efficiently estimate G[p * β (z|x)] at β (step 8). We then use the G[p * β (z|x)] value as the new β and do it again (step 7 in the next iteration). At convergence, we will reach the phase transition point given by G[p * β (z|x)] = β (Theorem 1). After convergence as measured by patience parameter K, we slightly increase β by δ (step 13), so that the algorithm can discover the subsequent phase transitions.

6. EMPIRICAL STUDY:

We quantitatively and qualitatively test the ability of our theory and Algorithm 1 to provide good predictions for IB phase transitions. We first verify them in fully categorical settings, where X, Y, Z are all discrete, and we show that the phase transitions can correspond to learning new classes as we increase β. We then test our algorithm on versions of the MNIST and CIFAR10 datasets with added label noise.

Algorithm 1 Phase transitions discovery for IB Require (X, Y ): the dataset Require f θ : a neural net with the same encoder architecture as the (variational) IB Require K: patience Require δ: precision floor Require R: maximum ratio between β (th) and β. // First stage: fit p(y|x) using neural net f θ : 1: p(y|x) ← fitting (X, Y ) using f θ via maximum likelihood. 2: p(x) ← 1 N // Second stage: coordinate descent using G[p(z|x)] and IB algorithm:

3: β c 0 ← β (th) (1) 4: B ← {β c 0 } //B is a set collecting the phase transition points 5: (β (new) , β, count) ← (β c 0 , 1, 0) 6: while β (new) β < R do: 7: β ← β (new) 8: β (new) ← β (th) (β) 9: if β (new) -β < δ do: 10: count ← count + 1 11: if count > K do: 12: B ← B ∪ {β (new) } 13: β (new) ← β (new) + δ 14:

end if 15:

else: count ← 0 16:

end if 17: end while 18: return B subroutine β (th) (β): s1: Compute p * β (z|x) using the IB algorithm (Tishby et al., 2000) . s2: Furthermore, in Appendix G we show that the phase transitions correspond to the onset of separation of p(z|x) for subsets of X that correspond to different classes. This supports our conjecture from Section 4.1 that there are at most |Y| -1 phase transitions in classification problems.

EQUATION

6.2. MNIST DATASET:

For continuous X, how does our algorithm perform, and will it reveal aspects of the dataset? We first test our algorithm in a 4-class MNIST with noisy labels3 , whose confusion matrix and experimental settings are given in Appendix H. Fig. 3 (a) shows the path Alg. 1 takes. We see again that in each phase Alg. 1 converges to the phase transition points within a few iterations, and it discovers in total 3 phase transition points. Similar to the categorical case, we expect that each phase transition corresponds to the onset of learning a new class, and that the last class is much harder to learn due to a larger separation of β. Therefore, this class should have a much larger label noise so that it is hard to capture this component of maximum correlation between X and Y , as analyzed in representational maximum correlation (Section 4.2). Fig. 3 (b) plots the per-class accuracy with increasing β for running the Conditional Entropy Bottleneck (Fischer, 2018) (another variational bound on IB). We see that the first two predicted phase transition points β c 0 , β c 1 closely match the observed onset of learning class 3 and class 0. Class 1 is observed to learn earlier than expected, possibly due to the gap between the variational IB objective and the true IB objective in continuous settings. By looking at the confusion matrix for the label noise (Fig. 7 ), we see that the ordering of onset of learning: class 2, 3, 0, 1, corresponds exactly to the decreasing diagonal element p(ỹ = 1|y = 1) (increasing noise) of the classes, and as predicted, class 1 has a much smaller diagonal element p(ỹ = 1|y = 1) than the other three classes, which makes it much more difficult to learn. This ordering of classes by difficulty is what our representational maximum correlation predicts. β with the same sets of points identified. The most interesting region is right before β = 2, where accuracy decreases with β. Alg. 1 identifies both sides of that region, as well as points at or near all of the early obvious phase transitions. It also seems to miss later transitions, possibly due to the gap between the variational IB objective and the true IB objective in continuous settings.

6.3. CIFAR10 DATASET:

Finally, we investigate the CIFAR10 experiment from Section 1. The details of the experimental setup are described in Appendix I. This experiment stretches the current limits of our discrete approximation to the underlying continuous representation being learned by the models. Nevertheless, we can see in Fig. 4 that many of the visible empirical phase transitions are tightly identified by Alg. 1. Particularly, the onset of learning is predicted quite accurately; the large interval between the predicted β 3 = 1.21 and β 4 = 1.61 corresponds well to the continuous increase of I(X; Z) and I(Y ; Z) at the same interval. And Alg. 1 is able to identify many dense transitions not obviously seen by just looking at I(Y ; Z) vs. β curve alone. Alg. 1 predicts 9 phase transitions, exactly equal to |Y| -1 for CIFAR10.

7. CONCLUSION:

In this work, we observe and study the phase transitions in IB as we vary β. We introduce the definition for IB phase transitions, and based on it derive a formula that gives a practical condition for IB phase transitions. We further understand the formula via Jensen's inequality and representational maximum correlation. We reveal the close interplay between the IB objective, the dataset and the learned representation, as each phase transition is learning a nonlinear maximum correlation component in the orthogonal space of the learned representation. We present an algorithm for finding the phase transitions, and show that it gives tight matches with observed phase transitions in categorical datasets, predicts onset of learning new classes and class difficulty in MNIST, and predicts prominent transitions in CIFAR10 experiments. This work is a first theoretical step towards a deeper understanding of the phenomenon of phase transitions in the Information Bottleneck. We believe our approach will be applicable to other "trade-off" objectives, like β-VAE (Higgins et al., 2017) and InfoDropout (Achille & Soatto, 2018a) , where the model's ability to predict is balanced against a measure of complexity. Proof. Suppose that we perform a relative perturbation r(z|x) on p(z|x) such that the perturbed conditional probability is p (z|x) = p(z|x) (1 + • r(z|x)), then we have

EQUATION

p (z) = p(x)p (z|x)dx = dxp(x)p(z|x) (1 + • r(z|x)) = p(z) + • dxp(x)p(z|x)r(z|x)

Therefore, we can denote the corresponding relative perturbation r(z) on p(z) as

r(z) ≡ 1 p (z) -p(z) p(z) = 1 p(z) dxp(x)p(z|x)r(z|x) = E x∼p(x|z) [r(z|x)]

Similarly, we have Published as a conference paper at ICLR 2020

p (z|y) = p (y, z) p(y) = 1 p(y) dxp(x, y)p(z|x) (1 + • r(z|x)) = p(z|y) + • 1 p(y) dxp(x,

We have

IB β [p (z|x)] = IB β [p(z|x)(1 + • r(z|x))] = dxdzp(x)p (z|x)log p (z|x) p (z) -β • dydzp(y)p (z|y)log p (z|y) p (z) = dxdzp(x)p(z|x)(1 + • r(z|x))log p(z|x)(1 + • r(z|x)) p(z)(1 + • r(z)) -β • dydzp(y)p(z|y)(1 + • r(z|y))log p(z|y)(1 + • r(z|y)) p(z)(1 + • r(z)) = dxdzp(x)p(z|x)(1 + • r(z|x)) log p(z|x) p(z) + log (1 + • r(z|x)) -log (1 + • r(z)) -β • dydzp(y)p(z|y)(1 + • r(z|y)) log p(z|y) p(z) + log (1 + • r(z|y)) -log (1 + • r(z)) = dxdzp(x)p(z|x)(1 + • r(z|x)) log p(z|x) p(z) + ∞ n=1 (-1) n-1 n n (r(z|x) -r(z)) -β • dydzp(y)p(z|y)(1 + • r(z|y)) log p(z|y) p(z) + ∞ n=1 (-1) n-1 n n (r(z|y) -r(z))

The 0 th -order term is simply IB β [p(z|x)]. The first order term is

δIB β [p(z|x)] = • E x,z∼p(x,z) r(z|x)log p(z|x) p(z) -β • E y,z∼p(y,z) r(z|y)log p(z|y) p(z)

The n th -order term for n ≥ 2 is

δ n IB β [p(z|x)] =(-1) n n dxdzp(x)p(z|x) - 1 n [r n (z|x) -r n (z)] + r(z|x) 1 n -1 r n-1 (z|x) -r n (z) -β • (-1) n n dydzp(y)p(z|y) - 1 n [r n (z|y) -r n (z)] + r(z|y) 1 n -1 r n-1 (z|y) -r n (z) = (-1) n n n(n -1) E x,z∼p(x,z) [r n (z|x)] -nE x,z∼p(x,z) [r(z|x)r n-1 (z)] -(n -1)E z∼p(z) [r n (z)] -β • (-1) n n n(n -1) E y,z∼p(y,z) [r n (z|y)] -nE y,z∼p(y,z) [r(z|y)r n-1 (z)] -(n -1)E z∼p(z) [r n (z)] = (-1) n n n(n -1) {(E[r n (z|x)] -E[r n (z)]) -β • (E[r n (z|y)] -E[r n (z)])}

In the last equality we have used

E x,z∼p(x,z) [r(z|x)r n-1 (z)] = E z∼p(z) [r n-1 (z)E x∼p(x|z) [r(z|x)]] = E z∼p(z) [r n-1 (z)r(z)] = E z∼p(z) [r n (z)]

Combining the terms with all orders, we have

IB β [p(z|x)(1 + • r(z|x))] =IB β [p(z|x)] + • E x,z∼p(x,z) r(z|x)log p(z|x) p(z) -β • E y,z∼p(y,z) r(z|y)log p(z|y) p(z) + ∞ n=2 (-1) n n n(n -1) {(E[r n (z|x)] -E[r n (z)]) -β • (E[r n (z|y)] -E[r n (z)])}

As a side note, the KL-divergence between p (z|x) = p(z|x)(1 + • r(z|x)) and p(z|x) is

KL (p (z|x)||p(z|x)) = dzp(z|x)(1 + • r(z|x))log p(z|x)(1 + • r(z|x)) p(z|x) = dzp(z|x)(1 + • r(z|x)) • r(z|x) - 2 2 • r 2 (z|x)) + O( 3 ) = • dzp(z|x)r(z|x) + 2 2 dzp(z|x)r 2 (z|x) + O( 3 ) = 2 2 E z∼p(z|x) [r 2 (z|x)] + O( 3 )

Therefore, to the second order, we have

EQUATION

Similarly, we have z|x) ] up to the second order. Using similar procedure, we have up to the second-order,

E x∼p(x) [KL (p(z|x)||p (z|x))] = 2 2 E[r 2 (

E y∼p(y) [KL (p (z|y)||p(z|y))] = E y∼p(y) [KL (p(z|y)||p (z|y))] = 2 2 E[r 2 (z|y)] KL (p (z)||p(z)) = KL (p(z)||p (z)) = 2 2 E[r 2 (z)]

B PROOF OF LEMMA 0.1

Proof. From Lemma 2.1, we have

EQUATION

The condition of ∀r(z|x

EQUATION

Using Jensen's inequality and the convexity of the square function, we have

E[r 2 (z|y)] = E y,z∼p(y,z) E x∼p(x|y,z) [r(z|x)] 2 = E z∼p(z) E y∼p(y|z) E x∼p(x|y,z) [r(z|x)] 2 ≥ E z∼p(z) E y∼p(y|z) E x∼p(x|y,z) [r(z|x)] 2 = E z∼p(z) E x∼p(x|z) [r(z|x)] 2 = E[r 2 (z)]

The equality holds iff r(z|y) = E x∼p(x|y,z) [r(z|x)] is constant w.r.t. y, for any z.

Using Jensen's inequality on E[r 2 (z)], we have

E[r 2 (z)] = E z∼p(z) E x∼p(x|z) [r(z|x)] 2 ≤ E z∼p(z) E x∼p(x|z) [r 2 (z|x)] = E[r 2 (z|x)],

where the equality holds iff r(z|x) is constant w.r.t. x for any z.

When E[r 2 (z|y)] -E[r 2 (z)] > 0, we have that the condition Eq. ( 16) is equivalent to ∀r If E[r 2 (z|y)] -E[r 2 (z)] = 0, substituting into Eq. ( 16), we have

EQUATION

EQUATION

which is always true due to that E[r 2 (z|x)] ≥ E[r 2 (z)], and will be a looser condition than Eq. ( 17) above. Above all, we have Eq. ( 17).

Empirical estimate of G[p(z|x)]:

To empirically estimate G[p(z|x)] from a minibatch of {(x i , y i )}, i = 1, 2, ...N and the encoder p(z|x), we can make the following Monte Carlo importance sampling estimation, where we use the samples {x j } ∼ p(x) and also get samples of {z i } ∼ p(z) = p(x)p(z|x), and have:

E x,z∼p(x,z) [r 2 (z|x)] = dxdzp(x)p(z) p(x, z) p(x)p(z) r 2 (z|x) 1 N 2 N i=1 N j=1 p(x j , z i ) p(x j )p(z i ) r 2 (z i |x j ) E z∼p(z) [r 2 (z)] =E z∼p(z) E x∼p(x|z) [r(z|x)] 2 1 N N i=1 dxp(x|z i )r(z i |x) 2 = 1 N N i=1 dxp(x) p(z i |x) p(z i ) r(z i |x) 2 1 N N i=1   1 N N j=1 p(z i |x j ) p(z i ) r(z i |x j )   2 1 N N i=1   1 N N j=1 p(z i |x j ) 1 N N k=1 p(z i |x k ) r(z i |x j )   2 = 1 N N i=1 N j=1 p(z i |x j )r(z i |x j ) N j=1 p(z i |x j ) 2 E y,z∼p(y,z) [r 2 (z|y)] =E y,z∼p(y,z) E x∼p(x|y,z) [r(z|x)] 2 1 N N i=1 dxp(x|y i , z i )r(z i |x) 2 = 1 N N i=1 1 p(y i , z i ) dxp(y i )p(x|y i )p(z i |x)r(z i |x) 2 = 1 N N i=1 dxp(y i )p(x|y i )p(z i |x)r(z i |x) dxp(y i )p(x|y i )p(z i |x) 2 1 N N i=1 xj ∈Ωx(yi) p(z i |x j )r(z i |x j ) xj ∈Ωx(yi) p(z i |x j ) 2 = 1 N N i=1 N j=1 p(z i |x j )r(z i |x j )1 [y i = y j ] N j=1 p(z i |x j )1 [y i = y j ] 2

Here Ω x (y i ) denotes the set of x examples that has label of y i , and 1[•] is an indicator function that takes value 1 if its argument is true, 0 otherwise.

EQUATION

for any x j .

Combining all terms, we have that the empirical Ĝ[p(z|x)] is given by Ĝ

EQUATION

where {z i } ∼ p(z) and {x i } ∼ p(x). It is also possible to use different distributions for importance sampling, which will results in different formulas for empirical estimation of

G[p(z|x)]. C G Θ [p θ (z|x)] FOR PARAMETERIZED DISTRIBUTION p θ (z|x)

Proof. For the parameterized 4 p θ (z|x) with θ ∈ Θ, after θ ← θ + ∆θ, where5 ∆θ ∈ Θ is an infinitesimal perturbation on θ, we have that the distribution changes from p θ (z|x) to p θ+∆θ (z|x), 4 In this paper, θ = (θ1, θ2, ...θ k ) T and ∂p θ (z|x)

∂θ = ∂p θ (z|x) ∂θ 1 , ∂p θ (z|x) ∂θ 2 , ... ∂p θ (z|x) ∂θ k T are all column vectors. ∂ 2 p θ (z|x) ∂θ 2 is a k × k matrix with (i, j) element of ∂ 2 p θ (z|x) ∂θ i ∂θ j .

and thus the relative perturbation on p θ (z|x) is

• r(z|x) = p θ+∆θ (z|x) -p θ (z|x) p θ (z|x) = 1 p θ (z|x) p θ (z|x) + ∆θ T ∂p θ (z|x) ∂θ + 1 2 ∆θ T ∂ 2 p θ (z|x) ∂θ 2 ∆θ + O( ∆θ 3 ) -p θ (z|x) ∆θ T ∂ ∂θ logp θ (z|x) + 1 2 ∆θ T 1 p θ (z|x) ∂ 2 p θ (z|x) ∂θ 2 ∆θ + O( ∆θ 3 )

where ∆θ is the norm of ∆θ in the parameter field Θ.

Similarly, we have

• r(z|y) = ∆θ T ∂ ∂θ log p θ (z|y) + 1 2 ∆θ T 1 p θ (z|y) ∂ 2 p θ (z|y) ∂θ 2 ∆θ + O( ∆θ 3 ) • r(z) = ∆θ T ∂ ∂θ log p θ (z) + 1 2 ∆θ T 1 p θ (z) ∂ 2 p θ (z) ∂θ 2 ∆θ + O( ∆θ 3 )

Substituting the above expressions into the expansion of IB β [p(z|x)] in Eq. ( 12), and preserving to the second order ∆θ 2 , we have

IB β [p θ (z|x)(1 + • r(z|x))] =IB β [p θ (z|x)] + • E x,z∼p θ (x,z) r(z|x)log p θ (z|x) p θ (z) -β • E y,z∼p θ (y,z) r(z|y)log p θ (z|y) p θ (z) + 2 1 • 2 E x,z∼p θ (x,z) [r 2 (z|x)] -E z∼p θ (z) [r 2 (z)] -β • E y,z∼p θ (y,z) [r 2 (z|y)] -E z∼p θ (z) [r 2 (z)] =IB β [p θ (z|x)] + E x,z∼p θ (x,z) ∆θ T ∂ ∂θ log p θ (z|x) + 1 2 ∆θ T 1 p θ (z|x) ∂ 2 p θ (z|x) ∂θ 2 ∆θ log p θ (z|x) p θ (z) -β • E y,z∼p θ (y,z) ∆θ T ∂ ∂θ log p θ (z|y) + 1 2 ∆θ T 1 p θ (z|y) ∂ 2 p θ (z|y) ∂θ 2 ∆θ log p θ (z|y) p θ (z) + 1 2 E x,z∼p θ (x,z) ∆θ T ∂ ∂θ log p θ (z|x) 2 -E z∼p θ (z) ∆θ T ∂ ∂θ log p θ (z) 2 - β 2 E y,z∼p θ (y,z) ∆θ T ∂ ∂θ log p θ (z|y) 2 -E z∼p θ (z) ∆θ T ∂ ∂θ log p θ (z) 2 =IB β [p θ (z|x)] + ∆θ T E x,z∼p θ (x,z) log p θ (z|x) p θ (z) ∂ ∂θ logp θ (z|x) -β • E x,z∼p θ (x,z) log p θ (z|x) p θ (z) ∂ ∂θ logp θ (z|x) + 1 2 ∆θ T I Z|X (θ) -I Z (θ) -β I Z|X (θ) -I Z (θ) ∆θ

In the last equality we have used are the conditional Fisher information matrix (Zegers, 2015) of θ for Z conditioned on X and Y , respectively.

E x,z∼p θ (x,z) [ 1 p θ (z|x) ∂ 2 p θ (z|x) ∂θ 2 ] = dxp(x) ∂ 2 ∂θ 2 dzp θ (z|x) = dxp(x) ∂ 2 ∂θ 2 1 = 0,

Let us look at

EQUATION

Firstly, note that δ 2 IB β [p θ (z|x)] is a quadratic function of ∆θ, and the scale of ∆θ does not change the sign of δ 2 IB β [p θ (z|x)], so the condition of ∀∆θ ∈ Θ, δ 2 IB β [p θ (z|x)] ≥ 0 is invariant to the scale of ∆θ, and is describing the "curvature" in the infinitesimal neighborhood of θ. Therefore, ∆θ can explore any value in Θ. Secondly, we see that Eq. ( 21) is a special case of Eq. ( 14) with • r(z|x) = ∆θ T ∂ ∂θ log p θ (z|x). Therefore, The inequalities due to Jensen still hold:

EQUATION

If ∆θ T I Z|Y (θ) -I Z (θ) ∆θ = 0, we have that Eq. ( 21) always holds, which is a looser condition than Eq. ( 22). Above all, we have that the condition of ∀∆θ ∈ Θ,

δ 2 IB β [p θ (z|x)] is equivalent to β ≤ G Θ [p θ (z|x)].

Moreover, (G Θ [p θ (z|x)]) -1 given by Eq. ( 22 

D PROOF OF THEOREM 1:

Proof. Define

EQUATION

where E β [•] denotes taking expectation w.r.t. the optimal solution p * β (x, y, z) = p(x, y)p * β (z|x) at β. Using Lemma 2.1, we have that the IB phase transition as defined in Definition 3 corresponds to satisfying the following two equations:

EQUATION

Now we prove that T β (β ) is continuous at

β = β, i.e. ∀ε > 0, ∃δ > 0 s.t. ∀β ∈ (β -δ, β + δ), we have |T β (β ) -T β (β)| < .

From Eq. ( 23), we have

T β (β ) -T β (β) = -(β -β) • E β [r 2 (z|y)] -E β [r 2 (z)] . Since r(z|x) is bounded, i.e. ∃M > 0 s.t. ∀z ∈ Z, x ∈ X , |r(z|x)| ≤ M , we have E β r 2 (z|y) = E β E x∼p(x|y,z) [r(z|x)] 2 ≤ E β E x∼p(x|y,z) [M ] 2 = M 2

Similarly, we have

E β r 2 (z) = E β E x∼p(x|z) [r(z|x)] 2 ≤ E β E x∼p(x|z) [M ] 2 = M 2 Hence, |T β (β ) -T β (β)| = |β -β| E β [r 2 (z|y)] -E β [r 2 (z)] ≤ 2|β -β|M 2 .

To prove that T β (β ) is continuous at

β = β, we have ∀ε > 0, ∃δ = ε 2M 2 > 0, s.t. ∀β ∈ (β -δ, β + δ), we have |T β (β ) -T β (β)| ≤ 2|β -β|M 2 < 2δM 2 = 2 ε 2M 2 M 2 = ε Hence T β (β ) is continuous at β = β.

Combining the continuity of T β (β ) at β = β, and Eq. ( 24) and ( 25 The numerator of G[r(z|x); p(z|x)] is then

E x,z∼p(x,z) (r (z|x)) 2 -E z∼p(z) (r (z)) 2 =E x,z∼p(x,z) (r(z|x) + s(z)) 2 -E z∼p(z) (r(z) + s(z)) 2 = E x,z∼p(x,z) r 2 (z|x) + 2E x,z∼p(x,z) [r(z|x)s(z)] + E x,z∼p(x,z) s 2 (z) -E z∼p(z) r 2 (z) + 2E z∼p(z) [r(z)s(z)] + E z∼p(z) s 2 (z) = E x,z∼p(x,z) r 2 (z|x) + 2E z∼p(z) s(z)E x∼p(x|z) [r(z|x)] + E z∼p(z) s 2 (z) -E z∼p(z) r 2 (z) + 2E z∼p(z) [r(z)s(z)] + E z∼p(z) s 2 (z) = E x,z∼p(x,z) r 2 (z|x) + 2E z∼p(z) [r(z)s(z)] + E z∼p(z) s 2 (z) -E z∼p(z) r 2 (z) + 2E z∼p(z) [r(z)s(z)] + E z∼p(z) s 2 (z) =E x,z∼p(x,z) r 2 (z|x) -E z∼p(z) r 2 (z)

Symmetrically, we have z)] is invariant to r (z|x) ← r(z|x) + s(z).

E y,z∼p(y,z) (r (z|y)) 2 -E z∼p(z) (r (z)) 2 = E y,z∼p(y,z) r 2 (z|y) -E z∼p(z) r 2 (z) Therefore, G[r(z|x); p(z|x)] = E x,z∼p(x,z) [r 2 (z|x)]-E z∼p(z)[ r 2 (z)] E y,z∼p(y,z) [r 2 (z|y)]-E z∼p(z) [r 2 (

F PROOF OF THEOREM 2:

Proof. Using the condition of the theorem, we have that ∀r(z|x) ∈ Q 0 Z|X , there exists r 1 (z|x) ∈ Q Z|X and s(z) ∈ {s : Z → R|s bounded} s.t. r(z|x) = r 1 (z|x) + s(z). Note that the only difference between Q Z|X and Q

(0) Z|X is that Q Z|X requires E p(z|x) [r 1 (z|x)] = 0. Using Lemma 2.2, we have inf r(z|x)∈Q (0) Z|X G[r(z|x); p(z|x)] = inf r1(z|x)∈Q Z|X G[r 1 (z|x); p(z|x)] = G[p(z|x)]

where r(z|x) doesn't have the constraint of E p(z|x) [•] = 0.

After dropping the constraint of E z∼p(z|x) [r(z|x)] = 0, again using Lemma 2.2, we can let r(z) = E x∼p(x|z) [r(z|x)] = 0 (since we can perform the transformation r (z|x) ← r(z|x) -r(z), so that the new r (z) ≡ 0). Now we get a simpler formula for G[p(z|x)], as follows:

EQUATION

where Q

(1)

Z|X := {r : X × Z → R E x∼p(x|z) [r(z|x)] = 0, r bounded}.

From Eq. ( 26), we can further require that E x,z∼p(x,z) [r 2 (z|x)] = 1. Define 

EQUATION

G[p(z|x)] = 1 ρ 2 s (X, Y ; Z) (i)

We only have to prove that ρ s (X, Y ; Z) = ρ r (X, Y ; Z), where ρ r (X, Y ; Z) is defined in Definition 4.

We have

E[f (X, Z)g(Y, Z)] = dxdydzp(x, y, z)f (x, z)g(y, z) = dydzp(y, z)g(y, z) dxp(x|y, z)f (x, z) ≡ dydzp(y, z)g(y, z)F (y, z) ≤ dydzp(y, z)g 2 (y, z) • dydzp(y, z)F 2 (y, z)

where F (y, z) := dxp(x|y, z)f (x, z). We have used Cauchy-Schwarz inequality, where the equality holds when g(y, z) = αF (y, z) for some α. Since E[g 2 (y, z)] = 1, we have α 2 E[F 2 (y, z)] = 1.

Taking the supremum of (E[f (X, Z)g(Y, Z)]) 2 w.r.t. f and g, we have

ρ 2 r (X, Y ; Z) = sup (f (X,Z),g(Y,Z))∈S1 (E[f (X, Z)g(Y, Z)]) 2 = sup (f (x,z),g(y,z))∈S1 dydzp(y, z)g 2 (y, z) • dydzp(y, z)F 2 (y, z) = sup f (x,z)∈Q (2) Z|X dydzp(y, z)F 2 (y, z) = sup f (x,z)∈Q (2) Z|X dydzp(y, z) dxp(x|y, z)f (x, z) 2 = sup f (X,Z)∈Q (2) Z|X E[(E[f (X, Z)|Y, Z]) 2 ] ≡ρ 2 s (X, Y ; Z)

Here S 1 is defined in Definition 4. By definition both ρ r (X, Y ; Z) and ρ s (X, Y ; Z) take non-negative values. Therefore,

EQUATION

(ii) Using the definition of ρ r (X, Y ; Z), we have

ρ 2 r (X, Y ; Z) ≡ sup f (x,z)∈Q (2) Z|X dydzp(y, z) dxp(x|y, z)f (x, z) 2 = sup f (x,z)∈Q (2) Z|X dzp(z) dyp(y|z) dxp(x|y, z)f (x, z) 2 ≡ sup f (x,z)∈Q (2) Z|X dzp(z)W [f (x, z)]

where , z) ] is equivalent to the following two-stage supremum:

W [f (x, z)] := dyp(y|z) dxp(x|y, z)f (x, z) 2 . Denote c(z) := p(z)E x∼p(x|z) [f 2 (x, z)], we have c(z)dz = E x,z∼p(x,z) [f 2 (x, z)] = 1. Then the supremum ρ 2 r (X, Y ; Z) = sup f (x,z)∈Q (2) Z|X dzp(z)W [f (x

EQUATION

where Q

(3)

Z|X := {X × Z → R E x∼p(x|z) [f 2 (x, z)] = c(z) p(z) , E x∼p(x|z) [f (x, z)] = 0, f bounded} We can think of the inner supremum sup f (x,z)∈Q (3) Z|X W [f (x, z)] as only w.r.t.

x, for some given z. Now let's consider another supremum:

EQUATION

where

Q (h) X := {h : X → R E p(x|z) [h(x)] = 0, E p(x|z) [h 2 (x)] = 1, h bounded}.

Using similar technique in (ii), it is easy to prove that it equals ρ 2 m (X, Y |Z) as defined in Definition 4.

Comparing Eq. ( 30) and the supremum:

sup f (x,z)∈Q (3) Z|X W [f (x, z)]

we see that the only difference is that in the latter

E x∼p(x|z) [f 2 (x, z)] equals c(z) p(z) instead of 1. Since W [f (x, z)] is a quadratic functional of f (x, z), we have sup f (x,z)∈Q (3) Z|X W [f (x, z)] = c(z) p(z) ρ 2 m (X, Y |Z) Therefore, ρ r (X, Y ; Z) = sup c(z): c(z)dz=1 dz p(z) sup f (x,z)∈Q (3) Z|X W [f (x, z)] = sup c(z): c(z)dz=1 dz p(z) c(z) p(z) ρ 2 m (X, Y |Z) = sup c(z): c(z)dz=1 dz c(z)ρ 2 m (X, Y |Z = z) = sup Z∈Z ρ 2 m (X, Y |Z)

where in the last equality we have let c(z) have "mass" only on the place where ρ

2 m (X, Y |Z = z) attains supremum w.r.t. z. (iii) When Z is a continuous variable, let f (x, z) = f X (x) δ(z-z0) p(z) , where δ(•) is the Dirac-delta function, z 0 is a parameter, f X (x) ∈ Q (f ) X |Z , with Q (f ) X |Z := {f X : X → R f X bounded; ∀Z ∈ Z : E X∼p(X|Z) [f X (x)] = 0, E X∼p(X|Z) [f 2 X (x)] = 1}. We have E x∼p(x|z) [f (x, z)] = p(x|z)f (x, z)dx = δ(z -z 0 ) p(z) p(x|z)f X (x)dx = δ(z -z 0 ) p(z) • 0 = 0 And E[f 2 (X, Z)] = p(x, z)f 2 (x, z)dxdz = p(x, z)f 2 X (x) δ(z -z 0 ) p(z) dxdz = dzδ(z -z 0 ) dxp(x|z)f 2 X (x)dx = dzδ(z -z 0 ) • 1 =1 Therefore, such constructed f (x, z) = f X (x) δ(z-z0) p(z) ∈ Q (2)

Z|X , satisfying the requirement for ρ s (X, Y ; Z) (which equals ρ r (X, Y ; Z) by Eq. 28).

Substituting in the special form of f (x, z) into the expression of ρ s (X, Y ; Z) in Eq. ( 27), we have

sup f (x,z):f (x,z)=f X (x) δ(z-z 0 ) p(z) ,f X (x)∈Q (f ) X |Z dzp(z) dyp(y|z) dxp(x|y, z)f (x, z) 2 = sup f X (x)∈Q (f ) X |Z ,z0∈Z dzp(z) dyp(y|z) dxp(x|y, z)f X (x) δ(z -z 0 ) p(z) 2 = sup z0∈Z dzp(z) δ(z -z 0 ) p(z) sup f X (x)∈Q (f ) X |Z dyp(y|z) dxp(x|y, z)f X (x) 2 = sup z0∈Z dzδ(z -z 0 ) sup f X (X)∈Q (f ) X |Z E[(E[f X (X)|Y, Z = z]) 2 |Z = z] = sup z0∈Z dzδ(z -z 0 )ρ 2 m (X, Y |Z = z) = sup z0∈Z ρ 2 m (X, Y |Z = z 0 ) = sup Z∈Z ρ 2 m (X, Y |Z) We can identify sup f X (X)∈Q (f ) X |Z E[(E[f X (X)|Y, Z = z]) 2 |Z = z] with ρ 2 m (X, Y |Z = z) because f X (x) satisfies the requirement for conditional maximum correlation that E p(x|z) [f X (x)] = 0 and E p(x|z) [f 2 X (x)] = 1

, for any z, and using the same technique in (i), it is straightforward to prove that sup f X (X)∈Q (f )

X |Z E[(E[f X (X)|Y, Z = z]) 2 |Z = z]

equals the conditional maximum correlation as defined in Definition 4.

Since the conditional maximum correlation can be viewed as the maximum correlation between X and Y , where X, Y ∼ p(X, Y |Z), using the equality of (β 0 [h(x)])

-1 = ρ 2 m (X; Y ) (Eq. 7 in Wu et al. ( 2019)), we can identify the h(x) in β 0 [h(x)] with the f X (X) here, and an optimal f * X (X) that maximizes ρ 2 m (X, Y |Z) is also an optimal h * (x) that minimizes β 0 [h(x)]. 

G SUBSET SEPARATION AT PHASE TRANSITIONS:

In this section we study the behavior of p(z|x) on the phase transitions. We use the same categorical dataset (where |X| = |Y | = |Z| = 3 and p(x) is uniform, and p(y|x) is given in Fig. 5 ). In Fig. 6 we show the p(z|x) on the simplex before and after each phase transition. We see that the first phase transition corresponds to the separation of x = 2 (belonging to y = 2) w.r.t. x ∈ {0, 1} (belonging to classes y ∈ {0, 1}), on the p(z|x) simplex. The second phase transition corresponds to the separation of x = 0 with x = 1. Therefore, each phase transition corresponds to the ability to distinguish subset of examples, and learning of new classes.

H MNIST EXPERIMENT DETAILS:

We use the MNIST training examples with class 0, 1, 2, 3, with a hidden label-noise matrix as given in Fig. 7 , based on which at each minibatch we dynamically sample the observed label. We use conditional entropy bottleneck (CEB) (Fischer, 2018) as the variational IB objective, and run multiple independent instances with different the target β. We jump start learning by started training at β = 100 for 100 epochs, annealing β from 100 down to the target β over 600 epochs, and continue to train at the target epoch for another 800 epochs. The encoder is a three-layer neural net, where each hidden layer has 512 neurons and leakyReLU activation, and the last layer has linear activation. The classifier p(y|z) is a 2-layer neural net with a 128-neuron ReLU hidden layer. The backward encoder p(z|y) is also a 2-layer neural net with a 128-neuron ReLU hidden layer. We trained with Adam (Kingma & Welling, 2013) at learning rate of 10 -3 , and anneal down with factor 1/(1 + 0.01 • epoch). For Alg. 1, for the f θ we use the same architecture as the encoder of CEB, and use |Z| = 50 in Alg. 1.

I CIFAR10 EXPERIMENT DETAILS:

We use the same CIFAR10 class confusion matrix provided in Wu et al. (2019) to generate noisy labels with about 20% label noise on average (reproduced in Table 1 ). We trained 28 × 1 Wide ResNet (He et al., 2016; Zagoruyko & Komodakis, 2016) models using the open source implementation from Cubuk et al. (2018) as encoders for the Variational Information Bottleneck (VIB) (Alemi et al., 2016) . The 10 dimensional output of the encoder parameterized a mean-field Gaussian with unit covariance. Samples from the encoder were passed to the classifier, a 2 layer MLP. The marginal distributions were mixtures of 500 fully covariate 10-dimensional Gaussians, all parameters of which are trained.

With this standard model, we trained 251 different models at β from 1.0 to 6.0 with step size of 0.02. As in Wu et al. (2019) , we jump-start learning by annealing β from 100 down to the target β. We do this over the first 4000 steps of training. The models continued to train for another 56,000 gradient steps after that, a total of 600 epochs. We trained with Adam (Kingma & Ba, 2015) at a base learning rate of 10 -3 , and reduced the learning rate by a factor of 0.5 at 300, 400, and 500 epochs. The models converged to essentially their final accuracy within 40,000 gradient steps, and then remained stable. 

:

For example, in classification, it reduces to cross-entropy loss.

d 2 F [f (x)] d 2is a quadratic functional of • h(x), and is called the second variation, denoted as δ 2 F [f (x)].

We use

classes since it is simpler than the full 10 classes, but still potentially possesses phase transitions. We use noisy label to mimic realistic settings where the data may be noisy and also to have controllable difficulty for different classes.

Note that since Θ is a field, it is closed under subtraction, we have ∆θ ∈ Θ.

Here all the Fisher information matrices are real symmetric, thus Hermitian.

In the definition of ρr(X, Y ; Z), we have used an equivalent format f (x, z) instead of r(z|x).

Figures and Tables:

Figure 1: CIFAR10 plots (a) showing the information plane, as well as β vs (b) I(X; Z) and I(Y ; Z), and (c) accuracy, all on the training set with 20% label noise. The arrows point to empirically-observed phase transitions. The vertical lines correspond to phase transitions found with Alg. 1.

information matrix of θ for Z, I Z|X (θ) := dxdzp(x)p θ (z|x) ∂logp θ (z|x) ∂θ ∂logp θ (z|x) ∂θ T , I Z|Y (θ) := dydzp(y)p θ (z|y) ∂logp θ (z|y) ∂θ ∂logp θ (z|y) ∂θ T

For discrete X, Z such that the cardinality |Z| ≥ |X |, this is generally true since in this scenario, h(x, z) and s(z) have |X ||Z| + |Z| unknown variables, but the condition has only |X ||Z| + |X | linear equations. The difference between Q Z|X and Q (0) Z|X is that Q (0) Z|X does not have the requirement of E p(z|x) [r(z|x)] = 0. Combined with Lemma 2.2, this condition allows us to replace r(z|x) ∈ Q Z|X by r(z|x) ∈ Q (0) Z|X in Eq. (2).

z|x)] using SVD (Eq. 8 and 11). s3: return β(new)   6.1 CATEGORICAL DATASET For categorical datasets, X and Y are discrete, and p(X) and p(Y |X) are given. To test Theorem 1, we use the Blahut-Arimoto IB algorithm to compute the optimal p * β (z|x) for each β. I(Y ; Z * ) vs. β is plotted in Fig.2 (a). There are two phase transitions at β c 0 and β c 1 . For each β and the corresponding p * β (z|x), we use the SVD formula (Theorem 2) to compute G[p * β (z|x)], shown in Fig. 2 (b). We see that G[p * β (z|x)] = β at exactly the observed phase transition points β c 0 and β c 1 . Moreover, starting at β = 1, Alg. 1 converges to each phase transition points within few iterations. Our other experiments with random categorical datasets show similarly tight matches.

Figure 2: (a) I(Y ; Z * ) vs. β for a categorical dataset with |X| = |Y | = |Z| = 3, where Z * is given by p * β (z|x), and the vertical lines are the experimentally discovered phase transition points β c 0 and β c 1 . (b) G[p * β (z|x)] vs. β for the same dataset, and the path for Alg. 1, with β c 0 and β c 1 in (a) also plotted. The dataset is given in Fig. 5.

Figure 4: (a) Accumulated G[p * β (z|x)] vs. β by running Alg. 1 with varying starting β (blue dots). Also plotted are predicted phase transition points. (b) I(X; Z) and I(Y ; Z) vs. β. The manuallyidentified phase transition points are labelled with arrows. The vertical black lines are the phase transitions identified by Alg. 1, denoted as β c 0 , β c 1 ,... β c 8 , from left to right. (c) Accuracy vs.β with the same sets of points identified. The most interesting region is right before β = 2, where accuracy decreases with β. Alg. 1 identifies both sides of that region, as well as points at or near all of the early obvious phase transitions. It also seems to miss later transitions, possibly due to the gap between the variational IB objective and the true IB objective in continuous settings.

where r(z|y) = E x∼p(x|y,z) [r(z|x)] and r(z) = E x∼p(x|z) [r(z|x)]. The expectations in the equations are all w.r.t. all variables. For example E[r 2 (z|x)] = E x,z∼p(x,z) [r 2 (z|x)].

y)p(z|x)r(z|x) And we can denote the corresponding relative perturbation r(z|y) on p(z|y) as r(z|y) ≡ 1 p (z|y) -p(z|y) p(z|y) = 1 p(z|y)p(y) dxp(x, y)p(z|x)r(z|x) = E x∼p(x|y,z) [r(z|x)] Since IB β [p(z|x)] = I(X; Z) -β • I(Y ; Z) = dxdzp(x, z)log p(z|x) p(z) -β • dydzp(y, z)log p(z|y) p(z)

) where r(z|y) = E x∼p(x|y,z) [r(z|x)] and r(z) = E x∼p(x|z) [r(z|x)].

and similarly E y,z∼p θ (y,z) [ 1 p θ (z|y) ∂ 2 p θ (z|y) ∂θ 2 In other words, the ∆θ 2 terms in the first-order variation δIB β [p θ (z|x)] vanish, and the remaining ∆θ 2 are all in δ 2 IB β [p θ (z|x)]. Also in the last expression, I Z (θ) ≡ dzp θ (z) ∂logp θ (z) ∂θ ∂logp θ (z) ∂θ T is the Fisher information matrix of θ for Z, I Z|X (θ) ≡ dxdzp(x)p θ (z|x) ∂logp θ (z|x) ∂θ ∂logp θ (z|x) ∂θ T , I Z|Y (θ) ≡ dydzp(y)p θ (z|y) ∂logp θ (z|y) ∂θ ∂logp θ (z|y) ∂θ T

) has the format of a generalized Rayleigh quotient R(A, B; x) ≡ ∆θ T A∆θ ∆θ T B∆θ where A = I Z|Y (θ)-I Z (θ) and B = I Z|X (θ)-I Z (θ) are both Hermitian matrices 6 , which can be reduced to Rayleigh quotient R(D, C T ∆θ) = (C T ∆θ) T D(C T ∆θ) (C T ∆θ) T (C T ∆θ) , with the transformation D = C -1 A(C T ) -1 where CC T is the Cholesky decomposition of B = I Z|X (θ) -I Z (θ). Moreover, we have that when G Θ [p θ (z|x)] attains its minimum value, the Reyleigh quotient R(D, C T ∆θ) attains its maximum value of λ max with C T ∆θ = v max , i.e. ∆θ = (C T ) -1 v max , where λ max is the largest eigenvalue of D and v max the corresponding eigenvector.

), we have T β (β) = 0, which is equivalent to G[p * β (z|x)] = β after simple manipulation. E INVARIANCE OF G[r(z|x); p(z|x)] TO ADDITION OF A GLOBAL REPRESENTATION Here we prove the following lemma: Lemma 2.2. G[r(z|x); p(z|x)] defined in Lemma 0.1 is invariant to the transformation r (z|x) ← r(z|x) + s(z). Proof. When we r(z|x) is shifted by a global transformation r (z|x) ← r(z|x) + s(z), we have r (z) ← E x∼p(x|z) [r(z|x) + s(z)] = E x∼p(x|z) [r(z|x)] + s(z)E x∼p(x|z) [1] = r(z) + s(z), and similarly r (z|y) ← r(z|y) + s(z).

2) Z|X E y,z∼p(y,z) E x∼p(x|y,z) [f (x, z)] {r : X × Z → R E x∼p(x|z) [r(z|x)] = 0, E x,z∼p(x,z) [r 2 (z|x)] = 1, r bounded}.Comparing with Eq. (26), it immediately follows that

iv) For discrete X, Y and Z and a given Z = z, let Q X,Y |Z := prove that its second largest singular value isρ 2 m (X, Y |Z) = sup (f,g)∈S2 E x,y∼p(x,y|z) [f (x)g(y)] (S 2 is defined in Definition 4).Let column vectors u 1 = p(x|z) and v 1 = p(y|z) (note that z is given and fixed). Also let u 2 = f (x) p(x|z) and v 2 = g(y) p(y|z). Denote inner product u, v ≡ i u i v i , and the length of a vector as ||u|| = u, u . We have||u 1 || = ||v 1 || = 1 due to the normalization of probability, ||u 2 || = ||v 2 || = 1 due to E x∼p(x|z) [f 2 (x)] = E y∼p(y|z) [g 2 (y)] = 1, and u 1 , u 2 = v 1 , v 2 = 0 due to E x∼p(x|z) [f (x)] = E y∼p(y|z) [g(y)] = 0. Furthermore, we have sup (f,g)∈S2 E x,y∼p(x,y|z) [f (x)g(y)] = max u,v u T Q X,Y |Z vwhich is exactly the second largest singular value σ 2 (Z) of the matrix Q X,Y |Z . Using the result in (ii), we have that ρ r (X, Y ; Z) = max Z∈Z σ 2 (Z).

Figure 5: p(y|x) for the categorical dataset in Fig. 2 and Fig. 6. The value in i th row and j th column denotes p(y = j|x = i). p(x) is uniform.

off, have been observed in many different learning domains, for multiple learning objectives. In Rezende & Viola (2018), the authors observe phase transitions in the latent representation of β-VAE for varying β. Strouse & Schwab (2017b) utilize the kink angle of the phase transitions in the Deterministic Information Bottleneck (DIB) (Strouse & Schwab, 2017a) to determine the optimal number of clusters for geometric clustering. Tegmark & Wu (2019) explicitly considers critical points in binary classification tasks using a discrete information bottleneck with a non-convex

Table content: <table/>

Peer review (iclr2020):
In this paper, the authors studied the phase transition in the information bottleneck, which is defined as the point where the IB loss landscape changes. The authors give a theorem showing the practical condition for IB phase transitions which is related to Fisher information. Then an algorithm finding IB phase transitions are proposed and applied to MNIST and CIFAR10 dataset. The paper is well-organized. The problem formulation, definition, theorem, algorithm and application form a whole story about finding IB phase transition. However, I found that the motivation of this paper is not very clear. Why do we want to find the IB phase transition accurately by running such a complicated algorithm (Algorithm 1). Usually, the scalar _x0008_eta is considered as a hyperparameter. By tuning _x0008_eta, one can easily produce Figure 1 and see some phase transition points. One can also tune _x0008_eta to learn different classed as shown in Figure 3(b). So designing a complicated coordinate-descent based algorithm to find a parameter seems overkilling the problem. Moreover, the mathematical part of this paper can be made more accurate. For example, in definition 1, the authors write *eta << _x0008_eta* and *|epsilon| << 1*, where the quantities eta and epsilon are used many times in the proof. However, *<<* is not a rigorous mathematical symbol. The authors should consider rewrite their proof by using limit. Therefore, I think the motivation and the mathematical quality of this paper can be further improved before getting accepted.