wolffd@0: % SGS p118
wolffd@0: % Try learning the structure using an oracle for the cond indep tests
wolffd@0: 
wolffd@0: n = 5;
wolffd@0: 
wolffd@0: A = 1; B = 2; C = 3; D = 4; E = 5;
wolffd@0: 
wolffd@0: G = zeros(n);
wolffd@0: G(A,B)=1;
wolffd@0: G(B,[C D]) = 1;
wolffd@0: G(C,E)=1;
wolffd@0: G(D,E)=1;
wolffd@0: 
wolffd@0: k = 2;
wolffd@0: 
wolffd@0: pdag = learn_struct_pdag_pc('dsep', n, k, G)
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: if 0
wolffd@0: N = 4; 
wolffd@0: dag = zeros(N,N);
wolffd@0: C = 1; S = 2; R = 3; W = 4;
wolffd@0: dag(C,[R S]) = 1;
wolffd@0: dag(R,W) = 1;
wolffd@0: dag(S,W)=1;
wolffd@0: 
wolffd@0: pdag = learn_struct_pdag_pc('dsep', N, 2, dag)
wolffd@0: end