Chapter 1.6 : Information Theory
Chapter 1.6 : Information Theory
Christopher M. Bishop, PRML, Chapter 1 Introdcution
1. Information h(x)
Given a random variable and we ask how much information is received when we observe a specific value for this variable.
- The amount of information can be viewed as the “degree of surprise” on learning the value of .
- information : where the negative sign ensures that information is positive or zero.
- the units of :
- using logarithms to the base of 2: the units of are bits (‘binary digits’).
- using logarithms to the base of , i.e., natural logarithms: the units of are nats.
2. Entropy H(x): average amount of information
2.1 Entropy H(x)
Firstly we interpret the concept of entropy in terms of the average amount of information needed to specify the state of a random variable.
Now suppose that a sender wishes to transmit the value of a random variable to a receiver. The average amount of information that they transmit in the process is obtained by taking the expectation of (1.92) with respect to the distribution and is given as
- discrete entropy for discrete random variable by
- or differential/continuous entropy for continuous random variable by
- Note that and so we shall take whenever we encounter a value for such that .
- The nonuniform distribution has a smaller entropy than the uniform one.
2.2 Noiseless coding theorem (Shannon, 1948)
The noiseless coding theorem states that the entropy is a lower bound on the number of bits needed to transmit the state of a random variable.
2.3 Alternative view of entropy H(x)
Secondly, let us introduces the concept of entropy in physics in the context of equilibrium thermodynamics and later given a deeper interpretation as a measure of disorder through developments in statistical mechanics.
Consider a set of identical objects that are to be divided amongst a set of bins, such that there are objects in the bin. Consider the number of different ways of allocating the objects to the bins.
- There are ways to choose the first object, ways to choose the second object, and so on, leading to a total of ways to allocate all objects to the bins.
- However, we don’t wish to distinguish between rearrangements of objects within each bin. In the bin there are ways of reordering the objects, and so the total number of ways of allocating the objects to the bins is given by which is called the multiplicity.
- The entropy is then defined as the logarithm of the multiplicity scaled by an appropriate constant
- We now consider the limit , in which the fractions are held fixed, and apply Stirling’s approximation
- which gives
i.e.,
- where we have used , and is the probability of an object being assigned to the bin.
- microstate: In physics terminology, the specific arrangements of objects in the bins is called a microstate,
- macrostate: the overall distribution of occupation numbers, expressed through the ratios , is called macrostate.
- The multiplicity is also known as the weight of the macrostate.
- We can interpret the bins as the states of a discrete random variable , where . The entropy of the random variable X is then
2.4 Comparison between discrete entropy and continuous entropy
H(x) | Discrete Distribution X | Continuous Distribution X |
---|---|---|
Maximum | Uniform X | Gaussian X |
+/- | could be negative |
- Maximum entropy H(x) :
- In the case of discrete distributions, the maximum entropy configuration corresponded to an equal distribution of probabilities across the possible states of the variable.
- For a continuous variable, the distribution that maximizes the differential entropy is the Gaussian [see Page 54 in PRML].
- Is Negative or Positive?
- The discrete entropy in (1.93) is always , because . It will equal its minimum value of when one of the and all other .
- The differential entropy can be negative, because in (1.110) for .
If we evaluate the differential entropy of the Gaussian, we obtain This result also shows that the differential entropy, unlike the discrete entropy, can be negative, because in (1.110) for .
2.5 Conditional entropy H(y|x)
- Conditional entropy:
Suppose we have a joint distribution from which we draw pairs of values of and . If a value of is already known, then the additional information needed to specify the corresponding value of is given by . Thus the average additional information needed to specify y can be written as
which is called the conditional entropy of y given x. - It is easily seen, using the product rule, that the conditional entropy satisfies the relation where is the differential entropy (i.e., continuous entropy) of , and is the differential entropy of the marginal distribution .
- From (1.112) we get to know that
the information needed to describe and is given by the sum of the information needed to describe alone plus the additional information required to specify given .
3. Relative entropy or KL divergence
3.1 Relative entropy or KL divergence
Problem: How to relate the notion of entropy to pattern recognition?
Description: Consider some unknown distribution , and suppose that we have modeled this using an approximating distribution .
-
Relative entropy or Kullback-Leibler divergence, or KL divergence between the distributions and : If we use to construct a coding scheme for the purpose of transmitting values of to a receiver, then the average additional amount of information (in nats) required to specify the value of (assuming we choose an efficient coding scheme) as a result of using instead of the true distribution is given by
and
It can be rewritten as [see Ref-1]
-
Cross Entropy: where is called the cross entropy,
- Understanding of Cross Entropy: One can show that the cross entropy is the average number of bits (or nats) needed to encode data coming from a source with distribution when we use model to define our codebook.
- Understanding of (“Regular”) Entropy: Hence the “regular” entropy , is the expected number of bits if we use the true model.
- Understanding of Relative Entropy: So the KL divergence is the difference between these (shown in 2.111). In other words, the KL divergence is the average number of extra bits (or nats) needed to encode the data, due to the fact that we used approximation distribution to encode the data instead of the true distribution .
Asymmetric: Note that KL divergence is not a symmetrical quantity, that is to say .
- KL divergence is a way to measure the dissimilarity of two probability distributions, and [see Ref-1].
3.2 Information inequality [see Ref-1]
The “extra number of bits” interpretation should make it clear that , and that the KL is only equal to zero i.f.f. . We now give a proof of this important result.
Proof:
- 1) Convex functions: To do this we first introduce the concept of convex functions. A function is said to be convex if it has the property that every chord lies on or above the function, as shown in Figure 1.31.
- Convexity then implies
- 2) Jensen’s inequality:
- Using the technique of proof by induction(数学归纳法), we can show from (1.114) that a convex function satisfies where and , for any set of points . The result (1.115) is known as Jensen’s inequality.
- If we interpret the as the probability distribution over a discrete variable taking the values , then (1.115) can be written For continuous variables, Jensen’s inequality takes the form
- 3) Apply Jensen’s inequality in the form (1.117) to the KL divergence (1.113) to give where we have used the fact that is a convex function (In fact, is a strictly convex function, so the equality will hold if, and only if, for all ), together with the normalization condition .
- 4) Similarly, Let be the support of , and apply Jensen’s inequality in the form (1.115) to the discrete form KL divergence (2.110) to get [see Ref-1] where the first inequality follows from Jensen’s. Since is a strictly concave (i.e., the inverse of convex) function, we have equality in Equation (2.115) iff for some . We have equality in Equation (2.116) iff , which implies .
- 5) Hence iff for all .
3.3 How to use KL divergence
Note that:
- we can interpret the KL divergence as a measure of the dissimilarity of the two distributions and .
- If we use a distribution that is different from the true one, then we must necessarily have a less efficient coding, and on average the additional information that must be transmitted is (at least) equal to the Kullback-Leibler divergence between the two distributions.
Problem description:
- Suppose that data is being generated from an unknown distribution that we wish to model.
- We can try to approximate this distribution using some parametric distribution , governed by a set of adjustable parameters , for example a multivariate Gaussian.
- One way to determine is to minimize the KL divergence between and with respect to .
- We cannot do this directly because we don’t know . Suppose, however, that we have observed a finite set of training points , for , drawn from . Then the expectation with respect to can be approximated by a finite sum over these points, using (1.35), so that (???don’t know how to derive it)
- the first term is the negative log likelihood function for under the distribution evaluated using the training set.
- Thus we see that minimizing this KL divergence is equivalent to maximizing the likelihood function.
3.4 Mutual information
Now consider the joint distribution between two sets of variables and given by .
Mutual information between the variables and :
- If and are independent, p(x, y) = p(x)p(y).
- If and are not independent, we can gain some idea of whether they are “close” to being independent by considering the KL divergence between the joint distribution and the product of the marginals, given by which is called the mutual information between the variables and .
- Using the sum and product rules of probability, we see that the mutual information is related to the conditional entropy
through
Understanding of Mutual information:
- Thus we can view the mutual information as the reduction in the uncertainty about by virtue of being told the value of (or vice versa).
- From a Bayesian perspective, we can view as the prior distribution for and as the posterior distribution after we have observed new data . The mutual information therefore represents the reduction in uncertainty about as a consequence of the new observation .
Reference
[1]: Section 2.8.2, Page 57, Kevin P. Murphy. 2012. Machine Learning: A Probabilistic Perspective. The MIT Press.
CCJ PRML Study Note - Chapter 1.6 : Information Theory的更多相关文章
-
信息熵 Information Theory
信息论(Information Theory)是概率论与数理统计的一个分枝.用于信息处理.信息熵.通信系统.数据传输.率失真理论.密码学.信噪比.数据压缩和相关课题.本文主要罗列一些基于熵的概念及其意 ...
-
Tree - Information Theory
This will be a series of post about Tree model and relevant ensemble method, including but not limit ...
-
information entropy as a measure of the uncertainty in a message while essentially inventing the field of information theory
https://en.wikipedia.org/wiki/Claude_Shannon In 1948, the promised memorandum appeared as "A Ma ...
-
Better intuition for information theory
Better intuition for information theory 2019-12-01 21:21:33 Source: https://www.blackhc.net/blog/201 ...
-
信息论 | information theory | 信息度量 | information measures | R代码(一)
这个时代已经是多学科相互渗透的时代,纯粹的传统学科在没落,新兴的交叉学科在不断兴起. life science neurosciences statistics computer science in ...
-
Beginning Scala study note(5) Pattern Matching
The basic functional cornerstones of Scala: immutable data types, passing of functions as parameters ...
-
TIJ——Chapter Fourteen:Type Information
Runtime type information(RTTI) allows you to discover and use type information while a program is ru ...
-
Beginning Scala study note(9) Scala and Java Interoperability
1. Translating Java Classes to Scala Classes Example 1: # a class declaration in Java public class B ...
-
Beginning Scala study note(8) Scala Type System
1. Unified Type System Scala has a unified type system, enclosed by the type Any at the top of the h ...
随机推荐
-
APP性能测试
方法一: 本地安装安卓模拟器,用LR选择模拟器录制方式录制 方法二: 手机真机需要root,可以在电脑上下载一键root工具(如卓大师),然后手机和电脑用数据线连接,然后root. 在手机上运行 Mo ...
-
2014 UESTC暑前集训搜索专题解题报告
A.解救小Q BFS.每次到达一个状态时看是否是在传送阵的一点上,是则传送到另一点即可. 代码: #include <iostream> #include <cstdio> # ...
-
AndroidStudio创建新项目报错
创建新项目自动执行时报错: Failed to import new Gradle project: failed to find Build Tools revision 17.0.0 Consul ...
-
dip的全称
dip的全称是Density-independent pixel,密度无关像素.很多地方误认为是device independent pixel,即设备无关像素.这是错误的. 因为dip也就是dp只能 ...
-
18SpringMvc_在业务控制方法中收集数组参数
这篇文章我们要解决的问题的多选框选中,并批量删除. 比如:
-
JS学习笔记(一) 概述
参考资料: 1. http://www.w3school.com.cn/js/ 2. http://www.liaoxuefeng.com/wiki/001434446689867b27157e896 ...
-
UVA 10896 Sending Email
这个题目真是伤透脑筋了,一直RE,连着改了好几个版本,又是spfa,又是单调队列dijkstra+单调队列,总是不过,后来发现M开小了,双向边应该开m的两倍,悲剧啊!!!以后不管怎样,数组一定要尽量开 ...
-
SQL点滴25—T-SQL面试语句,练练手
原文:SQL点滴25-T-SQL面试语句,练练手 1. 用一条SQL语句查询出每门课都大于80分的学生姓名 name kecheng fenshu 张三 语文 81张三 ...
-
C++—动态内存管理之深入探究new和delete
C++中程序存储空间除栈空间和静态区外,每个程序还拥有一个内存池,这部分内存被称为*空间(free store)或堆(heap).程序用堆来存储动态分配的对象,即,那些程序运行时分配的对象.动态对象 ...
-
配置rpm包安装的jdk环境变量
最近在搭建james邮件服务的时候,由于这个服务是用Java开发的,之前这台服务器跑过tomcat服务,故有Java环境,就没在意有无配置环境变量,但在启动james的时候报没有配置环境变量: 那么问 ...