1 - d2.cdf(xx), err_msg='cdf' + d1.name) assert_almost_equal(d1.pdf(nxx), d2.pdf(xx)) assert_almost_equal(d1.sf(nxx), 1 - d2.sf(xx)) assert_almost_equal(d1.ppf(ppfq), -d2.ppf(ppfq)[::-1]) assert_almost_equal(d1.isf(ppfq), -d2.isf(ppfq)[::-1]) assert_almost_equal(d1.moment(3), -d2.moment(3)) ch2oddneg = [ v * (-1)**(i + 1) for i, v in enumerate(d2.stats(moments='mvsk')) ] assert_almost_equal(d1.stats(moments='mvsk'), ch2oddneg, err_msg='stats ' + d1.name + d2.name) if __name__ == '__main__': tt = Test_Transf2() tt.test_equivalent() tt.test_equivalent_negsq() debug = 0 if debug: print(negsquarenormalg.ppf([0.1, 0.5, 0.9])) print(stats.chi2.ppf([0.1, 0.5, 0.9], 1)) print(negsquarenormalg.a) print(negsquarenormalg.b) print(absnormalg.stats(moments='mvsk')) print(stats.foldnorm(1e-10).stats(moments='mvsk')) print(stats.halfnorm.stats(moments='mvsk'))
xx, nxx, ppfq = self.xx, self.nxx, self.ppfq d1,d2 = (negsquarenormalg, stats.chi2(1)) #print d1.name assert_almost_equal(d1.cdf(nxx), 1-d2.cdf(xx), err_msg='cdf'+d1.name) assert_almost_equal(d1.pdf(nxx), d2.pdf(xx)) assert_almost_equal(d1.sf(nxx), 1-d2.sf(xx)) assert_almost_equal(d1.ppf(ppfq), -d2.ppf(ppfq)[::-1]) assert_almost_equal(d1.isf(ppfq), -d2.isf(ppfq)[::-1]) assert_almost_equal(d1.moment(3), -d2.moment(3)) ch2oddneg = [v*(-1)**(i+1) for i,v in enumerate(d2.stats(moments='mvsk'))] assert_almost_equal(d1.stats(moments='mvsk'), ch2oddneg, err_msg='stats '+d1.name+d2.name) if __name__ == '__main__': tt = Test_Transf2() tt.test_equivalent() tt.test_equivalent_negsq() debug = 0 if debug: print(negsquarenormalg.ppf([0.1,0.5,0.9])) print(stats.chi2.ppf([0.1,0.5,0.9],1)) print(negsquarenormalg.a) print(negsquarenormalg.b) print(absnormalg.stats( moments='mvsk')) print(stats.foldnorm(1e-10).stats( moments='mvsk')) print(stats.halfnorm.stats( moments='mvsk'))