def testSigma(self):
unadjusted = sl.stouffer_liptak(self.pvals)
sigma = np.eye(3)
sigma[0,1] = sigma[1,0] = 0.9
sigma[2,1] = sigma[1,2] = 0.7
sigma[2,2] = sigma[1,2] = 0.5
sigma[2,0] = sigma[0,2] = 0.6
adjusted = sl.stouffer_liptak(self.pvals, sigma)
self.assert_(adjusted['p'] > unadjusted['p'])
def testSigma(self):
unadjusted = sl.stouffer_liptak(self.pvals)
sigma = np.eye(3)
sigma[0, 1] = sigma[1, 0] = 0.9
sigma[2, 1] = sigma[1, 2] = 0.7
sigma[2, 2] = sigma[1, 2] = 0.5
sigma[2, 0] = sigma[0, 2] = 0.6
adjusted = sl.stouffer_liptak(self.pvals, sigma)
self.assert_(adjusted['p'] > unadjusted['p'])
def testSLOK(self):
res = sl.stouffer_liptak(self.pvals)
self.assert_(all(k in res for k in "OK p C".split()))
def testSLOK(self):
res = sl.stouffer_liptak(self.pvals)
self.assert_(all(k in res for k in "OK p C".split()))