# HG changeset patch # User samer # Date 1328739846 0 # Node ID 917d5d11aa48fcb7f6043d3de73ea657b1c8df99 # Parent 0f16c069423bec2eb36bb0cbfc8b6db18539d757 mt_show_transmat now computes information measures if not given. diff -r 0f16c069423b -r 917d5d11aa48 mt_show_transmat.m --- a/mt_show_transmat.m Mon Feb 06 15:55:08 2012 +0000 +++ b/mt_show_transmat.m Wed Feb 08 22:24:06 2012 +0000 @@ -9,7 +9,6 @@ figure(Id); imagesc(T,[0,1]); axis xy; - if nargin>2 - title(sprintf('H=%.2f, R=%.2f, PI=%.2f',I(1),I(2),I(3))); - end + if nargin<3, I=mc_global_info(T); + title(sprintf('H=%.2f, R=%.2f, PI=%.2f',I(1),I(2),I(3))); end diff -r 0f16c069423b -r 917d5d11aa48 private/mc_global_info.m --- a/private/mc_global_info.m Mon Feb 06 15:55:08 2012 +0000 +++ b/private/mc_global_info.m Wed Feb 08 22:24:06 2012 +0000 @@ -12,6 +12,7 @@ % Entropy rate H(X|Z) % Redundancy I(X,Z) % Predictive information rate I(X,Y|Z) +% All in NATS, not bits. n=size(T,1);