The decoder described in http://edwardlib.org/tutorials/decoder

should be generating variance as well as the mean of the reconstructed X.

However, the code has just

`x = Bernoulli(logits=tf.layers.dense(hidden, 28 * 28))`

which is only the mean. How would you capture variance of the reconstruction as well?