@statistic('diff_sdev', pairwise=True, is_distribution=True, is_estimate=True)
    def get_diff_sdev(self, key):
        sigmas_1 = self.mcmcs[key].trace('sigma_1')[:]
        sigmas_2 = self.mcmcs[key].trace('sigma_2')[:]
        return sigmas_2 ** 0.5 - sigmas_1 ** 0.5
Compute difference of standard deviations from sampler
dataset/ETHPy150Open HHammond/kcbo/kcbo/statistical_tests/t_test.py/TTest.get_diff_sdev