function h = cond_entropy(pjoint)
% COND_ENTROPY Compute H(X|Y)
% h = cond_entropy(pjoint)
%
% pjoint(Y,X) : child is last dim.

pcond = mk_stochastic(pjoint);
h = -sum(sum(pjoint .* log2(pcond+eps)));
