def runTest(self): data = CompactFactor(read_csv(open('tetrad_asia.csv')),domain=Domain()) ci = PCCI(G2Separator(data)) g = ICPattern(ci) self.assertEquals(g.shd(self._asia_pdag),5) self.assertEquals(self._tetrad_pdag.shd(self._asia_pdag),4) # I think tetrad is wrong (in terms of implementation) self.assertEquals(g.shd(self._tetrad_pdag),1)