From 7cc31810d53176e805532b2789955f4eedbce6bb Mon Sep 17 00:00:00 2001
From: ziejd2
Date: Thu, 14 Sep 2017 16:17:47 -0500
Subject: Add files via upload

Adding the *.m files used in parameter learning.---
 .../Predictmultipleintrvention.m                   | 73 ++++++++++++++++++++++
 1 file changed, 73 insertions(+)
 create mode 100644 BNW_parameter_learning/Predictmultipleintrvention.m

(limited to 'BNW_parameter_learning/Predictmultipleintrvention.m')

diff --git a/BNW_parameter_learning/Predictmultipleintrvention.m b/BNW_parameter_learning/Predictmultipleintrvention.m
new file mode 100644
index 00000000..4675b846
--- /dev/null
+++ b/BNW_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
-- 
cgit 1.4.1