Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/KPMtools/marg_table.m @ 0:e9a9cd732c1e tip
first hg version after svn
| author | wolffd |
|---|---|
| date | Tue, 10 Feb 2015 15:05:51 +0000 |
| parents | |
| children |
| rev | line source |
|---|---|
| wolffd@0 | 1 function smallT = marg_table(bigT, bigdom, bigsz, onto, maximize) |
| wolffd@0 | 2 % MARG_TABLE Marginalize a table |
| wolffd@0 | 3 % smallT = marg_table(bigT, bigdom, bigsz, onto, maximize) |
| wolffd@0 | 4 |
| wolffd@0 | 5 if nargin < 5, maximize = 0; end |
| wolffd@0 | 6 |
| wolffd@0 | 7 |
| wolffd@0 | 8 smallT = myreshape(bigT, bigsz); % make sure it is a multi-dim array |
| wolffd@0 | 9 sum_over = mysetdiff(bigdom, onto); |
| wolffd@0 | 10 ndx = find_equiv_posns(sum_over, bigdom); |
| wolffd@0 | 11 if maximize |
| wolffd@0 | 12 for i=1:length(ndx) |
| wolffd@0 | 13 smallT = max(smallT, [], ndx(i)); |
| wolffd@0 | 14 end |
| wolffd@0 | 15 else |
| wolffd@0 | 16 for i=1:length(ndx) |
| wolffd@0 | 17 smallT = sum(smallT, ndx(i)); |
| wolffd@0 | 18 end |
| wolffd@0 | 19 end |
| wolffd@0 | 20 |
| wolffd@0 | 21 |
| wolffd@0 | 22 ns = zeros(1, max(bigdom)); |
| wolffd@0 | 23 %ns(bigdom) = mysize(bigT); % ignores trailing dimensions of size 1 |
| wolffd@0 | 24 ns(bigdom) = bigsz; |
| wolffd@0 | 25 |
| wolffd@0 | 26 smallT = squeeze(smallT); % remove all dimensions of size 1 |
| wolffd@0 | 27 smallT = myreshape(smallT, ns(onto)); % put back relevant dims of size 1 |