clear getgo who Your variables are: ds5 ds7 pg1 pg3 ds6 p pg2 pg4 pg4 pg4 = 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 1 1 0 1 0 cpts = learnNet(ds5,p); cpts cpts = Columns 1 through 7 0.4333 0.5667 -1.0000 -1.0000 -1.0000 -1.0000 -1.0000 0.8654 0.1346 0.5147 0.4853 -1.0000 -1.0000 -1.0000 0.5778 0.4222 1.0000 0 0.4286 0.5714 0.2121 0.3902 0.6098 0.3590 0.6410 0.3571 0.6429 0.3077 0.7619 0.2381 0.3529 0.6471 0.7273 0.2727 0.3953 Column 8 -1.0000 -1.0000 0.7879 0.6923 0.6047 likeNet(p,cpts,ds5) ans = -374.7600 diary off % parents matrices: % (defined in getgo.m) pg1 pg1 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 pg2 pg2 = 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 pg3 pg3 = 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 0 1 0 pg4 pg4 = 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 1 1 0 1 0 EDUğ cpt1 = learnNet(ds6,pg1); EDUğ cpt1 cpt1 = 0.3083 0.6917 0.6083 0.3917 0.6583 0.3417 0.6333 0.3667 0.3833 0.6167 EDUğ likeNet(pg1,cpt1,ds7) ans = -388.0465 EDUğ lik(1) = likeNet(pg1,cpt1,ds7) lik = -388.0465 EDUğ cpt2 = learnNet(ds6,pg2); EDUğ lik(2) = likeNet(pg2,cpt2,ds7) lik = -388.0465 -347.9791 EDUğ cpt3 = learnNet(ds6,pg3); EDUğ lik(3) = likeNet(pg3,cpt3,ds7) lik = -388.0465 -347.9791 -348.8306 EDUğ cpt4 = learnNet(ds6,pg4); EDUğ lik(4) = likeNet(pg4,cpt4,ds7) lik = -388.0465 -347.9791 -348.8306 -351.9809 EDUğ trnlik(1) = likeNet(pg1,cpt1,ds6); EDUğ trnlik(2) = likeNet(pg2,cpt2,ds6) trnlik = -390.2656 -339.6508 EDUğ trnlik(3) = likeNet(pg3,cpt3,ds6) trnlik = -390.2656 -339.6508 -338.7481 EDUğ trnlik(4) = likeNet(pg4,cpt4,ds6) trnlik = -390.2656 -339.6508 -338.7481 -335.9058 EDUğ lik lik = -388.0465 -347.9791 -348.8306 -351.9809 EDUğ trnlik trnlik = -390.2656 -339.6508 -338.7481 -335.9058 EDUğ plot(lik) EDUğ plot(lik,'.') EDUğ bar(lik) EDUğ plot(lik,'.-') EDUğ plot(trnlik,'r.-') EDUğ hold on EDUğ plot(lik,'.-') EDUğ title('log likelihood of various graph structures'); EDUğ xlabel('Graph #') EDUğ set(gca,'XTick',[1 2 3 4]) EDUğ ylabel('Log lik') EDUğ ylabel('Log likelihood') EDUğ legend('Train','Test')