#%% plt.imshow(gc + 0.5*bsel, cmap='bwr') plt.tight_layout() #%% plt.imshow(bsel + gsel) #%% plt.imshow(bsel + rsel) # coef = ogre.calc_coloc_fromlabs(blabs, bc, gc, gmask) #%% from scipy.stats import pearsonr plive = pearsonr(bcslice.flatten(), gcslice.flatten())[0] pdead = pearsonr(bcslice.flatten(), rcslice.flatten())[0] 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)