about summary refs log tree commit diff
path: root/sourcecodes/parameter_learning/Predictmultipleintrvention.m
diff options
context:
space:
mode:
Diffstat (limited to 'sourcecodes/parameter_learning/Predictmultipleintrvention.m')
-rw-r--r--sourcecodes/parameter_learning/Predictmultipleintrvention.m73
1 files changed, 73 insertions, 0 deletions
diff --git a/sourcecodes/parameter_learning/Predictmultipleintrvention.m b/sourcecodes/parameter_learning/Predictmultipleintrvention.m
new file mode 100644
index 00000000..d3d509cb
--- /dev/null
+++ b/sourcecodes/parameter_learning/Predictmultipleintrvention.m
@@ -0,0 +1,73 @@
+function Predictmultipleintrvention(pre)
+dfile=strcat(pre,'structure_input.txt');
+sfile=dfile;
+dfile=strcat(pre,'continuous_input.txt');
+nnodefile=strcat(pre,'nnode.txt');
+
+fnnode = fopen(nnodefile,'r');
+nnodes = fscanf(fnnode,'%d');
+
+fvarnamefile=strcat(pre,'varname.txt');
+
+varfile = fopen(fvarnamefile,'r');
+
+%nnodes=5;
+Std_flag=true;
+[labels,cases,bnet]=readInput(dfile,sfile,nnodes,Std_flag);
+[bnet]=parameterLearning(bnet,cases);
+
+%[predict_mean,predict_sd,q_sq]=looCrossValid(bnet,cases);
+
+
+fvarfile=strcat(pre,'var.txt');
+fvar = fopen(fvarfile,'r');
+                            
+select_var_new = fscanf(fvar,'%d');
+
+nm = numel(select_var_new);
+
+varlabels = cell(1,nm);
+varbuffer = fgetl(varfile);    %get header line as a string
+for j=1:nm
+    [varnext,varbuffer] = strtok(varbuffer);
+    varlabels{j} = varnext;
+    for i=1:nnodes    
+        if strcmp(varlabels{j},labels{i})
+            select_var_new(j)=i;
+        end
+     end    
+    
+end
+
+
+
+
+fvardfile=strcat(pre,'vardata.txt');
+
+fvard = fopen(fvardfile,'r');
+
+select_var_data_new = fscanf(fvard,'%f');
+
+filename=strcat(pre,'net_figure_new.txt');
+
+drawFigureM(nnodes,bnet,labels,filename,cases,select_var_new,select_var_data_new);
+
+%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);
+end
\ No newline at end of file