Skip to content

Commit dc330c3

Browse files
authored
Merge pull request #102 from JoramSoch/master
added 4 definitions and 3 proofs
2 parents 170c4d2 + de5683e commit dc330c3

8 files changed

Lines changed: 466 additions & 11 deletions

File tree

D/cvlme.md

Lines changed: 63 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,63 @@
1+
---
2+
layout: definition
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2020-11-19 04:55:00
9+
10+
title: "Cross-validated log model evidence"
11+
chapter: "Model Selection"
12+
section: "Bayesian model selection"
13+
topic: "Log model evidence"
14+
definition: "Cross-validated log model evidence"
15+
16+
sources:
17+
- authors: "Soch J, Allefeld C, Haynes JD"
18+
year: 2016
19+
title: "How to avoid mismodelling in GLM-based fMRI data analysis: cross-validated Bayesian model selection"
20+
in: "NeuroImage"
21+
pages: "vol. 141, pp. 469-489, eqs. 13-15"
22+
url: "https://www.sciencedirect.com/science/article/pii/S1053811916303615"
23+
doi: "10.1016/j.neuroimage.2016.07.047"
24+
- authors: "Soch J, Meyer AP, Allefeld C, Haynes JD"
25+
year: 2017
26+
title: "How to improve parameter estimates in GLM-based fMRI data analysis: cross-validated Bayesian model averaging"
27+
in: "NeuroImage"
28+
pages: "vol. 158, pp. 186-195, eq. 6"
29+
url: "https://www.sciencedirect.com/science/article/pii/S105381191730527X"
30+
doi: "10.1016/j.neuroimage.2017.06.056"
31+
- authors: "Soch J, Allefeld C"
32+
year: 2018
33+
title: "MACS – a new SPM toolbox for model assessment, comparison and selection"
34+
in: "Journal of Neuroscience Methods"
35+
pages: "vol. 306, pp. 19-31, eqs. 14-15"
36+
url: "https://www.sciencedirect.com/science/article/pii/S0165027018301468"
37+
doi: "10.1016/j.jneumeth.2018.05.017"
38+
- authors: "Soch J"
39+
year: 2018
40+
title: "cvBMS and cvBMA: filling in the gaps"
41+
in: "arXiv stat.ME"
42+
pages: "arXiv:1807.01585"
43+
url: "https://arxiv.org/abs/1807.01585"
44+
45+
def_id: "D111"
46+
shortcut: "cvlme"
47+
username: "JoramSoch"
48+
---
49+
50+
51+
**Definition:** Let there be a [data set](/D/data) $y$ with mutually exclusive and collectively exhaustive subsets $y_1, \ldots, y_S$. Assume a [generative model](/D/gm) $m$ with model parameters $\theta$ implying a [likelihood function](/D/lf) $p(y \vert \theta, m)$ and a [non-informative](/D/prior-inf) [prior density](/D/prior) $p(\theta \vert m)$.
52+
53+
Then, the cross-validated log model evidence of $m$ is given by
54+
55+
$$ \label{eq:cvLME}
56+
\mathrm{cvLME}(m) = \sum_{i=1}^{S} \log \int p( y_i \vert \theta, m ) \, p( \theta \vert y_{\neg i}, m ) \, \mathrm{d}\theta
57+
$$
58+
59+
where $y_{\neg i} = \bigcup_{j \neq i} y_j$ is the union of all data subsets except $y_i$ and $p( \theta \vert y_{\neg i}, m )$ is the [posterior distribution](/D/post) obtained from $y_{\neg i}$ when using the [prior distribution](/D/prior) $p(\theta \vert m)$:
60+
61+
$$ \label{eq:post}
62+
p( \theta \vert y_{\neg i}, m ) = \frac{p( y_{\neg i} \vert \theta, m ) \, p(\theta \vert m)}{p( y_{\neg i} \vert m )} \; .
63+
$$

D/ind-cond.md

Lines changed: 54 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,54 @@
1+
---
2+
layout: definition
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2020-11-19 05:40:00
9+
10+
title: "Conditional independence"
11+
chapter: "General Theorems"
12+
section: "Probability theory"
13+
topic: "Probability"
14+
definition: "Conditional independence"
15+
16+
sources:
17+
- authors: "Wikipedia"
18+
year: 2020
19+
title: "Conditional independence"
20+
in: "Wikipedia, the free encyclopedia"
21+
pages: "retrieved on 2020-11-19"
22+
url: "https://en.wikipedia.org/wiki/Conditional_independence#Conditional_independence_of_random_variables"
23+
24+
def_id: "D112"
25+
shortcut: "ind-cond"
26+
username: "JoramSoch"
27+
---
28+
29+
30+
**Definition:** Generally speaking, [random variables](/D/rvar) are conditionally independent given another random variable, if they are [statistically independent](/D/ind) in their [conditional probability distributions](/D/dist-cond) given this random variable.
31+
32+
<br>
33+
1) A set of [discrete random variables](/D/rvar-disc) $X_1, \ldots, X_n$ with possible values $\mathcal{X}_1, \ldots, \mathcal{X}_n$ is called conditionally independent given the random variable $Y$ with possible values $\mathcal{Y}$, if
34+
35+
$$ \label{eq:disc-ind}
36+
p(X_1 = x_1, \ldots, X_n = x_n|Y = y) = \prod_{i=1}^{n} p(X_i = x_i|Y = y) \quad \text{for all} \; x_i \in \mathcal{X}_i \quad \text{and all} \; y \in \mathcal{Y}
37+
$$
38+
39+
where $p(x_1, \ldots, x_n \vert y)$ are the [joint (conditional) probabilities](/D/prob-joint) of $X_1, \ldots, X_n$ given $Y$ and $p(x_i)$ are the [marginal (conditional) probabilities](/D/prob-marg) of $X_i$ given $Y$.
40+
41+
<br>
42+
2) A set of [random variables](/D/rvar) $X_1, \ldots, X_n$ with possible values $\mathcal{X}_1, \ldots, \mathcal{X}_n$ is called conditionally independent given the random variable $Y$ with possible values $\mathcal{Y}$, if
43+
44+
$$ \label{eq:cond-ind-F}
45+
F_{X_1,\ldots,X_n|Y=y}(x_1,\ldots,x_n) = \prod_{i=1}^{n} F_{X_i|Y=y}(x_i) \quad \text{for all} \; x_i \in \mathcal{X}_i \quad \text{and all} \; y \in \mathcal{Y}
46+
$$
47+
48+
or equivalently, if the [probability densities](/D/pdf) exist, if
49+
50+
$$ \label{eq:cont-ind-f}
51+
f_{X_1,\ldots,X_n|Y=y}(x_1,\ldots,x_n) = \prod_{i=1}^{n} f_{X_i|Y=y}(x_i) \quad \text{for all} \; x_i \in \mathcal{X}_i \quad \text{and all} \; y \in \mathcal{Y}
52+
$$
53+
54+
where $F$ are the [joint (conditional)](/D/dist-joint) or [marginal (conditional)](/D/dist-marg) [cumulative distribution functions](/D/cdf) and $f$ are the respective [probability density functions](/D/pdf).

D/reve.md

Lines changed: 36 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,36 @@
1+
---
2+
layout: definition
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2020-11-19 04:33:00
9+
10+
title: "Random event"
11+
chapter: "General Theorems"
12+
section: "Probability theory"
13+
topic: "Random variables"
14+
definition: "Random event"
15+
16+
sources:
17+
- authors: "Wikipedia"
18+
year: 2020
19+
title: "Event (probability theory)"
20+
in: "Wikipedia, the free encyclopedia"
21+
pages: "retrieved on 2020-11-19"
22+
url: "https://en.wikipedia.org/wiki/Event_(probability_theory)"
23+
24+
def_id: "D110"
25+
shortcut: "reve"
26+
username: "JoramSoch"
27+
---
28+
29+
30+
**Definition:** A random event $E$ is the outcome of a [random experiment](/D/rexp) which can be described by a statement that is either true or false.
31+
32+
* If the statement is true, the event is said to take place, denoted as $E$.
33+
34+
* If the statement is false, the complement of $E$ occurs, denoted as $\overline{E}$.
35+
36+
In other words, a random event is a [random variable](/D/rvar) with two possible values (true and false, or 1 and 0). A [random experiment](/D/rexp) with two possible outcomes is called a [Bernoulli trial](/D/bern).

D/rexp.md

Lines changed: 38 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,38 @@
1+
---
2+
layout: definition
3+
mathjax: true
4+
5+
author: "Joram Soch"
6+
affiliation: "BCCN Berlin"
7+
e_mail: "joram.soch@bccn-berlin.de"
8+
date: 2020-11-19 04:10:00
9+
10+
title: "Random experiment"
11+
chapter: "General Theorems"
12+
section: "Probability theory"
13+
topic: "Random variables"
14+
definition: "Random experiment"
15+
16+
sources:
17+
- authors: "Wikipedia"
18+
year: 2020
19+
title: "Experiment (probability theory)"
20+
in: "Wikipedia, the free encyclopedia"
21+
pages: "retrieved on 2020-11-19"
22+
url: "https://en.wikipedia.org/wiki/Experiment_(probability_theory)"
23+
24+
def_id: "D109"
25+
shortcut: "rexp"
26+
username: "JoramSoch"
27+
---
28+
29+
30+
**Definition:** A random experiment is any repeatable procedure that [results in one](/D/rvar) out of a well-defined set of possible outcomes.
31+
32+
* The set of possible outcomes is called sample space.
33+
34+
* A set of zero or more outcomes is called a [random event](/D/reve).
35+
36+
* A function that maps from events to probabilities is called a [probability function](/D/pmf).
37+
38+
Together, sample space, event space and probability function characterize a random experiment.

I/Table_of_Contents.md

Lines changed: 18 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -17,12 +17,14 @@ title: "Table of Contents"
1717
1. Probability theory
1818

1919
1.1. Random variables <br>
20-
&emsp;&ensp; 1.1.1. *[Random variable](/D/rvar)* <br>
21-
&emsp;&ensp; 1.1.2. *[Random vector](/D/rvec)* <br>
22-
&emsp;&ensp; 1.1.3. *[Random matrix](/D/rmat)* <br>
23-
&emsp;&ensp; 1.1.4. *[Constant](/D/const)* <br>
24-
&emsp;&ensp; 1.1.5. *[Discrete vs. continuous](/D/rvar-disc)* <br>
25-
&emsp;&ensp; 1.1.6. *[Univariate vs. multivariate](/D/rvar-uni)* <br>
20+
&emsp;&ensp; 1.1.1. *[Random experiment](/D/rexp)* <br>
21+
&emsp;&ensp; 1.1.2. *[Random event](/D/reve)* <br>
22+
&emsp;&ensp; 1.1.3. *[Random variable](/D/rvar)* <br>
23+
&emsp;&ensp; 1.1.4. *[Random vector](/D/rvec)* <br>
24+
&emsp;&ensp; 1.1.5. *[Random matrix](/D/rmat)* <br>
25+
&emsp;&ensp; 1.1.6. *[Constant](/D/const)* <br>
26+
&emsp;&ensp; 1.1.7. *[Discrete vs. continuous](/D/rvar-disc)* <br>
27+
&emsp;&ensp; 1.1.8. *[Univariate vs. multivariate](/D/rvar-uni)* <br>
2628

2729
1.2. Probability <br>
2830
&emsp;&ensp; 1.2.1. *[Probability](/D/prob)* <br>
@@ -31,6 +33,7 @@ title: "Table of Contents"
3133
&emsp;&ensp; 1.2.4. *[Conditional probability](/D/prob-cond)* <br>
3234
&emsp;&ensp; 1.2.5. *[Exceedance probability](/D/prob-exc)* <br>
3335
&emsp;&ensp; 1.2.6. *[Statistical independence](/D/ind)* <br>
36+
&emsp;&ensp; 1.2.7. *[Conditional independence](/D/ind-cond)* <br>
3437

3538
1.3. Probability distributions <br>
3639
&emsp;&ensp; 1.3.1. *[Probability distribution](/D/dist)* <br>
@@ -265,6 +268,7 @@ title: "Table of Contents"
265268
&emsp;&ensp; 3.2.14. **[Variance](/P/norm-var)** <br>
266269
&emsp;&ensp; 3.2.15. **[Full width at half maximum](/P/norm-fwhm)** <br>
267270
&emsp;&ensp; 3.2.16. **[Differential entropy](/P/norm-dent)** <br>
271+
&emsp;&ensp; 3.2.17. **[Kullback-Leibler divergence](/P/norm-kl)** <br>
268272

269273
3.3. Gamma distribution <br>
270274
&emsp;&ensp; 3.3.1. *[Definition](/D/gam)* <br>
@@ -273,11 +277,12 @@ title: "Table of Contents"
273277
&emsp;&ensp; 3.3.4. **[Relation to standard gamma distribution](/P/gam-sgam2)** (2) <br>
274278
&emsp;&ensp; 3.3.5. **[Probability density function](/P/gam-pdf)** <br>
275279
&emsp;&ensp; 3.3.6. **[Cumulative distribution function](/P/gam-cdf)** <br>
276-
&emsp;&ensp; 3.3.7. **[Mean](/P/gam-mean)** <br>
277-
&emsp;&ensp; 3.3.8. **[Variance](/P/gam-var)** <br>
278-
&emsp;&ensp; 3.3.9. **[Logarithmic expectation](/P/gam-logmean)** <br>
279-
&emsp;&ensp; 3.3.10. **[Expectation of x ln x](/P/gam-xlogx)** <br>
280-
&emsp;&ensp; 3.3.11. **[Kullback-Leibler divergence](/P/gam-kl)** <br>
280+
&emsp;&ensp; 3.3.7. **[Quantile function](/P/gam-qf)** <br>
281+
&emsp;&ensp; 3.3.8. **[Mean](/P/gam-mean)** <br>
282+
&emsp;&ensp; 3.3.9. **[Variance](/P/gam-var)** <br>
283+
&emsp;&ensp; 3.3.10. **[Logarithmic expectation](/P/gam-logmean)** <br>
284+
&emsp;&ensp; 3.3.11. **[Expectation of x ln x](/P/gam-xlogx)** <br>
285+
&emsp;&ensp; 3.3.12. **[Kullback-Leibler divergence](/P/gam-kl)** <br>
281286

282287
3.4. Exponential distribution <br>
283288
&emsp;&ensp; 3.4.1. *[Definition](/D/exp)* <br>
@@ -297,6 +302,7 @@ title: "Table of Contents"
297302
3.6. Beta distribution <br>
298303
&emsp;&ensp; 3.6.1. *[Definition](/D/beta)* <br>
299304
&emsp;&ensp; 3.6.2. **[Probability density function](/P/beta-pdf)** <br>
305+
&emsp;&ensp; 3.6.3. **[Cumulative distribution function](/P/beta-cdf)** <br>
300306

301307
3.7. Wald distribution <br>
302308
&emsp;&ensp; 3.7.1. *[Definition](/D/wald)* <br>
@@ -466,6 +472,7 @@ title: "Table of Contents"
466472
&emsp;&ensp; 3.1.1. *[Definition](/D/lme)* <br>
467473
&emsp;&ensp; 3.1.2. **[Derivation](/P/lme-der)** <br>
468474
&emsp;&ensp; 3.1.3. **[Partition into accuracy and complexity](/P/lme-anc)** <br>
475+
&emsp;&ensp; 3.1.4. *[Cross-validated log model evidence](/D/cvlme)* <br>
469476

470477
3.2. Log family evidence <br>
471478
&emsp;&ensp; 3.2.1. *[Definition](/D/lfe)* <br>

P/beta-cdf.md

Lines changed: 71 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,71 @@
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: 2020-11-19 08:01:00
9+
10+
title: "Cumulative distribution function of the beta distribution"
11+
chapter: "Probability Distributions"
12+
section: "Univariate continuous distributions"
13+
topic: "Beta distribution"
14+
theorem: "Cumulative distribution function"
15+
16+
sources:
17+
- authors: "Wikipedia"
18+
year: 2020
19+
title: "Beta function"
20+
in: "Wikipedia, the free encyclopedia"
21+
pages: "retrieved on 2020-11-19"
22+
url: "https://en.wikipedia.org/wiki/Beta_function#Incomplete_beta_function"
23+
24+
proof_id: "P195"
25+
shortcut: "beta-cdf"
26+
username: "JoramSoch"
27+
---
28+
29+
30+
**Theorem:** Let $X$ be a positive [random variable](/D/rvar) following a [beta distribution](/D/gam):
31+
32+
$$ \label{eq:beta}
33+
X \sim \mathrm{Bet}(\alpha, \beta) \; .
34+
$$
35+
36+
Then, the [cumulative distribution function](/D/cdf) of $X$ is
37+
38+
$$ \label{eq:beta-cdf}
39+
F_X(x) = \frac{B(x; \alpha, \beta)}{B(\alpha, \beta)}
40+
$$
41+
42+
where $B(a,b)$ is the beta function and $B(x;a,b)$ is the incomplete gamma function.
43+
44+
45+
**Proof:** The [probability density function of the beta distribution](/P/beta-pdf) is:
46+
47+
$$ \label{eq:beta-pdf}
48+
f_X(x) = \frac{1}{\mathrm{B}(\alpha, \beta)} \, x^{\alpha-1} \, (1-x)^{\beta-1} \; .
49+
$$
50+
51+
Thus, the [cumulative distribution function](/D/cdf) is:
52+
53+
$$ \label{eq:beta-cdf-app}
54+
\begin{split}
55+
F_X(x) &= \int_{0}^{x} \mathrm{Bet}(z; \alpha, \beta) \, \mathrm{d}z \\
56+
&= \int_{0}^{x} \frac{1}{\mathrm{B}(\alpha, \beta)} \, z^{\alpha-1} \, (1-z)^{\beta-1} \, \mathrm{d}z \\
57+
&= \frac{1}{B(x;a,b)} \int_{0}^{x} z^{\alpha-1} \, (1-z)^{\beta-1} \, \mathrm{d}z \; .
58+
\end{split}
59+
$$
60+
61+
With the definition of the incomplete beta function
62+
63+
$$ \label{eq:inc-beta-fct}
64+
B(x;a,b) = \int_{0}^{x} t^{a-1} \, (1-t)^{b-1} \, \mathrm{d}t \; ,
65+
$$
66+
67+
we arrive at the final result given by equation \eqref{eq:beta-cdf}:
68+
69+
$$ \label{eq:beta-cdf-qed}
70+
F_X(x) = \frac{B(x; \alpha, \beta)}{B(\alpha, \beta)} \; .
71+
$$

0 commit comments

Comments
 (0)