For the following:
_sample_proba_nn1 = approx.sample_node({bsnn.out1.distribution.p, bsnn.out2.distribution.p},
size = n1,
more_replacements = {bsnn['bnn_input1']: x1,
bsnn['bnn_input2']: x2})
I get:
TypeError: ('output must be a theano Variable or Out instance (or list of them)', {sigmoid.0, sigmoid.0})