Statistical Inference

Suppose we have a variable XX (may or may not be a random variable) that represents the state of nature. We observe a variable YY which is obtained by some model of the world PYXP_{Y|X}.

Suppose we know that XπX\sim \pi where π\pi is a probability distribution. If we observe Y=yY=y, then the a posteriori estimate of XX is given by Bayes Rule

Pr{X=xY=y}=PYX(yx)π(x)x~PYX(yx~)π(x~)PYX(yx)π(x).\text{Pr}\left\{X=x | Y=y\right\} = \frac{P_{Y|X}(y|x)\pi(x)}{\sum_{\tilde{x}}P_{Y|X}(y|\tilde{x})\pi(\tilde{x})} \propto P_{Y|X}(y|x)\pi(x).

Since the estimate is only dependent on the model and the prior, we don’t actually need to compute the probabilities to figure out the most likely XX.

Definition 77

The Maximum A Posteriori (MAP) estimate is given by

X^MAP(y)=argmaxxPYX(yx)π(x)\hat{X}_{MAP}(y) = \text{argmax}_x P_{Y|X}(y|x)\pi(x)

If we have no prior information on XX, then we can assume π\pi is uniform, reducing Definition 77 to only optimize over the model.

Definition 78

The Maximum Likelihood (ML) estimate is given by

X^ML(y)=argmaxxPYX(yx)\hat{X}_{ML}(y) = \text{argmax}_x P_{Y|X}(y|x)

Binary Hypothesis Testing

Definition 79

A Binary Hypothesis Test is a type of statistical inference where the unknown variable X{0,1}X\in\{ 0, 1 \}.

Since there are only two possible values of XX in a binary test, there are two “hypotheses” that we have, and we want to accept the more likely one.

Definition 80

The Null Hypothesis H0H_0 says that YPYX=0Y\sim P_{Y|X=0}

Definition 81

The Alternate Hypothesis H1H_1 says that YPYX=1Y\sim P_{Y|X=1}

With two possible hypotheses, there are two kinds of errors we can make.

Definition 82

A Type I error (false positive) is when we incorrectly reject the null hypothesis. The Type I error probability is then

Pr{X^(Y)=1X=0}\text{Pr}\left\{\hat{X}(Y) = 1 | X = 0\right\}

Definition 83

A Type II error (false negative) is when we incorrectly accept the null hypothesis. The Type II error probability is then

Pr{X^(Y)=0X=1}\text{Pr}\left\{\hat{X}(Y) = 0 | X = 1\right\}

Our goal is to create a decision rule X^:Y{0,1}\hat{X}: \mathcal{Y} \to \{0, 1\} that we can use to predict XX. Based on what the decision rule is used for, there will be requirements on how large the probability of Type I and Type II errors can be. We can formulate the search for a hypothesis test as an optimization. For some β[0,1]\beta \in [0, 1], we want to find

X^β(Y)=argminPr{X^(Y)=0X=1}:Pr{X^(Y)=1X=0}β.(1)\hat{X}_\beta(Y) = \text{argmin} \text{Pr}\left\{\hat{X}(Y)=0 | X=1\right\} \quad : \quad \text{Pr}\left\{\hat{X}(Y)=1|X=0\right\} \leq \beta. \qquad (1)

Intuitively, our test should depend on pYX(y1)p_{Y|X}(y|1) and pYX(y0)p_{Y|X}(y|0) since these quantities give us how likely we are to get our observations if we knew the ground truth. We can define a ratio that formally compares these two quantities.

Definition 84

The likelihood ratio is given by

L(y)=pYX(y1)pYX(y0)L(y) = \frac{p_{Y|X}(y|1)}{p_{Y|X}(y|0)}

Notice that we can write MLE as a threshold on the likelihood ratio since if L(y)1L(y) \geq 1, then we say X=1X=1, and vice versa. However, there is no particular reason that 11 must always be the number at which we threshold our likelihood ratio, and so we can generalize this idea to form different forms of tests.

Definition 85

For some threshold cc and randomization probability γ\gamma, a threshold test is of the form

X^(y)={1 if L(y)>c0 if L(y)<c Bernoulli(γ) if L(y)=c.\hat{X}(y) = \begin{cases} 1 & \text{ if } L(y) > c\\ 0 & \text{ if } L(y) < c\\ \text{ Bernoulli}(\gamma) & \text { if } L(y) = c. \end{cases}

MAP fits into the framework of a threshold test since we can write

X^MAP={1 if L(y)π0π10 if L(y)<π0π1\hat{X}_{MAP} = \begin{cases} 1 & \text{ if } L(y) \geq \frac{\pi_0}{\pi_1}\\ 0 & \text{ if } L(y) < \frac{\pi_0}{\pi_1} \end{cases}

It turns out that threshold tests are optimal with respect to solving Equation 1.

Theorem 44 (Neyman Pearson Lemma)

Given β[0,1]\beta\in[0, 1], the optimal decision rule to

X^β(Y)=argminPr{X^(Y)=0X=1}:Pr{X^(Y)=1X=0}β\hat{X}_\beta(Y) = \text{argmin} \text{Pr}\left\{\hat{X}(Y)=0 | X=1\right\} \quad : \quad \text{Pr}\left\{\hat{X}(Y)=1|X=0\right\} \leq \beta

is a threshold test.

When L(y)L(y) is monotonically increasing or decreasing, we can make the decision rule even simpler since it can be turned into a threshold on yy. For example, if L(y)L(y) is monotonically inreasing, then an optimal decision rule might be

X^(y)={1 if y>c0 if y<cBernoulli(γ) if y=c.\hat{X}(y) = \begin{cases} 1 & \text{ if } y > c\\ 0 & \text{ if } y < c\\ \text{Bernoulli}(\gamma) & \text{ if } y = c. \end{cases}

Last updated