This page is a collection of all theorems taught in EECS126: Probability and Random Processes, Spring 2021. A good reference!
Link to PDF version here.
Probability Basics
Conditional Probability
P ( A ∣ B ) = P ( A ∩ B ) P ( B ) , P ( B ) > 0 P(A | B) = \frac{P(A \cap B)}{P(B)}, P(B) > 0 P ( A ∣ B ) = P ( B ) P ( A ∩ B ) , P ( B ) > 0
Total Probability Theorem
P ( B ) = ∑ i = 1 n P ( A i ) P ( B ∣ A i ) P(B) = \sum_{i=1}^{n} P(A_i)P(B | A_i) P ( B ) = i = 1 ∑ n P ( A i ) P ( B ∣ A i )
Bayes Rules
P ( A i ∣ B ) = P ( A i ) P ( B ∣ A i ) P ( B ) P(A_i | B) = \frac{P(A_i)P(B|A_i)}{P(B)} P ( A i ∣ B ) = P ( B ) P ( A i ) P ( B ∣ A i )
Union Bound
P ( ⋃ i = 1 ∞ A i ) ≤ ∑ i = 1 ∞ P ( A i ) P(\bigcup_{i=1}^{\infty} A_i) \leq \sum_{i=1}^{\infty}P(A_i) P ( i = 1 ⋃ ∞ A i ) ≤ i = 1 ∑ ∞ P ( A i )
Independence
P ( A ∣ B ) = P ( A ) ⟺ A and B are independent P(A | B) = P(A) \iff \text{A and B are independent} P ( A ∣ B ) = P ( A ) ⟺ A and B are independent
Conditional Independence
P ( A ∩ B ∣ C ) = P ( A ∣ C ) P ( B ∣ C ) ⟹ A and B conditionally independent P(A \cap B | C) = P(A | C) P(B | C) \implies \text{A and B conditionally independent} P ( A ∩ B ∣ C ) = P ( A ∣ C ) P ( B ∣ C ) ⟹ A and B conditionally independent
Independence of Several Events
P ( ∩ i ∈ S A i ) = ∏ i ∈ S P ( A i ) P(\cap_{i \in S} A_i) = \prod_{i \in S} P(A_i) P ( ∩ i ∈ S A i ) = i ∈ S ∏ P ( A i )
Counting Permutations of Size k k k in n n n Objects
n P k = n ! ( n − k ) ! ^nP_k = \frac{n!}{(n-k)!} n P k = ( n − k )! n !
Counting Ways to Choose k k k Objects in n n n Objects
( n k ) = n ! k ! ( n − k ) ! {n\choose k} = \frac{n!}{k!(n-k)!} ( k n ) = k ! ( n − k )! n !
Counting Ways To Partition n n n Objects into n i n^i n i Groups
( n n 1 , n 2 . . . n k ) = n ! n 1 ! n 2 ! . . . n k ! {n\choose n_1, n_2... n_k} = \frac{n!}{n_1! n_2!...n_k!} ( n 1 , n 2 ... n k n ) = n 1 ! n 2 ! ... n k ! n !
Discrete Random variables
Bernoulli Random Variable
P ( X = k ) = { 1 with probability p 0 with probability 1 − p
\begin{aligned}
P(X = k) = \begin{cases}
1 & \text{ with probability } p \newline
0 & \text{ with probability } 1-p
\end{cases}
\end{aligned}
P ( X = k ) = { 1 0 with probability p with probability 1 − p
E ( X ) = p \mathbb{E}(X) = p E ( X ) = p
v a r ( X ) = p ( 1 − p ) var(X) = p(1-p) v a r ( X ) = p ( 1 − p )
Binomial Random Variable
P ( X = k ) = ( n k ) p k ( 1 − p ) n − k P(X = k) = {n\choose k} p^k (1-p)^{n-k} P ( X = k ) = ( k n ) p k ( 1 − p ) n − k
E ( X ) = n p \mathbb{E}(X) = np E ( X ) = n p
v a r ( X ) = n p ( 1 − p ) var(X) = np(1-p) v a r ( X ) = n p ( 1 − p )
Geometric Random Variable
P ( X = k ) = ( 1 − p ) k − 1 p P(X = k) = (1-p)^{k-1} p P ( X = k ) = ( 1 − p ) k − 1 p
E ( X ) = 1 p \mathbb{E}(X) = \frac{1}{p} E ( X ) = p 1
v a r ( X ) = 1 − p p 2 var(X) = \frac{1-p}{p^2} v a r ( X ) = p 2 1 − p
Poisson Random Variable
P ( X λ = k ) = e − λ λ k k ! P(X_\lambda = k) = \frac{e^{-\lambda}\lambda^k}{k!} P ( X λ = k ) = k ! e − λ λ k
E ( X ) = λ \mathbb{E}(X) = \lambda E ( X ) = λ
v a r ( X ) = λ var(X) = \lambda v a r ( X ) = λ
Linearity of a Poisson RV
P o i s s o n ( λ ) + P o i s s o n ( μ ) ∼ P o i s s o n ( λ + μ ) Poisson(\lambda) + Poisson(\mu) \sim Poisson(\lambda + \mu) P o i sso n ( λ ) + P o i sso n ( μ ) ∼ P o i sso n ( λ + μ )
Uniform Random Variable
P ( X = k ) = { 1 b − a + 1 k ∈ [ a , b ] 0 otherwise \begin{equation*}
P(X = k) = \begin{cases}
\frac{1}{b-a+1} & \ k \in [a,b] \\
0 & \text{ otherwise}
\end{cases}
\end{equation*} P ( X = k ) = { b − a + 1 1 0 k ∈ [ a , b ] otherwise
Joint PMFs
P X , Y ( x , y ) = Pr ( X = x , Y = y ) P_{X, Y} (x, y) = \Pr(X = x, Y = y) P X , Y ( x , y ) = Pr ( X = x , Y = y )
P X ( x ) = ∑ y P X , Y ( x , y ) and vice versa P_X(x) = \sum_{y} P_{X, Y}(x, y) \text{ and vice versa} P X ( x ) = y ∑ P X , Y ( x , y ) and vice versa
Conditional PMFs
P X ∣ A ( X = x ∣ A ) = P ( { X = x } ∩ A ) P ( A ) P_{X|A}(X = x|A) = \frac{P(\{X = x\} \cap A)}{P(A)} P X ∣ A ( X = x ∣ A ) = P ( A ) P ({ X = x } ∩ A )
P X Y ( x ∣ y ) = P X , Y ( x , y ) P Y ( y ) P_{X_Y}(x|y) = \frac{P_{X,Y}(x,y)}{P_Y(y)} P X Y ( x ∣ y ) = P Y ( y ) P X , Y ( x , y )
Expectation, Variance and Covariance
Expectation
E ( X ) = ∑ x x P ( X = x ) \mathbb{E}(X) = \sum_{x} xP(X = x) E ( X ) = x ∑ x P ( X = x )
Law of The Unconscious Statistician
E ( g ( X ) ] = ∑ x g ( x ) P ( X = x ) \mathbb{E}(g(X)] = \sum_{x} g(x)P(X = x) E ( g ( X )] = x ∑ g ( x ) P ( X = x )
Variance
v a r ( X ) = E [ ( X − E ( X ) ) 2 ] ≥ 0 var(X) = \mathbb{E}[(X - \mathbb{E}(X))^2] \geq 0 v a r ( X ) = E [( X − E ( X ) ) 2 ] ≥ 0
Standard Deviation
σ = v a r \sigma = \sqrt{var} σ = v a r
Linearity of Expectation
E ( a X + b Y ) = a E ( X ) + b E ( Y ) \mathbb{E}(aX + bY) = a\mathbb{E}(X) + b\mathbb{E}(Y) E ( a X + bY ) = a E ( X ) + b E ( Y )
Expectation of Joint Distribution
E ( g ( X , Y ) ) = ∑ x ∑ y g ( x , y ) P X , Y ( x , y ) \mathbb{E}(g(X, Y)) = \sum_{x}\sum_{y} g(x, y)P_{X, Y}(x, y) E ( g ( X , Y )) = x ∑ y ∑ g ( x , y ) P X , Y ( x , y )
Variance of a Sum of Random Variables
v a r ( X + Y ) = v a r ( X ) + v a r ( Y ) + 2 c o v ( X , Y ) var(X+Y) = var(X) + var(Y) + 2cov(X, Y) v a r ( X + Y ) = v a r ( X ) + v a r ( Y ) + 2 co v ( X , Y )
Conditional Expectation
E ( X ∣ Y = y ) = ∑ x x P X ∣ Y ( x ∣ y ) \mathbb{E}(X | Y = y) = \sum_{x} x\ P_{X|Y}(x|y) E ( X ∣ Y = y ) = x ∑ x P X ∣ Y ( x ∣ y )
Total Expectation Theorem
E ( X ) = ∑ y P Y ( y ) E ( X ∣ Y = y ) \mathbb{E}(X) = \sum_{y} P_Y(y)\mathbb{E}(X|Y=y) E ( X ) = y ∑ P Y ( y ) E ( X ∣ Y = y )
Iterated Expectation
E ( X ) = E ( E ( X ∣ Y ) ) \mathbb{E}(X) = \mathbb{E}(\mathbb{E}(X|Y)) E ( X ) = E ( E ( X ∣ Y ))
Tower Property
E [ E [ X ∣ Y ] g ( Y ) ] = E [ X g ( Y ) ] \mathbb{E}[\mathbb{E}[X|Y]g(Y)] = \mathbb{E}[Xg(Y)] E [ E [ X ∣ Y ] g ( Y )] = E [ X g ( Y )]
Expectation of Independent Variables
E ( X Y ) = E ( X ) E ( Y ) if X, Y independent \mathbb{E}(XY) = \mathbb{E}(X)\mathbb{E}(Y) \text{ if X, Y independent} E ( X Y ) = E ( X ) E ( Y ) if X, Y independent
Covariance
c o v ( X , Y ) = E ( X Y ) − E ( X ) E ( Y ) cov(X,Y) = \mathbb{E}(XY) - \mathbb{E}(X)\mathbb{E}(Y) co v ( X , Y ) = E ( X Y ) − E ( X ) E ( Y )
Correlation Coefficient
ρ ( X , Y ) = c o v ( X , Y ) V a r ( X ) V a r ( Y ) \rho(X, Y) = \frac{cov(X,Y)}{\sqrt{Var(X)Var(Y)}} ρ ( X , Y ) = V a r ( X ) V a r ( Y ) co v ( X , Y )
∣ ρ ∣ ≤ 1 |\rho| \leq 1 ∣ ρ ∣ ≤ 1
Variance of Two Independent Variables
V a r [ X Y ] = E [ X 2 ] E [ Y 2 ] − E [ X ] 2 E [ Y ] 2 Var[XY] = \mathbb{E}[X^2]\mathbb{E}[Y^2] - \mathbb{E}[X]^2\mathbb{E}[Y]^2 V a r [ X Y ] = E [ X 2 ] E [ Y 2 ] − E [ X ] 2 E [ Y ] 2
Law of Total Variance
v a r ( X ) = V a r ( E ( X ∣ Y ) ) + E ( v a r ( X ∣ Y ) ) var(X) = Var(\mathbb{E}(X|Y)) + \mathbb{E}(var(X|Y)) v a r ( X ) = V a r ( E ( X ∣ Y )) + E ( v a r ( X ∣ Y ))
Continuous Random Variables
Probability Density Functions
P ( X ∈ [ a , b ] ) = ∫ a b f X ( x ) d x P(X \in [a, b]) = \int_{a}^{b} f_X(x)dx P ( X ∈ [ a , b ]) = ∫ a b f X ( x ) d x
Cumulative Ditribution Function
F X ( x ) = ∫ − ∞ x f ( t ) d t F_X(x) = \int_{-\infty}^{x} f(t)dt F X ( x ) = ∫ − ∞ x f ( t ) d t
Uniform Distribution
f X ( x ) = 1 b − a , a < x < b E ( X ) = a + b 2 v a r ( X ) = ( b − a ) 2 12
\begin{aligned}
f_X(x) &= \frac{1}{b-a},\ a<x<b \newline
\mathbb{E}(X) &= \frac{a+b}{2} \newline
var(X) &= \frac{(b-a)^2}{12}
\end{aligned}
f X ( x ) E ( X ) v a r ( X ) = b − a 1 , a < x < b = 2 a + b = 12 ( b − a ) 2
Exponential Distribution
f X ( x ) = λ e − λ x , x > 0 F X ( x ) = 1 − e − λ x E ( X ) = 1 λ v a r ( X ) = 1 λ 2
\begin{aligned}
f_X(x) &= \lambda e^{-\lambda x},\ x>0 \newline
F_X(x) &= 1 - e^{-\lambda x} \newline
\mathbb{E}(X) &= \frac{1}{\lambda} \newline
var(X) &= \frac{1}{\lambda^2}
\end{aligned}
f X ( x ) F X ( x ) E ( X ) v a r ( X ) = λ e − λ x , x > 0 = 1 − e − λ x = λ 1 = λ 2 1
Gaussian Distribution
f X ( X ) = 1 2 π σ 2 e − ( x − μ ) 2 / 2 σ 2 f_X(X) = \frac{1}{\sqrt{2\pi\sigma^2}}e^{-(x-\mu)^2/2\sigma^2} f X ( X ) = 2 π σ 2 1 e − ( x − μ ) 2 /2 σ 2
Sum of Two Gaussian Variables
a N ( μ 1 , σ 1 2 ) + b N ( μ 2 , σ 2 2 ) ∼ N ( a μ 1 + b μ 2 , a 2 σ 1 2 + b 2 σ 2 2 ) aN(\mu_1, \sigma^2_1) + bN(\mu_2, \sigma^2_2) \sim N(a\mu_1 + b\mu_2, a^2\sigma^2_1 + b^2\sigma_2^2) a N ( μ 1 , σ 1 2 ) + b N ( μ 2 , σ 2 2 ) ∼ N ( a μ 1 + b μ 2 , a 2 σ 1 2 + b 2 σ 2 2 )
Joint PDFs
f X ∣ Y ( x ∣ y ) = f X , y ( x , y ) f Y ( y ) f_{X|Y}(x|y) = \frac{f_{X, y}(x, y)}{f_Y(y)} f X ∣ Y ( x ∣ y ) = f Y ( y ) f X , y ( x , y )
Independence of Continuous Variables
f X , Y ( x , y ) = f x ( x ) f Y ( y ) f_{X, Y}(x, y) = f_x(x)f_Y(y) f X , Y ( x , y ) = f x ( x ) f Y ( y )
Order Statistics
Smallest RV in a set of RVs
Let Y = min 1 ≤ k ≤ n X k , iid with CDF F X \text{Let } Y = \min_{1 \leq k \leq n} X_k \text{ , iid with CDF $F_X$} Let Y = 1 ≤ k ≤ n min X k , iid with CDF F X
F Y ( y ) = 1 − ( 1 − F X ( y ) ) n F_Y(y) = 1 - (1 - F_X(y))^n F Y ( y ) = 1 − ( 1 − F X ( y ) ) n
Largest RV in a set of RVs
Let Y = max 1 ≤ k ≤ n X k , iid with CDF F X \text{Let } Y = \max_{1 \leq k \leq n} X_k \text{ , iid with CDF $F_X$} Let Y = 1 ≤ k ≤ n max X k , iid with CDF F X
F Y ( y ) = ( F X ( y ) ) n F_Y(y) = (F_X(y))^n F Y ( y ) = ( F X ( y ) ) n
Convolution
Discrete Convolution
p Z ( z ) = P ( X + Y = z ) = ∑ x P ( X = x , Y = z − x ) = ∑ x P x ( x ) P Y ( z − x ) if X, Y independent
\begin{aligned}
p_Z(z) &= P(X+Y=z) = \sum_{x}P(X=x, Y=z-x) \newline
&= \sum_{x}P_x(x)P_Y(z-x) \text{ if X, Y independent}
\end{aligned}
p Z ( z ) = P ( X + Y = z ) = x ∑ P ( X = x , Y = z − x ) = x ∑ P x ( x ) P Y ( z − x ) if X, Y independent
Continuous Convolution
f Z ( z ) = ∫ − ∞ ∞ f X ( x ) f Y ( z − x ) d x
f_Z(z) = \int_{-\infty}^{\infty} f_X(x)f_Y(z-x)dx
f Z ( z ) = ∫ − ∞ ∞ f X ( x ) f Y ( z − x ) d x
Moment Generating Function
MGF for a RV
M x ( s ) = E [ e s x ] = ∫ − ∞ ∞ e s x f X ( x ) d x
\begin{aligned}
M_x(s) &= \mathbb{E}[e^{sx}] \newline
&=\int_{-\infty}^{\infty} e^{sx} f_X(x)dx
\end{aligned}
M x ( s ) = E [ e s x ] = ∫ − ∞ ∞ e s x f X ( x ) d x
Derivative of an MGF
d n M ( s ) d s n ∣ s = 0 = ∫ x n f ( x ) d x = E [ X n ]
\frac{d^nM(s)}{ds^n} |_{s=0} = \int x^nf(x)dx = \mathbb{E}[X^n]
d s n d n M ( s ) ∣ s = 0 = ∫ x n f ( x ) d x = E [ X n ]
MGF of a Poisson RV
M ( s ) = e λ ( e s − 1 ) M(s) = e^{\lambda (e^s-1)} M ( s ) = e λ ( e s − 1 )
MGF of a Exponential RV
M ( s ) = λ λ − s , s < λ M(s) = \frac{\lambda}{\lambda - s}\text{ , $s < \lambda$} M ( s ) = λ − s λ , s < λ
MGF of the Standard Normal Gaussian RV
M ( s ) = e s 2 / 2 M(s) = e^{s^2/2} M ( s ) = e s 2 /2
Moments of Standard Normal RV
E ( X m ) = { 0 , m odd 2 − m / 2 m ! ( m / 2 ) ! , m even
\mathbb{E}(X^m) = \begin{cases}
0 & \text{ , m odd} \newline
2^{-m/2}\frac{m!}{(m/2)!} & \text{ , m even}
\end{cases}
E ( X m ) = { 0 2 − m /2 ( m /2 )! m ! , m odd , m even
MGF of a Geometric RV
M ( s ) = p e s 1 − ( 1 − p ) e s M(s) = \frac{pe^s}{1- (1-p)e^s} M ( s ) = 1 − ( 1 − p ) e s p e s
MGF of a Bernoulli RV
M ( s ) = 1 − p + p e s M(s) = 1 - p + pe^s M ( s ) = 1 − p + p e s
MGF of a Binomial RV
M ( s ) = ( 1 − p + p e s ) n M(s) = (1 - p + pe^s)^n M ( s ) = ( 1 − p + p e s ) n
MGF of a Uniform RV
M ( s ) = { e b s − e a s s ( b − a ) s ≠ 0 1 s = 0
\begin{equation*}
M(s) = \begin{cases}
\frac{e^{bs} - e^{as}}{s(b-a)} &\ s \neq 0 \newline
1 &\ s = 0
\end{cases}
\end{equation*}
M ( s ) = { s ( b − a ) e b s − e a s 1 s = 0 s = 0
MGF of a Sum of RVs
Let Z = ∑ X i M Z ( s ) = ∏ M X i ( s )
\begin{aligned}
\text{Let } Z &= \sum X_i \newline
M_Z(s) &= \prod M_{X_i}(s)
\end{aligned}
Let Z M Z ( s ) = ∑ X i = ∏ M X i ( s )
MGF of a Y = a T X Y = a^TX Y = a T X , X is Gaussian Vector
M Y ( s ) = M X ( s a ) = exp ( s ( a T μ x ) + 1 2 s 2 a T Σ a ) M_Y(s) = M_X(sa) = \exp{(s(a^T\mu_x) +\frac{1}{2}s^2a^T\Sigma a)} M Y ( s ) = M X ( s a ) = exp ( s ( a T μ x ) + 2 1 s 2 a T Σ a )
Bounds
Markov Inequality
P ( X ≥ a ) ≤ E [ X ] a P(X \geq a) \leq \frac{\mathbb{E}[X]}{a} P ( X ≥ a ) ≤ a E [ X ]
Chebyshev’s Inequality
P ( ∣ X − μ ∣ ≥ c ) ≤ σ 2 c 2 P(|X - \mu| \geq c) \leq \frac{\sigma^2}{c^2} P ( ∣ X − μ ∣ ≥ c ) ≤ c 2 σ 2
Chernoff Bound
P ( X ≥ a ) ≤ E [ e s x ] e s a , s > 0 P(X \geq a) \leq \frac{\mathbb{E}[e^{sx}]}{e^{sa}} \text{ , $s > 0$} P ( X ≥ a ) ≤ e s a E [ e s x ] , s > 0
P ( X ≤ a ) ≤ M ( s ) e s a , s ≤ 0 P(X \leq a) \leq \frac{M(s)}{e^{sa}} \text{ , $s \leq 0$} P ( X ≤ a ) ≤ e s a M ( s ) , s ≤ 0
Jensen Inequality
f ( E ( x ) ) ≤ E [ f ( x ) ] , f is convex, f ′ ′ ( x ) > 0 f(\mathbb{E}(x)) \leq \mathbb{E}[f(x)] \text{ , f is convex, $f''(x) > 0$} f ( E ( x )) ≤ E [ f ( x )] , f is convex, f ′′ ( x ) > 0
Weak Law of Large Numbers
lim n → ∞ P ( ∣ 1 n ∑ i = 1 n X i − E [ X ] ∣ ≥ ϵ ) = 0 \lim_{n \xrightarrow{} \infty} P(|\frac{1}{n}\sum_{i=1}^{n}X_i - \mathbb{E}[X]| \geq \epsilon) = 0 n ∞ lim P ( ∣ n 1 i = 1 ∑ n X i − E [ X ] ∣ ≥ ϵ ) = 0
Strong Law of Large Numbers
P ( lim n → ∞ M n = μ ) = 1 P(\lim_{n\xrightarrow{} \infty} M_n = \mu) = 1 P ( n ∞ lim M n = μ ) = 1
Central Limit Theorem
Define Z = S n − n μ n σ F Z ( z ) → ϕ ( z )
\begin{aligned}
\text{Define } Z &= \frac{S_n - n\mu}{\sqrt{n}\sigma} \newline
F_Z(z) &\xrightarrow{} \phi(z)
\end{aligned}
Define Z F Z ( z ) = n σ S n − n μ ϕ ( z )
Convergences
Almost Sure Convergence
P ( lim n → ∞ X n = X ) = 1 P(\lim_{n \xrightarrow{} \infty} X_n = X) = 1 P ( n ∞ lim X n = X ) = 1
Convergence in Probability
lim n → ∞ P ( ∣ X n − X ∣ ≥ ϵ ) = 0 \lim_{n \xrightarrow{} \infty} P(|X_n - X| \geq \epsilon) = 0 n ∞ lim P ( ∣ X n − X ∣ ≥ ϵ ) = 0
Convergence in Distribution
lim n → ∞ F X n ( x ) = F X ( x ) ∀ x \lim_{n \xrightarrow{} \infty} F_{X_n}(x) = F_X(x)\ \ \forall x n ∞ lim F X n ( x ) = F X ( x ) ∀ x
Entropy
Entropy
H ( X ) = − ∑ i = 1 n p i ln ( p i ) H(X) = - \sum_{i=1}^{n} p_i \ln(p_i) H ( X ) = − i = 1 ∑ n p i ln ( p i )
Chain Rule of Entropy
H ( X , Y ) = H ( Y ) + H ( X ∣ Y ) = H ( X ) + H ( Y ∣ X ) H(X, Y) = H(Y) + H(X|Y) = H(X) + H(Y|X) H ( X , Y ) = H ( Y ) + H ( X ∣ Y ) = H ( X ) + H ( Y ∣ X )
Convergence of Joint Entropy
− 1 n log p ( x 1 , x 2 . . . x n ) → p H ( X ) -\frac{1}{n} \log p(x_1, x_2... x_n) \xrightarrow{p} H(X) − n 1 log p ( x 1 , x 2 ... x n ) p H ( X )
Source Coding Theorem
As n → ∞ n \xrightarrow{} \infty n ∞ , consider N iid RVs with entropy H ( X ) H(X) H ( X ) . You can compress this into no more and no less than N H ( X ) NH(X) N H ( X ) bits without sending over extra bits or losing information.
Channel Coding Theorem
Define channel capacity as the C = # of message input bits / # of bits transmitted. Any sequence of codes with error probability p → 0 p \xrightarrow{} 0 p 0 has a rate R < C R < C R < C .
Capacity of a BEC
C = 1 − p C = 1 - p C = 1 − p
Capacity of a BSC
C = 1 − H ( p ) C = 1 - H(p) C = 1 − H ( p )
Average Number of Bits Transmitted
E [ n u m b e r b i t s ] ≤ n ( H ( X ) + ϵ ) \mathbb{E}[number\ bits] \leq n(H(X) + \epsilon) E [ n u mb er bi t s ] ≤ n ( H ( X ) + ϵ )
Mutual Information
I ( X ; Y ) = ∑ p X Y ( x , y ) log P X Y ( x , y ) P X ( x ) P y ( Y ) I(X;Y) = \sum p_{XY}(x, y)\log\frac{P_{XY}(x, y)}{P_X(x)P_y(Y)} I ( X ; Y ) = ∑ p X Y ( x , y ) log P X ( x ) P y ( Y ) P X Y ( x , y )
Mutual Information and Entropy
I ( X ; Y ) = H ( X ) + H ( Y ) − H ( X , Y ) I(X;Y) = H(X) + H(Y) - H(X, Y) I ( X ; Y ) = H ( X ) + H ( Y ) − H ( X , Y )
Capacity of A Channel
C = max p x I ( X ; Y ) C = \max_{p_x} I(X;Y) C = p x max I ( X ; Y )
Upper Bound on Probability of Error in BEC
P ( error ) = 2 − n ( 1 − p ) + L ( n ) P(\text{error}) = 2^{-n(1-p) + L(n)} P ( error ) = 2 − n ( 1 − p ) + L ( n ) where n = # bits of bits sent and L = # of bits in message
Discrete Time Markov Chains
Markov Property
P ( X n + 1 ∣ X n . . . X 1 ) = P ( X n + 1 ∣ X n ) P(X_{n+1} | X_n...X_1) = P(X_{n+1} | X_n) P ( X n + 1 ∣ X n ... X 1 ) = P ( X n + 1 ∣ X n )
Chapman Komogorov Equations
P i j n = [ P n ] i j P_{ij}^n = [P^n]_{ij} P ij n = [ P n ] ij
Periodicity
d ( i ) = gcd { n ≥ 1 : P i i n > 0 } d(i) = \gcd\{n \geq 1: P_{ii}^n > 0\} d ( i ) = g cd{ n ≥ 1 : P ii n > 0 }
Stationary Distribution
π P = π \pi P = \pi π P = π
Hitting Time
β ( i ) = { 1 + ∑ j p i j β j i ∉ A 0 i ∈ A \begin{equation*}
\beta(i) = \begin{cases}
1 + \sum_j p_{ij}\beta_j & i \notin A \\
0 & i \in A
\end{cases}
\end{equation*} β ( i ) = { 1 + ∑ j p ij β j 0 i ∈ / A i ∈ A
Detailed Balance Equations
π i P j i = π i P i j , i , j ∈ S \pi_i P_{ji} = \pi_i P_{ij},\ \ i, j \in S π i P j i = π i P ij , i , j ∈ S
Stationary Distribution of an Undirected Graph
π ( i ) = d ( i ) ∑ j d ( j ) = d e g r e e ( i ) 2 E \pi(i) = \frac{d(i)}{\sum_{j}d(j)} = \frac{degree(i)}{2E} π ( i ) = ∑ j d ( j ) d ( i ) = 2 E d e g r ee ( i )
Poisson Processes
Number of arrivals within t t t
P ( N t = n ) ∼ P o i s s o n ( λ t ) = e − λ t ( λ t ) n n ! P(N_t = n)\sim Poisson(\lambda t) = \frac{e^{-\lambda t}(\lambda t)^n}{n!} P ( N t = n ) ∼ P o i sso n ( λ t ) = n ! e − λ t ( λ t ) n
Inter-arrival Time
S i ∼ E x p ( λ ) S_i \sim Exp(\lambda) S i ∼ E x p ( λ )
Sum of Inter-arrival Times: Erlang Distribution
f T n ( s ) = λ e − λ s ( λ s ) n − 1 ( n − 1 ) ! f_{T_n}(s) = \frac{\lambda e^{-\lambda s}(\lambda s)^{n-1}}{(n-1)!} f T n ( s ) = ( n − 1 )! λ e − λ s ( λ s ) n − 1
Memoryless Property
N T i − N T i − 1 ∼ P o i s s o n ( λ ( t i − t i − 1 ) ) N_{T_i} - N_{T_{i-1}} \sim Poisson(\lambda (t_i - t_{i-1})) N T i − N T i − 1 ∼ P o i sso n ( λ ( t i − t i − 1 ))
Poisson Merging
P P ( λ 1 ) + P P ( λ 2 ) ∼ P P ( λ 1 + λ 2 ) PP(\lambda_1) + PP(\lambda_2) \sim PP(\lambda_1 +
\lambda_2) P P ( λ 1 ) + P P ( λ 2 ) ∼ P P ( λ 1 + λ 2 )
Poisson Splitting
P ( min { T a , T b } = T a ) = λ a λ a + λ b P(\min\{T_a, T_b\} = T_a) = \frac{\lambda_a}{\lambda_a + \lambda_b} P ( min { T a , T b } = T a ) = λ a + λ b λ a
Random Incidence Paradox
L ∼ E r l a n g ( 2 , k ) L \sim Erlang(2, k) L ∼ E r l an g ( 2 , k )
Continuous Time Markov Chains
Temporal Homogeneity
P ( X t + τ ∣ X t = i , X s = i s ∀ 0 ≤ s < t ) = P ( X τ = j ∣ X 0 = i ) P(X_{t+\tau}\ |\ X_t = i, X_s = i_s \forall\ 0 \leq s < t) = P(X_\tau = j | X_0 = i) P ( X t + τ ∣ X t = i , X s = i s ∀ 0 ≤ s < t ) = P ( X τ = j ∣ X 0 = i )
Rate of Self-Transition
Q ( i , i ) = − ∑ j ≠ i Q ( i , j ) Q(i, i) = -\sum_{j\neq i} Q(i, j) Q ( i , i ) = − j = i ∑ Q ( i , j )
Balance Equations
∑ i ≠ j π i Q ( i , j ) = π j ∑ k ≠ j Q ( j , k ) \sum_{i\neq j} \pi_i Q(i, j) = \pi_j \sum_{k \neq j} Q(j, k) i = j ∑ π i Q ( i , j ) = π j k = j ∑ Q ( j , k )
Uniformization (Simulated DTMC)
Let q = sup q ( i ) , strongest self-loop R = I + 1 q Q
\begin{aligned}
\text{Let } q &= \text{sup}\ q(i) \text{ , strongest self-loop} \newline
R &= I + \frac{1}{q}Q
\end{aligned}
Let q R = sup q ( i ) , strongest self-loop = I + q 1 Q
Hitting Time
β ( i ) = { 1 q ( i ) + ∑ j ≠ i Q ( i , j ) q ( i ) β ( j ) i ∉ A 0 i ∈ A
\begin{equation*}
\beta(i) = \begin{cases}
\frac{1}{q(i)} + \sum_{j \neq i} \frac{Q(i, j)}{q(i)} \beta(j) & i \notin A \newline
0 & i \in A
\end{cases}
\end{equation*}
β ( i ) = { q ( i ) 1 + ∑ j = i q ( i ) Q ( i , j ) β ( j ) 0 i ∈ / A i ∈ A
Random Graph
Probability of a Random Graph Being Given Graph
P ( G = G 0 ) ∼ B i n o m i a l ( ( n 2 ) , p ) P(G = G_0) \sim Binomial({n\choose 2}, p) P ( G = G 0 ) ∼ B in o mia l ( ( 2 n ) , p )
Distribution of Degree of Vertex in Random Graph
P ( D = d ) ∼ B i n o m i a l ( n − 1 , p ) → n → ∞ P o i s s o n ( ( n − 1 ) p ) P(D = d) \sim Binomial(n-1, p) \xrightarrow{n \xrightarrow{} \infty} Poisson((n-1)p) P ( D = d ) ∼ B in o mia l ( n − 1 , p ) n ∞ P o i sso n (( n − 1 ) p )
Erdos Renyi Theorem
L e t p ( n ) = λ ln ( n ) n P ( G is connected ) → n → ∞ 0 , λ < 1 P ( G is connected ) → n → ∞ 1 , λ > 1
\begin{aligned}
Let\ p(n) = \lambda\frac{\ln(n)}{n} \newline
P(G \text{ is connected}) \xrightarrow{n \xrightarrow{} \infty} 0 &\ ,\ \lambda < 1 \newline
P(G \text{ is connected})\xrightarrow{n \xrightarrow{} \infty} 1 &\ ,\ \lambda > 1
\end{aligned}
L e t p ( n ) = λ n ln ( n ) P ( G is connected ) n ∞ 0 P ( G is connected ) n ∞ 1 , λ < 1 , λ > 1
Combining Graphs
P ( e ∈ G = G 1 ∪ G 2 ∣ e ∈ G 1 ∪ e ∈ G 2 ) = p 1 + p 2 − p 1 p 2 P(e \in G = G_1 \cup G_2 | e \in G_1 \cup e \in G2) = p_1 + p_2 - p_1p_2 P ( e ∈ G = G 1 ∪ G 2 ∣ e ∈ G 1 ∪ e ∈ G 2 ) = p 1 + p 2 − p 1 p 2
Statistical Inference
Bayes Rule Redux
P ( X = x ∣ Y = y ) = P Y ∣ X ( y ∣ x ) π ( x ) ∑ i P Y ∣ X ( y ∣ i ) π ( i ) P(X = x | Y = y) = \frac{P_{Y|X}(y|x) \pi(x)}{\sum_{i} P_{Y|X}(y|i) \pi(i)} P ( X = x ∣ Y = y ) = ∑ i P Y ∣ X ( y ∣ i ) π ( i ) P Y ∣ X ( y ∣ x ) π ( x )
Maximum A-Posteriori Estimation (MAP)
MAP ( X ∣ Y = y ) = max x P X ∣ Y ( x ∣ y ) = max x P Y ∣ X ( y ∣ x ) π ( x ) \text{MAP}(X|Y=y) = \max_{x} P_{X|Y}(x|y) = \max_{x} P_{Y|X}(y|x)\pi(x) MAP ( X ∣ Y = y ) = x max P X ∣ Y ( x ∣ y ) = x max P Y ∣ X ( y ∣ x ) π ( x )
Maximum Likelihood Estimation (MLE)
M L E ( X ∣ Y = y ) = max x P Y ∣ X ( y ∣ x ) MLE(X|Y=y) = \max_x P_{Y|X}(y|x) M L E ( X ∣ Y = y ) = x max P Y ∣ X ( y ∣ x )
Likelihood Ratio
L ( y ) = P Y ∣ X ( y ∣ 1 ) P Y ∣ X ( y ∣ 0 ) L(y) = \frac{P_{Y|X}(y|1)}{P_{Y|X}(y|0)} L ( y ) = P Y ∣ X ( y ∣0 ) P Y ∣ X ( y ∣1 )
MLE of a BSC
M L E ( X ∣ Y = y ) = { y if p ≤ 1 / 2 1 − y if p > 1 / 2
\begin{equation*}
MLE(X|Y=y) = \begin{cases}
y & \text{if }p \leq 1/2 \newline
1-y & \text{if }p > 1/2
\end{cases}
\end{equation*}
M L E ( X ∣ Y = y ) = { y 1 − y if p ≤ 1/2 if p > 1/2
M L E ( X ∣ Y = y ) = { 1 if L ( y ) ≥ 1 0 if L ( y ) < 1
\begin{equation*}
MLE(X|Y=y) = \begin{cases}
1 & \text{if } L(y) \geq 1 \newline
0 & \text{if } L(y) < 1
\end{cases}
\end{equation*}
M L E ( X ∣ Y = y ) = { 1 0 if L ( y ) ≥ 1 if L ( y ) < 1
MAP of a BSC
MAP ( X ∣ Y = y ) = { 0 if L ( y ) < π 0 π 1 1 if L ( y ) ≥ π 0 π 1
\begin{equation*}
\text{MAP}(X | Y = y) =
\begin{cases}
0 &\ \text{if } L(y) < \frac{\pi_0}{\pi_1} \newline
1 &\ \text{if } L(y) \geq \frac{\pi_0}{\pi_1}
\end{cases}
\end{equation*}
MAP ( X ∣ Y = y ) = { 0 1 if L ( y ) < π 1 π 0 if L ( y ) ≥ π 1 π 0
Likelihood Ratio for X ∈ { 0 , 1 } X\in \{0, 1\} X ∈ { 0 , 1 } with Gaussian Noise
L ( y ) = exp [ y σ 2 − 1 2 σ 2 ] L(y) = \exp{[\frac{y}{\sigma^2} - \frac{1}{2\sigma^2}]} L ( y ) = exp [ σ 2 y − 2 σ 2 1 ]
MAP for X ∈ { 0 , 1 } X\in \{0, 1\} X ∈ { 0 , 1 } with Gaussian Noise
MAP ( X ∣ Y = y ) = { 0 if L ( y ) < π 0 π 1 = y ≥ 1 2 + σ 2 l o g ( π 0 π 1 ) 1 if L ( y ) ≥ π 0 π 1
\text{MAP}(X | Y = y) =
\begin{cases}
0 &\ \text{if } L(y) < \frac{\pi_0}{\pi_1} = y \geq \frac{1}{2} + \sigma^2 log(\frac{\pi_0}{\pi_1}) \newline
1 &\ \text{if } L(y) \geq \frac{\pi_0}{\pi_1}
\end{cases}
MAP ( X ∣ Y = y ) = { 0 1 if L ( y ) < π 1 π 0 = y ≥ 2 1 + σ 2 l o g ( π 1 π 0 ) if L ( y ) ≥ π 1 π 0
MLE for X ∈ { 0 , 1 } X\in \{0, 1\} X ∈ { 0 , 1 } with Gaussian Noise
M L E ( X ∣ Y = y ) = { 1 if L ( y ) ≥ 1 = y ≥ 1 2 0 if L ( y ) < 1
\begin{equation*}
MLE(X|Y=y) = \begin{cases}
1 & \text{if } L(y) \geq 1 = y \geq \frac{1}{2} \newline
0 & \text{if } L(y) < 1 \newline
\end{cases}
\end{equation*}
M L E ( X ∣ Y = y ) = { 1 0 if L ( y ) ≥ 1 = y ≥ 2 1 if L ( y ) < 1
Binary Error Testing
Neyman-Pearson Lemma
Minimizes P(false negatives) with P(false positive) ≤ β X ^ = { 1 L ( y ) > λ 0 L ( y ) < λ B e r n ( γ ) L ( y ) = λ Setting P ( X ^ = 1 ∣ X = 0 ) = β
\begin{aligned}
&\text{Minimizes P(false negatives) with P(false positive) $\leq \beta$} \newline
&\hat{X} = \begin{cases}
1\ & L(y) > \lambda \newline
0\ & L(y) < \lambda \newline
Bern(\gamma)\ & L(y) = \lambda
\end{cases} \newline
&\text{Setting } P(\hat{X} = 1 | X = 0) = \beta
\end{aligned}
Minimizes P(false negatives) with P(false positive) ≤ β X ^ = ⎩ ⎨ ⎧ 1 0 B er n ( γ ) L ( y ) > λ L ( y ) < λ L ( y ) = λ Setting P ( X ^ = 1∣ X = 0 ) = β
Estimations
Mean Square Error (MSE)
E [ ( X − X ^ ( Y ) ) 2 ] \mathbb{E}[(X - \hat{X}(Y))^2] E [( X − X ^ ( Y ) ) 2 ]
Minimum Mean-Squared Estimation (MMSE)
MMSE ( X ∣ Y ) = argmin X ^ E [ ( X − X ^ ( Y ) ) 2 ] = E ( X ∣ Y ) \text{MMSE}(X|Y) = \text{argmin}_{\hat{X}}\mathbb{E}[(X - \hat{X}(Y))^2] = \mathbb{E}(X|Y) MMSE ( X ∣ Y ) = argmin X ^ E [( X − X ^ ( Y ) ) 2 ] = E ( X ∣ Y )
MMSE Theorem
E [ ( X − g ( Y ) ) f ( Y ) ] = 0 ∀ f ⟹ g ( Y ) = MMSE E[(X-g(Y))f(Y)] = 0 \ \forall f \implies g(Y) = \text{ MMSE} E [( X − g ( Y )) f ( Y )] = 0 ∀ f ⟹ g ( Y ) = MMSE
Linear Least Squares Estimation
L [ X ∣ Y ] = min linear X ^ E [ ∣ X − X ^ ( Y ) ∣ 2 ] = min a , b 1 . . . b n E [ ∣ X − ( a + ∑ b i Y i ) ∣ 2 ] Let Y be a vector of all observations Y i Define Σ X Y = E [ ( X − μ x ) ( Y − μ y ) T ] Σ Y = E [ ( Y − μ y ) ( Y − μ y ) T ] L [ X ∣ Y ] = μ x + Σ X Y Σ Y − 1 ( Y − μ y ) L [ X ∣ Y ] = μ x + c o v ( X , Y ) v a r ( Y ) ( Y − μ y )
\begin{aligned}
&\mathbb{L}[X|Y] = \min_{\text{linear} \hat{X}} \mathbb{E}[|X - \hat{X}(Y)|^2] = \min_{a, b_1...b_n} \mathbb{E}[|X - (a+\sum b_iY_i)|^2] \newline
&\text{Let Y be a vector of all observations $Y_i$} \newline
&\text{Define } \Sigma_{XY} = \mathbb{E}[(X-\mu_x)(Y-\mu_y)^T] \newline
&\Sigma_{Y} = \mathbb{E}[(Y-\mu_y)(Y-\mu_y)^T] \newline
&\mathbb{L}[X|Y] = \mu_x + \Sigma_{XY} \Sigma_{Y}\ ^{-1} (Y - \mu_y) \newline
&\mathbb{L}[X|Y] = \mu_x + \frac{cov(X, Y)}{var(Y)} (Y - \mu_y)
\end{aligned}
L [ X ∣ Y ] = linear X ^ min E [ ∣ X − X ^ ( Y ) ∣ 2 ] = a , b 1 ... b n min E [ ∣ X − ( a + ∑ b i Y i ) ∣ 2 ] Let Y be a vector of all observations Y i Define Σ X Y = E [( X − μ x ) ( Y − μ y ) T ] Σ Y = E [( Y − μ y ) ( Y − μ y ) T ] L [ X ∣ Y ] = μ x + Σ X Y Σ Y − 1 ( Y − μ y ) L [ X ∣ Y ] = μ x + v a r ( Y ) co v ( X , Y ) ( Y − μ y )
Linear Least Squared Error
L L S E = v a r ( X ) − Σ X Y Σ Y − 1 Σ Y X LLSE = var(X) - \Sigma_{XY}\Sigma_{Y}\ ^{-1}\Sigma_{YX} LL S E = v a r ( X ) − Σ X Y Σ Y − 1 Σ Y X
Hilbert Spaces
Hilbert Projection Theorem
∀ v ∈ H , U ⊆ H , ∃ min u ∈ U ∣ ∣ u − v ∣ ∣ : u is unique < u − v , u ′ > = 0 ∀ u ′ ∈ U
\forall v \in H, U \subseteq H,\ \exists\ \text{$\min_{u \in U}$} ||u-v||:\text{ u is unique} \\
<u-v, u'> = 0 \ \forall\ u' \in U
∀ v ∈ H , U ⊆ H , ∃ min u ∈ U ∣∣ u − v ∣∣ : u is unique < u − v , u ′ >= 0 ∀ u ′ ∈ U
Hilbert Random Variable Theorem
< X , Y > = E [ X Y ] <X, Y> = \mathbb{E}[XY] < X , Y >= E [ X Y ]
LLSE in Hilbert Spaces
< L [ X ∣ Y ] − X , u > = E [ ( L [ X ∣ Y ] − X ) u ] = 0 ∀ u
<\mathbb{L}[X|Y] - X, u> = \mathbb{E}[(\mathbb{L}[X|Y] - X)u] = 0\ \forall\ u
< L [ X ∣ Y ] − X , u >= E [( L [ X ∣ Y ] − X ) u ] = 0 ∀ u
Orthogonality Principle
E ( L [ X ∣ Y ] ) = E [ X ] E [ ( L [ X ∣ Y ] − X ) Y i ] = 0 E [ ( L [ X ∣ Y ] ⋅ Y T ] = E [ X Y T ]
\begin{aligned}
&\mathbb{E}(\mathbb{L}[X|Y]) = \mathbb{E}[X] \newline
&\mathbb{E}[(\mathbb{L}[X|Y] - X)Y_i] = 0 \newline
&\mathbb{E}[(\mathbb{L}[X|Y] \cdot Y^T] = \mathbb{E}[XY^T]
\end{aligned}
E ( L [ X ∣ Y ]) = E [ X ] E [( L [ X ∣ Y ] − X ) Y i ] = 0 E [( L [ X ∣ Y ] ⋅ Y T ] = E [ X Y T ]
Magnitude
∣ ∣ X ∣ ∣ = < X , X > = E ( ∣ X ∣ 2 ) ||X|| = \sqrt{<X, X>} = \sqrt{\mathbb{E}(|X|^2)} ∣∣ X ∣∣ = < X , X > = E ( ∣ X ∣ 2 )
Zero-Mean Multiple RVs
Let X, Y, Z zero-mean L [ X ∣ Y , Z ] = L [ X ∣ Y ] − L [ X ∣ Z − L [ Z ∣ Y ] ] L [ X ∣ Y , Z ] = L [ X ∣ Y ] − L [ X ∣ Z ] if Y, Z uncorrelated
\begin{aligned}
&\text{Let X, Y, Z zero-mean} \newline
&L[X|Y, Z] = L[X|Y] - L[X|Z - L[Z|Y]] \newline
&L[X|Y, Z] = L[X|Y] - L[X|Z] \text { if Y, Z uncorrelated}
\end{aligned}
Let X, Y, Z zero-mean L [ X ∣ Y , Z ] = L [ X ∣ Y ] − L [ X ∣ Z − L [ Z ∣ Y ]] L [ X ∣ Y , Z ] = L [ X ∣ Y ] − L [ X ∣ Z ] if Y, Z uncorrelated