test_Bayesian.py 文件源码

def test_fisher_matrix(self):
        """
        """
        new = self.bayes.update(simulations={'simple1': [self.sim1, self.exp1],
                                             'simple2': [self.sim1, self.exp1]},
                                models={'simp':
                                        self.model.update_dof([1, 0.5])})

        initial_data = new.get_data()
        sens_matrix = new._get_sens(initial_data=initial_data)

        sigma = inv(new.simulations['simple1']['exp'].get_sigma())

        fisher_1 = np.dot(np.dot(sens_matrix['simple1'].T,
                                 sigma),
                          sens_matrix['simple1'])

        fisher_2 = np.dot(np.dot(sens_matrix['simple2'].T,
                                 sigma),
                          sens_matrix['simple2'])

        npt.assert_array_almost_equal(self.exp1.get_fisher_matrix(
                                      sens_matrix=sens_matrix['simple1']),
                                      fisher_1)
        npt.assert_array_almost_equal(self.exp1.get_fisher_matrix(
                                      sens_matrix=sens_matrix['simple2']),
                                      fisher_2)