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