More info: when calculated through passing the stats_funcs parameter of the az.stats.summary(), it results in the same value of 7.339, so something happens to the trace on the way, which is not done when calling az.plots.plot_utils.calculate_point_estimate('mode', trace['tau'], bw=4.5) directly.