view toolboxes/FullBNT-1.0.7/KPMstats/sample.m @ 0:cc4b1211e677 tip

initial commit to HG from
Changeset:
646 (e263d8a21543) added further path and more save "camirversion.m"

function x = sample(p, n)
% SAMPLE    Sample from categorical distribution.
% Returns a row vector of integers, sampled according to the probability
% distribution p.
% Uses the stick-breaking algorithm.
% Much faster algorithms are also possible.

if nargin < 2
  n = 1;
end

cdf = cumsum(p(:));
for i = 1:n
  x(i) = sum(cdf < rand) + 1;
end