@statistic('diff_means', pairwise=True, is_distribution=True, is_estimate=True)
    def get_diff_means(self, key):
        mus_1 = self.mcmcs[key].trace('mu_1')[:]
        mus_2 = self.mcmcs[key].trace('mu_2')[:]
        return mus_2 - mus_1
Compute difference of means from sampler
dataset/ETHPy150Open HHammond/kcbo/kcbo/statistical_tests/t_test.py/TTest.get_diff_means