Teoria informacji/TI Wykład 5

Entropia zależna i informacja wzajemna

Definicja [Entropia zmiennej losowej]

Jeśli $X:S\to {𝒳}$ jest zmienną losową, określamy jej entropię jako

$H_r(X)=\sum _{t\in {𝒳}}p(X=t)\cdot {\log }_r\frac{1}{p(X=t)}$

Innymi słowy, $H_r(X)$ jest równe wartości oczekiwanej

$H_r(X)=E\left({\log }_r\frac{1}{p(X)}\right)$

gdzie p(X) jest zmienną losową na S zdefiniowaną jako $p(X):s\mapsto p(X=X(s))$

$\sum _{t\in {𝒳}}p(X=t)\cdot {\log }_r\frac{1}{p(X=t)}=\sum _{t\in {𝒳}}\sum _{X(s)=t}p(s)\cdot \frac{1}{p(X=t)}=\sum _{s\in S}p(s)\cdot \frac{1}{p(X=X(s))}$

Umowa notacyjna Jeśli zmienne losowe o których mowa będą wynikały z kontekstu, często będziemy omijać zapis $X=a$ i pisać po prostu a. Przykładowo będziemy pisać p(x|y) zamiast $p(X=x|Y=y)$, $p(x\wedge y)$ zamiast $p((X=x)\wedge (Y=y))$ itp.

Definicja [Entropia zależna]

Niech $A:S\to {𝒜}$,$B:S\to {ℬ}$ będą dwiema zmiennymi losowymi. Dla $b\in {ℬ}$ określamy

$H_r(A|b)=\sum _{a\in {𝒜}}p(a|b)\cdot {\log }_r\frac{1}{p(a|b)}$

i ogólnie

$H_r(A|B)=\sum _{b\in {ℬ}}p(b)H_r(A|b)$.

Powyższą wartość nazywamy entopią zależną A od B

Zauważmy że jeśli A i B są niezależne, to w powyższej formule $p(a|b)=a$ a więc $H_r(A|B)=A$. Z drugiej strony $H_r(A|A)=0$. Ogólnie dla dowolnej funkcji $\phi :{𝒜}\to {ℬ}$ mamy

$H_r(\phi (A)|A)=0$

Rzeczywiście, jeśli $p(A=a)>0$ to $p(\phi (A)=\phi (a)|A=a)=1$, i w konsekwencji ${\log }_r\frac{1}{p(\phi (A)=\phi (a)|A=a)}=0$.

Entropia produktowa Będziemy również rozważać pary (A,B) jako jedną zmienną losową $(A,B):S\to {𝒜}\times {ℬ}$,

$(A,B)(s)=(A(s),B(s))$

Prawdopodobieństwo że ta zmienna przyjmie wartość (a,b) wynosi $p((A,B)=(a,b))=p((A=a)\wedge (B=b))$, co zapisujemy w skrócie jako $p(a\wedge b)$. To prawdopodobieństwo w ogólności jest inne niż $p(a)\cdot p(b)$. Jeśli dla dowolnych $a\in {𝒜},b\in {ℬ}$

$p(a\wedge b)=p(a)\cdot p(b)$,

zmienne losowe A i B są niezależne.

Entropia $H_r(A,B)$ wprost z definicji wynosi

$H_r(A,B)=\sum _{a\in {𝒜},b\in {ℬ}}p(a\wedge b)\cdot {\log }_r\frac{1}{p(a\wedge b)}$

Jeśli A i B są niezależne, to

${\log }_r\frac{1}{p(A,B)}={\log }_r\frac{1}{p(A)}+{\log }_r\frac{1}{p(B)}$

Z liniowości wartości oczekiwanej dostajemy wtedy

$H_r(A,B)=H_r(A)+H_r(B)$

W ogólnym przypadku możemy udowodnić:

Twierdzenie

$H_r(A,B)\le H_r(A)+H_r(B)$

i równość zachodzi jedynie gdy A i B są niezależne.

Dowód Rozpiszemy prawą stronę tak żebyśmy mogli użyć Złotego Lematu. Użyjemy w tym celu oczywistych równości $p(a)=\sum _{b\in {ℬ}}p(a\wedge b)$ i $p(b)=\sum _{a\in {𝒜}}p(a\wedge b)$. $H_r(A)+H_r(B)=\sum _{a\in {𝒜}}p(a){\log }_r\frac{1}{p(a)}+\sum _{b\in {ℬ}}p(b){\log }_r\frac{1}{p(b)}$

$=\sum _{a\in {𝒜}}\sum _{b\in {ℬ}}p(a\wedge b){\log }_r\frac{1}{p(a)}+\sum _{b\in {ℬ}}\sum _{a\in {𝒜}}p(a\wedge b){\log }_r\frac{1}{p(b)}$

$=\sum _{a\in {𝒜},b\in {ℬ}}p(a\wedge b){\log }_r\frac{1}{p(a)p(b)}$

Ważne że powyższe wyrażenie jest dobrze zdefiniowane, bo gdy $p(a)=0$ lub $p(b)=0$, to również $p(a\wedge b)=0$.

Oznaczmy chwilowo

$({𝒜}\times {ℬ}{)}^{+}=\{(a,b):p(a)>0\text{ i }p(b)>0\}$

Mamy wtedy

$\sum _{(a,b)\in ({𝒜}\times {ℬ}{)}^{+}}p(a\wedge b)=\sum _{(a,b)\in ({𝒜}\times {ℬ}{)}^{+}}p(a)\cdot p(b)=1$.

Używając Złotego Lematu dla $x=p(a\wedge b)$, $y=p(a)\cdot p(b)$ dla wszystkich $(a,b)\in ({𝒜}\times {ℬ}{)}^{+}$ otrzymujemy $H_r(A,B)=\sum _{(a,b)\in ({𝒜}\times {ℬ}{)}^{+}}p(a\wedge b){\log }_r\frac{1}{p(a\wedge b)}$

$\le \sum _{(a,b)\in ({𝒜}\times {ℬ}{)}^{+}}p(a\wedge b){\log }_r\frac{1}{p(a)p(b)}$

$=H_r(A)+H_r(B)$

Dodatkowo równość zachodzi wyłącznie gdy $p(a\wedge b)=p(a)\cdot p(b)$ dla wszystkich $(a,b)\in ({𝒜}\times {ℬ}{)}^{+}$ (czyli w ogóle dla wszystkich $a\in {𝒜},b\in {ℬ}$. W drugą stronę, wiemy już że niezależność A i B implikuje tutaj równość. QED.

Definicja [Informacja]

Wartość

$I(A;B)=H_r(A)+H_r(B)-H_r(A,B)$

nazywamy informacją wzajemną zmiennych A i B.

Komentarz Powyższę definicję łatwo zrozumieć w odniesieniu do Gry w 20 pytań. Przypuścmy że mamy zidentyfikować obiekt który jest parą (a,b) gdzie a i b są wartościami zmiennych losowych A i B. Jeśli A i B są niezależne, najlepsze co możemy zrobić to zidentyfikować niezależnie a i b. Tym samym gramy w dwie niezależne gry „pytania o a” i „pytania o b” (co odpowiada równości $H_r(A,B)=H_r(A)+H_r(B)$). Jeśli jednak A i B są zależne, możemy wykorzystać tę wzajemną informację do zmniejszenia liczby pytań.

Dla zwiększenia czytelności tekstu, od tej poty będziemy zwykle omijać dolny indeks r, pisząc H, I, itp. Wszędzie tam gdzie nie napisano inaczej, wszystkie twierdzenia odnoszą się do przypadku dowolnego $r>1$. Bez utraty ogólności czytelnik może założyć r=2.

Komentarz Przekształcając definicję informacji analogicznie jak w ostatnim dowodzie, otrzymujemy:

$I(A;B)=\sum _{a\in {𝒜},b\in {ℬ}}p(a\wedge b)(\log \frac{1}{p(a)p(b)}-\log \frac{1}{p(a\wedge b)})$

W takiej postaci widać że informacja jest pewną miarą odległości pomiędzy faktycznym rozkładem zmiennej (A;B), a jej rozkładem gdyby A i B były niezależne.

Warto zauważyć że powyższa suma jest nieujemna, choć niektóre składniki $(\log \frac{1}{p(a)p(b)}-\log \frac{1}{p(a\wedge b)})$ mogą być ujemne.

Istnieje odpowiednik równości $H(A,B)=H(A)+H(B)$, który stosuje się do zmiennych zależnych:

Fakt [Zasada łańcuchowa]

$H(A,B)=H(A|B)+H(B)$

Dowód Obliczamy:

$H(A,B)=\sum _{a\in {𝒜},b\in {ℬ}}p(a\wedge b)\cdot \log \frac{1}{p(a\wedge b)}$

$=\sum _{a\in {𝒜}}\sum _{b\in {ℬ}}p(a|b)p(b)\cdot \log \frac{1}{p(a|b)p(b)}$

$=\sum _{a\in {𝒜}}\sum _{b\in {ℬ}}p(a|b)p(b)\cdot (\log \frac{1}{p(a|b)}+\log \frac{1}{p(b)})$

$=\sum _{b\in {ℬ}}p(b)\cdot \sum _{a\in {𝒜}}p(a|b)\cdot \log \frac{1}{p(a|b)}+\sum _{b\in {ℬ}}p(b)\log \frac{1}{p(b)}\cdot \sum _{a\in {𝒜}}p(a|b)$

$=H(A|B)+H(B)$

QED.

Używając zasady łańcuchowej, możemy wyliczać informację na różne sposoby: $I(A;B)=H(A)-H(A|B)=H(B)-H(B|A)$

Kolejną rzeczą jaką możemy zauważyć to $I(A;B)\le \min \{H(A),H(B)\}$

Łatwo możemy też uogólnić zasadę łańcuchową na przypadek $n\ge 2$ zmiennych $A_1,A_2,\dots ,A_n$

$H(A_1,\dots ,A_n)=H(A_1|A_2,\dots ,A_n)+H(A_2,\dots ,A_n)$

$=H(A_1|A_2,\dots ,A_n)+H(A_2|A_3,\dots ,A_n)+H(A_3,\dots ,A_n)$

$=\sum _{i=1}^{n}H(A_i|A_{i+1},\dots ,A_n)$

(przyjmujemy konwencję $H(A|\mathrm{\varnothing })=A$)

Bardziej wyrafinowane uogólnienie możemy uzyskać stosując entropię zależną:

Fakt [Zależna zasada łańcuchowa]

$H(A,B|C)=H(A|B,C)+H(B|C)$

Dowód Mamy:

$H(A,B|c)=\sum _{a\in {𝒜},b\in {ℬ}}p(a\wedge b|c)\cdot \log \frac{1}{p(a\wedge b|c)}$

$=\sum _{a,b}p(a|b\wedge c)\cdot p(b|c)\cdot (\log \frac{1}{p(a|b\wedge c)}+\log \frac{1}{p(b|c)})$

$=\sum _{b}p(b|c)\cdot \sum _{a}p(a|b\wedge c)\cdot \log \frac{1}{p(a|b\wedge c)}+\sum _{b}p(b|c)\cdot \log \frac{1}{p(b|c)}\cdot \underset{1}{\underbrace{\sum _{a}p(a|b\wedge c)}}$

W powyższym wyliczeniu sumy po a i b obejmują te wartości, dla których odpowiednie prawdopodobieństwa zależne są zdefiniowane (p(x|y) nie jest określone jeśli p(y)=0). Używamy tu łatwego faktu, że jeśli $p(a\wedge b|c)>0$, to

$p(a\wedge b|c)=\frac{p(a\wedge b\wedge c)}{p(c)}=\frac{p(a\wedge b\wedge c)}{p(b\wedge c)}\cdot \frac{p(b\wedge c)}{p(c)}=p(a|b\wedge c)\cdot p(b|c)$

Uśredniając po p(c) dostajemy:

$H(A,B|C)=\sum _{c\in {𝒞}}p(c)\cdot H(A,B|c)$

$=\sum _{c}p(c)\cdot \sum _{b}p(b|c)\cdot \sum _{a}p(a|b\wedge c)\cdot \log \frac{1}{p(a|b\wedge c)}+\sum _{c}p(c)\cdot \sum _{b}p(b|c)\cdot \log \frac{1}{p(b|c)}$

$=\underset{H(A|B,C)}{\underbrace{\sum _{b,c}p(b\wedge c)\cdot \sum _{a}p(a|b\wedge c)\cdot \log \frac{1}{p(a|b\wedge c)}}}+\underset{H(B|C)}{\underbrace{\sum _{c}p(c)\cdot \sum _{b}p(b|c)\cdot \log \frac{1}{p(b|c)}}}$

QED.

Czytelnik może teraz łatwo pokazać że:

$H(A,B|C)\le H(A|C)+H(B|C)$

i równość zachodzi wtedy i tylko wtedy gdy A i B są niezależne w odniesieniu do C, czyli

$p(A=a\wedge B=b|C=c)=p(A=a|C=c)\cdot p(B=b|C=c)$

(dowód na ćwiczeniach)

Definicja [Informacja zależna]

Definiujemy informację wzajemną A i B w odniesieniu do C jako

$I(A;B|C)=H(A|C)+H(B|C)-\underset{H(A|B,C)+H(B|C)}{\underbrace{H(A,B|C)}}$

$=H(A|C)-H(A|B,C)$

I wreszcie, informacją wzajemną A, B i C definiujemy jako:

$R(A;B;C)=I(A;B)-I(A;B)|C)$

Łatwo sprawdzimy że ta definicja jest rzeczywiście symetryczna, tzn nie zależy od kolejności A, B i C:

$I(A;C)-I(A;C|B)=H(A)-H(A|C)-(H(A|B)-H(A|B,C))$

$=\underset{I(A;B)}{\underbrace{H(A)-H(A|B)}}-\underset{I(A;B|C)}{\underbrace{H(A|C)-H(A|B,C)}}$

Należy jednak pamiętać że w przeciwieństwie do I(A;B) i I(A;B|C), zdefiniowana powyżej R(A;B;C) może mieć ujemną wartość.

Zależności pomiędzy wartościami H(X), H(Y), H(Z), H(X,Y), H(X,Y|Z), I(X;Y), I(X;Y|Z), R(X;Y;Z) itd. można przedstawić w postaci diagramu: