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.
?-
mc_sample(dm,1000,Prob)
. ?-mc_mh_sample(dm,e,10000,Prob)
.*/