Daniel@0: % SGS p141 (female orgasm data set) Daniel@0: Daniel@0: C = eye(7,7); Daniel@0: C(2,1:1) = [-0.132]; Daniel@0: C(3,1:2) = [0.009 -0.136]; Daniel@0: C(4,1:3) = [0.22 -0.166 0.403]; Daniel@0: C(5,1:4) = [-0.008 0.008 0.598 0.282]; Daniel@0: C(6,1:5) = [0.119 -0.076 0.264 0.514 0.176]; Daniel@0: C(7,1:6) = [0.118 -0.137 0.368 0.414 0.336 0.338]; Daniel@0: Daniel@0: n = 7; Daniel@0: for i=1:n Daniel@0: for j=i+1:n Daniel@0: C(i,j)=C(j,i); Daniel@0: end Daniel@0: end Daniel@0: Daniel@0: max_fan_in = 4; Daniel@0: nsamples = 281; Daniel@0: alpha = 0.05; Daniel@0: pdag = learn_struct_pdag_pc('cond_indep_fisher_z', n, max_fan_in, C, nsamples, alpha)