Welcome to Randal Douc's wiki

A collaborative site on maths but not only!

User Tools

Site Tools


world:kullback

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
world:kullback [2023/11/10 10:10]
rdouc
world:kullback [2023/11/10 10:19] (current)
rdouc [Proof]
Line 20: Line 20:
     &= -\frac{p}{2} ​ + \frac{(\mu_0-\mu_1)^T \Sigma_1^{-1} (\mu_0-\mu_1) }{2} + \frac 1 2 \PE\lrb{U_0^T \Sigma_0^{1/​2} \Sigma_1^{-1} \Sigma_0^{1/​2} U_0} - \frac{1}{2} \log \frac{\mathrm{det} \Sigma_0}{\mathrm{det} \Sigma_1}\\ ​     &= -\frac{p}{2} ​ + \frac{(\mu_0-\mu_1)^T \Sigma_1^{-1} (\mu_0-\mu_1) }{2} + \frac 1 2 \PE\lrb{U_0^T \Sigma_0^{1/​2} \Sigma_1^{-1} \Sigma_0^{1/​2} U_0} - \frac{1}{2} \log \frac{\mathrm{det} \Sigma_0}{\mathrm{det} \Sigma_1}\\ ​
     &= -\frac{p}{2} ​ + \frac{(\mu_0-\mu_1)^T \Sigma_1^{-1} (\mu_0-\mu_1) }{2} + \frac{Tr\lr{\Sigma_1^{-1} \Sigma_0}}{2} - \frac{1}{2} \log \frac{\mathrm{det} \Sigma_0}{\mathrm{det} \Sigma_1}     &= -\frac{p}{2} ​ + \frac{(\mu_0-\mu_1)^T \Sigma_1^{-1} (\mu_0-\mu_1) }{2} + \frac{Tr\lr{\Sigma_1^{-1} \Sigma_0}}{2} - \frac{1}{2} \log \frac{\mathrm{det} \Sigma_0}{\mathrm{det} \Sigma_1}
 +\end{align*}
 +where in the last line, we have used that 
 +
 +\begin{align*}
 +    \PE\lrb{U_0^T \Sigma_0^{1/​2} \Sigma_1^{-1} \Sigma_0^{1/​2} U_0}&​=\lr{\PE\lrb{Tr\lr{U_0^T \Sigma_0^{1/​2} \Sigma_1^{-1} \Sigma_0^{1/​2} U_0}}}=Tr\lr{\PE\lrb{\Sigma_1^{-1} \Sigma_0^{1/​2} U_0 U_0^T \Sigma_0^{1/​2} }} \\ 
 +    &​=Tr\lr{\Sigma_1^{-1} \Sigma_0^{1/​2} \PE\lrb{U_0 U_0^T} \Sigma_0^{1/​2} }=Tr\lr{\Sigma_1^{-1} \Sigma_0}
 \end{align*} \end{align*}
world/kullback.1699607448.txt.gz · Last modified: 2023/11/10 10:10 by rdouc