wolffd@0: function oracle = mk_vstruct_bnet()
wolffd@0: % MK_VSTRUCT_BNET Make a simple V-structured 3-node noisy-AND Bayes net
wolffd@0: % oracle = mk_vstruct_bnet()
wolffd@0: 
wolffd@0: N = 3;
wolffd@0: dag = zeros(N);
wolffd@0: A = 1; B = 2; C = 3;
wolffd@0: dag(A,C)=1;
wolffd@0: dag(B,C)=1;
wolffd@0: ns = 2*ones(1,N);
wolffd@0: 
wolffd@0: oracle = mk_bnet(dag, ns);
wolffd@0: oracle.CPD{1} = tabular_CPD(oracle, 1, [0.5 0.5]);
wolffd@0: oracle.CPD{2} = tabular_CPD(oracle, 2, [0.5 0.5]);
wolffd@0: pnoise = 0.1; % degree of noise
wolffd@0: oracle.CPD{3} = boolean_CPD(oracle, 3, 'named', 'all', pnoise);