{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"# EXERCISES\n",
"\n",
"### Bayesian Machine Learning and Information Processing (5SSD0)"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Prelude \n",
"\n",
"- In this notebook, we provide a set of exercises that should help you prepare for the exam. There are two sets of exercises: \n",
" 1. (**5SSB0**). The first set contains questions from previous exams of class 5SSB0 (Adaptive Information Processing). 5SSB0 was last taught at TU/e in the academic year 2018/19. Starting with the academic year 2019/20, the current class (5SSD0) is offered instead. In comparison to 5SSB0, the emphasis of the current class is more on the Bayesian framework, rather than maximum likelihood. Below, we present a selection of exercises that were used in the 5SSB0 exams. While these exercises may not represent the emphasis on Bayesian methods, they do represent the \"style of questions\" and the \"level of difficulty\" that you may expect in upcoming 5SSD0 exams. Also, there is nothing in the questions below that is considered outside the scope of 5SSD0. \n",
" 2. (**Rehearsal**). A second set of exercises are categorized by lesson headers, e.g., \"Probability Theory\" or \"Generative Classification\". These exercises intend to test you on the contents of the corresponding lesson. Of course, for some exercises you may need some contents of some other (usually earlier) lessons or background materials. A perfect categorization is not feasible, but we've tried to link each question to a lesson so as to make it easier to test yourself after studying a specific lesson. \n",
" \n",
"- For some of these exercises you will be able to find solutions quickly on the internet. Try to resist this route to solving the problems. You will not be graded for these exercises and solutions will be made available in a separate notebook. **Your ability to solve these exercises without external help provides an excellent indicator of your readiness to pass the exam**. \n",
"\n",
"- This notebook is still under construction. We will be adding more exercises to help you prepare. Also, not all solutions are provided yet. We are working on this too. If you absolutely need the solution for a problem that doesn't have a solution yet, please contact one of the TAs (see http://bmlip.nl for contact info.)\n",
"\n",
"- As 2019/20 is the first time we teach this class, please be alert to errors (and let us know if you find any!) or more generally, let us know if things are unclear to you or if a question can be improved. \n",
"\n",
"- We are aware there are some rendering problems in some browsers. We are trying to fix that as well. \n",
"\n",
"- Good luck!"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Cheatsheet\n",
"\n",
"- You are not allowed to bring books or notes to the exam. Instead, feel free to make use of the following cheatsheet as we will provide this or a similar cheatsheet in an appendix of the exam papers.\n",
"\n",
"- Some Matrix Calculus, see also Bishop, appendix C. \n",
"$$\\begin{align*}\n",
"|A^{-1}|&=|A|^{-1} \\\\\n",
"\\nabla_A \\log |A| &= (A^{T})^{-1} = (A^{-1})^T \\\\\n",
"\\mathrm{Tr}[ABC]&= \\mathrm{Tr}[CAB] = \\mathrm{Tr}[BCA] \\\\\n",
"\\nabla_A \\mathrm{Tr}[AB] &=\\nabla_A \\mathrm{Tr}[BA]= B^T \\\\\n",
"\\nabla_A \\mathrm{Tr}[ABA^T] &= A(B+B^T)\\\\\n",
" \\nabla_x x^TAx &= (A+A^T)x\\\\\n",
"\\nabla_X a^TXb &= \\nabla_X \\mathrm{Tr}[ba^TX] = ab^T\n",
"\\end{align*}$$\n",
"\n",
"- Definition of the Multivariate Gaussian Distribution (MVG)\n",
"$$\n",
"\\mathcal{N}(x|\\,\\mu,\\Sigma) = |2 \\pi \\Sigma|^{-\\frac{1}{2}} \\exp\\left\\{-\\frac{1}{2}(x-\\mu)^T\n",
"\\Sigma^{-1} (x-\\mu) \\right\\}\n",
"$$\n",
"\n",
"- A **linear transformation** $z=Ax+b$ of a Gaussian variable $\\mathcal{N}(x|\\mu,\\Sigma)$ is Gaussian distributed as\n",
"$$\n",
"p(z) = \\mathcal{N} \\left(z \\,|\\, A\\mu+b, A\\Sigma A^T \\right) \n",
"$$\n",
"\n",
"- **Multiplication** of 2 Gaussian distributions\n",
"$$ \n",
" \\mathcal{N}(x|\\mu_a,\\Sigma_a) \\cdot \\mathcal{N}(x|\\mu_b,\\Sigma_b) = \\alpha \\cdot \\mathcal{N}(x|\\mu_c,\\Sigma_c)\n",
"$$\n",
"with\n",
"$$\\begin{align*}\n",
"\\Sigma_c^{-1} &= \\Sigma_a^{-1} + \\Sigma_b^{-1} \\\\\n",
"\\Sigma_c^{-1}\\mu_c &= \\Sigma_a^{-1}\\mu_a + \\Sigma_b^{-1}\\mu_b \\\\\n",
"\\alpha &= \\mathcal{N}(\\mu_a | \\mu_b, \\Sigma_a + \\Sigma_b)\n",
"\\end{align*}$$\n",
"\n",
"- **Conditioning** and **marginalization** of Gaussians. Let $z = \\begin{bmatrix} x \\\\ y \\end{bmatrix}$ be jointly normal distributed as\n",
"$$\\begin{align*}\n",
"p(z) &= \\mathcal{N}(z | \\mu, \\Sigma) \n",
" =\\mathcal{N} \\left( \\begin{bmatrix} x \\\\ y \\end{bmatrix} \\,\\left|\\, \\begin{bmatrix} \\mu_x \\\\ \\mu_y \\end{bmatrix}, \n",
" \\begin{bmatrix} \\Sigma_x & \\Sigma_{xy} \\\\ \\Sigma_{yx} & \\Sigma_y \\end{bmatrix} \\right. \\right)\\,,\n",
"\\end{align*}$$\n",
"then $p(z) = p(y|x)\\cdot p(x)$, with \n",
"$$\\begin{align*}\n",
"p(y|x) &= \\mathcal{N}\\left(y\\,|\\,\\mu_y + \\Sigma_{yx}\\Sigma_x^{-1}(x-\\mu_x),\\, \\Sigma_y - \\Sigma_{yx}\\Sigma_x^{-1}\\Sigma_{xy} \\right) \\\\\n",
"p(x) &= \\mathcal{N}\\left( x\\,|\\,\\mu_x, \\Sigma_x \\right)\n",
"\\end{align*}$$\n",
"\n",
"- For a binary variable $x \\in \\{0,1\\}$, the **Bernoulli** distribution is given by\n",
"$$ \n",
"p(x|\\mu) = \\mu^{x}(1-\\mu)^{1-x}\n",
"$$\n",
"\n",
"- The conjugate prior for $\\mu$ is the **Beta** distribution, given by\n",
"$$\n",
"p(\\mu) = \\mathcal{B}(\\mu|\\alpha,\\beta) = \\frac{\\Gamma(\\alpha+\\beta)}{\\Gamma(\\alpha)\\Gamma(\\beta)} \\mu^{\\alpha-1}(1-\\mu)^{\\beta-1}\n",
"$$\n",
"where $\\alpha$ and $\\beta$ are \"hyperparameters\" that you can set to reflect your prior beliefs about $\\mu$. "
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"# Selected exercises from previous exams \"5SSB0\"\n",
"\n",
"\n",
"- **[1]** Answer shortly (max. 3 sentences): What is the difference between supervised and unsupervised learning? \n",
" \n",
"- **[2]** Which of the following statements is true (or justified)? \n",
" (a) If $X$ and $Y$ are independent Gaussian distributed variables, then $Z = 3X+Y$ is also a Gaussian distributed variable. \n",
" (b) The sum of two Gaussian functions is always also a Gaussian function. \n",
" (c) Discriminative classification is more similar to regression than to density estimation. \n",
" (d) Density estimation is more similar to generative classification than to discriminative classification. \n",
" (e) Clustering is more similar to supervised classification than to unsupervised classification. \n",
"\n",
"- **[3]** Consider a binary classification problem with two classes $\\{y_1,y_2\\}$ and input vector $x$. Outputs $y_k$ are recorded by a one-hot encoding scheme. We are given a data set to train the parameters $\\theta$ for a likelihood model of the form\n",
"$$\n",
"p(y_k=1|x,\\theta) = \\frac{1}{1 + e^{-\\theta_k^T x}}\n",
"$$\n",
"There a two fundamentally different ways to train $\\theta$, namely through a generative model or by discriminative training. \n",
" (a) Explain shortly how we train $\\theta$ through a generative model. No need to work out all equations for Gaussian models, but explain the strategy in probabilistic modeling terms. \n",
" (b) Explain shortly how we train $\\theta$ through a discriminative\n",
"approach. \n",
"\n",
"- **[4]** What is the difference between supervised and unsupervised learning? Express the goals of these two learning methods in terms of a probability distribution. (I'm looking here for a statement such as: \" Given $\\ldots$, the goals of supervised/unsupervised learning is to estimate $p(\\cdot|\\cdot)$\".) \n",
"\n",
"- **[5]** In a particular model with hidden variables, the log-likelihood can be worked out to the following expression:\n",
"$$\n",
" L(\\theta) = \\sum_n \\log \\left(\\sum_k \\pi_k\\,\\mathcal{N}(x_n|\\mu_k,\\Sigma_k)\\right)\n",
"$$\n",
"Do you prefer a gradient descent or EM algorithm to estimate maximum likelihood values for the parameters? Explain your answer. (No need to work out the equations.)\n",
"\n",
"- **[6]** Consider a data set $D = \\{x_1,x_2,\\ldots,x_N\\}$ where we assume that each sample $x_n$ is IID distributed by a multivariate Gaussian (MVG) $\\mathcal{N}(x_n|\\,\\mu,\\Sigma)$. \n",
"Proof that the maximum likelihood estimate (MLE) of the mean value of this distribution is given by\n",
"$$\\begin{equation*}\n",
"\\hat \\mu = \\frac{1}{N}\\sum_n x_n \n",
"\\end{equation*}$$ \n",
"(Note the list of [matrix calculus formulas](#matrix-calculus) above).\n",
"\n",
"- **[7]** Consider a data set $D = \\{(x_1,y_1), (x_2,y_2),\\dots,(x_N,y_N)\\}$ with 1-of-$K$ notation for the discrete classes, i.e.,\n",
"\n",
"\\begin{equation*} y_{nk} = \\begin{cases} 1 & \\text{if $y_n$ in $k$th class} \\\\\n",
" 0 & \\text{otherwise} \n",
" \\end{cases}\n",
"\\end{equation*}\n",
"\n",
"together with class-conditional distribution $p(x_n| y_{nk}=1,\\theta) = \\mathcal{N}(x_n|\\mu_k,\\Sigma)$ and multinomial prior $p(y_{nk}=1) = \\pi_k$. \n",
" (a) Proof that the joint log-likelihood is given by\n",
"$$\\begin{equation*}\n",
"\\log p(D|\\theta) = \\sum_{n,k} y_{nk} \\log \\mathcal{N}(x_n|\\mu_k,\\Sigma) + \\sum_{n,k} y_{nk} \\log \\pi_k\n",
"\\end{equation*}$$ \n",
" (b) Show now that the MLE of the *class-conditional* mean is given by\n",
"$$\\begin{equation*}\n",
" \\hat \\mu_k = \\frac{\\sum_n y_{nk} x_n}{\\sum_n y_{nk}} \n",
"\\end{equation*}\n",
"$$\n",
"\n",
"- **[8]** Consider an IID data set $D=\\{(x_1,y_1),\\ldots,(x_N,y_N)\\}$. We will model this data set by a model $$y_n =\\theta^T f(x_n) + e_n\\,,$$ where $f(x_n)$ is an $M$-dimensional feature vector of input $x_n$; $y_n$ is a scalar output and $e_n \\sim \\mathcal{N}(0,\\sigma^2)$. \n",
" (a) Rewrite the model in matrix form by lumping input features in a matrix $F=[f(x_1),\\ldots,f(x_N)]^T$, outputs and noise in the vectors $y=[y_1,\\ldots,y_N]^T$ and $e=[e_1,\\ldots,e_N]^T$, respectively. \n",
"\n",
" (b) Now derive an expression for the log-likelihood $\\log p(y|\\,F,\\theta,\\sigma^2)$. \n",
"\n",
" (c) Proof that the maximum likelihood estimate for the parameters is given by \n",
"$$\\hat\\theta_{\\text{ml}} = (F^TF)^{-1}F^Ty$$ \n",
"(Note the list of [matrix calculus formulas](#matrix-calculus) above). \n",
"\n",
" (d) What is the predicted output value $y_\\bullet$, given an observation $x_\\bullet$ and the maximum likelihood parameters $\\hat \\theta_{\\text{ml}}$. Work this expression out in terms of $F$, $y$ and $f(x_\\bullet)$. \n",
"\n",
" (e) Suppose that, before the data set $D$ was observed, we had reason to assume a prior distribution $p(\\theta)=\\mathcal{N}(0,\\sigma_0^2)$. Derive the Maximum a posteriori (MAP) estimate $\\hat \\theta_{\\text{map}}$.(hint: work this out in the $\\log$ domain.) \n",
" \n",
"- **[9]** For each of the following sub-questions, provide a *short but essential* answer. \n",
" (a) The joint distribution for feature vector $x$ and membership of $k$th class is given by (we use one-hot encoding of classes)\n",
"$$p(x,y_k=1) = \\pi_k \\cdot \\mathcal{N}(x |\\mu_k,\\Sigma)$$ \n",
"Write down an expression for the posterior class probability $p(y_k=1|x)$ (No derivations are needed, just a proper expression)? \n",
"\n",
" (b) Why does maximum likelihood estimation become a better approximation to Bayesian learning as you collect more data? \n",
"\n",
" (c) Given is a model \n",
"$$\\begin{align*}\n",
"p(x|z) &= \\mathcal{N}(x | W z,\\Psi) \\\\\n",
"p(z) &= \\mathcal{N}(z|0,I)\n",
"\\end{align*}$$ \n",
"Work out an expression for the marginal distribution $p(x)$. \n",
"\n",
"- **[10]** Explain shortly how Bayes rule relates to machine learning in the context of an observed data set $D$ and a model $M$ with parameters $\\theta$. Your answer must contain the expression for Bayes rule. \n",
"\n",
"- **[11]** Consider the following state-space model:\n",
"$$\\begin{align*}\n",
"z_k &= A z_{k-1} + w_k \\\\\n",
"x_k &= C z_k + v_k \n",
"\\end{align*}$$\n",
"where $k=1,2,\\ldots,n$ is the time step counter; $z_k$ is an *unobserved* state sequence; $x_k$ is an *observed* sequence; $w_k \\sim \\mathcal{N}(0,\\Sigma_w)$ and $v_k \\sim \\mathcal{N}(0,\\Sigma_v)$ are (unobserved) state and observation noise sequences respectively; $z_0 \\sim \\mathcal{N}(0,\\Sigma_0)$ is the initial state and $A$, $C$, $\\Sigma_v$,$\\Sigma_w$ and $\\Sigma_0$ are known parameters. The Forney-style factor graph (FFG) for one time step is depicted here: \n",
" \n",
" (a) Rewrite the state-space equations as a set of conditional probability distributions. \n",
" $$\\begin{align*}\n",
" p(z_k|z_{k-1},A,\\Sigma_w) &= \\ldots \\\\\n",
" p(x_k|z_k,C,\\Sigma_v) &= \\ldots \\\\\n",
" p(z_0|\\Sigma_0) &= \\ldots\n",
"\\end{align*}$$ \n",
" (b) Define $z^n \\triangleq (z_0,z_1,\\ldots,z_n)$, $x^n \\triangleq (x_1,\\ldots,x_n)$ and $\\theta=\\{A,C,\\Sigma_w,\\Sigma_v\\}$. Now write out the generative model $p(x^n,z^n|\\theta)$ as a product of factors. \n",
" (c) We are interested in estimating $z_k$ from a given estimate for $z_{k-1}$ and the current observation $x_k$, i.e., we are interested in computing $p(z_k|z_{k-1},x_k,\\theta)$. Can $p(z_k|z_{k-1},x_k,\\theta)$ be expressed as a Gaussian distribution? Explain why or why not in one sentence. \n",
" (d) Copy the graph onto your exam paper and draw the message passing schedule for computing $p(z_k|z_{k-1},x_k,\\theta)$ by drawing arrows in the factor graph. Indicate the order of the messages by assigning numbers to the arrows. \n",
" (e) Now assume that our belief about parameter $\\Sigma_v$ is instead given by a distribution $p(\\Sigma_v)$ (rather than a known value). Adapt the factor graph drawing of the previous answer to reflects our belief about $\\Sigma_v$. \n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"# Rehearsal Exercises\n",
"\n",
"- Below you will find more questions that test your knowledge of each lesson. My perception of the difficulty level of each question may differ from yours, but I ll indicate my ratings with three levels, (#) for easy, (##) for intermediate level and (###) for (very) challenging. While all questions are in principle within the scope of the lessons, there should be enough questions of levels (#) and (##) in the exams to pass the class. "
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Machine Learning Overview\n",
"\n",
"- **[1]** (##) Pick three applications from the [\"Some Machine Learning Applications\"](https://nbviewer.jupyter.org/github/bertdv/BMLIP/blob/master/lessons/notebooks/Machine-Learning-Overview.ipynb#some-ml-apps)-slide and (shortly) describe for each application how (a combination of) clustering, dimensionality reduction, regression classification or reinforcement learning could accomplish the task."
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Probability Theory Review\n",
"\n",
"\n",
"- **[1]** (a) (#) Proof that the \"elementary\" sum rule $p(A) + p(\\bar{A}) = 1$ follows from the (general) sum rule $$p(A+B) = p(A) + p(B) - p(A,B)\\,.$$ \n",
" (b) (###) Conversely, derive the general sum rule\n",
" $p(A + B) = p(A) + p(B) - p(A,B)$\n",
"from the elementary sum rule $p(A) + p(\\bar A) = 1$ and the product rule. Here, you may make use of the (Boolean logic) fact that $A + B = \\overline {\\bar A \\bar B }$. \n",
" \n",
"- **[2]** Box 1 contains 8 apples and 4 oranges. Box 2 contains 10 apples and 2 oranges. Boxes are chosen with equal probability. \n",
" (a) (#) What is the probability of choosing an apple? \n",
" (b) (##) If an apple is chosen, what is the probability that it came from box 1?\n",
"\n",
"- **[3]** (###) The inhabitants of an island tell the truth one third of the time. They lie with probability $2/3$. On an occasion, after one of them made a statement, you ask another \"was that statement true?\" and he says \"yes\". What is the probability that the statement was indeed true?\n",
"\n",
"- **[4]** (##) A bag contains one ball, known to be either white or black. A white ball is put in, the bag is shaken, and a ball is drawn out, which proves to be white. What is now the chance of drawing a white ball? (Note that the state of the bag, after the operations, is exactly identical to its state before.)\n",
"\n",
"- **[5]** A dark bag contains five red balls and seven green ones. \n",
" (a) (#) What is the probability of drawing a red ball on the first draw? \n",
" (b) (##) Balls are not returned to the bag after each draw. If you know that on the second draw the ball was a green one, what is now the probability of drawing a red ball on the first draw? \n",
"\n",
"- **[6]** (#) Is it more correct to speak about the likelihood of a _model_ (or model parameters) than about the likelihood of an _observed data set_. And why? \n",
"\n",
"- **[7]** (##) Is a speech signal a 'probabilistic' (random) or a deterministic signal?\n",
" \n",
"- **[8]** (##) Proof that, for any distribution of $x$ and $y$ and $z=x+y$\n",
"$$\\begin{align*}\n",
" \\mathrm{E}[z] &= \\mathrm{E}[x] + \\mathrm{E}[y] \\\\\n",
" \\mathrm{var}[z] &= \\mathrm{var}[x] + \\mathrm{var}[y] + 2\\mathrm{cov}[x,y] \n",
"\\end{align*}$$\n",
"You may make use of the more general theorem that the mean and variance of any distribution $p(x)$ is processed by a linear tranformation as\n",
" $$\\begin{align*}\n",
"\\mathrm{E}[Ax +b] &= A\\mathrm{E}[x] + b \\\\\n",
"\\mathrm{var}[Ax +b] &= A\\,\\mathrm{var}[x]\\,A^T \n",
"\\end{align*}$$\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Bayesian Machine Learning\n",
"\n",
"\n",
"- **[1]** (#) (a) Explain shortly the relation between machine learning and Bayes rule. \n",
" (b) How are Maximum a Posteriori (MAP) and Maximum Likelihood (ML) estimation related to Bayes rule and machine learning?\n",
"\n",
"- **[2]** (#) What are the four stages of the Bayesian design approach?\n",
"\n",
"- **[3]** (##) The Bayes estimate is a summary of a posterior distribution by a delta distribution on its mean, i.e., \n",
"$$\n",
"\\hat \\theta_{bayes} = \\int \\theta \\, p\\left( \\theta |D \\right)\n",
"\\,\\mathrm{d}{\\theta}\n",
"$$\n",
"Proof that the Bayes estimate minimizes the expected mean-square error, i.e., proof that\n",
"$$\n",
"\\hat \\theta_{bayes} = \\arg\\min_{\\hat \\theta} \\int_\\theta (\\hat \\theta -\\theta)^2 p \\left( \\theta |D \\right) \\,\\mathrm{d}{\\theta}\n",
"$$\n",
"\n",
"- **[4]** (##) We make $N$ IID observations $D=\\{x_1 \\dots x_N\\}$ and assume the following model\n",
"$$\n",
"x_k = A + \\epsilon_k \n",
"$$\n",
" where $\\epsilon_k = \\mathcal{N}(\\epsilon_k | 0,\\sigma^2)$ with known $\\sigma^2=1$. We are interested in deriving an estimator for $A$. \n",
" (a) Make a reasonable assumption for a prior on $A$ and derive a Bayesian (posterior) estimate. \n",
" (b) Derive the Maximum Likelihood estimate for $A$. \n",
" (c) Derive the MAP estimates for $A$. \n",
" (d) Now assume that we do not know the variance of the noise term? Describe the procedure for Bayesian estimation of both $A$ and $\\sigma^2$ (No need to fully work out to closed-form estimates). \n",
"\n",
" \n",
"- **[5]** (##) We consider the coin toss example from the notebook and use a conjugate prior for a Bernoulli likelihood function. \n",
" (a) Derive the Maximum Likelihood estimate. \n",
" (b) Derive the MAP estimate. \n",
" (c) Do these two estimates ever coincide (if so under what circumstances)? \n",
"\n",
"- **[6]** (###) Given a single observation $x_0$ from a uniform distribution $\\mathrm{Unif}[0,1/\\theta]$, where $\\theta > 0$. \n",
" (a) Show that $\\mathbb{E}[g(x_0)] = \\theta$ if and only if $\\int_0^{1/\\theta} g(u) du =1$. \n",
" (b) Show that there is no function $g$ that satisfies the condition for all $\\theta > 0$."
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Continuous Data and the Gaussian Distribution\n",
"\n",
"\n",
"- **[1]** (##) We are given an IID data set $D = \\{x_1,x_2,\\ldots,x_N\\}$, where $x_n \\in \\mathbb{R}^M$. Let's assume that the data were drawn from a multivariate Gaussian (MVG),\n",
"$$\\begin{align*}\n",
"p(x_n|\\theta) = \\mathcal{N}(x_n|\\,\\mu,\\Sigma) = |2 \\pi \\Sigma|^{-\\frac{1}{2}} \\exp\\left\\{-\\frac{1}{2}(x_n-\\mu)^T\n",
"\\Sigma^{-1} (x_n-\\mu) \\right\\}\n",
"\\end{align*}$$ \n",
" (a) Derive the log-likelihood of the parameters for these data. \n",
" (b) Derive the maximum likelihood estimates for the mean $\\mu$ and variance $\\Sigma$ by setting the derivative of the log-likelihood to zero.\n",
"\n",
"- **[2]** (#) Shortly explain why the Gaussian distribution is often preferred over other distributions with the same support?\n",
"\n",
"- **[3]** (###) Proof that the Gaussian distribution is the maximum entropy distribution over the reals with specified mean and variance.\n",
"\n",
"- **[4]** (##) Proof that a linear transformation $z=Ax+b$ of a Gaussian variable $\\mathcal{N}(x|\\mu,\\Sigma)$ is Gaussian distributed as\n",
"$$\n",
"p(z) = \\mathcal{N} \\left(z \\,|\\, A\\mu+b, A\\Sigma A^T \\right) \n",
"$$\n",
"\n",
"- **[5]** (#) Given independent variables\n",
"$x \\sim \\mathcal{N}(\\mu_x,\\sigma_x^2)$ and $y \\sim \\mathcal{N}(\\mu_y,\\sigma_y^2)$, what is the PDF for $z = A\\cdot(x -y) + b$? \n",
"\n",
"- **[6]** (##) Compute\n",
"\n",
"\\begin{equation*}\n",
" \\int_{-\\infty}^{\\infty} \\exp(-x^2)\\mathrm{d}x \\,.\n",
" \\end{equation*}\n",
" \n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Discrete Data and the Multinomial Distribution\n",
"\n",
"- **[1]** (##) We consider IID data $D = \\{x_1,x_2,\\ldots,x_N\\}$ obtained from tossing a $K$-sided die. We use a *binary selection variable*\n",
"$$x_{nk} \\equiv \\begin{cases} 1 & \\text{if $x_n$ lands on $k$th face}\\\\\n",
" 0 & \\text{otherwise}\n",
"\\end{cases}\n",
"$$\n",
"with probabilities $p(x_{nk} = 1)=\\theta_k$. \n",
" (a) Write down the probability for the $n$th observation $p(x_n|\\theta)$ and derive the log-likelihood $\\log p(D|\\theta)$. \n",
" (b) Derive the maximum likelihood estimate for $\\theta$.\n",
"\n",
"- **[2]** (#) In the notebook, Laplace's generalized rule of succession (the probability that we throw the $k$th face at the next toss) was derived as \n",
"$$\\begin{align*}\n",
"p(x_{\\bullet,k}=1|D) = \\frac{m_k + \\alpha_k }{ N+ \\sum_k \\alpha_k}\n",
"\\end{align*}$$\n",
"Provide an interpretation of the variables $m_k,N,\\alpha_k,\\sum_k$ and $\\alpha_k$.\n",
"\n",
"- **[3]** (##) Show that Laplace's generalized rule of succession can be worked out to a prediction that is composed of a prior prediction and data-based correction term. \n",
"\n",
"- **[4]** (#) Verify that \n",
" (a) the categorial distribution is a special case of the multinomial for $N=1$. \n",
" (b) the Bernoulli is a special case of the categorial distribution for $K=2$. \n",
" (c) the binomial is a special case of the multinomial for $K=2$. \n",
"\n",
"- **[5]** (###) Determine the mean, variance and mode of a Beta distribution.\n",
"\n",
"- **[6]** (###) Consider a data set of binary variables $D=\\{x_1,x_2,\\ldots,x_N\\}$ with a Bernoulli distribution $\\mathrm{Ber}(x_k|\\mu)$ as data generating distribution and a Beta prior for $\\mu$. Assume that you make $n$ observations with $x=1$ and $N-n$ observations with $x=0$. Now consider a new draw $x_\\bullet$. We are interested in computing $p(x_\\bullet|D)$. Show that the mean value for $p(x_\\bullet|D)$ lies in between the prior mean and Maximum Likelihood estimate.\n",
"\n",
""
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Regression\n",
"\n",
"\n",
"- **[1]** (#) (a) Write down the generative model for Bayesian linear ordinary regression (i.e., write the likelihood and prior). \n",
" (b) State the inference task for the weight parameter in the model. \n",
" (c) Why do we call this problem linear?\n",
"\n",
"- **[2]** (##) Consider a linear regression problem \n",
"$$\\begin{align*}\n",
"p(y\\,|\\,\\mathbf{X},w,\\beta) &= \\mathcal{N}(y\\,|\\,\\mathbf{X} w,\\beta^{-1} \\mathbf{I}) \\\\\n",
" &= \\prod_n \\mathcal{N}(y_n\\,|\\,w^T x_n,\\beta^{-1}) \n",
"\\end{align*}$$\n",
"with $y, X$ and $w$ as defined in the notebook. \n",
" (a) Work out the maximum likelihood solution for linear regression by solving\n",
"$$\n",
"\\nabla_{w} \\log p(y|X,w) = 0 \\,.\n",
"$$ \n",
" (b) Work out the MAP solution. How does it relate to the ML solution?\n",
"\n",
"- **[3]** (###) Show that the variance of the predictive distribution for linear regression decreases as more data becomes available.\n",
"\n",
"- **[4]** (#) Assume a given data set $D=\\{(x_1,y_1),(x_2,y_2),\\ldots,(x_N,y_N)\\}$ with $x \\in \\mathbb{R}^M$ and $y \\in \\mathbb{R}$. We propose a model given by the following data generating distribution and weight prior functions: \n",
"$$\\begin{equation*} p(y|x,w)\\cdot p(w)\\,. \\end{equation*}$$\n",
" (a) Write down Bayes rule for generating the posterior $p(w|D)$ from a prior and likelihood. \n",
" (b) Work out how to compute a distribution for the predicted value $y_\\bullet$, given a new input $x_\\bullet$. \n",
" \n",
"\n",
"- **[5]** (#) In the class we use the following prior for the weights:\n",
"$$\\begin{equation*}\n",
"p(w|\\alpha) = \\mathcal{N}\\left(w | 0, \\alpha^{-1} I \\right)\n",
"\\end{equation*}$$\n",
" (a) Give some considerations for choosing a Gaussian prior for the weights. \n",
" (b) We could have chosen a prior with full (not diagonal) covariance matrix $p(w|\\alpha) = \\mathcal{N}\\left(w | 0, \\Sigma \\right)$. Would that be better? Give your thoughts on that issue. \n",
" (c) Generally we choose $\\alpha$ as a small positive number. Give your thoughts on that choice as opposed to choosing a large positive value. How about choosing a negative value for $\\alpha$?\n",
"\n",
"\n",
""
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Generative Classification\n",
"\n",
"- **[1]** You have a machine that measures property $x$, the \"orangeness\" of liquids. You wish to discriminate between $C_1 = \\text{`Fanta'}$ and $C_2 = \\text{`Orangina'}$. It is known that\n",
"$$\\begin{align*}\n",
"p(x|C_1) &= \\begin{cases} 10 & 1.0 \\leq x \\leq 1.1 \\\\\n",
" 0 & \\text{otherwise}\n",
" \\end{cases} \\\\\n",
"p(x|C_2) &= \\begin{cases} 200(x - 1) & 1.0 \\leq x \\leq 1.1 \\\\\n",
"0 & \\text{otherwise}\n",
"\\end{cases}\n",
"\\end{align*}$$\n",
"The prior probabilities $p(C_1) = 0.6$ and $p(C_2) = 0.4$ are also known from experience. \n",
" (a) (##) A \"Bayes Classifier\" is given by \n",
" \n",
"$$ \n",
"\\text{Decision} = \\begin{cases} C_1 & \\text{if } p(C_1|x)>p(C_2|x) \\\\\n",
" C_2 & \\text{otherwise}\n",
" \\end{cases}\n",
"$$\n",
" Derive the optimal Bayes classifier. \n",
" \n",
" (b) (###) The probability of making the wrong decision, given $x$, is\n",
" \n",
"$$\n",
"p(\\text{error}|x)= \\begin{cases} p(C_1|x) & \\text{if we decide $C_2$} \\\\\n",
" p(C_2|x) & \\text{if we decide $C_1$}\n",
"\\end{cases}\n",
"$$\n",
" Compute the **total** error probability $p(\\text{error})$ for the Bayes classifier in this example.\n",
" \n",
"- **[2]** (#) (see Bishop exercise 4.8): Using (4.57) and (4.58) (from Bishop's book), derive the result (4.65) for the posterior class probability in the two-class generative model with Gaussian densities, and verify the results (4.66) and (4.67) for the parameters $w$ and $w0$. \n",
"\n",
"- **[3]** (###) (see Bishop exercise 4.9). \n",
"\n",
"- **[4]** (##) (see Bishop exercise 4.10). \n",
"\n",
""
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Discriminative Classification\n",
"\n",
"- **[1]** Given a data set $D=\\{(x_1,y_1),\\ldots,(x_N,y_N)\\}$, where $x_n \\in \\mathbb{R}^M$ and $y_n \\in \\{0,1\\}$. The probabilistic classification method known as *logistic regression* attempts to model these data as\n",
"$$p(y_n=1|x_n) = \\sigma(\\theta^T x_n + b)$$\n",
"where $\\sigma(x) = 1/(1+e^{-x})$ is the *logistic function*. Let's introduce shorthand notation $\\mu_n=\\sigma(\\theta^T x_n + b)$. So, for every input $x_n$, we have a model output $\\mu_n$ and an actual data output $y_n$. \n",
" (a) Express $p(y_n|x_n)$ as a Bernoulli distribution in terms of $\\mu_n$ and $y_n$. \n",
" (b) If furthermore is given that the data set is IID, show that the log-likelihood is given by\n",
"$$\n",
"L(\\theta) \\triangleq \\log p(D|\\theta) = \\sum_n \\left\\{y_n \\log \\mu_n + (1-y_n)\\log(1-\\mu_n)\\right\\}\n",
"$$ \n",
" (c) Prove that the derivative of the logistic function is given by\n",
"$$\n",
"\\sigma^\\prime(\\xi) = \\sigma(\\xi)\\cdot\\left(1-\\sigma(\\xi)\\right)\n",
"$$ \n",
" (d) Show that the derivative of the log-likelihood is\n",
"$$\n",
"\\nabla_\\theta L(\\theta) = \\sum_{n=1}^N \\left( y_n - \\sigma(\\theta^T x_n +b)\\right)x_n\n",
"$$ \n",
" (e) Design a gradient-ascent algorithm for maximizing $L(\\theta)$ with respect to $\\theta$. \n",
"\n",
"\n",
"- **[2]** Describe shortly the similarities and differences between the discriminative and generative approach to classification.\n",
" \n",
"\n",
"- **[3]** (Bishop ex.4.7) (#) Show that the logistic sigmoid function $\\sigma(a) = \\frac{1}{1+\\exp(-a)}$ satisfies the property $\\sigma(-a) = 1-\\sigma(a)$ and that its inverse is given by $\\sigma^{-1}(y) = \\log\\{y/(1-y)\\}$.\n",
" \n",
" \n",
"- **[4]** (Bishop ex.4.16) (###) Consider a binary classification problem in which each observation $x_n$ is known to belong to one of two classes, corresponding to $y_n = 0$ and $y_n = 1$. Suppose that the procedure for collecting training data is imperfect, so that training points are sometimes mislabelled. For every data point $x_n$, instead of having a value $y_n$ for the class label, we have instead a value $\\pi_n$ representing the probability that $y_n = 1$. Given a probabilistic model $p(y_n = 1|x_n,\\theta)$, write down the log-likelihood function appropriate to such a data set.\n",
" \n",
"\n",
"\n",
"\n",
"- **[5]** (###) Let $X$ be a real valued random variable with probability density\n",
"$$\n",
"p_X(x) = \\frac{e^{-x^2/2}}{\\sqrt{2\\pi}},\\quad\\text{for all $x$}.\n",
"$$\n",
"Also $Y$ is a real valued random variable with conditional density\n",
"$$\n",
"p_{Y|X}(y|x) = \\frac{e^{-(y-x)^2/2}}{\\sqrt{2\\pi}},\\quad\\text{for all $x$ and $y$}. \n",
"$$\n",
" (a) Give an (integral) expression for $p_Y(y)$. Do not try to evaluate the integral. \n",
" (b) Approximate $p_Y(y)$ using the Laplace approximation.\n",
" Give the detailed derivation, not just the answer.\n",
"Hint: You may use the following results.\n",
"Let \n",
"$$g(x) = \\frac{e^{-x^2/2}}{\\sqrt{2\\pi}}$$\n",
"and\n",
"$$\n",
"h(x) = \\frac{e^{-(y-x)^2/2}}{\\sqrt{2\\pi}}$$\n",
"for some real value $y$. Then:\n",
"$$\\begin{align*}\n",
"\\frac{\\partial}{\\partial x} g(x) &= -xg(x) \\\\\n",
"\\frac{\\partial^2}{\\partial x^2} g(x) &= (x^2-1)g(x) \\\\\n",
"\\frac{\\partial}{\\partial x} h(x) &= (y-x)h(x) \\\\\n",
"\\frac{\\partial^2}{\\partial x^2} h(x) &= ((y-x)^2-1)h(x) \n",
"\\end{align*}$$ \n",
"\n",
"\n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"\n",
"### Latent Variable Models and Variational Bayes\n",
"\n",
"\n",
"- **[1]** (##) For a Gaussian mixture model, given by generative equations\n",
"\n",
"$$\n",
"p(x,z) = \\prod_{k=1}^K (\\underbrace{\\pi_k \\cdot \\mathcal{N}\\left( x | \\mu_k, \\Sigma_k\\right) }_{p(x,z_{k}=1)})^{z_{k}} \n",
"$$\n",
"\n",
"proof that the marginal distribution for observations $x_n$ evaluates to \n",
"\n",
"$$\n",
"p(x) = \\sum_{j=1}^K \\pi_k \\cdot \\mathcal{N}\\left( x | \\mu_j, \\Sigma_j \\right) \n",
"$$\n",
"\n",
"\n",
"- **[2]** (#) Given the free energy functional $F[q] = \\sum_z q(z) \\log \\frac{q(z)}{p(x,z)}$, proof the [EE, DE and AC decompositions](https://nbviewer.jupyter.org/github/bertdv/BMLIP/blob/master/lessons/notebooks/Latent-Variable-Models-and-VB.ipynb#fe-decompositions). \n",
" \n",
"\n",
"- **[3]** (#) The Free energy functional $\\mathrm{F}[q] = -\\sum_z q(z) \\log p(x,z) - \\sum_z q(z) \\log \\frac{1}{q(z)}$ decomposes into \"Energy minus Entropy\". So apparently the entropy of the posterior $q(z)$ is maximized. This entropy maximization may seem puzzling at first because inference should intuitively lead to *more* informed posteriors, i.e., posterior distributions whose entropy is smaller than the entropy of the prior. Explain why entropy maximization is still a reasonable objective. \n",
"\n",
"- **[4]** (#) Explain the following update rule for the mean of the Gaussian cluster-conditional data distribution (from the example about mean-field updating of a Gaussian mixture model):\n",
"\n",
"$$\n",
"m_k = \\frac{1}{\\beta_k} \\left( \\beta_0 m_0 + N_k \\bar{x}_k \\right) \\tag{B-10.61} \n",
"$$\n",
" \n",
"- **[5]** (##) Consider a model $p(x,z|\\theta)$, where $D=\\{x_1,x_2,\\ldots,x_N\\}$ is observed, $z$ are unobserved variables and $\\theta$ are parameters. The EM algorithm estimates the parameters by iterating over the following two equations ($i$ is the iteration index):\n",
"\n",
"$$\\begin{align*}\n",
"q^{(i)}(z) &= p(z|D,\\theta^{(i-1)}) \\\\\n",
"\\theta^{(i)} &= \\arg\\max_\\theta \\sum_z q^{(i)}(z) \\cdot \\log p(D,z|\\theta)\n",
"\\end{align*}$$\n",
"\n",
"Proof that this algorithm minimizes the Free Energy functional \n",
"$$\\begin{align*}\n",
"F[q,\\theta] = \\sum_z q(z) \\log \\frac{q(z)}{p(D,z|\\theta)} \n",
"\\end{align*}$$\n",
"\n",
"- **[6]** (###) Consult the internet on what *overfitting* and *underfitting* is and then explain how FE minimization finds a balance between these two (unwanted) extremes.\n",
" \n",
"\n",
"- **[7]** (##) Consider a model $p(x,z|\\theta) = p(x|z,\\theta) p(z|\\theta)$ where $x$ and $z$ relate to observed and unobserved variables, respectively. Also available is an observed data set $D=\\left\\{x_1,x_2,\\ldots,x_N\\right\\}$. One iteration of the EM-algorithm for estimating the parameters $\\theta$ is described by ($m$ is the iteration counter)\n",
"$$\n",
"\\hat{\\theta}^{(m+1)} := \\arg \\max_\\theta \\left(\\sum_z p(z|x=D,\\hat{\\theta}^{(m)}) \\log p(x=D,z|\\theta) \\right) \\,.\n",
"$$\n",
"\n",
" (a) Apparently, in order to execute EM, we need to work out an expression for the 'responsibility' $p(z|x=D,\\hat{\\theta}^{(m)})$. Use Bayes rule to show how we can compute the responsibility that allows us to execute an EM step. \n",
" (b) Why do we need multiple iterations in the EM algorithm? \n",
" (c) Why can't we just use simple maximum log-likelihood to estimate parameters, as described by \n",
"$$\n",
"\\hat{\\theta} := \\arg \\max_\\theta \\log p(x=D,z|\\theta) \\,?\n",
"$$ \n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Dynamic Models \n",
"\n",
"- **[1]** (##) Given the Markov property\n",
"\\begin{equation*}\n",
"p(x_n|x_{n-1},x_{n-2},\\ldots,x_1) = p(x_n|x_{n-1}) \\tag{A1}\n",
"\\end{equation*}\n",
"proof that, for any $n$,\n",
"\\begin{align*}\n",
"p(x_n,x_{n-1},&\\ldots,x_{k+1},x_{k-1},\\ldots,x_1|x_k) = \\\\\n",
"&p(x_n,x_{n-1},\\ldots,x_{k+1}|x_k) \\cdot p(x_{k-1},x_{k-2},\\ldots,x_1|x_k) \\tag{A2}\\,.\n",
"\\end{align*}\n",
"In other words, proof that, if the Markov property A1 holds, then, given the \"present\" ($x_k$), the \"future\" $(x_n,x_{n-1},\\ldots,x_{k+1})$ is _independent_ of the \"past\" $(x_{k-1},x_{k-2},\\ldots,x_1)$.\n",
"\n",
"- **[2]** (#) \n",
" (a) What's the difference between a hidden Markov model and a linear Dynamical system? \n",
" (b) For the same number of state variables, which of these two models has a larger memory capacity, and why? \n",
" \n",
"- **[3]** (#) \n",
"(a) What is the 1st-order Markov assumption? \n",
"(b) Derive the joint probability distribution $p(x_{1:T},z_{0:T})$ (where $x_t$ and $z_t$ are observed and latent variables respectively) for the state-space model with transition and observation models $p(z_t|z_{t-1})$ and $p(x_t|z_t)$. \n",
"(c) What is a Hidden Markov Model (HMM)? \n",
"(d) What is a Linear Dynamical System (LDS)? \n",
"(e) What is a Kalman Filter? \n",
"(f) How does the Kalman Filter relate to the LDS? \n",
"(g) Explain the popularity of Kalman filtering and HMMs? \n",
"(h) How relates a HMM to a GMM? \n",
"\n",
" \n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"### Intelligent Agents and Active Inference\n",
"\n",
"- **[1]** (##) I asked you to watch a video segment (https://www.vibby.com/watch?vib=71iPtUJxd) where Karl Friston talks about two main approaches to goal-directed acting by agents: (1) choosing actions that maximize (the expectation of) a value function $V(s)$ of the state ($s$) of the environment; or (2) choosing actions that minimize a functional ($F[q(s)]$) of *beliefs* ($q(s)$) over environmental states ($s$). Discuss the advantage of the latter appraoch. \n",
"\n",
"- **[2]** (#) The *good regulator theorem* states that a \"successful and efficient\" controller of the world must contain a model of the world. But it's hard to imagine how just learning a model of the world leads to goal-directed behavior, like learning how to read or drive a car. Which other ingredient do we need to get learning agents to behave as goal-directed agents? \n",
"\n",
"- **[3]** (##) The figure below reflects the state of a factor graph realization of an active inference agent after having pushed action $a_t$ onto the environment and having received observation $x_t$. In this graph, the variables $x_\\bullet$, $u_\\bullet$ and $s_\\bullet$ correspond to observations, and unobserved control and internal states respectively. Copy the figure onto your sheet and draw a message passing schedule to infer a posterior belief (i.e. after observing $x_t$) over the next control state $u_{t+1}$. \n",
"\n",
"\n",
"\n",
"\n",
"- **[4]** (##) The Free Energy Principle (FEP) is a theory about biological self-organization, in particular about how brains develop through interactions with their environment. Which of the following statements is not consistent with FEP (and explain your answer): \n",
" (a) We act to fullfil our predictions about future sensory inputs. \n",
" (b) Perception is inference about the environmental causes of our sensations. \n",
" (c) Our actions aim to reduce the complexity of our model of the environment. \n",
" \n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
""
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"slideshow": {
"slide_type": "skip"
}
},
"outputs": [
{
"data": {
"text/html": [
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n",
"\r\n"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"open(\"../../styles/aipstyle.html\") do f\n",
" display(\"text/html\", read(f,String))\n",
"end"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"@webio": {
"lastCommId": null,
"lastKernelId": null
},
"anaconda-cloud": {},
"celltoolbar": "Slideshow",
"kernelspec": {
"display_name": "Julia 1.5.2",
"language": "julia",
"name": "julia-1.5"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 4
}