import pylab import numpy as np import matplotlib as mpl import matplotlib.pyplot as pyplot import pylab import scipy #omr stats tp = 12 ''' ePreOMRstats = aba.getOMRScoreStatsMulti(e_omr, stateRange=[3,3], timePoint=tp) eEndOMRstats = aba.getOMRScoreStatsMulti(e_omr, stateRange=[8,8], timePoint=tp) ePostOMRstats_2a= aba.getOMRScoreStatsMulti(e_omr_post, stateRange=[3,3], timePoint=tp) ePostOMRstats_2b = aba.getOMRScoreStatsMulti(e_omr_post, stateRange=[8,8], timePoint=tp) ''' ePreOMRstatsSame_1a= aba.getOMRScoreStatsMulti(e_omr_same_pre, stateRange=[3,3], timePoint=tp) ePreOMRstatsSame_1b = aba.getOMRScoreStatsMulti(e_omr_same_pre, stateRange=[8,8], timePoint=tp) ePreOMRstatsSame = aba.getOMRScoreStatsMulti(e_omr_same, stateRange=[3,3], timePoint=tp) #eEndOMRstatsSame = aba.getOMRScoreStatsMulti(e_omr_same, stateRange=[8,8], timePoint=tp) eEndOMRstatsSame = aba.getOMRScoreStatsMulti(e_omr_same_last, stateRange=[8,8], timePoint=tp) ePostOMRstatsSame_2a= aba.getOMRScoreStatsMulti(e_omr_same_post, stateRange=[3,3], timePoint=tp) ePostOMRstatsSame_2b = aba.getOMRScoreStatsMulti(e_omr_same_post, stateRange=[8,8], timePoint=tp) cPreOMRstats_1a= aba.getOMRScoreStatsMulti(c_omr_pre, stateRange=[3,3], timePoint=tp) cPreOMRstats_1b = aba.getOMRScoreStatsMulti(c_omr_pre, stateRange=[8,8], timePoint=tp) cPreOMRstats = aba.getOMRScoreStatsMulti(c_omr, stateRange=[3,3], timePoint=tp) #cEndOMRstats = aba.getOMRScoreStatsMulti(c_omr, stateRange=[8,8], timePoint=tp) cEndOMRstats = aba.getOMRScoreStatsMulti(c_omr_last, stateRange=[8,8], timePoint=tp) cPostOMRstats_2a= aba.getOMRScoreStatsMulti(c_omr_post, stateRange=[3,3], timePoint=tp) cPostOMRstats_2b = aba.getOMRScoreStatsMulti(c_omr_post, stateRange=[8,8], timePoint=tp)
'/home/vburns/Dropbox/ConchisData/2013-09-30/f00686/f00686_2013-09-30-19-25-24.json', '/home/vburns/Dropbox/ConchisData/2013-09-30/f00687/f00687_2013-09-30-19-25-26.json', '/home/vburns/Dropbox/ConchisData/2013-09-30/f00688/f00688_2013-09-30-19-25-29.json' ] c_omr_10mm_slow = aba.loadMultipleDataFiles(c_omr_10mm_slow) c_1 = c_omr_10mm_fast cTitle = '10mm fast' e_1 = c_omr_5mm eTitle = '5mm' #omr analysis omr_tp = 12 eBaseResults = aba.getOMRScoreStatsMulti(e_1, stateRange=[3,3],timePoint=omr_tp) eEndResults = aba.getOMRScoreStatsMulti(e_1, stateRange=[8,8],timePoint=omr_tp) cBaseResults = aba.getOMRScoreStatsMulti(c_1, stateRange=[3,3],timePoint=omr_tp) cEndResults = aba.getOMRScoreStatsMulti(c_1, stateRange=[8,8],timePoint=omr_tp) def plotOMRMetric(cond1base, cond1end, cond2base, cond2end, title, baseTitle, endTitle, cond1label, cond2label): ax = pyplot.gca() ax.plot([0,1],[np.mean(cond1base), np.mean(cond1end)],'o-k',lw=3, label=cond1label) sem2 = (2*scipy.stats.sem(cond1base), 2*scipy.stats.sem(cond1end)) pyplot.errorbar([0,1],[np.mean(cond1base), np.mean(cond1end)],fmt='ok',yerr=sem2, lw=3) ax.plot([0,1],[np.mean(cond2base), np.mean(cond2end)],'o-m',lw=3, label=cond2label) sem2 = (2*scipy.stats.sem(cond2base), 2*scipy.stats.sem(cond2end)) pyplot.errorbar([0,1],[np.mean(cond2base), np.mean(cond2end)],fmt='om',yerr=sem2, lw=3) handles, labels= ax.get_legend_handles_labels()
"~/Dropbox/ConchisData/2013-10-03/f00767/f00767_2013-10-03-16-33-03.json", "~/Dropbox/ConchisData/2013-10-03/f00768/f00768_2013-10-03-16-33-16.json", "~/Dropbox/ConchisData/2013-10-04/f00769/f00769_2013-10-04-11-11-05.json", "~/Dropbox/ConchisData/2013-10-04/f00770/f00770_2013-10-04-11-11-14.json", "~/Dropbox/ConchisData/2013-10-04/f00771/f00771_2013-10-04-11-11-29.json", "~/Dropbox/ConchisData/2013-10-04/f00772/f00772_2013-10-04-11-11-35.json", # '~/Dropbox/ConchisData/2013-10-04/f00773/f00773_2013-10-04-11-12-15.json', "~/Dropbox/ConchisData/2013-10-04/f00774/f00774_2013-10-04-11-12-22.json", "~/Dropbox/ConchisData/2013-10-04/f00775/f00775_2013-10-04-11-12-28.json", # '~/Dropbox/ConchisData/2013-10-04/f00776/f00776_2013-10-04-11-12-36.json' ] e_omr_post = aba.loadMultipleDataFiles(e_omr_post) # omr stats tp = 12 cPreOMRstats = aba.getOMRScoreStatsMulti(c_omr, stateRange=[3, 3], timePoint=tp) cPostOMRstats = aba.getOMRScoreStatsMulti(c_omr, stateRange=[8, 8], timePoint=tp) ePreOMRstats = aba.getOMRScoreStatsMulti(e_omr, stateRange=[3, 3], timePoint=tp) ePostOMRstats = aba.getOMRScoreStatsMulti(e_omr, stateRange=[8, 8], timePoint=tp) esafePreOMRstats = aba.getOMRScoreStatsMulti(e_omr_post, stateRange=[3, 3], timePoint=tp) esafePostOMRstats = aba.getOMRScoreStatsMulti(e_omr_post, stateRange=[8, 8], timePoint=tp) # velocity stats sm = 15 cPreVel = aba.getMedianVelMulti(c_omr, stateRange=[3, 3], smoothWinLen=sm) cPostVel = aba.getMedianVelMulti(c_omr, stateRange=[8, 8], smoothWinLen=sm) ePreVel = aba.getMedianVelMulti(e_omr, stateRange=[3, 3], smoothWinLen=sm) ePostVel = aba.getMedianVelMulti(e_omr, stateRange=[8, 8], smoothWinLen=sm) esafePreVel = aba.getMedianVelMulti(e_omr_post, stateRange=[3, 3], smoothWinLen=sm) esafePostVel = aba.getMedianVelMulti(e_omr_post, stateRange=[8, 8], smoothWinLen=sm)