print 'plive = {0}, pdead = {1}'.format(plive, pdead) plive1 = ogre.pearsonr(bcslice.flatten(), gcslice.flatten())[0] pdead1 = ogre.pearsonr(bcslice.flatten(), rcslice.flatten())[0] print 'plive1 = {0}, pdead1 = {1}'.format(plive1, pdead1) plive2 = ogre.pearsonr2(bcslice, gcslice) pdead2 = ogre.pearsonr2(bcslice, rcslice) print 'plive2 = {0}, pdead2 = {1}'.format(plive2, pdead2) mlive = ogre.mandersr(bcslice, gcslice) mdead = ogre.mandersr(bcslice, rcslice) print 'mlive = {0}, mdead = {1}'.format(mlive, mdead) klive = ogre.overlapk(bcslice, gcslice) kdead = ogre.overlapk(bcslice, rcslice) print 'klive = {0}, kdead = {1}'.format(klive, kdead) '''need to figure out how the masks should work for bmask, gmask, and rmask''' #Mlive = ogre.overlapM(bcslice, gcslice, bbox_slice(bc*(blabs == sel), bprops[sel].bbox), bbox_slice(gc*(blabs == sel), bprops[sel].bbox)) #Mdead = ogre.overlapM(bcslice, gcslice, bmask, rmask) #print 'Mlive = {0}, Mdead = {1}'.format(Mlive, Mdead) # pearson = ogre.calc_PCC_fromlabs(blabs, bc, gc) # print pearson #%%