Skip to content

Commit ae31f1f

Browse files
authored
Merge pull request #176 from JoramSoch/master
added 5 proofs
2 parents dc5ac2e + aeeb3ca commit ae31f1f

6 files changed

Lines changed: 388 additions & 12 deletions

File tree

I/ToC.md

Lines changed: 17 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -134,9 +134,11 @@ title: "Table of Contents"
134134
&emsp;&ensp; 1.9.7. *[Covariance matrix](/D/covmat)* <br>
135135
&emsp;&ensp; 1.9.8. *[Sample covariance matrix](/D/covmat-samp)* <br>
136136
&emsp;&ensp; 1.9.9. **[Covariance matrix and expected values](/P/covmat-mean)** <br>
137-
&emsp;&ensp; 1.9.10. **[Covariance matrix and correlation matrix](/P/covmat-corrmat)** <br>
138-
&emsp;&ensp; 1.9.11. *[Precision matrix](/D/precmat)* <br>
139-
&emsp;&ensp; 1.9.12. **[Precision matrix and correlation matrix](/P/precmat-corrmat)** <br>
137+
&emsp;&ensp; 1.9.10. **[Invariance under addition of vector](/P/covmat-inv)** <br>
138+
&emsp;&ensp; 1.9.11. **[Scaling upon multiplication with matrix](/P/covmat-scal)** <br>
139+
&emsp;&ensp; 1.9.12. **[Covariance matrix and correlation matrix](/P/covmat-corrmat)** <br>
140+
&emsp;&ensp; 1.9.13. *[Precision matrix](/D/precmat)* <br>
141+
&emsp;&ensp; 1.9.14. **[Precision matrix and correlation matrix](/P/precmat-corrmat)** <br>
140142

141143
1.10. Correlation <br>
142144
&emsp;&ensp; 1.10.1. *[Definition](/D/corr)* <br>
@@ -468,10 +470,12 @@ title: "Table of Contents"
468470
&emsp;&ensp; 4.3.2. **[Special case of normal-Wishart distribution](/P/ng-nw)** <br>
469471
&emsp;&ensp; 4.3.3. **[Probability density function](/P/ng-pdf)** <br>
470472
&emsp;&ensp; 4.3.4. **[Mean](/P/ng-mean)** <br>
471-
&emsp;&ensp; 4.3.5. **[Differential entropy](/P/ng-dent)** <br>
472-
&emsp;&ensp; 4.3.6. **[Kullback-Leibler divergence](/P/ng-kl)** <br>
473-
&emsp;&ensp; 4.3.7. **[Marginal distributions](/P/ng-marg)** <br>
474-
&emsp;&ensp; 4.3.8. **[Conditional distributions](/P/ng-cond)** <br>
473+
&emsp;&ensp; 4.3.5. **[Covariance](/P/ng-cov)** <br>
474+
&emsp;&ensp; 4.3.6. **[Differential entropy](/P/ng-dent)** <br>
475+
&emsp;&ensp; 4.3.7. **[Kullback-Leibler divergence](/P/ng-kl)** <br>
476+
&emsp;&ensp; 4.3.8. **[Marginal distributions](/P/ng-marg)** <br>
477+
&emsp;&ensp; 4.3.9. **[Conditional distributions](/P/ng-cond)** <br>
478+
&emsp;&ensp; 4.3.10. **[Drawing samples](/P/ng-samp)** <br>
475479

476480
4.4. Dirichlet distribution <br>
477481
&emsp;&ensp; 4.4.1. *[Definition](/D/dir)* <br>
@@ -487,11 +491,12 @@ title: "Table of Contents"
487491
&emsp;&ensp; 5.1.3. **[Probability density function](/P/matn-pdf)** <br>
488492
&emsp;&ensp; 5.1.4. **[Mean](/P/matn-mean)** <br>
489493
&emsp;&ensp; 5.1.5. **[Covariance](/P/matn-cov)** <br>
490-
&emsp;&ensp; 5.1.6. **[Kullback-Leibler divergence](/P/matn-kl)** <br>
491-
&emsp;&ensp; 5.1.7. **[Transposition](/P/matn-trans)** <br>
492-
&emsp;&ensp; 5.1.8. **[Linear transformation](/P/matn-ltt)** <br>
493-
&emsp;&ensp; 5.1.9. **[Marginal distributions](/P/matn-marg)** <br>
494-
&emsp;&ensp; 5.1.10. **[Drawing samples](/P/matn-samp)** <br>
494+
&emsp;&ensp; 5.1.6. **[Differential entropy](/P/matn-dent)** <br>
495+
&emsp;&ensp; 5.1.7. **[Kullback-Leibler divergence](/P/matn-kl)** <br>
496+
&emsp;&ensp; 5.1.8. **[Transposition](/P/matn-trans)** <br>
497+
&emsp;&ensp; 5.1.9. **[Linear transformation](/P/matn-ltt)** <br>
498+
&emsp;&ensp; 5.1.10. **[Marginal distributions](/P/matn-marg)** <br>
499+
&emsp;&ensp; 5.1.11. **[Drawing samples](/P/matn-samp)** <br>
495500

496501
5.2. Wishart distribution <br>
497502
&emsp;&ensp; 5.2.1. *[Definition](/D/wish)* <br>

P/covmat-inv.md

Lines changed: 52 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,52 @@
1+
---
2+
layout: proof
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2022-09-22 11:29:00
9+
10+
title: "Invariance of the covariance matrix under addition of constant vector"
11+
chapter: "General Theorems"
12+
section: "Probability theory"
13+
topic: "Covariance"
14+
theorem: "Invariance under addition of vector"
15+
16+
sources:
17+
- authors: "Wikipedia"
18+
year: 2022
19+
title: "Covariance matrix"
20+
in: "Wikipedia, the free encyclopedia"
21+
pages: "retrieved on 2022-09-22"
22+
url: "https://en.wikipedia.org/wiki/Covariance_matrix#Basic_properties"
23+
24+
proof_id: "P347"
25+
shortcut: "covmat-inv"
26+
username: "JoramSoch"
27+
---
28+
29+
30+
**Theorem:** The [covariance matrix](/D/covmat) $\Sigma_{XX}$ of a [random vector](/D/rvec) $X$ is invariant under addition of a [constant vector](/D/const) $a$:
31+
32+
$$ \label{eq:covmat-inv}
33+
\Sigma(X+a) = \Sigma(X) \; .
34+
$$
35+
36+
37+
**Proof:** The [covariance matrix](/D/covmat) of $X$ [can be expressed](/P/covmat-mean) in terms of [expected values](/D/mean) as follows:
38+
39+
$$ \label{eq:covmat}
40+
\Sigma_{XX} = \Sigma(X) = \mathrm{E}\left[ (X-\mathrm{E}[X]) (X-\mathrm{E}[X])^\mathrm{T} \right] \; .
41+
$$
42+
43+
Using this and the [linearity of the expected value](/P/mean-lin), we can derive \eqref{eq:covmat-inv} as follows:
44+
45+
$$ \label{eq:covmat-inv-qed}
46+
\begin{split}
47+
\Sigma(X+a) &\overset{\eqref{eq:covmat}}{=} \mathrm{E}\left[ ([X+a]-\mathrm{E}[X+a]) ([X+a]-\mathrm{E}[X+a])^\mathrm{T} \right] \\
48+
&= \mathrm{E}\left[ (X + a - \mathrm{E}[X] - a) (X + a - \mathrm{E}[X] - a)^\mathrm{T} \right] \\
49+
&= \mathrm{E}\left[ (X - \mathrm{E}[X]) (X - \mathrm{E}[X])^\mathrm{T} \right] \\
50+
&\overset{\eqref{eq:covmat}}{=} \Sigma(X) \; .
51+
\end{split}
52+
$$

P/covmat-scal.md

Lines changed: 53 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,53 @@
1+
---
2+
layout: proof
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2022-09-22 11:45:00
9+
10+
title: "Scaling of the covariance matrix upon multiplication with constant matrix"
11+
chapter: "General Theorems"
12+
section: "Probability theory"
13+
topic: "Covariance"
14+
theorem: "Scaling upon multiplication with matrix"
15+
16+
sources:
17+
- authors: "Wikipedia"
18+
year: 2022
19+
title: "Covariance matrix"
20+
in: "Wikipedia, the free encyclopedia"
21+
pages: "retrieved on 2022-09-22"
22+
url: "https://en.wikipedia.org/wiki/Covariance_matrix#Basic_properties"
23+
24+
proof_id: "P348"
25+
shortcut: "covmat-scal"
26+
username: "JoramSoch"
27+
---
28+
29+
30+
**Theorem:** The [covariance matrix](/D/covmat) $\Sigma_{XX}$ of a [random vector](/D/rvec) $X$ scales upon multiplication with a [constant matrix](/D/const) $A$:
31+
32+
$$ \label{eq:covmat-scal}
33+
\Sigma(AX) = A \, \Sigma(X) A^\mathrm{T} \; .
34+
$$
35+
36+
37+
**Proof:** The [covariance matrix](/D/covmat) of $X$ [can be expressed](/P/covmat-mean) in terms of [expected values](/D/mean) as follows:
38+
39+
$$ \label{eq:covmat}
40+
\Sigma_{XX} = \Sigma(X) = \mathrm{E}\left[ (X-\mathrm{E}[X]) (X-\mathrm{E}[X])^\mathrm{T} \right] \; .
41+
$$
42+
43+
Using this and the [linearity of the expected value](/P/mean-lin), we can derive \eqref{eq:covmat-scal} as follows:
44+
45+
$$ \label{eq:covmat-scal-qed}
46+
\begin{split}
47+
\Sigma(AX) &\overset{\eqref{eq:covmat}}{=} \mathrm{E}\left[ ([AX]-\mathrm{E}[AX]) ([AX]-\mathrm{E}[AX])^\mathrm{T} \right] \\
48+
&= \mathrm{E}\left[ (A[X-\mathrm{E}[X]]) (A[X-\mathrm{E}[X]])^\mathrm{T} \right] \\
49+
&= \mathrm{E}\left[ A (X-\mathrm{E}[X]) (X-\mathrm{E}[X])^\mathrm{T} A^\mathrm{T} \right] \\
50+
&= A \, \mathrm{E}\left[ (X-\mathrm{E}[X]) (X-\mathrm{E}[X])^\mathrm{T} \right] A^\mathrm{T} \\
51+
&\overset{\eqref{eq:covmat}}{=} A \, \Sigma(X) A^\mathrm{T} \; .
52+
\end{split}
53+
$$

P/matn-dent.md

Lines changed: 70 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,70 @@
1+
---
2+
layout: proof
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2022-09-22 08:39:00
9+
10+
title: "Differential entropy for the matrix-normal distribution"
11+
chapter: "Probability Distributions"
12+
section: "Matrix-variate continuous distributions"
13+
topic: "Matrix-normal distribution"
14+
theorem: "Differential entropy"
15+
16+
sources:
17+
18+
proof_id: "P344"
19+
shortcut: "matn-dent"
20+
username: "JoramSoch"
21+
---
22+
23+
24+
**Theorem:** Let $X$ be an $n \times p$ [random matrix](/D/rmat) following a [matrix-normal distribution](/D/matn)
25+
26+
$$ \label{eq:matn}
27+
X \sim \mathcal{MN}(M, U, V) \; .
28+
$$
29+
30+
Then, the [differential entropy](/D/dent) of $X$ in nats is
31+
32+
$$ \label{eq:matn-dent}
33+
\mathrm{h}(X) = \frac{np}{2} \ln(2\pi) + \frac{n}{2} \ln|V| + \frac{p}{2} \ln|U| + \frac{np}{2} \; .
34+
$$
35+
36+
37+
**Proof:** The [matrix-normal distribution is equivalent to the multivariate normal distribution](/P/matn-mvn),
38+
39+
$$ \label{eq:matn-mvn}
40+
X \sim \mathcal{MN}(M, U, V) \quad \Leftrightarrow \quad \mathrm{vec}(X) \sim \mathcal{N}(\mathrm{vec}(M), V \otimes U) \; ,
41+
$$
42+
43+
and the [differential entropy for the multivariate normal distribution](/P/mvn-dent) in nats is
44+
45+
$$ \label{eq:mvn-dent}
46+
X \sim \mathcal{N}(\mu, \Sigma) \quad \Rightarrow \quad \mathrm{h}(X) = \frac{n}{2} \ln(2\pi) + \frac{1}{2} \ln|\Sigma| + \frac{1}{2} n
47+
$$
48+
49+
where $X$ is an $n \times 1$ [random vector](/D/rvec).
50+
51+
Thus, we can plug the distribution parameters from \eqref{eq:matn} into the differential entropy in \eqref{eq:mvn-dent} using the relationship given by \eqref{eq:matn-mvn}
52+
53+
$$ \label{eq:matn-dent-s1}
54+
\mathrm{h}(X) = \frac{np}{2} \ln(2\pi) + \frac{1}{2} \ln|V \otimes U| + \frac{1}{2} np \; .
55+
$$
56+
57+
Using the Kronecker product property
58+
59+
$$ \label{eq:kron-det}
60+
|A \otimes B| = |A|^m \, |B|^n \quad \text{where} \quad A \in \mathbb{R}^{n \times n} \quad \text{and} \quad B \in \mathbb{R}^{m \times m} \; ,
61+
$$
62+
63+
the differential entropy from \eqref{eq:matn-dent-s1} becomes:
64+
65+
$$ \label{eq:matn-dent-s2}
66+
\begin{split}
67+
\mathrm{h}(X) &= \frac{np}{2} \ln(2\pi) + \frac{1}{2} \ln\left(|V|^n |U|^p\right) + \frac{1}{2} np \\
68+
&= \frac{np}{2} \ln(2\pi) + \frac{n}{2} \ln|V| + \frac{p}{2} \ln|U| + \frac{np}{2} \; .
69+
\end{split}
70+
$$

P/ng-cov.md

Lines changed: 105 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,105 @@
1+
---
2+
layout: proof
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2022-09-22 09:17:00
9+
10+
title: "Covariance and variance of the normal-gamma distribution"
11+
chapter: "Probability Distributions"
12+
section: "Multivariate continuous distributions"
13+
topic: "Normal-gamma distribution"
14+
theorem: "Covariance"
15+
16+
sources:
17+
18+
proof_id: "P345"
19+
shortcut: "ng-cov"
20+
username: "JoramSoch"
21+
---
22+
23+
24+
**Theorem:** Let $x \in \mathbb{R}^n$ and $y > 0$ follow a [normal-gamma distribution](/D/ng):
25+
26+
$$ \label{eq:ng}
27+
x,y \sim \mathrm{NG}(\mu, \Lambda, a, b) \; .
28+
$$
29+
30+
Then,
31+
32+
1) the [covariance](/D/cov) of $x$, [conditional](/D/dist-cond) on $y$ is
33+
34+
$$ \label{eq:ng-cov-cond}
35+
\mathrm{Cov}[x|y] = \frac{1}{y} \Lambda^{-1} \; ;
36+
$$
37+
38+
2) the [covariance](/D/cov) of $x$, [unconditional](/D/dist-marg) on $y$ is
39+
40+
$$ \label{eq:ng-cov-x}
41+
\mathrm{Cov}[x] = \frac{b}{a-1} \Lambda^{-1} \; ;
42+
$$
43+
44+
3) the [variance](/D/var) of $y$ is
45+
46+
$$ \label{eq:ng-var-y}
47+
\mathrm{Var}[y] = \frac{a}{b^2} \; .
48+
$$
49+
50+
51+
**Proof:**
52+
53+
1) According to the [definition of the normal-gamma distribution](/D/ng), the distribution of $x$ given $y$ is a [multivariate normal distribution](/D/mvn):
54+
55+
$$ \label{eq:ng-mvn}
56+
x \vert y \sim \mathcal{N}(\mu, (y \Lambda)^{-1}) \; .
57+
$$
58+
59+
The [covariance of the multivariate normal distribution](/P/mvn-cov) is
60+
61+
$$ \label{eq:mvn-cov}
62+
x \sim \mathcal{N}(\mu, \Sigma) \quad \Rightarrow \quad \mathrm{Cov}(x) = \Sigma \; ,
63+
$$
64+
65+
such that we have:
66+
67+
$$ \label{eq:ng-cov-cond-qed}
68+
\mathrm{Cov}[x|y] = (y \Lambda)^{-1} = \frac{1}{y} \Lambda^{-1} \; .
69+
$$
70+
71+
2) The [marginal distribution of the normal-gamma distribution](/P/ng-marg) with respect to $x$ is a [multivariate t-distribution](/D/mvt):
72+
73+
$$ \label{eq:ng-marg-x}
74+
x \sim t\left( \mu, \left(\frac{a}{b} \Lambda \right)^{-1}, 2a \right) \; .
75+
$$
76+
77+
The [covariance of the multivariate t-distribution](/P/mvt-cov) is
78+
79+
$$ \label{eq:mvt-cov}
80+
x \sim t(\mu, \Sigma, \nu) \quad \Rightarrow \quad \mathrm{Cov}(x) = \frac{\nu}{\nu-2} \Sigma \; ,
81+
$$
82+
83+
such that we have:
84+
85+
$$ \label{eq:ng-cov-x-qed}
86+
\mathrm{Cov}[x] = \frac{2a}{2a-2} \left(\frac{a}{b} \Lambda \right)^{-1} = \frac{a}{a-1} \, \frac{b}{a} \, \Lambda^{-1} = \frac{b}{a-1} \Lambda^{-1} \; .
87+
$$
88+
89+
3) The [marginal distribution of the normal-gamma distribution](/P/ng-marg) with respect to $y$ is a [univariate gamma distribution](/D/gam):
90+
91+
$$ \label{eq:ng-marg-y}
92+
y \sim \mathrm{Gam}(a, b) \; .
93+
$$
94+
95+
The [variance of the gamma distribution](/P/gam-var) is
96+
97+
$$ \label{eq:gam-var}
98+
x \sim \mathrm{Gam}(a, b) \quad \Rightarrow \quad \mathrm{Var}(x) = \frac{a}{b^2} \; ,
99+
$$
100+
101+
such that we have:
102+
103+
$$ \label{eq:ng-var-y-qed}
104+
\mathrm{Var}[y] = \frac{a}{b^2} \; .
105+
$$

0 commit comments

Comments
 (0)