PumpkinBook/docs/chapter9/chapter9.md

9.5

$$ JC=\frac{a}{a+b+c} $$

[解析]：给定两个集合$A$和$B$，则Jaccard系数定义为如下公式

$$ JC=\frac{|A\bigcap B|}{|A\bigcup B|}=\frac{|A\bigcap B|}{|A|+|B|-|A\bigcap B|} $$ Jaccard系数可以用来描述两个集合的相似程度。

推论：假设全集$U$共有$n$个元素，且$A\subseteq U$，$B\subseteq U$，则每一个元素的位置共有四种情况：

1、元素同时在集合$A$和$B$中，这样的元素个数记为$M_{11}$；

2、元素出现在集合$A$中，但没有出现在集合$B$中，这样的元素个数记为$M_{10}$；

3、元素没有出现在集合$A$中，但出现在集合$B$中，这样的元素个数记为$M_{01}$；

4、元素既没有出现在集合$A$中，也没有出现在集合$B$中，这样的元素个数记为$M_{00}$。

根据Jaccard系数的定义，此时的Jaccard系数为如下公式 $$ JC=\frac{M{11}}{M{11}+M{10}+M{01}} $$ 由于聚类属于无监督学习，事先并不知道聚类后样本所属类别的类别标记所代表的意义，即便参考模型的类别标记意义是已知的，我们也无法知道聚类后的类别标记与参考模型的类别标记是如何对应的，况且聚类后的类别总数与参考模型的类别总数还可能不一样，因此只用单个样本无法衡量聚类性能的好坏。

由于外部指标的基本思想就是以参考模型的类别划分为参照，因此如果某一个样本对中的两个样本在聚类结果中同属于一个类，在参考模型中也同属于一个类，或者这两个样本在聚类结果中不同属于一个类，在参考模型中也不同属于一个类，那么对于这两个样本来说这是一个好的聚类结果。

总的来说所有样本对中的两个样本共存在四种情况：
1、样本对中的两个样本在聚类结果中属于同一个类，在参考模型中也属于同一个类；
2、样本对中的两个样本在聚类结果中属于同一个类，在参考模型中不属于同一个类；
3、样本对中的两个样本在聚类结果中不属于同一个类，在参考模型中属于同一个类；
4、样本对中的两个样本在聚类结果中不属于同一个类，在参考模型中也不属于同一个类。

综上所述，即所有样本对存在着书中公式(9.1)-(9.4)的四种情况，现在假设集合$A$中存放着两个样本都同属于聚类结果的同一个类的样本对，即$A=SS\bigcup SD$，集合$B$中存放着两个样本都同属于参考模型的同一个类的样本对，即$B=SS\bigcup DS$，那么根据Jaccard系数的定义有： $$ JC=\frac{|A\bigcap B|}{|A\bigcup B|}=\frac{|SS|}{|SS\bigcup SD\bigcup DS|}=\frac{a}{a+b+c} $$ 也可直接将书中公式(9.1)-(9.4)的四种情况类比推论，即$M{11}=a$，$M{10}=b$，$M{01}=c$，所以 $$ JC=\frac{M{11}}{M{11}+M{10}+M_{01}}=\frac{a}{a+b+c} $$

9.6

$$ FMI=\sqrt{\frac{a}{a+b}\cdot \frac{a}{a+c}} $$

[解析]：其中$\frac{a}{a+b}$和$\frac{a}{a+c}$为Wallace提出的两个非对称指标，$a$代表两个样本在聚类结果和参考模型中均属于同一类的样本对的个数，$a+b$代表两个样本在聚类结果中属于同一类的样本对的个数，$a+c$代表两个样本在参考模型中属于同一类的样本对的个数，这两个非对称指标均可理解为样本对中的两个样本在聚类结果和参考模型中均属于同一类的概率。由于指标的非对称性，这两个概率值往往不一样，因此Fowlkes和Mallows提出利用几何平均数将这两个非对称指标转化为一个对称指标，即Fowlkes and Mallows Index, FMI。

9.7

$$ RI=\frac{2(a+d)}{m(m-1)} $$ [解析]：Rand Index定义如下： $$ RI=\frac{a+d}{a+b+c+d}=\frac{a+d}{m(m-1)/2}=\frac{2(a+d)}{m(m-1)} $$ 即可以理解为两个样本都属于聚类结果和参考模型中的同一类的样本对的个数与两个样本都分别不属于聚类结果和参考模型中的同一类的样本对的个数的总和在所有样本对中出现的频率，可以简单理解为聚类结果与参考模型的一致性。

参看 https://en.wikipedia.org/wiki/Rand_index

9.33

$$ \sum{j=1}^m \frac{\alpha{i}\cdot p\left(\boldsymbol{x_{j}}|\boldsymbol\mu {i},\boldsymbol\Sigma{i}\right)}{\sum{l=1}^k \alpha{l}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{l},\boldsymbol\Sigma{l})}(\boldsymbol{x{j}-\mu_{i}})=0 $$

[推导]：根据公式(9.28)可知： $$ p(\boldsymbol{x{j}|\boldsymbol\mu{i},\boldsymbol\Sigma{i}})=\frac{1}{(2\pi)^\frac{n}{2}\left| \boldsymbol\Sigma{i}\right |^\frac{1}{2}}e^{-\frac{1}{2}(\boldsymbol{x{j}}-\boldsymbol\mu{i})^T\boldsymbol\Sigma{i}^{-1}(\boldsymbol{x{j}-\mu_{i}})} $$

又根据公式(9.32)，由 $$ \frac {\partial LL(D)}{\partial \boldsymbol\mu{i}}=0 $$ 可得 $$\begin{aligned} \frac {\partial LL(D)}{\partial\boldsymbol\mu{i}}&=\frac {\partial}{\partial \boldsymbol\mu{i}}\sum{j=1}^mln\Bigg(\sum{i=1}^k \alpha{i}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{i},\boldsymbol\Sigma{i})\Bigg) \ &=\sum{j=1}^m\frac{\partial}{\partial\boldsymbol\mu{i}}ln\Bigg(\sum{i=1}^k \alpha{i}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{i},\boldsymbol\Sigma{i})\Bigg) \ &=\sum{j=1}^m\frac{\alpha{i}\cdot \frac{\partial}{\partial\boldsymbol{\mu{i}}}(p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}}))}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})} \ &=\sum{j=1}^m\frac{\alpha{i}\cdot \frac{1}{(2\pi)^\frac{n}{2}\left| \boldsymbol\Sigma{i}\right |^\frac{1}{2}}e^{-\frac{1}{2}(\boldsymbol{x{j}}-\boldsymbol\mu{i})^T\boldsymbol\Sigma{i}^{-1}(\boldsymbol{x{j}-\mu{i}})}}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}\frac{\partial}{\partial \boldsymbol\mu{i}}\left(-\frac{1}{2}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \ &=\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}(-\frac{1}{2})\left(\left(\boldsymbol\Sigma{i}^{-1}+\left(\boldsymbol\Sigma{i}^{-1}\right)^T\right)\cdot\left(\boldsymbol{x{j}-\mu{i}}\right)\cdot(-1)\right) \ &=\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}(-\frac{1}{2})\left(-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)-\left(\boldsymbol\Sigma{i}^{-1}\right)^T\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \ &=\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}(-\frac{1}{2})\left(-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)-\left(\boldsymbol\Sigma{i}^T\right)^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \ &=\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}(-\frac{1}{2})\left(-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \ &=\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}(-\frac{1}{2})\left(-2\cdot\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \ &=\sum{j=1}^m \frac{\alpha{i}\cdot p\left(\boldsymbol{x{j}}|\boldsymbol\mu {i},\boldsymbol\Sigma{i}\right)}{\sum{l=1}^k \alpha{l}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{l},\boldsymbol\Sigma{l})}\boldsymbol\Sigma{i}^{-1}(\boldsymbol{x{j}-\mu{i}}) \ &=\sum{j=1}^m \frac{\alpha{i}\cdot p\left(\boldsymbol{x_{j}}|\boldsymbol\mu {i},\boldsymbol\Sigma{i}\right)}{\sum{l=1}^k \alpha{l}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{l},\boldsymbol\Sigma{l})}(\boldsymbol{x{j}-\mu_{i}})=0 \end{aligned}$$

9.35

$$ \boldsymbol\Sigma{i}=\frac{\sum{j=1}^m\gamma{ji}(\boldsymbol{x{j}-\mu{i}})(\boldsymbol{x{j}-\mu{i}})^T}{\sum{j=1}^m}\gamma_{ji} $$

[推导]：根据公式(9.28)可知： $$ p(\boldsymbol{x{j}|\boldsymbol\mu{i},\boldsymbol\Sigma{i}})=\frac{1}{(2\pi)^\frac{n}{2}\left| \boldsymbol\Sigma{i}\right |^\frac{1}{2}}e^{-\frac{1}{2}(\boldsymbol{x{j}}-\boldsymbol\mu{i})^T\boldsymbol\Sigma{i}^{-1}(\boldsymbol{x{j}-\mu{i}})} $$ 又根据公式(9.32)，由 $$ \frac {\partial LL(D)}{\partial \boldsymbol\Sigma{i}}=0 $$ 可得 $$\begin{aligned} \frac {\partial LL(D)}{\partial\boldsymbol\Sigma{i}}&=\frac {\partial}{\partial \boldsymbol\Sigma{i}}\sum{j=1}^mln\Bigg(\sum{i=1}^k \alpha{i}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{i},\boldsymbol\Sigma{i})\Bigg) \ &=\sum{j=1}^m\frac{\partial}{\partial\boldsymbol\Sigma{i}}ln\Bigg(\sum{i=1}^k \alpha{i}\cdot p(\boldsymbol{x{j}}|\boldsymbol\mu{i},\boldsymbol\Sigma{i})\Bigg) \ &=\sum{j=1}^m \frac{\alpha{i}\cdot \frac{\partial}{\partial\boldsymbol\Sigma{i}}p(\boldsymbol x{j}|\boldsymbol \mu{i},\boldsymbol\Sigma{i})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})} \ &=\sum{j=1}^m \frac{\alpha{i}\cdot \frac{\partial}{\partial\boldsymbol\Sigma{i}}\frac{1}{(2\pi)^\frac{n}{2}\left| \boldsymbol\Sigma{i}\right |^\frac{1}{2}}e^{-\frac{1}{2}(\boldsymbol{x{j}}-\boldsymbol\mu{i})^T\boldsymbol\Sigma{i}^{-1}(\boldsymbol{x{j}-\mu{i}})}}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})}\ &=\sum{j=1}^m \frac{\alpha{i}\cdot \frac{\partial}{\partial\boldsymbol\Sigma{i}}e^{ln\left(\frac{1}{(2\pi)^\frac{n}{2}\left| \boldsymbol\Sigma{i}\right |^\frac{1}{2}}e^{-\frac{1}{2}(\boldsymbol{x{j}}-\boldsymbol\mu{i})^T\boldsymbol\Sigma{i}^{-1}(\boldsymbol{x{j}-\mu{i}})}\right)}}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})} \ &=\sum{j=1}^m \frac{\alpha{i}\cdot \frac{\partial}{\partial\boldsymbol\Sigma{i}}e^{-\frac{n}{2}ln\left(2\pi\right)-\frac{1}{2}ln\left(|\boldsymbol\Sigma{i}|\right)-\frac{1}{2}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)}}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})} \ &=\sum{j=1}^m \frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol \mu{i},\boldsymbol\Sigma{i})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})}\frac{\partial}{\partial\boldsymbol\Sigma{i}}\left(-\frac{n}{2}ln\left(2\pi\right)-\frac{1}{2}ln\left(|\boldsymbol\Sigma{i}|\right)-\frac{1}{2}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \ &=\sum{j=1}^m \frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol \mu{i},\boldsymbol\Sigma{i})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})}\left(-\frac{1}{2}\left(\boldsymbol\Sigma{i}^{-1}\right)^T-\frac{1}{2}\frac{\partial}{\partial\boldsymbol\Sigma{i}}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right) \end{aligned}$$

为求得 $$ \frac{\partial}{\partial\boldsymbol\Sigma{i}}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right) $$

首先分析对$\boldsymbol \Sigma{i}$中单一元素的求导，用$r$代表矩阵$\boldsymbol\Sigma{i}$的行索引，$c$代表矩阵$\boldsymbol\Sigma{i}$的列索引，则 $$\begin{aligned} \frac{\partial}{\partial\Sigma{i{rc}}}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)&=\left(\boldsymbol{x{j}-\mu{i}}\right)^T\frac{\partial\boldsymbol\Sigma{i}^{-1}}{\partial\Sigma{i{rc}}}\left(\boldsymbol{x{j}-\mu{i}}\right) \ &=-\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\frac{\partial\boldsymbol\Sigma{i}}{\partial\Sigma{i{rc}}}\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right) \end{aligned}$$ 设$B=\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)$，则 $$\begin{aligned} B^T&=\left(\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right)^T \ &=\left(\boldsymbol{x{j}-\mu{i}}\right)^T\left(\boldsymbol\Sigma{i}^{-1}\right)^T \ &=\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1} \end{aligned}$$ 所以 $$\begin{aligned} \frac{\partial}{\partial\Sigma{i{rc}}}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)=-B^T\frac{\partial\boldsymbol\Sigma{i}}{\partial\Sigma{i{rc}}}B\end{aligned}$$ 其中$B$为$n\times1$阶矩阵，$\frac{\partial\boldsymbol\Sigma{i}}{\partial\Sigma{i{rc}}}$为$n$阶方阵，且$\frac{\partial\boldsymbol\Sigma{i}}{\partial\Sigma{i{rc}}}$仅在$\left(r,c\right)$位置处的元素值为1，其它位置处的元素值均为$0$，所以 $$\begin{aligned} \frac{\partial}{\partial\Sigma{i{rc}}}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)=-B^T\frac{\partial\boldsymbol\Sigma{i}}{\partial\Sigma{i{rc}}}B=-B{r}\cdot B{c}=-\left(B\cdot B^T\right){rc}=\left(-B\cdot B^T\right){rc}\end{aligned}$$ 即对$\boldsymbol\Sigma{i}$中特定位置的元素的求导结果对应于$\left(-B\cdot B^T\right)$中相同位置的元素值，所以 $$\begin{aligned} \frac{\partial}{\partial\boldsymbol\Sigma{i}}\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)&=-B\cdot B^T\ &=-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\left(\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\right)^T\ &=-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma_{i}^{-1} \end{aligned}$$

因此最终结果为 $$ \frac {\partial LL(D)}{\partial \boldsymbol\Sigma{i}}=\sum{j=1}^m \frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol \mu{i},\boldsymbol\Sigma{i})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j},|\boldsymbol \mu{l},\boldsymbol\Sigma{l})}\left( -\frac{1}{2}\left(\boldsymbol\Sigma{i}^{-1}-\boldsymbol\Sigma{i}^{-1}\left(\boldsymbol{x{j}-\mu{i}}\right)\left(\boldsymbol{x{j}-\mu{i}}\right)^T\boldsymbol\Sigma{i}^{-1}\right)\right)=0 $$

整理可得 $$ \boldsymbol\Sigma{i}=\frac{\sum{j=1}^m\gamma{ji}(\boldsymbol{x{j}-\mu{i}})(\boldsymbol{x{j}-\mu{i}})^T}{\sum{j=1}^m}\gamma_{ji} $$

9.38

$$ \alpha{i}=\frac{1}{m}\sum{j=1}^m\gamma_{ji} $$

[推导]：基于公式(9.37)进行恒等变形： $$ \sum{j=1}^m\frac{p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma_{l}})}+\lambda=0 $$

$$ \Rightarrow\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}+\alpha_{i}\lambda=0 $$

对所有混合成分进行求和： $$ \Rightarrow\sum{i=1}^k\left(\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}+\alpha{i}\lambda\right)=0 $$

$$ \Rightarrow\sum{i=1}^k\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}+\sum{i=1}^k\alpha_{i}\lambda=0 $$

$$ \Rightarrow\lambda=-\sum{i=1}^k\sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma_{l}})}=-m $$

又 $$ \sum{j=1}^m\frac{\alpha{i}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{i},\Sigma{i}})}{\sum{l=1}^k\alpha{l}\cdot p(\boldsymbol x{j}|\boldsymbol{\mu{l},\Sigma{l}})}+\alpha_{i}\lambda=0 $$

$$ \Rightarrow\sum{j=1}^m\gamma{ji}+\alpha_{i}\lambda=0 $$

$$ \Rightarrow\alpha{i}=-\frac{\sum{j=1}^m\gamma{ji}}{\lambda}=\frac{1}{m}\sum{j=1}^m\gamma_{ji} $$

附录

参考公式 $$ \frac{\partial\boldsymbol x^TB\boldsymbol x}{\partial\boldsymbol x}=\left(B+B^T\right)\boldsymbol x $$ $$ \frac{\partial}{\partial A}ln|A|=\left(A^{-1}\right)^T $$ $$ \frac{\partial}{\partial x}\left(A^{-1}\right)=-A^{-1}\frac{\partial A}{\partial x}A^{-1} $$ 参考资料
[1] Meilă, Marina. "Comparing clusterings—an information based distance." Journal of multivariate analysis 98.5 (2007): 873-895.
[2] Halkidi, Maria, Yannis Batistakis, and Michalis Vazirgiannis. "On clustering validation techniques." Journal of intelligent information systems 17.2-3 (2001): 107-145.
[3] Petersen, K. B. & Pedersen, M. S. The Matrix Cookbook.
[4] Bishop, C. M. (2006). Pattern Recognition and Machine Learning. Springer.