如何解决scipy.stats 连续分布的互信息
看到 scipy.stats
中的连续分布使用 scipy.integrate.quad 计算 .entropy
,它解决了一个单变量积分,下面的函数或类似的函数如何用于计算互信息使用 scipy.integrate.dblquad(双积分)的 2 个连续分布?换句话说,在一个工作示例中,它如何与 scipy.stats
中的连续分布兼容?
import numpy as np
from scipy.integrate import quad,dblquad
def mutual_information(self):
def mi_integrand(x,y):
return np.exp(self.logpdf_xy(x,y)) * \
(self.logpdf_xy(x,y) - self.logpdf_x(x) - self.logpdf_y(y))
return dblquad(
lambda y,x: mi_integrand(x,y),self.D[0],self.D[1],self._lower_y,self._upper_y)
Source 的代码
不确定子函数 mi_integrand
中写的是什么公式或运算符 * \
的作用。看起来并不熟悉。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。