annotate toolboxes/FullBNT-1.0.7/KPMtools/compute_counts.m @ 0:cc4b1211e677 tip

initial commit to HG from Changeset: 646 (e263d8a21543) added further path and more save "camirversion.m"
author Daniel Wolff
date Fri, 19 Aug 2016 13:07:06 +0200
parents
children
rev   line source
Daniel@0 1 function count = compute_counts(data, sz)
Daniel@0 2 % COMPUTE_COUNTS Count the number of times each combination of discrete assignments occurs
Daniel@0 3 % count = compute_counts(data, sz)
Daniel@0 4 %
Daniel@0 5 % data(i,t) is the value of variable i in case t
Daniel@0 6 % sz(i) : values for variable i are assumed to be in [1:sz(i)]
Daniel@0 7 %
Daniel@0 8 % Example: to compute a transition matrix for an HMM from a sequence of labeled states:
Daniel@0 9 % transmat = mk_stochastic(compute_counts([seq(1:end-1); seq(2:end)], [nstates nstates]));
Daniel@0 10
Daniel@0 11 assert(length(sz) == size(data, 1));
Daniel@0 12 P = prod(sz);
Daniel@0 13 indices = subv2ind(sz, data'); % each row of data' is a case
Daniel@0 14 %count = histc(indices, 1:P);
Daniel@0 15 count = hist(indices, 1:P);
Daniel@0 16 count = myreshape(count, sz);
Daniel@0 17