?- mc_sample(dm,1000,Prob). ?- mc_mh_sample(dm,e,10000,Prob).

*/

    8:- use_module(library(mcintyre)).    9
   10:- if(current_predicate(use_rendering/1)).   11:- use_rendering(c3).   12:- use_rendering(graphviz).   13:- endif.   14
   15:- mc.   16:- begin_lpad.   17
   18predisposition(average):0.698;predisposition(moderate):0.227;predisposition(high):0.075.
   19dm:0.054:-
   20  predisposition(average).
   21dm:0.131:-
   22  predisposition(moderate).
   23dm:0.266:-
   24  predisposition(high).
   25
   26gluc_if_dm(G):gaussian(G,7.5,3.8).
   27gluc_if_not_dm(G):gaussian(G,5.79,0.98).
   28
   29noise_if_dm(N):gaussian(N,0.0,3.3).
   30noise_if_not_dm(N):gaussian(N,0.0,0.3).
   31
   32hba1c(H):-
   33    dm,
   34    gluc_if_dm(G),
   35    noise_if_dm(N),
   36    {H=:=1.4+0.92*G+N}.
   37
   38hba1c(H):-
   39    \+ dm,
   40    gluc_if_not_dm(G),
   41    noise_if_not_dm(N),
   42    {H=:=0.6+0.9*G+N}.
   43
   44e:- hba1c(H),{H>7.2}.
   45:- end_lpad.