about summary refs log tree commit diff
path: root/BNW_parameter_learning/runBN_initial.m
diff options
context:
space:
mode:
authorziejd22018-03-14 23:19:16 -0500
committerziejd22018-03-14 23:19:16 -0500
commitc80226899f5cdd9f11c163817d59445213f5bef0 (patch)
treee0fd79d2e32fd2aedda2eadaed0f19af3514c520 /BNW_parameter_learning/runBN_initial.m
parent324ebc8ddacab8e154047b8518afd5cbb5bb2fa4 (diff)
downloadBNW-c80226899f5cdd9f11c163817d59445213f5bef0.tar.gz
Separating Octave and php calculations
Diffstat (limited to 'BNW_parameter_learning/runBN_initial.m')
-rw-r--r--BNW_parameter_learning/runBN_initial.m87
1 files changed, 23 insertions, 64 deletions
diff --git a/BNW_parameter_learning/runBN_initial.m b/BNW_parameter_learning/runBN_initial.m
index c3f2a34b..43caf402 100644
--- a/BNW_parameter_learning/runBN_initial.m
+++ b/BNW_parameter_learning/runBN_initial.m
@@ -15,84 +15,43 @@ mapfile = fopen(mapfilename,'w');
 mapval = fopen(mapvalfilename,'w');

 

 

-%nnodes=5;

 Std_flag=true;

 [labels,cases,bnet,node_sizes,data,labelsold]=readInput(dfile,sfile,nnodes,Std_flag);

-s = std(data,0,1);

+s=std(data,0,1);

 m=mean(data);

 

-

 for i=1:nnodes

   fprintf(mapval,'%s\t%d\t%f\t%f\n',labelsold{i},node_sizes(i),s(i),m(i));

 end

 

-

-% for j=1:nnodes

-%     [next,buffer] = strtok(buffer);

-%     name{j}=next;

-%     for i=1:nnodes    

-%         if strcmp(name{j},labels{i})

-%             map{j}=i;

-%             fprintf(mapfile,'%d\t',i);

-%         end

-%      end

-% end

-%name

-%labels

-%map

 fprintf(mapfile,'%s',labels{1});

 for i=2:nnodes

   fprintf(mapfile,'\t%s',labels{i});

 end

 fprintf(mapfile,'\n');

+fclose(mapval);

+fclose(mapfile);

+

+%Need to rearrange the means and stdevs to match the new labeling.

+means = cell(1,nnodes);

+stdevs = cell(1,nnodes);

+for i = 1:nnodes

+    for j = 1:nnodes

+       if strcmp(labels{i},labelsold{j})

+          means{i} = m(j);

+          stdevs{i} = s(j);

+          break

+       end

+    end

+end

+

 

 [bnet]=parameterLearning(bnet,cases);

-%[predict_mean,predict_sd,q_sq]=looCrossValid(bnet,cases);

-%engine=jtree_inf_engine(bnet);

-%evidence=cell(1,nnodes);

-

-%varfile='var.txt';

-%fvar = fopen(varfile,'r');

-%select_var = fscanf(fvar,'%d');

-%select_var=map{select_var};

-%varfiled='vardata.txt';

-%fvard = fopen(varfiled,'r');

-%select_var_data = fscanf(fvard,'%f');

-

-%evidence{select_var}=select_var_data;

-%[engine,loglik]=enter_evidence(engine,evidence);

-

-%outdata='prediction.txt';

-%fout = fopen(outdata,'w');

-

-%for ii = 1:nnodes 

- % i=map{ii};

- % data=marginal_nodes(engine,i);

- % fprintf(fout,'%d\t%d\t%f\t%f\t%f\n',ii,data.domain,data.T,data.mu,data.Sigma);

-  %fprintf(1,'%d\n',i);

-% end

 

 filename=strcat(pre,'net_figure.txt');

-drawFigure(nnodes,bnet,labels,filename,cases);

-

-%quit force;

-%marginal_nodes(engine,2)

-%marginal_nodes(engine,3)

-%marginal_nodes(engine,4)

-%marginal_nodes(engine,5)

-%evidence{1}=2;

-%[engine,loglik]=enter_evidence(engine,evidence)

-%marginal_nodes(engine,1)

-%marginal_nodes(engine,2)

-%marginal_nodes(engine,3)

-%marginal_nodes(engine,4)

-%marginal_nodes(engine,5)

-%evidence{2}=0.6;

-%evidence{1}=[];

-%[engine,loglik]=enter_evidence(engine,evidence);

-%marginal_nodes(engine,3);

-%marginal_nodes(engine,4);

-%marginal_nodes(engine,5);

-fclose(mapval);

-fclose(mapfile);

-end
\ No newline at end of file
+

+drawFigure(nnodes,bnet,labels,filename,cases,stdevs,means);

+

+writeParameters(pre,nnodes,bnet,labels,cases,labelsold,s,m);

+

+end