diff options
Diffstat (limited to 'sourcecodes/parameter_learning/Predictmultipleintrvention.m')
| -rw-r--r-- | sourcecodes/parameter_learning/Predictmultipleintrvention.m | 73 |
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 |
