VI node sampling with multiple parameters

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})