def refixLaunch(exp,corr=False, trim=True, bins=8):
"""
"""
src = 'selected_dm_%s_WITH_drift_corr_onlyControl_True.csv' % exp
dm = CsvReader(src).dataMatrix()
lat1 = "saccLat1"
lat2 = "saccLat2"
if corr:
fix1 = "endX1CorrNormToHandle"
fix2 = "endX2CorrNormToHandle"
else:
fix1 = "endX1NormToHandle"
fix2 = "endX2NormToHandle"
for var in [lat1, lat2, fix1, fix2]:
dm = dm.select("%s != ''" % var)
# Trim the data
if trim:
dm = dm.selectByStdDev(["file"],lat2,
verbose=False)
dm = dm.removeField("__dummyCond__")
dm = dm.removeField("__stdOutlier__")
dm = dm.selectByStdDev(["file"], fix2,
verbose=False)
dm = dm.addField('saccLat2_perc', dtype=float)
dm = dm.calcPerc(lat2, 'saccLat2_perc', nBin=bins)
cmLat1 = dm.collapse(['saccLat2_perc'], lat1)
cmX = cmLat2 = dm.collapse(['saccLat2_perc'], lat2)
cmFix1 = dm.collapse(['saccLat2_perc'], fix1)
cmFix2 = dm.collapse(['saccLat2_perc'], fix2)
fig = plt.figure(figsize = (3,6))
plt.subplots_adjust(left=.2, bottom=.15)
colList = [orange[1], blue[1]]
nRows = 3
nCols = 1
nPlot = 0
lTitles = ["Landing pos 2", "Landing pos 1"]
lTitles.reverse()
for cmY in [cmFix2, cmFix1]:#, cmLat1]:
#nPlot +=1
#plt.subplot(nRows, nCols, nPlot)
color = colList.pop()
plt.plot(cmX['mean'], cmY['mean'], marker = 'o', color=color, \
markerfacecolor='white', markeredgecolor=color, \
markeredgewidth=1)
plt.xlabel("Sacc lat 2")
plt.legend(lTitles, frameon=False, loc='best')
plt.axhline(0, linestyle = "--", color = gray[3])
plt.savefig("Launch_site_refixations_%s_corr_%s.png" % (exp, corr))
def ovp(exp,corr=False, trim=True, bins=8):
"""
"""
src = 'selected_dm_%s_WITH_drift_corr_onlyControl_True.csv' % exp
dm = CsvReader(src).dataMatrix()
lat1 = "saccLat1"
lat2 = "saccLat2"
prob = "saccCount"
dur1 = "durationFix1"
dur2 = "durationFix2"
total = "gazeDur"
rt = "rtFromStim"
if corr:
fix1 = "endX1CorrNormToHandle"
fix2 = "endX2CorrNormToHandle"
else:
fix1 = "endX1NormToHandle"
fix2 = "endX2NormToHandle"
for var in [lat1, lat2, fix1, fix2, dur1, dur2, total, rt]:
dm = dm.select("%s != ''" % var)
# Trim the data
if trim:
dm = dm.selectByStdDev(["file"],lat2,
verbose=False)
dm = dm.removeField("__dummyCond__")
dm = dm.removeField("__stdOutlier__")
dm = dm.selectByStdDev(["file"], fix2,
verbose=False)
dm = dm.addField('land2_perc', dtype=float)
dm = dm.calcPerc(fix2, 'land2_perc', nBin=bins)
cmLat1 = dm.collapse(['land2_perc'], lat1)
cmLat2 = dm.collapse(['land2_perc'], lat2)
cmFix1 = dm.collapse(['land2_perc'], fix1)
cmX = cmFix2 = dm.collapse(['land2_perc'], fix2)
cmProb = dm.collapse(['land2_perc'], prob)
cmDur1 = dm.collapse(['land2_perc'], dur1)
cmDur2 = dm.collapse(['land2_perc'], dur2)
cmTotal = dm.collapse(['land2_perc'], total)
cmRt = dm.collapse(['land2_perc'], rt)
fig = plt.figure(figsize = (3,10))
plt.subplots_adjust(left=.2, bottom=.15, hspace = .2)
lTitles = ["dur 1", "prob", "land 1", "rt", "total", "lat 1", "lat 2", "dur 2"]
lTitles.reverse()
nRows = len(lTitles)
nCols = 1
nPlot = 0
for cmY in [cmDur1, cmProb, cmFix1, cmRt, cmTotal, \
cmLat1, cmLat2, cmDur2]:
nPlot +=1
ax = plt.subplot(nRows, nCols, nPlot)
#color = colList.pop()
color = blue[1]
plt.plot(cmX['mean'], cmY['mean'], marker = 'o', color=color, \
markerfacecolor='white', markeredgecolor=color, \
markeredgewidth=1)
plt.ylabel(lTitles.pop())
plt.axvline(0, linestyle = "--", color = gray[3])
if nPlot == nRows:
plt.xlabel("Second landing position")
else:
ax.xaxis.set_ticklabels([])
plt.savefig("ovp2_%s_corr_%s.png" % (exp, corr))
