A self-contained package for computing mutual information, joint/conditional probability, entropy
function h = mutualinfo(vec1,vec2)
%This is a prog in the MutualInfo 0.9 package written by
% Hanchuan Peng.
%Disclaimer: The author of program is Hanchuan Peng
% at <email@example.com> and <firstname.lastname@example.org>.
%The CopyRight is reserved by the author.
%Last modification: April/19/2002
% h = mutualinfo(vec1,vec2)
% calculate the mutual information of two vectors
% By Hanchuan Peng, April/2002
[p12, p1, p2] = estpab(vec1,vec2);
h = estmutualinfo(p12,p1,p2);