[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)):