Line 107: | Line 107: | ||
From Eq. (--------) we assign <math> x \in \omega_1 </math> | From Eq. (--------) we assign <math> x \in \omega_1 </math> | ||
if | if | ||
− | |||
<center> | <center> | ||
Line 139: | Line 138: | ||
which has the form | which has the form | ||
+ | <center> | ||
<math>a x^2 + b x + c > 0</math> | <math>a x^2 + b x + c > 0</math> | ||
+ | </center> | ||
+ | |||
+ | where | ||
+ | |||
+ | <center> | ||
+ | <math>a = \frac{1}{2} \left( \frac{1}{\sigma_2^2} - | ||
+ | \frac{1}{\sigma_1^2} \right),</math> | ||
+ | </center> | ||
+ | |||
+ | <center> | ||
+ | <math>b = \frac{\mu_1}{\sigma_1^2} - | ||
+ | \frac{\mu_2}{\sigma_2^2},</math> | ||
+ | </center> | ||
+ | |||
+ | and | ||
+ | |||
+ | <center> | ||
+ | <math>c = \frac{1}{2} \left( \frac{\mu_2^2}{\sigma_2^2} | ||
+ | - \frac{\mu_1^2}{\sigma_1^2} \right) | ||
+ | + \text{ln} \left( \frac{\sigma_2}{\sigma_1} \right) | ||
+ | + \text{ln} \left( \frac{P_1 c_{21}}{P_2 c_{12}} \right).</math> | ||
+ | </center> | ||
+ | |||
+ | This form is similar to Bayes rule for minimizing error, | ||
+ | except for the factor of <math> \text{ln} \left( \frac{P_1 | ||
+ | c_{21}}{P_2 c_{12}} \right) </math>, which shifts the decision | ||
+ | thresholds. An equivalent formulation is to decide <math> | ||
+ | x \in \omega_1 </math> if | ||
+ | |||
+ | <center> | ||
+ | <math>\frac{p_1(x)}{p_2(x)} > \frac{P_2 c_{12}}{P_1 c_{21}}</math> | ||
+ | </center> | ||
+ | |||
+ | or | ||
+ | |||
+ | <center> | ||
+ | <math>\frac{p_1(x)}{p_2(x)} > \lambda</math> | ||
+ | </center> | ||
+ | |||
+ | where | ||
+ | |||
+ | <center> | ||
+ | <math>\lambda = \frac{P_2 c_{12}}{P_1 c_{21}}.</math> | ||
+ | </center> | ||
+ | We can interpret the decision rule as a modification of the | ||
+ | Neyman-Pearson criterion that takes into account the priors | ||
+ | and the cost. | ||
Revision as of 16:45, 12 April 2014
Bayes Error for Minimizing Risk
Contents
Introduction
In class we discussed Bayes rule for minimizing the probability of error. Our goal is to generalize this rule to minimize risk instead of probability of error. For simplicity we deal with the two class case. Then we provide examples for the cases of 1D and 2D features, and we derive Bayes rule for minimizing risk in these cases.
We briefly motivate this topic by considering the following scenario. Suppose that a doctor has to classify a tumor as cancerous or benign. We consider the cost of misclassifying the tumor as benign to be high, since we would like the patient to be treated quickly if cancer is present. Therefore, we would assign a high cost to misclassifying the tumor as benign. We shall show how to incorporate this cost into Bayes rule in this article.
Bayes rule for minimizing risk
Let $ x \in \mathbb R^n $ be a feature vector. Let $ q_i(x) $ be the posterior probability of class $ i $ (denoted as $ \omega_i $) given $ x $, and let $ P_i $ be the prior probability for $ omega_i $. Let $ p_i(x) $ be the class conditional density for class $ i $. Denote $ c_{ij} $ as the cost of deciding $ x \in \omega_i $ with $ \omega_j $ as the true class. The conditional cost of assigning $ x \in \omega_i $ given $ x $ is
$ r_i(x) = c_{i1} q_1(x) + c_{i2} q_2(x) $
where $ i = 1, 2 $.
We assign $ x $ so that cost is minimum:
$ r_1(x) \lessgtr^{\omega_1}_{\omega_2} r_2(x) $
which says to decide $ \omega_1 $ if $ r_1(x) < r_2(x) $ and $ \omega_2 $ otherwise. If we make the decision this way, the total cost becomes
$ \begin{align} E[r(x)] &= \int \text{min}[ r_1(x), r_2(x) ] p(x) dx \\ &= \int \text{min} [c_{11} q_1(x) + c_{12} q_2(x), c_{21} q_1(x) + c_{22} q_2(x)] p(x) dx \\ &= \int \text{min} [c_{11} P_1 p_1(x) + c_{12} P_2 p_2(x), c_{21} P_1 p_1(x) + c_{22} P_2 p_2(x) ] dx \\ &= \int_{R_1} c_{11} P_1 p_1(x) + c_{12} P_2 p_2(x) dx + \int_{R_2} c_{21} P_1 p_1(x) + c_{22} P_2 p_2(x) dx \\ &= (c_{21} P_1 + c_{22} P_2) + \int_{R_1} (c_{11} - c_{21}) P_1 p_1(x) + (c_{12} - c_{22}) P_2 p_2(x) dx \end{align} $
where $ R_1 $ and $ R_2 $ are partitions of $ \mathbb R^n $, $ R_1 $ and $ R_2 $ are determined by the decision rule from Eq. (-----------), and we use $ \int_{R_2} p_1(x) dx = 1 - \int_{R_1} p_1(x) dx $ in the last line. To minimize Eq. (---------), we make the integrand as negative as possible:
$ (c_{12} - c_{22}) P_2 p_2(x) \lessgtr^{\omega_1}_{\omega_2} (c_{21} - c_{11}) P_1 p_1(x) $
which is equivalent to assigning $ x \in \omega_1 $ if $ x $ makes the integrand negative, and $ x \in \omega_2 $ otherwise. Bayes test for minimum cost can now be stated as
$ \frac{p_1(x)}{p_2(x)} \gtrless^{\omega_1}_{\omega_2} \frac{(c_{12} - c_{22})P_2}{(c_{21} - c_{11})P_1}. $
Example 1: 1D features
Consider two classes and $ x \in \mathbb R $. Let $ p_1(x) = N(\mu_1, \sigma_1) $ and $ p_2(x) = N(\mu_2, \sigma_2) $. For simplicity, let $ c_{11} = c_{22} = 0 $. From Eq. (--------) we assign $ x \in \omega_1 $ if
$ \begin{align} c_{21} P_1 p_1(x) &> c_{12} P_2 p_2(x) \\ \iff \text{ln}(c_{21} P_1) + \text{ln}(p_1(x)) &> \text{ln}(c_{12} P_2) + \text{ln}(p_2(x)) \\ \iff \text{ln}(c_{21} P_1) - \text{ln}(\sigma_1) - \frac{1}{2}\text{ln}(2 \pi) - \frac{(x - \mu_1)^2}{2 \sigma_1^2} &> \text{ln}(c_{12} P_2) - \text{ln}(\sigma_2) - \frac{1}{2}\text{ln}(2 \pi) - \frac{(x - \mu_2)^2}{2 \sigma_2^2} \end{align} $
so the discriminant function becomes
$ \frac{1}{2} \left( \frac{1}{\sigma_2^2} - \frac{1}{\sigma_1^2} \right) x^2 + \left( \frac{\mu_1}{\sigma_1^2} - \frac{\mu_2}{\sigma_2^2} \right) x + \frac{1}{2} \left( \frac{\mu_2^2}{\sigma_2^2} - \frac{\mu_1^2}{\sigma_1^2} \right) + \text{ln} \left( \frac{\sigma_2}{\sigma_1} \right) + \text{ln} \left( \frac{P_1 c_{21}}{P_2 c_{12}} \right) > 0 $
which has the form
$ a x^2 + b x + c > 0 $
where
$ a = \frac{1}{2} \left( \frac{1}{\sigma_2^2} - \frac{1}{\sigma_1^2} \right), $
$ b = \frac{\mu_1}{\sigma_1^2} - \frac{\mu_2}{\sigma_2^2}, $
and
$ c = \frac{1}{2} \left( \frac{\mu_2^2}{\sigma_2^2} - \frac{\mu_1^2}{\sigma_1^2} \right) + \text{ln} \left( \frac{\sigma_2}{\sigma_1} \right) + \text{ln} \left( \frac{P_1 c_{21}}{P_2 c_{12}} \right). $
This form is similar to Bayes rule for minimizing error, except for the factor of $ \text{ln} \left( \frac{P_1 c_{21}}{P_2 c_{12}} \right) $, which shifts the decision thresholds. An equivalent formulation is to decide $ x \in \omega_1 $ if
$ \frac{p_1(x)}{p_2(x)} > \frac{P_2 c_{12}}{P_1 c_{21}} $
or
$ \frac{p_1(x)}{p_2(x)} > \lambda $
where
$ \lambda = \frac{P_2 c_{12}}{P_1 c_{21}}. $
We can interpret the decision rule as a modification of the Neyman-Pearson criterion that takes into account the priors and the cost.
Example 2: 2D features
![](/rhea/images/thumb/a/a7/C12_c21.png/600px-C12_c21.png)
Summary and Conclusions
In this lecture we have shown that the probability of error ($Prob \left[ Error \right] $) when using Bayes error, is upper bounded by the Chernoff Bound. Therefore,
for $ \beta \in \left[ 0, 1 \right] $.
When $ \beta =\frac{1}{2} $ then $ \varepsilon_{\frac{1}{2}} $ in known as the Bhattacharyya bound.
References
[1]. Duda, Richard O. and Hart, Peter E. and Stork, David G., "Pattern Classication (2nd Edition)," Wiley-Interscience, 2000.
[2]. Mireille Boutin, "ECE662: Statistical Pattern Recognition and Decision Making Processes," Purdue University, Spring 2014.
Questions and comments
If you have any questions, comments, etc. please post them On this page.