diff 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"

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/toolboxes/FullBNT-1.0.7/KPMstats/sample.m	Fri Aug 19 13:07:06 2016 +0200
@@ -0,0 +1,15 @@
+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