[t, time_diff_five] = scipy.stats.ttest_ind(np.mean(ef_frac,1), np.mean(c_frac,1))
[t, time_diff_old] = scipy.stats.ttest_ind(np.mean(et_frac,1), np.mean(c_frac,1))
[t, dist_diff_five] = scipy.stats.ttest_ind(np.mean(ef_dist,1), np.mean(c_dist,1))
[t, dist_diff_old] = scipy.stats.ttest_ind(np.mean(et_dist,1), np.mean(c_dist,1))
print 'BY ANIMAL 5 DPF Avoidance Post LH_5V (frac, distance): ', ef_frac_stat, ef_dist_stat
print 'BY ANIMAL 13 DPF Avoidance Post LH_5V (frac, distance): ', et_frac_stat, et_dist_stat
print 'BY ANIMAL Controls diff from Exper (time, dist) at 5 days: ', time_diff_five, dist_diff_five
print 'BY ANIMAL Controls diff from Exper (time, dist) at 13 days: ', time_diff_five, dist_diff_five
for n in range(len(f_fish)):
pylab.figure(n+10);
ax1 = pylab.subplot2grid((4,1),(0,0),rowspan=2);
aba.plotFishXPosition(f_fish[n])
ax = pylab.subplot2grid((4,1),(2,0),sharex=ax1);
[sh,s,d,dsh,ds,st] = aba.getSidePreference(f_fish[n], cond=[3,4], refState='On')
pylab.bar(st,np.array(s)/np.array(d),width=d)
pylab.axhline(.5,color='k')
pylab.ylabel('% time side1')
ax.set_yticks((0,.5,1))
ax = pylab.subplot2grid((4,1),(3,0),sharex=ax1);
pylab.bar(st,np.array(ds),width=d);
pylab.ylim((0,48))
pylab.axhline(24,color='k')
ax.set_yticks((0,24,48))
pylab.ylabel('Avg dist from\nside 1 (mm)')
pylab.xlabel('Time (s)')
"""
pylab.figure(5)
pylab.subplot(2,2,1)
aba.plotFishXPosition(t_fish[0])
[tv, c_distRT_stat] = scipy.stats.ttest_1samp(np.mean(c_distRT, axis = 1), 24)
[t, time_diff] = scipy.stats.ttest_ind(np.mean(e_fracRT,axis = 1), np.mean(c_fracRT,1))
[t, dist_diff] = scipy.stats.ttest_ind(np.mean(e_distRT,axis = 1), np.mean(c_distRT,1))
print 'BY ANIMAL Avoidance Post LH_5V 14dpf (frac, distance): ', e_fracRT_stat, e_distRT_stat
print 'By ANIMAL Avoidance Post LH_5V 27dpf (frac, distance):', e_fracRT2_stat, e_distRT2_stat
print 'BY ANIMAL Avoidance Control: ', c_fracRT_stat, c_distRT_stat
print 'BY ANIMAL Controls diff from Exper (time, dist): ', time_diff, dist_diff
e_RTn = e_RT + e_RT2
for n in range(len(e_RTn)):
pylab.figure(n);
ax1 = pylab.subplot2grid((4,1),(0,0),rowspan=2)
aba.plotFishXPosition(e_RTn[n])
ax = pylab.subplot2grid((4,1),(2,0),sharex=ax1)
[sh,s,d,dsh,ds,st] = aba.getSidePreference(e_RTn[n], cond=[3,4], refState='On')
pylab.bar(st,np.array(s)/np.array(d),width=d)
pylab.axhline(.5,color='k')
pylab.ylabel('% time side1')
ax.set_yticks((0,.5,1))
ax = pylab.subplot2grid((4,1),(3,0),sharex=ax1);
pylab.bar(st,np.array(ds),width=d);
pylab.ylim((0,48))
pylab.axhline(24,color='k')
ax.set_yticks((0,24,48))
pylab.ylabel('Avg dist from\nside 1 (mm)')
pylab.xlabel('Time (s)')
pylab.show()
for y in range(len(c_RT)):
