Mercurial > hg > camir-ismir2012

diff toolboxes/FullBNT-1.0.7/graph/mk_rnd_dag_given_edge_prob.m @ 0:cc4b1211e677 tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/toolboxes/FullBNT-1.0.7/graph/mk_rnd_dag_given_edge_prob.m	Fri Aug 19 13:07:06 2016 +0200
@@ -0,0 +1,15 @@
+function dag = sample_dag(P)
+% SAMPLE_DAG Create a random directed acyclic graph with edge probabilities P(i,j)
+% dag = sample_dag(P)
+%
+% This uses rejection sampling to reject graphs with directed cycles.
+
+done = 0;
+directed = 1;
+iter = 1;
+while ~done
+  dag = binornd(1, P); % each edge is an indep Bernoulli (0/1) random variable
+  dag = setdiag(dag, 0);
+  done = acyclic(dag, directed);
+  iter = iter + 1
+end