@statistic('effect_size', pairwise=True, is_distribution=True, is_estimate=True)
    def get_effect_size(self, key):
        mus_1 = self.mcmcs[key].trace('mu_1')[:]
        mus_2 = self.mcmcs[key].trace('mu_2')[:]
        sigmas_1 = self.mcmcs[key].trace('sigma_1')[:]
        sigmas_2 = self.mcmcs[key].trace('sigma_2')[:]
        return (mus_2 - mus_1) / (np.sqrt((sigmas_2 + sigmas_1) / 2.0))
Compute effect size from sampler
dataset/ETHPy150Open HHammond/kcbo/kcbo/statistical_tests/t_test.py/TTest.get_effect_size