The generative model of the functional CHURCH program is contained in the CHURCH function "take-a-sample". The direction of inference in the conditional P[D | L=0, S=1) = P(Difficulty | strong_letter, SAT=high) is inverse to the direction of CPDs in the generative Bayesian network.
Here is a summary of the domains due to the modifications of Gordon et al. (2014):
Val(D) = <d0, d1> = <easy, hard>
Val(I) = <i0, i1> = <non smart, smart>
Val(G) = <g0, g1> = <A, B+C> = <excellent, good+average>
Val(S) = <s0, s1> = <low score, high score>
Val(L) = <l0, l1> = <strong_letter, weak_letter>
The number of samples was set to 20000 in this run. This number could in principle be increased to get a better precision of estimates. The value of the variable math_expect (= parameter theta in the HOEFFDING inequality) was computed by the exact 'variable-elimination' inference method in the probabilistic programming language FIGARO. The sampling method used in our CHURCH-program is the simple-to-understand 'forward sampling'. The screen-shot presented was generated by using the PlaySpace environment of WebCHURCH.
The inferred E(D | L=0, S=1) = P(D=1 | L=0, S=1) is near 0.34 and lower than its unconditional counterparts E(D) = P(D=1) = 0.40. So the verbal interpretation is "If you have a strong recommendation letter and an above average SAT score, then the course is more easy for you than for the average student".