|
| 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: 2024-10-25 11:59:37 |
| 9 | + |
| 10 | +title: "Linear combination of bivariate normal random variables" |
| 11 | +chapter: "Probability Distributions" |
| 12 | +section: "Multivariate continuous distributions" |
| 13 | +topic: "Multivariate normal distribution" |
| 14 | +theorem: "Linear combination of bivariate normal random variables" |
| 15 | + |
| 16 | +sources: |
| 17 | + - authors: "Wikipedia" |
| 18 | + year: 2024 |
| 19 | + title: "Misconceptions about the normal distribution" |
| 20 | + in: "Wikipedia, the free encyclopedia" |
| 21 | + pages: "retrieved on 2024-10-25" |
| 22 | + url: "https://en.wikipedia.org/wiki/Misconceptions_about_the_normal_distribution" |
| 23 | + |
| 24 | +proof_id: "P475" |
| 25 | +shortcut: "bvn-lincomb" |
| 26 | +username: "JoramSoch" |
| 27 | +--- |
| 28 | + |
| 29 | + |
| 30 | +**Theorem:** Let $X$ and $Y$ follow a [bivariate normal distribution](/D/bvn): |
| 31 | + |
| 32 | +$$ \label{eq:bvn} |
| 33 | +\left[ \begin{matrix} X \\ Y \end{matrix} \right] \sim |
| 34 | +\mathcal{N}\left( \left[ \begin{matrix} \mu_1 \\ \mu_2 \end{matrix} \right], \left[ \begin{matrix} \sigma_1^2 & \sigma_{12} \\ \sigma_{12} & \sigma_2^2 \end{matrix} \right] \right) \; . |
| 35 | +$$ |
| 36 | + |
| 37 | +Then, any linear combination of $X$ and $Y$ follows a [univariate normal distribution](/D/norm): |
| 38 | + |
| 39 | +$$ \label{eq:bvn-lincomb} |
| 40 | +Z = a X + b Y \sim |
| 41 | +\mathcal{N}\left( a \mu_1 + b \mu_2, a^2 \sigma_1^2 + 2ab \sigma_{12} + b^2 \sigma_2^2 \right) \; . |
| 42 | +$$ |
| 43 | + |
| 44 | + |
| 45 | +**Proof:** The [linear transformation theorem for the multivariate normal distribution](/P/mvn-ltt) states that |
| 46 | + |
| 47 | +$$ \label{eq:mvn-ltt} |
| 48 | +X \sim \mathcal{N}(\mu, \Sigma) \quad \Rightarrow \quad Y = AX + c \sim \mathcal{N}(A\mu + c, A \Sigma A^\mathrm{T}) |
| 49 | +$$ |
| 50 | + |
| 51 | +where $X \in \mathbb{R}^n$, $A \in \mathbb{R}^{n \times n}$ and $c \in \mathbb{R}^n$. In the present case, we have |
| 52 | + |
| 53 | +$$ \label{eq:X-A-a} |
| 54 | +X \in \mathbb{R}^2 |
| 55 | +\quad \text{and} \quad |
| 56 | +A = \left[ \begin{matrix} a & b \end{matrix} \right] \in \mathbb{R}^{1 \times 2} |
| 57 | +\quad \text{and} \quad |
| 58 | +c = 0 \in \mathbb{R} \; , |
| 59 | +$$ |
| 60 | + |
| 61 | +such that |
| 62 | + |
| 63 | +$$ \label{eq:Z-X-Y} |
| 64 | +Z |
| 65 | += A \left[ \begin{matrix} X \\ Y \end{matrix} \right] + c |
| 66 | += \left[ \begin{matrix} a & b \end{matrix} \right] \left[ \begin{matrix} X \\ Y \end{matrix} \right] + 0 |
| 67 | += a X + b Y \; . |
| 68 | +$$ |
| 69 | + |
| 70 | +Combining \eqref{eq:mvn-ltt}, \eqref{eq:bvn} and \eqref{eq:Z-X-Y}, it follows that |
| 71 | + |
| 72 | +$$ \label{eq:bvn-lincomb-qed} |
| 73 | +\begin{split} |
| 74 | +Z |
| 75 | +&\sim \mathcal{N}\left( \left[ \begin{matrix} a & b \end{matrix} \right] \left[ \begin{matrix} \mu_1 \\ \mu_2 \end{matrix} \right], \left[ \begin{matrix} a & b \end{matrix} \right] \left[ \begin{matrix} \sigma_1^2 & \sigma_{12} \\ \sigma_{12} & \sigma_2^2 \end{matrix} \right] \left[ \begin{matrix} a \\ b \end{matrix} \right] \right) \\ |
| 76 | +&\sim \mathcal{N}\left( a \mu_1 + b \mu_2, \left[ \begin{matrix} a \sigma_1^2 + b \sigma_{12} & a \sigma_{12} + b \sigma_2^2 \end{matrix} \right] \left[ \begin{matrix} a \\ b \end{matrix} \right] \right) \\ |
| 77 | +&\sim \mathcal{N}\left( a \mu_1 + b \mu_2, (a^2 \sigma_1^2 + ab \sigma_{12}) + (ab \sigma_{12} + b^2 \sigma_2^2) \right) \\ |
| 78 | +&\sim \mathcal{N}\left( a \mu_1 + b \mu_2, a^2 \sigma_1^2 + 2ab \sigma_{12} + b^2 \sigma_2^2 \right) \; . |
| 79 | +\end{split} |
| 80 | +$$ |
0 commit comments