In the paper, there seems like the R is a random variable. So the variance of Z shouldn't be calculate by (1-lambda)**2, the way calculated in the code(IBA-paper-code/attribution_bottleneck/bottleneck/per_sample_bottleneck.py). This way of calculation seems that the epsilon in the paper is a N(0,I) variable and R is a fixed variable.