A T HEORY OF U SABLE I NFORMATION UNDER C OMPUTATIONAL C ONSTRAINTS
Xu
Shengjia Zhao sjzhao@stanford.edu
Jiaming Song
Stanford University tsong@cs.stanford.edu
Russell Stewart russell.sb.nebel@gmail.com
Stefano Ermon
Stanford ermon@cs.stanford.edu
A BSTRACT
We propose a new framework for reasoning about information in complex systems.
Our foundation is based on a variational extension of Shannon’s information theory that takes into account the modeling power and computational constraints of the observer.
The resulting predictive V -information encompasses mutual information and other notions of informativeness such as the coefﬁcient of determination.
Unlike Shannon’s mutual information and in violation of the data processing inequality
, V -information can be created through computation.
This is consistent with deep neural networks extracting hierarchies of progressively more informative features in representation learning.
Additionally, we show that by incorporating computational constraints, V -information can be reliably estimated from data even in high dimensions with PAC-style guarantees.
Empirically, we demonstrate predictive V -information is more effective than mutual information for structure learning and fair representation learning
I NTRODUCTION
Extracting actionable information from noisy, possibly redundant, and high-dimensional data sources is a key computational and statistical challenge at the core of AI and machine learning.
Information theory, which lies at the foundation of AI and machine learning, provides a conceptual framework to characterize information in a mathematically rigorous sense (Shannon & Weaver, 1948; Cover &
Thomas, 1991).
However, important computational aspects are not considered in information theory.
To illustrate this, consider a dataset of encrypted messages intercepted from an opponent.
According to information theory, these encrypted messages have high mutual information with the opponent’s plans.
Indeed, with inﬁnite computation, the messages can be decrypted and the plans revealed.
Modern cryptography originated from this observation by Shannon that perfect secrecy is (essentially) impossible if the adversary is computationally unbounded (Shannon & Weaver, 1948).
This motivated cryptographers to consider restricted classes of adversaries that have access to limited computational resources (Pass & Shelat, 2010).
More generally, it is known that information theoretic quantities can be expressed in terms of betting games (Cover & Thomas, 1991).
For example, the (conditional) entropy of a random variable X is directly related to how predictable X is in a certain betting game, where an agent is rewarded for correct guesses.
Yet, the standard deﬁnition unrealistically assumes agents are computationally unbounded, i.e., they can employ arbitrarily complex prediction schemes.
Leveraging modern ideas from variational inference and learning (Ranganath et al., 2013; Kingma & Welling, 2013; LeCun et al., 2015), we propose an alternative formulation based on realistic computational constraints that is in many ways closer to our intuitive notion of information, which we term predictive V -information .
Without constraints, predictive V -information specializes to classic mutual information.
Under natural restrictions, V -information specializes to other well-known notions of predictiveness, such as the coefﬁcient of determination ( R 2 ).
A consequence of this new formulation is that computation can “create usable information” (e.g., by decrypting the intercepted messages), invalidating the famous data processing inequality.
This generalizes the idea that clever feature extraction enables prediction with extremely simple (e.g., linear) classiﬁers, a key notion in modern representation and deep learning (LeCun et al., 2015).
As an additional beneﬁt, we show that predictive V -information can be estimated with statistical guarantees using the Probably Approximately Correct framework (Valiant, 1984).
This is in sharp contrast with Shannon information, which is well known to be difﬁcult to estimate for high dimensional or continuous random variables (Battiti, 1994).
Theoretically we show that the statistical guarantees of estimating V information translate to statistical guarantees for a variant of the Chow-Liu algorithm for structure learning.
In practice, when the observer employs deep neural networks as a prediction scheme, V -information outperforms methods that approximate Shannon information in various applications, including Chow-Liu tree contruction in high dimension and gene regulatory network inference.
D EFINITIONS AND N OTATIONS
To formally deﬁne the predictive V -information, we begin with a formal model of a computationally bounded agent trying to predict the outcome of a real-valued random variable Y ; the agent is either provided another real-valued random variable X as side information, or provided no side information ∅ .
We use X and Y to denote the samples spaces of X and Y respectively (while assuming they are separable), and use P ( X ) to denote the set of all probability measures over the Borel algebra on X ( P ( Y ) similarly deﬁned for Y ).
Deﬁnition 1 (Predictive Family) .
1 Let Ω =
f : X ∪ { ∅ } → P ( Y ) } .
We say that V ⊆ Ω is a predictive family if it satisﬁes
[Equation]
A predictive family is a set of predictive models the agent is allowed to use, e.g., due to computational or statistical constraints.
We refer to the additional condition in Eq.(1) as optional ignorance .
Intuitively, it means that the agent can, in the context of the prediction game we deﬁne next, ignore the side information if she chooses to.
Deﬁnition 2 (Predictive conditional V -entropy) .
Let X, Y be two random variables taking values in
X × Y , and V be a predictive family.
Then the predictive conditional V -entropy is deﬁned as
[Equation]
We additionally call H V ( Y | ∅ ) the V -entropy, and also denote it as H V ( Y )
In our notation f is a function X ∪ { ∅ } → P ( Y ) , so f [ x ] ∈ P ( Y ) is a probability measure on Y chosen based on the received side information x (we use f [ · ] instead of the more conventional f ( · ) ); and f [ x ]( y ) ∈ R is the value of the density evaluated at y ∈ Y .
Intuitively, V (conditional) entropy is the smallest expected negative log-likelihood that can be achieved predicting Y given observation (side information) X
(or no side information ∅ ), using models from V .
Eq.(1) means that whenever the agent can use P to predict Y ’s outcomes, it has the option to ignore the input, and use P no matter whether X is observed or not.
Deﬁnition 2 generalizes several known deﬁnitions of uncertainty.
For example, as shown in proposition 2, if the V is the largest possible predictive family that includes all possible models, i.e. V = Ω , then Deﬁnition 2 reduces to Shannon entropy: H Ω ( Y | X )
=
H ( Y | X ) and
H V ( Y | ∅ )
=
H Ω ( Y ) =
H ( Y ) .
By choosing more restrictive families V , we recover several other notions of uncertainty such as trace of covariance, as will be shown in Proposition 1.
Shannon mutual information is a measure of changes in entropy when conditioning on new variables:
[Equation]
Here, we will use predictive V -entropy to deﬁne an analogous quantity, I V ( X → Y ) , to represent the change in predictability of an output variable Y when given side information X .
Deﬁnition 3 (Predictive V -information) .
Let X, Y be two random variables taking values in X × Y , and V be a predictive family.
The predictive V -information from X to Y is deﬁned as
[Equation]
2. I MPORTANT S PECIAL C ASES
Several important notions of uncertainty and predictiveness are special cases of our deﬁnition.
Note that when we are deﬁning V -entropy of a random variable Y in
sample space Y ∈ R d (without side information), out of convenience we can assume X is empty X = ∅
(this does not violate our requirement that ∅ (cid:54)∈ X .)
Proposition 1.
For V -entropy and V -information, we have 1.
Let Ω be as in Def. 1.
Then H Ω ( Y ) is the Shannon entropy, H Ω ( Y | X ) is the Shannon conditional entropy, and
I Ω ( Y → X ) is the Shannon mutual information.
Let Y =
R d and V = { f : { ∅ } → P µ | µ ∈ R d } , where P µ is the distribution with density
y (cid:55)→ 1 Z e −(cid:107) y − µ (cid:107) 2 where Z = (cid:82) e −(cid:107)
y − µ (cid:107) 2 dy , then the V -entropy of a random variable Y equals its mean absolute deviation, up to an additive constant.
Let Y = R d
and V = { f : { ∅ } → N ( µ, Σ)
| µ ∈ R d , Σ
= 1 / 2 I d × d } , then the V -entropy of a random variable Y equals the trace of its covariance tr
(Cov( Y )) , up to an additive constant.
Let V =
{ f : { ∅ } → Q t ,θ , θ ∈ Θ } , where Q t ,θ is a distribution in a minimal exponential family with sufﬁcient statistics t : Y → R d and set of natural parameters Θ .
For a random variable Y with expected sufﬁcient
statistics µ Y =
E [ t ( Y )] , the V -entropy of Y is the maximum Shannon entropy over all random variables ˆ Y with identical expected sufﬁcient statistics, i.e. E
[ t ( ˆ Y )] = µ Y .
Let Y =
R d , X be any vector space, and V = { f : x (cid:55)→ N ( φ ( x ) , Σ) , x ∈ X ; ∅ (cid:55)→ N ( µ, Σ) | µ ∈ R d ; Σ = 1 / 2 I d × d , φ ∈ Φ } , where Φ is the set of linear functions { φ : X → R d } , then V -information
I V ( X → Y ) equals the (unnormalized) maximum coefﬁcient of determination R 2 · tr (Cov( Y )) for linear regression.
The trace of covariance represents a natural notion of uncertainty – for example, a random variable with zero variance (when d = 1 , tr (Cov( Y ))
=
Var( Y ) )) is trivial to predict.
Proposition 1.
shows that the trace of covariance corresponds to a notion of surprise (in the Shannon sense) for an agent restricted to make predictions using certain Gaussian models.
More broadly, a similar analogy can be drawn for other exponential families of distributions.
In the same spirit, the coefﬁcient of determination, also known as the fraction of variance explained, represents a natural notion of informativeness for computationally bounded agents.
Also note that in the case of Proposition 1.4, the V -entropy is invariant if the expected sufﬁcient statistics remain the same.
P ROPERTIES OF V - INFORMATION 3. E LEMENTARY P ROPERTIES
We ﬁrst show several elementary properties of V -entropy and V -information.
In particular, Vinformation preserves many properties of Shannon information that are desirable in a machine learning context.
For example, mutual information (and V -information) should be non-negative as conditioning on additional side information X should not reduce an agent’s ability to predict Y .
Proposition 2.
Let Y and X be any random variables on Y and X , and V and U be any predictive families, then we have
Monotonicity :
If V ⊆ U , then H V ( Y ) ≥ H U ( Y ) , H V ( Y | X ) ≥ H U
( Y | X ) .
Non-Negativity : I V ( X → Y ) ≥ 0 .
Independence : If X is independent of Y , I V ( X → Y ) =
I V ( Y → X )
= 0 .
The optional ignorance requirement in Eq.(1) is a technical condition needed for these properties to hold.
Intuitively, it guarantees that conditioning on side information does not restrict the class of densities the agent can use to predict Y .
This property is satisﬁed by many existing machine learning models, often by setting some weights to zero so that an input is effectively ignored.
3. O N THE PRODUCTION OF INFORMATION THROUGH PREPROCESSING
The Data Processing Inequality guarantees that computing on data cannot increase its mutual information with other random variables.
Formally, letting t : X → X be any function, t ( X ) cannot have higher mutual information with Y than X : I ( t ( X ); Y ) ≤ I ( X ; Y ) .
But is this property desirable?
In analyzing optimal communication, yes - it demonstrates a fundamental limit to the number of bits that can be transmitted through a communication channel.
However, we argue that in machine learning settings this property is less appropriate.
Consider an RSA encryption scheme where the public key is known.
Given plain text and its corresponding encrypted text X , if we have inﬁnite computation, we can perfectly compute one from the other.
Therefore, the plain text and the encrypted text should have identical Shannon mutual information with respect to any label Y we want to predict.
However, to any human (or machine learning algorithm), it is certainly easier to predict the label from the plain text than the encrypted text.
In other words, decryption increases a human’s ability to predict the label: processing increases the “usable information”.
More formally, denoting t as the decryption algorithm and V as a class of natural language processing functions, we have that: I V ( t ( X ) → Y ) >
I V ( X → Y )
≈ 0 .
As another example, consider the mutual information between an image’s pixels and its label.
Due to data processing inequality, we cannot expect to use a function to map raw pixels to “features” that have higher mutual information with the label.
However, the fundamental principle of representation learning is precisely the ability to learn predictive features — functions of the raw inputs that enable predictions with higher accuracy.
Because of this key difference between V -information and
Shannon information, machine learning practices such as representation learning can be justiﬁed in the information theoretic context.
3. O N THE ASYMMETRY OF PREDICTIVE V -I NFORMATION
V -information also captures the intuition that sometimes, it is easy to predict Y from X but not vice versa.
In fact, modern cryptography is founded on the assumption that certain functions h : X → Y are one-way, meaning that there exists an polynomial algorithm to compute h ( x ) but no polynomial algorithm to compute h − 1 ( y ) .
This means that if V contains all polynomial-time computable functions, then I V ( X → h ( X ))
(cid:29) I V ( h ( X ) → X ) .
This property is also reasonable in the machine learning context.
For example, several important methods for causal discovery (Peters et al., 2017) rely on this asymmetry: if X causes Y , then usually it is easier to predict Y from X than vice versa; another commonly used assumption is that Y | X can be accurately modeled by a Gaussian distribution, while X | Y cannot (Pearl, 2000).
PAC G UARANTEES FOR V - INFORMATION E STIMATION
For many practical applications of mutual information (e.g., structure learning), we do not know the joint distribution of X, Y , so cannot directly compute the mutual information.
Instead we only have samples { ( x i , y i ) } Ni =1 ∼ X, Y and need to estimate mutual information from data.
Shannon information is notoriously difﬁcult to estimate for high dimensional random variables.
Although non-parametric estimators of mutual information exist (Kraskov et al., 2004; Darbellay &
Vajda, 1999; Gao et al., 2017), these estimators do not scale to high dimensions.
Several variational estimators for Shannon information have been recently proposed (van den Oord et al., 2018; Nguyen et al., 2010; Belghazi et al., 2018), but have two shortcomings: due to their variational assumptions, their bias/variance tradeoffs are poorly understood and they are still not efﬁcient enough for high dimensional problems.
For example, the CPC estimator suffers from large bias, since its estimates saturate at log N where N is the batch size (van den Oord et al., 2018; Poole et al., 2019); the NWJ estimator suffers from large variance that grows at least exponentially in the ground-truth mutual information (Song & Ermon, 2019).
Please see Appendix B for more details and proofs.
On the other hand, V -information is explicit about the assumptions (as a feature instead of a bug).
V -information is also easy to estimate with guarantees if we can bound the complexity of V (such as its Radamacher or covering number complexity)
As we will show, bounds on the complexity of V directly translate to PAC (Valiant, 1984) bounds for V -information estimation.
In practice, we can efﬁciently optimize over V , e.g., via gradient descent.
In this paper we will present the Rademacher complexity version; other complexity measures (such as covering number) can be derived similarly.
Deﬁnition 4 (Empirical V -information) .
Let X, Y be two random variables taking values in X , Y and
D = { ( x i , y i ) } Ni =1 ∼ X, Y denotes the set of samples drawn from the joint distribution over X and
Y .
V is a predictive family.
The empirical V -information (under D ) is the following V -information under the empirical distribution deﬁned via D :
[Equation]
Then we have the following PAC bound over the empirical V -information:
Theorem 1.
Assume ∀ f ∈ V , x ∈ X , y ∈ Y , log f
[ x ]( y ) ∈ [ − B, B ] .
Then for any δ ∈ (0 , 0 . 5) , with probability at least 1 − 2 δ , we have:
[Equation]
where we deﬁne the function family G V
[Equation] [Equation]
log f
[ x ]( y ) , f ∈ V} , and R N ( G ) denotes the Rademacher complexity of G with sample number N .
Typically, the Rademacher complexity term satisﬁes R |D| ( G V ) =
O ( |D| − 12 ) (Bartlett & Mendelson, 2001; Gao & Zhou, 2016).
It’s worth noticing that a complex function family V (i.e., with large
Rademacher complexity) could lead to overﬁtting.
On the other hand, an overly-simple V may not be expressive enough to capture the relationship between X and Y .
As an example of the theorem, we provide a concrete estimation bound when V is chosen to be linear functions mapping X to the mean of a Gaussian distribution.
This was shown in Proposition 1 to lead to the coefﬁcient of determination.
Corollary 1.1.
Assume X = { x ∈ R d x , (cid:107) x (cid:107) 2 ≤ k x } and Y = { y ∈ R d y , (cid:107) y (cid:107) 2 ≤ k y } .
If
V = { f : f [ x ] =
N
( Wx + b, I ) , f [ ∅ ] =
N ( c, I ) , W ∈ R d y × d x , b, c ∈ R d y , (cid:107) ( W, b ) (cid:107) 2 ≤ 1 }
Denote M =
( k x + k y ) 2 + log 2 π , then ∀ δ ∈ (0 , 0 . 5) , with probability at least 1 − 2 δ : (cid:12)(cid:12)(cid:12) I
V ( X → Y ) − ˆ
I V
( X → Y ; D ) (cid:12)(cid:12)(cid:12) ≤
M (cid:112) 4 |D| (cid:32)
+(cid:114) 2 log 1 δ (cid:33)
Similar results can be obtained using other classes of machine learning models with known (Rademacher) complexity.
S TRUCTURE LEARNING WITH V - INFORMATION
Among many possible applications of V -information, we show how to use it to perform structure learning with provable guarantees.
The goal of structure learning is to learn a directed graphical model (Bayesian network) or undirected graphical model (Markov network) that best captures the (conditional) independence structure of an underlying data generating process.
Structure learning is difﬁcult in general, but if we restrict ourselves to certain set of graphs G , there are efﬁcient algorithms.
In particular, the Chow-Liu algorithm (Chow & Liu, 1968) can efﬁciently learn tree graphs (i.e. G is the set of trees).
Chow & Liu (1968) show that the problem can be reduced to:
[Equation]
where I ( X i , X j ) is the Shannon mutual information between variables
X i and X j .
In other words, it sufﬁces to construct the maximal weighted spanning tree where the weight between two vertices is their Shannon mutual information.
Chow & Wagner (1973) show that the Chow-Liu algorithm is consistent,
i.e, it recovers the true solution as the dataset size goes to inﬁnity.
However, the ﬁnite sample behavior of the Chow-Liu algorithm for high dimensional problems is much less studied, due to the difﬁculty of estimating mutual information.
In fact, we show in our experiments that the empirical performance is often poor, even with state-of-the-art estimators.
Additionally, methods based on mutual information cannot take advantage of intrinsically asymmetric relationships, which are common for example in gene regulatory networks (Meyer et al., 2007).
To address these issues, we propose a new structure learning algorithm based on V -information instead of Shannon information.
The idea is that we can associate to each directed edge in G (i.e., each pair of variables) a suitable predictive family V i,j (cf.
Def 1).
The main challenge is that we cannot simply replace mutual information with V -information in Eq. 6 because V -information is asymmetric – we now have to optimize over directed trees:
[Equation]
where G d − tree is the set of directed trees, and t ( g ) :
N → N is the function mapping each non-root node of directed tree g to its parent, and V i,j is the predictive family for random variables
X i and X j .
After estimating V -information on each edge, we use the Chu-Liu algorithm (Chu & Liu, 1965) to construct the maximal directed spanning tree .
This allows us to solve (7) exactly, even though there is a combinatorially large number of trees to consider.
Pseudocode is summarized in Algorithm 1 in
Appendix.
Denote C ( g ) =
(cid:80) mi
=
I V t ( g )( i ) ,i ( X t ( g )( i ) →
X i ) , we show in the following theorem that unlike the original Chow-Liu algorithm, our algorithm has guarantees in the ﬁnite samples regime, even in continuous settings:
Theorem 2.
Let { X i } mi =1 be the set of m random variables, D i,j (resp.
D j ) be the set of samples drawn from P ( X i , X j ) (resp.
P ( X j ) ).
Denote the optimal directed tree with maximum expected edge weights sum C ( g ) as g ∗ and the optimal directed tree constructed on the dataset D as ˆ g .
Then with the assumption in theorem 1, for any δ ∈ (0 , 1 2 m ( m − 1) ) , with probability at least 1 − 2 m ( m − 1) δ , we have:
[Equation]
Theorem 2 shows that the total edge weights of the maximal directed spanning tree constructed by algorithm 1 would be close to the optimal total edge weights if the Rademacher term is small.
Although larger C ( g ) does not necessarily lead to better Chow-Liu trees, empirically we ﬁnd that the optimal tree in the sense of equation (7) is consistent with the optimal tree in equation (6) under commonly used V .
E XPERIMENTAL RESULTS 6. S TRUCTURE LEARNING WITH CONTINUOUS HIGH - DIMENSIONAL DATA
We generate synthetic data using various ground-truth tree structures g ∗ with between 7 andvariables, where each variable is 10-dimensional.
We use Gaussians, Exponentials, and Uniforms as ground truth edge-conditionals.
We use V -information(Gaussian) and V -information(Logistic) to denote Algorithm 1 with two different V families.
Please refer to Appendix D.1 for more details.
We compare with the original Chow-Liu algorithm equipped with state-of-the-art mutual information estimators: CPC (van den Oord et al., 2018), NWJ (Nguyen et al., 2010) and MINE (Belghazi et al., 2018), with the same neural network architecture as the V -families for fair comparison.
All the experiments are repeated for 10 times.
As a performance metric, we use the wrong-edges-ratio (the ratio of edges that are different from ground truth) as a function of the amount of training data.
We show two illustrative experiments in ﬁgure 1a; please refer to Appendix D.1 for all simulations.
We can see that although the two V -families used are misspeciﬁed with respect to the true underlying (conditional) distributions , the estimated Chow-Liu trees are much more accurate across all data regimes, with CPC (blue) being the best alternative.
Surprisingly, V -information(Gaussian) works consistently well in all cases and only requires about 100 samples to recover the ground-truth
Chow-Liu tree in simulation-A.
[Figure]
[Figure]
AU C information
[Figure]
KDE Partitioning Mixed
[Figure]
inference 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 Frame Distance     - i n f o r m a ti on stochastic dynamic deterministic dynamic
[Figure]
6. G ENE REGULATORY NETWORK INFERENCE
Mutual information between pairs of gene expressions is often used to construct gene regulatory networks.
We evaluate V -information on the in-silico dataset from the DREAM5 challenge (Marbach et al., 2012) and use the setup of Gao et al.
(2017), where 20 genes with 660 datapoints are utilized to evaluate all methods.
We compare with state-of-the-art non-parametric Shannon mutual information estimators in this low dimensional setting: KDE , the traditional kernel density estimator; the KSG estimator (Kraskov et al., 2004); the Mixed KSG estimator (Gao et al., 2017) and Partitioning , an adaptive partitioning estimator (Darbellay & Vajda, 1999) implemented by Szab´o (2014).
For fair comparison with these low dimensional estimators, we select V =
{ f : f [ x ] = N ( g ( x ) , 12 ) , x ∈
X ; f [ ∅ ] =
N ( µ, 12 ) | µ ∈ range( g ) } , where g is a 3 -rd order polynomial.
The task is to predict whether a directed edge between genes exists in the ground-truth gene network.
We use the estimated mutual information and V -information for gene pairs as the test statistic to obtain the AUC for various methods.
As shown in Figure 1b, our method outperforms all other methods in network inference under different fractions of data used for estimation.
The natural information measure in this task is asymmetry since the goal is to ﬁnd the pairs of genes ( A i , B i )
s in which A i regulates B i , thus V -information is more suitable for such case than mutual information.
6. R ECOVERING THE ORDER OF VIDEO FRAMES
Let X 1 , · · · , X 20 be random variables each representing a frame in videos from the Moving-MNIST dataset, which contains 10,000 sequences each of length 20 showing two digits moving with stochastic dynamics.
Can Algorithm 1 be used to recover the natural (causal) order of the frames?
Intuitively, predictability should be inversely related with frame distance, thus enabling structure learning.
Using a conditional PixelCNN++ (Salimans et al., 2017) as predictive family V , we shown in Figure 1c that predictive V -information does indeed decrease with frame distance, despite some ﬂuctuations when the frame distances are large.
Using Algorithm 1 to construct a Chow-Liu tree, we ﬁnd that the tree perfectly recovers the relative order of the frames .
We also generate a Deterministic-Moving-MNIST dataset, where digits move according to deterministic dynamics.
From the perspective of Shannon mutual information, every pair of frames has the same mutual information.
Hence, standard Chow-Liu tree learning algorithm would fail to discover the natural ordering of the frames (causal structure).
In contrast, once we constrain the observer to PixelCNN++ models, algorithm 1 with predictive V -information can still recover the order of different frames when the frame distances are relatively small (less than 9).
Compared to the stochastic dynamics case, V -information is more irregular with increasing frame distance, since the
PixelCNN++ tends to overﬁt.
6. I NFORMATION THEORETIC APPROACHES TO FAIRNESS
The goal of fair representation learning is to map inputs X ∈ X to a feature space
Z ∈
Z such that the mutual information between Z and some sensitive attribute U ∈ U (such as race or gender) is minimized.
The motivation is that using Z (instead of X ) as input we can no longer use the sensitive attributes U to make decisions, thus ensuring some notion of fairness.
Existing methods obtain fair representations by optimizing against an “adversarial” discriminator so that the discriminator cannot predict U from Z (Edwards & Storkey, 2015; Louizos et al., 2015; Madras et al., 2018; Song et al., 2018).
Under some assumptions on U and V , we show in Appendix D.2 that these works actually use
V -information minimization as part of their objective, where V depends on the functional form of the discriminator.
However, it is clear from the V -information perspective that features trained with V A -information minimization might not generalize to V B -information and vice versa.
To illustrate this, we use a function family V j as the attacker to extract information from features trained with I
V i
( Z → U ) minimization, where all the V s are neural nets.
On three datasets commonly used in the fairness literature (Adult, German, Heritage), previous methods work well at preventing information “leak” against the class of adversary they’ve been trained on, but fail when we consider different ones.
As shown in Figure 3b in Appendix, the diagonal elements in the matrix are usually the smallest in rows, indicating that the attacker function family V i extracts more information on featured trained with V j ( j (cid:54)
= i )
-information minimization.
This challenges the generalizability of fair representations in previous works.
Please refer to Appendix D.2 for details.
R ELATED WORK
Alternative deﬁnitions of Information Several alternative deﬁnitions of mutual information are available in the literature.
Renyi entropy and Renyi mutual information (Lenzi et al., 2000) extend
Shannon information by replacing KL divergence with f -divergences.
However, they have the same difﬁculty when applied to high dimensional problems as Shannon information.
The line of work most related to ours is the H entropy and H mutual information (DeGroot et al., 1962; Gr¨unwald et al., 2004), which associate a deﬁnition of entropy to every prediction loss.
However, there are two key differences.
First, literatures in H entropy only consider a few special types of prediction functions that serve unique theoretical purposes; for example, (Duchi et al., 2018) considers the set of all functions on a feature space to prove surrogate risk consistency, and (Gr¨unwald et al., 2004) only considers the H entropy to prove the duality between maximum entropy and worst-case loss minimization.
In contrast, our deﬁnition takes a completely different perspective — emphasizing bounded computation and intuitive properties of “usable” information.
Furthermore H entropy still suffers from difﬁculty of estimation in high dimension because the deﬁnitions do not restrict to functions with small complexity (e.g. Rademacher complexity).
Mutual information estimation The estimation of mutual information in the machine learning ﬁeld is often on the continuous underlying distribution.
For non-parametric mutual information estimators, many methods have exploited the 3 H principle to calculate the mutual information, such as the Kernel density estimator (Paninski & Yajima, 2008), k-Nearest-Neighbor estimator and the
KSG estimator (Kraskov et al., 2004).
However, these non-parametric estimators usually aren’t scalable to high dimension.
Recently, several works utilize the variational lower bounds of MI to design MI estimator based on deep neural network in order to estimate MI of high dimension continuous random variables (Nguyen et al., 2010; van den Oord et al., 2018; Belghazi et al., 2018).
C ONCLUSION
We deﬁned and investigated V -information, a variational extension to classic mutual information that incorporates computational constraints.
Unlike Shannon mutual information, V -information attempts to capture usable information and has very different properties, such as invalidating the data processing inequality.
In addition, V -information can be provably estimated, and can thus be more effective for structure learning and fair representation learning.
A CKNOWLEDGEMENTS
This research was supported by AFOSR (FA9550-19-1-0024), NSF (#1651565, #1522054, #1733686),
ONR and FLI.
A P ROOFS A. P ROOF OF P ROPOSITION
Proposition 1.
For V -entropy and V -information, we have
Let Ω be as in Def. 1.
Then H Ω ( Y ) is the Shannon entropy, H Ω ( Y | X ) is the Shannon conditional entropy, and
I Ω ( Y → X ) is the Shannon mutual information.
Let Y =
R d and V = { f : { ∅ } → P µ | µ ∈ R d } , where P µ is the distribution with density
y (cid:55)→ 1 Z e −(cid:107) y − µ (cid:107) 2 where Z = (cid:82) e −(cid:107)
y − µ (cid:107) 2 dy , then the V -entropy of a random variable Y equals its mean absolute deviation, up to an additive constant.
Let Y = R d
and V = { f : { ∅ } → N ( µ, Σ)
| µ ∈ R d , Σ
= 1 / 2 I d × d } , then the V -entropy of a random variable Y equals the trace of its covariance tr
(Cov( Y )) , up to an additive constant.
Let V =
{ f : { ∅ } → Q t ,θ , θ ∈ Θ } , where Q t ,θ is a distribution in a minimal exponential family with sufﬁcient statistics t : Y → R d and set of natural parameters Θ .
For a random variable Y with expected sufﬁcient statistics µ Y =
E [ t ( Y )] , the V -entropy of Y is the maximum Shannon entropy over all random variables ˆ Y with identical expected sufﬁcient statistics, i.e. E
[ t ( ˆ Y )] = µ Y .
Let Y =
R d , X be any vector space, and V = { f : x (cid:55)→ N ( φ ( x ) , Σ) , x ∈ X ; ∅ (cid:55)→ N ( µ, Σ) | µ ∈ R d ; Σ = 1 / 2 I d × d , φ ∈ Φ } , where Φ is the set of linear functions { φ : X → R d } , then V -information
I V ( X →
Y ) equals the (unnormalized) maximum coefﬁcient of determination R 2 · tr (Cov( Y )) for linear regression.
Proof.
Let P Y | x denote the density function of random variable Y conditioned on X
= x (we denote this random variable as Y | x ).
[Equation]
where inﬁmum is achieved for f where f [ x ] = P Y | x and H is the Shannon (conditional) entropy.
The same proof technique can be used to show that H Ω ( Y )
=
H ( Y ) , with the inﬁmum achieved by f where f [ ∅ ] = P Y .
Hence we have
[Equation]
where MAD denotes mean absolute deviation E y ∼ Y
[ (cid:107)
y − E
[ Y ] (cid:107) 2 ] .
[Equation]
of
[Equation]
(4) The density function of an exponential family distribution with sufﬁcient statistics t is y (cid:55)→ exp ( θ · t ( y )
− A
θ ))
where A ( θ ) is the partition function.
[Equation]
where A ∗ is the Fenchel dual of the log-partition function A ( θ ) .
Under mild conditions (Wainwright et al., 2008)
[Equation]
where P µ is the maximum entropy distribution out of all distributions satisfying E y ∼ P µ [ t ( y )]
= µ (Jaynes, 1982), and H ( · ) is the Shannon entropy.
(5) Assume random variable Y ∈ R d , V = { f : x (cid:55)→ N ( φ ( x ) , Σ) , x ∈ X ; ∅ (cid:55)→ N ( µ, Σ) | µ ∈
R d ; Σ
=
I
d
× d ; φ ∈ Φ } .
Then the V -information from X to Y is
I V
[Equation]
X → Y
[Equation]
A. P ROOF OF P ROPOSITION
Proposition 2.
Let Y and X be any random variables on Y and X , and V and U be any predictive families, then we have
Monotonicity :
If V ⊆ U , then H V ( Y ) ≥ H U ( Y ) , H V ( Y | X ) ≥ H U
( Y | X ) .
Non-Negativity : I V ( X → Y ) ≥ 0 .
Independence : If X is independent of Y , I V ( X → Y ) =
I V ( Y → X )
= 0 .
Proof.
[Equation] [Equation]
Denote V ∅ ⊂ V as the subset of f that satisfy f
[ x ] = f [ ∅ ] , ∀ x ∈ X .
[Equation]
Therefore
I V (
[Equation]
→ X
[Equation]
Denote V ∅ ⊂ V as the subset of f that satisfy f
[ x ] = f [ ∅ ] , ∀ x ∈ X .
[Equation]
Therefore I V ( Y → X ) =
H V ( Y ) − H V
( Y | X ) ≤ 0 .
Combined with the Proposition 2.2 that
I V ( X → Y ) must be non-negative, I V ( X → Y ) must be 0 .
A. P ROOF OF T HEOREM
Theorem 1.
Assume ∀ f ∈ V , x ∈ X , y ∈ Y , log f
[ x ]( y ) ∈ [ − B, B ] .
Then for any δ ∈ (0 , 0 . 5) , with probability at least 1 − 2 δ , we have:
[Equation]
where we deﬁne the function family G V
= { g | g ( x, y )
=
log f
[ x ]( y ) , f ∈ V} , and R N ( G ) denotes the Rademacher complexity of G with sample number N .
Before proving theorem 1, we introduce two lemmas.
Proofs for these Lemmas follow the same strategy as theorem 8 in Bartlett & Mendelson (2001):
Lemma 3.
Let X, Y be two random variables taking values in X , Y and D denotes the set of samples drawn from the joint distribution over X ×Y .
Assume ∀ f ∈ V , x ∈ X , y ∈ Y , log f
[ x ]( y ) ∈ [ − B, B ] .
Take ˆ f = arg min f ∈V
[Equation]
− log f
[ x i ]( y i ) , then ∀ δ ∈ (0 , 1) , with probability at least 1 − δ , we have:
[Equation]
1 δ |D| (16)
Proof.
We apply McDiarmid’s inequality to the function Φ deﬁned for any sample D by
[Equation]
Let D and D (cid:48) be two samples differing by exactly one point, then since the difference of suprema does not exceed the supremum of the difference and ∀ f ∈ V , x ∈ X , y ∈ Y , log f
[ x ]( y ) ∈ [ − B, B ] , we have:
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
≤
[Equation]
|D| then by McDiarmid’s inequality
[Equation]
for any δ ∈ (0 , 1)
[Equation]
with probability at least 1 − δ
[Equation]
the following holds:
[Equation]
Then we bound the E D
[Φ( D )] term:
[Equation] [Equation]
[Equation]
  sup f ∈V E D (cid:48) (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) 1 |D (cid:48) | (cid:88) x (cid:48)
i ,y (cid:48)
i ∈D (cid:48) log f
[ x (cid:48)
i ]( y (cid:48) i ) | − 1 |D| (cid:88) x
i ,y i ∈D log f
[ x i ]( y i ) (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) (21) ≤ E D , D (cid:48)   sup f ∈V (cid:12) (cid:12)(cid:12) (cid:12)(cid:12)(cid:12) |D (cid:48) | (cid:88) x (cid:48)
i ,y (cid:48)
i ∈D (cid:48) log f
[ x (cid:48)
i ]( y (cid:48) i ) | − 1 |D| (cid:88) x
i ,y i ∈D log f
[ x i ]( y i ) (cid:12) (cid:12)(cid:12) (cid:12)(cid:12)(cid:12) 
 (22) =
E D , D (cid:48)   sup f ∈V (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) |D| |D| (cid:88) i = (log f
[ x (cid:48)
i ]( y (cid:48) i )
− log f
[ x i ]( y i ))
(cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) (23) ≤ E D , D (cid:48) ,σ   sup f ∈V (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) |D| |D| (cid:88)
i =
σ i (log f
[ x (cid:48)
i ]( y (cid:48) i )
− log f
[ x i ]( y i ))
(cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) (24) ≤ E D ,σ   sup f ∈V (cid:12)(cid:12) (cid:12) (cid:12) (cid:12)(cid:12) |D| |D| (cid:88) i
=
σ i log f
[ x i ]( y i ) (cid:12)(cid:12) (cid:12) (cid:12) (cid:12)(cid:12) 
 + E D (cid:48) ,σ   sup f ∈V (cid:12)(cid:12) (cid:12) (cid:12) (cid:12)(cid:12) |D| |D| (cid:88) i =
σ i log f
[ x (cid:48)
i ]( y (cid:48) i ) (cid:12)(cid:12) (cid:12) (cid:12) (cid:12)(cid:12) 
 (25) = 2 E D ,σ   sup f ∈V (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) |D| |D| (cid:88) i
=
σ i log f
[ x i ]( y i ) (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12) (26)
[Equation]
where σ i s are Rademacher variables that is uniform in {− 1 , +1 } .
Inequality (22) follows from the convexity of sup , inequality (24) follows from the symmetrization argument for (cid:96) 1 norm for
Radermacher random variables (Ledoux & Talagrand (2013), Section 6.1), inequality (21) follows from the convexity of | x − c | .
(27) follows from the deﬁnition of G and Rademacher complexity.
Finally, combining inequality (18) and (27) yields for all f ∈ V , with probability at least 1 − δ
[Equation]
In particular, the inequality holds for ˆ f
[Equation]
x
[Equation] [Equation]
− log f
[ x i ]( y i ) and ˜ f = arg min f ∈V
E x,y ∼ X,Y [ − log f
[ x ]( y )] .
Then we have:
[Equation]
Hence the bound (16) holds.
Similar bounds can be derived for H V ( Y ) when we choose the domain of x to be X = { ∅ } :
Lemma 4.
Let Y be random variable taking values in Y and D denotes the set of samples drawn from the underlying distribution P ( Y ) .
Assume ∀ f ∈ V , y ∈ Y , log f [ ∅ ]( y ) ∈ [ − B, B ] .
Take ˆ f = arg min
[Equation]
− log f
[ ∅ ]( y i ) , then for any δ ∈ (0 , 1) , with probability at least 1 − δ , we have:
[Equation]
where G V ∅ = { g | g ( y ) = log f
[ ∅ ]( y ) , f ∈ V} .
Proof.
The ﬁrst inequality (29) can be derived similarly as Lemma 3.
Since V is a predictive family, hence there exits a function h : V → V , such that h ( f ) =
f (cid:48) and ∀ x ∈ X , f (cid:48)
[ x ] = f [ ∅ ] .
[Equation]
[Equation]
The inequality (31) holds because of h ( V ) ⊆ V .
Now we prove theorem 1:
Theorem 1.
Assume ∀ f ∈ V , x ∈ X , y ∈ Y , log f
[ x ]( y ) ∈ [ − B, B ] , for any δ ∈ (0 , 0 . 5) , with probability at least 1 − 2 δ , we have:
[Equation]
Proof.
Deﬁne ˆ f = arg
min
f ∈V
[Equation]
i
[Equation]
i
[Equation] [Equation]
and ˆ f
[Equation]
= arg
[Equation]
∈D − log f [ ∅ ]( y i ) .
Using the triangular inequality we have:
[Equation]
For simplicity let
[Equation]
and
[Equation]
With inequality (32), Lemma 3 and Lemma 4, we have:
[Equation] [Equation]
Union bound)
[Equation]
(Lemma 3 and Lemma 4)
Hence we have:
[Equation]
which completes the proof.
A. P ROOF OF C OROLLARY 1.
Corollary 1.1.
Assume X = { x ∈ R d x , (cid:107) x (cid:107) 2 ≤ k x } and Y = { y ∈ R d y , (cid:107) y (cid:107) 2 ≤ k y } .
If
V = { f : f [ x ] =
N
( Wx + b, I ) , f [ ∅ ] =
N ( c, I ) , W ∈ R d y × d x , b, c ∈ R d y , (cid:107) ( W, b ) (cid:107) 2 ≤ 1 }
Denote M =
( k x + k y ) 2 + log 2 π , then ∀ δ ∈ (0 , 0 . 5) , with probability at least 1 − 2 δ :
[Equation]
The proof is an adaptation of the proof for theorem 3 in Kakade et al.
Proof.
From theorem 1 we have:
[Equation]
In the following (cid:107) ( W, b ) (cid:107) 2 is the matrix 2-norm of ( W, b ) , then the Rademacher term can be bounded as follows:
[Equation]
(cid:12)(cid:12)(cid:12)(cid:12)(cid:12)(cid:12)
[Equation]
The second term in RHS can be bounded as follows:
[Equation] [Equation]
The ﬁrst term in RHS can be bounded as follows:
[Equation]
[Equation]
The inequalities (36) and (35) follow the same proof in (34).
Hence we have:
[Equation]
In this example
[Equation]
we can bound the upper bound of functions g ∈ G V by
[Equation]
Combining inequality (38) we arrive at the theorem.
A. P ROOF OF T HEOREM
Theorem 2.
Let { X i } mi =1 be the set of m random variables, D i,j (resp.
D j ) be the set of samples drawn from P ( X i , X j ) (resp.
P ( X j ) ).
Denote the optimal directed tree with maximum expected edge weights sum C ( g ) as g ∗ and the optimal directed tree constructed on the dataset D as ˆ g .
Then with the assumption in theorem 1, for any δ ∈ (0 , 1 2 m ( m − 1) ) , with probability at least 1 − 2 m ( m − 1) δ , we have:
[Equation]
Proof.
Let C D ( g ∗ ) be the estimated sum of edge weights on dataset D of the tree g ∗ ,
i
[Equation]
e.,
[Equation]
i
[Equation]
D ) .
where t ( g ) : N → N is the function mapping each non-root node of directed tree g to its parent.
The same notation for tree ˆ g .
Let
[Equation]
be the maximum absolute estimation error of single edge weight.
By the deﬁnition of (cid:15) we have ∀ g, | C (ˆ g ) − C D
(ˆ g ) | ≤ ( m − 1) (cid:15) , then:
[Equation] [Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
[Table]
[Equation]
∨ 
[Equation]
≤ m
[Equation]
m − 1)2 δ
(By lemma 3, 4 and union bound)
Hence
Pr
[Equation]
Then combining inequality (39) and (40) we arrive at the result.
B A NALYSIS OF APPROXIMATE ESTIMATORS FOR S HANNON INFORMATION
We consider two approximate estimators for Shannon information.
The ﬁrst is the CPC (or InfoNCE in Poole et al.
(2019)) estimator ( I CPC ) proposed by van den Oord et al.
[Equation]
where the expectation is over N independent samples form the joint distribution (cid:81)
i p ( x i , y i ) .
The second is the NWJ estimator ( I NWJ ) proposed by Nguyen et al.
[Equation]
In both cases, f θ is a parameterized function, and the objectives are to maximize these lower bounds parameterized by θ to approximate mutual information.
Ideally, with sufﬁciently ﬂexible models and data, we would be able recover the true mutual information.
However, these ideal cases does not carry over to practical scenarios.
For I CPC , van den Oord et al.
(2018) show that I CPC is no larger than log N , where N is the batch size.
This means that the I CPC estimator will incur large bias when I ( X ; Y ) ≥ log N .
We provide a proof for completeness as follows.
Proposition 3.
∀ f
θ : X × Y
→ R + ,
I CPC ≤ log N.
Proof.
We have:
[Equation]
which completes the proof.
For NWJ, we note that the I NWJ involves a term denoted as E x,y ∼ p ( x ) p ( y ) (cid:2)
e f θ ( x,y ) (cid:3) /e , which could be dominated by rare data-points that have high f θ values.
Intuitively, this would make it a poor mutual information estimator by optimizing θ .
The NWJ estimator may suffer from high variance when the estimator is optimal (Song & Ermon, 2019), this is also empirically observed in Poole et al.
We provide a proof for completeness as follows.
Proposition 4.
Assume that f θ achieves the optimum value for I NWJ .
Then the variance of the empirical NWJ estimator satisﬁes Var (cid:16) ˆ I NWJ (cid:17) ≥
e I ( X ; Y ) − 1 N , where
[Equation]
is the empirical NWJ estimator with N i.i.d.
samples { ( x i , y i ) } Ni =1 from p ( x, y ) and N i.i.d.
samples { (¯ x i , ¯ y i ) } Ni =1 from p ( x ) p ( y ) .
Proof.
Let us denote z i
= p ( x i ,y i ) p ( x i ) p ( y i ) .
Clearly E p ( x ) p ( y )
[ z i ] = 1 .
Then we have:
[Equation]
where we use Jensen’s inequality for log at the last step.
From Nguyen et al.
(2010), we have:
[Equation] [Equation]
for all x, y .
Since { ( x i , y i ) } Ni =1 (resp.
{ (¯ x i , ¯ y i ) } Ni =1 ) are N datapoints independently sampled from the distribution p ( x, y )
(resp.
p ( x ) p ( y ) ), we have
[Equation]
which completes the proof
[Equation]
Require:
D = { ˆ X i } mi =1 , with each ˆ X i being a set of datapoints sampled from the underlying distribution of random variable X i .
The set of function families {V i,j } mi,j =1 ,i (cid:54) =
j between all the nodes.
1: for i = 1 , . . .
, m do 2: for j = 1 , . . .
, m do 3: if i (cid:54) =
j then 4:
Calculate the edge weight: e i → j = ˆ I V i,j ( X i → X j ; { ˆ X i , ˆ X j } ) .
5: end if 6: end for 7: end for 8: Construct the fully connected graph G =
( V, E ) , with node set V =
( X 1 , . . .
, X m )
and edge set E = { e i → j } mi,j =1 ,i (cid:54) =
j .
9: Construct the maximal directed spanning tree g on G by Chow-Liu algorithm, where mutual information is replaced by V -information.
10: return g
C T HE NEW ALGORITHM FOR C HU -L IU TREE CONSTRUCTION
See Algorithm 1; ˆ I V i,j ( X i → X j ; { ˆ X i , ˆ X j } ) denotes the empirical V -information.
D D ETAILED E XPERIMENTS SETUP D. C HU -L IU TREE CONSTRUCTION
Figure 2 shows the Chu-Liu tree construction of Simulation-1 ∼ Simulation-6.
The Simulation-A and
Simulation-B in the main body correspond to Simulation-1 and Simulation-4.
Simulation-
[Equation]
Simulation
[Equation]
The ground-truth Chu-Liu tree is a star tree (i.e. all random variables are conditionally independent given X 1 ).
We conduct all experiments for 10 times, each time with random simulated orthogonal matrices { W i } 20 i =2 .
Simulation-1:
X 1 ∼ U (0 , 10) and X i | X 1 ∼ N ( W i X 1 , 6 I ) , (2 ≤ i ≤ 20) ;
Simulation-2: X 1 ∼ U (0 , 10) and X i | X 1 ∼ W i E
( X 1 + (cid:15) i ) , (2 ≤ i ≤ 20) , (cid:15) i ∼ E (0 . 1) ;
Simulation3 is a mixed version: X 1 ∼ U (0 , 10) , X i | X 1 ∼ 12 N ( W i X 1 , 6 I ) + 12 W i E
( X 1 + (cid:15) 1 ) , (2 ≤ i ≤ 20) .
Simulation
-4 ∼ Simulation-
The ground-truth Chu-Liu tree is a tree of depth two.
We conduct all experiments for 10 times, each time with random simulated orthogonal matrices { W i } 7 i =2 .
Simulation-4: X 1 ∼ U (0 , 10) , X i | X 1 ∼
N ( W i X 1 , 2 I )( i = 2 , 3) , X i | X 2 ∼ N
( W i X 2 , 2 I )
( i = 4 , 5) , X i | X 3 ∼ N
( W i X 3 , 2 I )( i = 6 , 7) ;
Simulation-5: X 1 ∼ U (0 , 10) , X i | X 1 ∼ E ( X 1 + (cid:15) i )( i = 2 , 3) , X i | X 2 ∼ W i E
( X 2 + (cid:15) i )
( i = 4 , 5) ,
X i | X 3 ∼ W i E
( X 3 + (cid:15) i )
( i = 6 , 7) , (cid:15) i ∼ E (0 . 1) ; Simulation-6 is a mixed version: X 1 ∼
U (0 , 10) , X i | X 1 ∼ W i E
( X 1 + (cid:15) i )
( i = 2 , 3) , X i | X 2 ∼ N
( W i X 2 , 2 I )
( i = 4 , 5) , X i | X 3 ∼
N ( W i X 3 , 2 I )( i = 6 , 7) , (cid:15) i ∼ E (0 . 1) .
D. F AIRNESS
We can adapt the V -information perspective to fairness.
Denote the random variable that represents sensitive data and the representation as U and Z respectively.
Assume U is discrete and V belongs to preditive family 1.
Then we have H V ( U )
=
H ( U ) as long as V has softmax on the top and belongs to predictive family.
In this case, minimizing I V ( Z → U ) equals to minimize − H V ( Y | X ) .
Let the joint distribution of Z and U be paramterized by φ .
Hence the ﬁnal objective is:
[Equation]
In Edwards & Storkey (2015); Madras et al.
; Louizos et al.
; Song et al.
(2018), functions in V are parameterized by a discriminator.
For the ( F i , F j ) elements described in the main body, please refer to ﬁgure 3b.
The three datasets are: the UCI Adult dataset 2 which has gender as the sensitive attribute; the UCI German credit dataset 3 which has age as the sensitive attribute and the Heritage Health dataset 4 which has theconﬁgurations of ages and gender as the sensitive attribute.
The models in the ﬁgure are:
[Equation]
D is the training set.
V B = { f : f [ z ] = softmax( g ( z )) } , where g is a two-layer MLP with Relu as the activation function.
V C =
{ f : f [ z ] = softmax( g ( z )) } , where g is a three-layer MLP with LeakyRelu as the activation function.
We further visualize a special case of the ( V A , V B ) pair in ﬁgure 3a, where the V i
= { f : Z →
P ( U ) | f
[ z ]( u ) =
[Equation]
∈ R } explicitly makes the features of different sensitivity attributes more evenly spread, and functions in V B is a simple two layers MLP with softmax at the top.
The leaned features by V A -information minimization appear more evenly spread as expected, however, the attacker functions in V B can still achieve a high AUC of 0 . 857 .
The ( i, j ) elements of tables in Figure 3b stand for using function family V i to attack features trained with V j -information minimization.
The diagonal elements in the matrix are usually the smallest in rows, indicating that the attacker function family V i extracts more information on featured trained with V j ( j (cid:54) =
i ) -information minimization.
2 https://archive.ics.uci.edu/ml/datasets/adulthttps://archive.ics.uci.edu/ml/datasetshttps://www.kaggle.com/c/hhp
E M INIMALITY OF P REDICTIVE F AMILY
Deﬁne V X→P
( Y ) =
{ g : X → P ( Y ) |∃ f ∈ V , ∀ x ∈ X , g [ x ] = f
[ x ] } .
Similarly deﬁne V ∅ →P
( Y ) =
{ g : ∅ → P ( Y ) |∃ f ∈ V , g [ ∅ ] = f [ ∅ ] } .
Intuitively, V X→P ( Y ) (resp.
V ∅ →P
( Y ) ) restricts the domain of functions in V to X (resp. ∅ ).
Non-Negativity As we demonstrated in Proposition 2, optional-ignorance guarantees that information will be non-negative for any X and Y .
Conversely, given any discrete X , Z , V ∅ →P
( Y ) , V X→P ( Y ) that does not satisfy optional-ignorance, there exists distribution X , Y such that I V ( X → Y ) < 0 .
Choose Y ∼ f ∗ [ ∅ ] where f ∗ is the function that has no corresponding g ∈ V
X→P ( Y ) that can ignore its inputs.
Pick X as the uniform distribution, and note that for all g ∈ G , there exists some measurable subset X (cid:48)
⊂ X on which g will produce a distribution unequal to f ∗ [ ∅ ] , and therefore having higher cross entropy.
The expected cross entropy expressed in H V
X→P ( Y )
( Y | X ) is thus higher than in H V ∅ →P
( Y ) ( Y ) , and I V ( X → Y ) < 0 .
Thus, if the function class does not satisfy optional ignorance, then the V -information could be negative.
Independence Given any discrete X , Y , V ∅ →P
( Y ) , V X→P ( Y ) that does not satisfy optionalignorance, there exists an independent X , Y such that I V ( X → Y ) > 0 .
Choose Y such that the distribution P Y can be expressed as g [ x ] for some x ∈ X,
g ∈ V
X→P ( Y ) , but cannot be expressed by any f ∈ V ∅ →P
( Y ) .
Let X be the distribution with all its mass on x ; note that the cross entropy of
P Y with g [ x ] will be zero, and is less than that of the function f
[ ∅ ] (because f [ ∅ ] and P Y differs on a measurable subset, the cross entropy will be positive).
Thus, if the function class does not satisfy optional ignorance, then the V -information does not take value
when the two distributions are independent.
F L IMITATIONS AND F UTURE W ORK
V -information is empirically useful, has several intuitive theoretical properties, but exhibits certain limitations.
For example, Shannon information can be manipulated with certain additive algebra (e.g. H ( X, Y ) =
H ( X ) + H ( Y | X ) ), while the same does not hold true for general V -Information.
However, this could be possible if we choose V to be a mathematically simple set, such as the set of polynomial time computable functions.
It would be interesting to ﬁnd special classes of
V -Information where additional theoretical development is possible.
Another interesting direction is better integration of V -Information with machine learning.
The production of usable information (representation learning), acquisition of usable information (active learning) and exploitation of usable information (classiﬁcation and reinforcement learning) could potentially be framed in a similar V -information-theoretic manner.
It is interesting to see whether fruitful theories can arise from these analyses.