Exemplos de ArenaBehaviorAnalysis.getOMRScoreStatsMulti em Python

Exemplo n.º 1

Arquivo: scratchpad_OMR_lh_8mm_safeandpost.py Projeto: nerduno/fishTrax

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)

Exemplo n.º 2

Arquivo: scratchpad_OMR_testing.py Projeto: nerduno/fishTrax

                   '/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()

Exemplo n.º 3

Arquivo: scratchpad_OMR_lh_8mm_post.py Projeto: nerduno/fishTrax

    "~/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)