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