this code doesn't produce the 5% quantile:

v = np.sort( posteriors[i] )[ int(0.05*N) ]

this would do it:

v = np.sort( posteriors[i] , axis=0)[ int(0.05*N) ]