我尝试使用mpmi软件包来计算两组连续变量之间的互信息。我对放在GutHub上的源代码感到困惑:
ans = ans + log(s12(i) / (s1(i) * s2(i)))
end do
ans = ans / lv + log(dble(lv))
s12看起来像p(x,y),而s1,s2看起来像p(x),p(y)。为什么不在log(s12(i) / (s1(i) * s2(i)))之前乘以p(x,y),因为计算MI的公式是p(x,y)*log(p(x,y)/(p(x)p(y)))
为什么在完成求和后会有ans = ans / lv + log(dble(lv))?