def calculate_power_hit_miss_locations(prm, data, channel):
hit_data, miss_data = vr_split_data.split_all_data_hit_miss(prm, data)
hit_data_ob, hit_data_rz, hit_data_hb = vr_track_location_analysis.split_locations(
prm, hit_data)
miss_data_ob, miss_data_rz, miss_data_hb = vr_track_location_analysis.split_locations(
prm, miss_data)
#plot power spectra
#hit_data_theta = vr_track_location_analysis.add_theta_and_gamma_signal(prm, hit_data, hit_data)
#miss_data_theta = vr_track_location_analysis.add_theta_and_gamma_signal(prm, miss_data, miss_data)
#vr_track_location_plots.plot_power_spectrum_all_hitmiss(prm, hit_data_theta,miss_data_theta, channel)
#vr_track_location_plots.plot_power_spectrum_all_hitmiss_gamma(prm, hit_data[:,5],miss_data[:,5], channel)
#average area under curve etc for all stops
hit_power_ob = get_gamma_per_window(prm, hit_data_ob[:, 5])
hit_power_rz = get_gamma_per_window(prm, hit_data_rz[:, 5])
hit_power_hb = get_gamma_per_window(prm, hit_data_hb[:, 5])
miss_power_ob = get_gamma_per_window(prm, miss_data_ob[:, 5])
miss_power_rz = get_gamma_per_window(prm, miss_data_rz[:, 5])
miss_power_hb = get_gamma_per_window(prm, miss_data_hb[:, 5])
vr_power_calculations_plots.plot_gamma_power_hit_miss_locations(
prm, hit_power_ob, miss_power_ob, hit_power_rz, miss_power_rz,
hit_power_hb, miss_power_hb, channel)
hit_power_ob = get_gamma_per_window(prm, hit_data_ob[:, 4])
hit_power_rz = get_gamma_per_window(prm, hit_data_rz[:, 4])
hit_power_hb = get_gamma_per_window(prm, hit_data_hb[:, 4])
miss_power_ob = get_gamma_per_window(prm, miss_data_ob[:, 4])
miss_power_rz = get_gamma_per_window(prm, miss_data_rz[:, 4])
miss_power_hb = get_gamma_per_window(prm, miss_data_hb[:, 4])
vr_power_calculations_plots.plot_theta_power_hit_miss_locations(
prm, hit_power_ob, miss_power_ob, hit_power_rz, miss_power_rz,
hit_power_hb, miss_power_hb, channel)
def process_movement_and_stationary_locations(prm, before_stop, after_stop,
channel):
#before/after power spectra at different locations
#split based on location
before_stop_outbound, before_stop_rewardzone, before_stop_homebound = vr_track_location_analysis.split_locations(
prm, before_stop)
after_stop_outbound, after_stop_rewardzone, after_stop_homebound = vr_track_location_analysis.split_locations(
prm, after_stop)
# make example power spectras
vr_track_location_analysis.stop_start_power_spectra_locations(
prm, before_stop_outbound, before_stop_rewardzone,
before_stop_homebound, after_stop_outbound, after_stop_rewardzone,
after_stop_homebound, channel) # FIGURE 2 B :
# pull data from whole recording
vr_power_calculations.calculate_and_plot_pooled_gamma_power_stop_start_locations(
prm, after_stop_outbound, after_stop_rewardzone, after_stop_homebound,
before_stop_outbound, before_stop_rewardzone, before_stop_homebound,
channel)
# plot example continuous data
vr_track_location_plots.make_stop_start_locations_continuous_plot(
prm, channel, after_stop_outbound, after_stop_rewardzone,
after_stop_homebound, before_stop_outbound, before_stop_rewardzone,
before_stop_homebound)
return before_stop_outbound, before_stop_rewardzone, before_stop_homebound, after_stop_outbound, after_stop_rewardzone, after_stop_homebound
def calculate_power_speed_hit_miss_locations(prm, data, channel):
hit_data, miss_data = vr_split_data.split_all_data_hit_miss(prm, data)
hit_data_ob, hit_data_rz, hit_data_hb = vr_track_location_analysis.split_locations(
prm, hit_data)
miss_data_ob, miss_data_rz, miss_data_hb = vr_track_location_analysis.split_locations(
prm, miss_data)
speed_hit_ob = get_average_speed(prm, hit_data_ob[:, 6])
speed_hit_rz = get_average_speed(prm, hit_data_rz[:, 6])
speed_hit_hb = get_average_speed(prm, hit_data_hb[:, 6])
speed_miss_ob = get_average_speed(prm, miss_data_ob[:, 6])
speed_miss_rz = get_average_speed(prm, miss_data_rz[:, 6])
speed_miss_hb = get_average_speed(prm, miss_data_hb[:, 6])
gamma_power_hit_ob = get_gamma_per_window(prm, hit_data_ob[:, 5])
gamma_power_hit_rz = get_gamma_per_window(prm, hit_data_rz[:, 5])
gamma_power_hit_hb = get_gamma_per_window(prm, hit_data_hb[:, 5])
gamma_power_miss_ob = get_gamma_per_window(prm, miss_data_ob[:, 5])
gamma_power_miss_rz = get_gamma_per_window(prm, miss_data_rz[:, 5])
gamma_power_miss_hb = get_gamma_per_window(prm, miss_data_hb[:, 5])
theta_power_hit_ob = get_gamma_per_window(prm, hit_data_ob[:, 4])
theta_power_hit_rz = get_gamma_per_window(prm, hit_data_rz[:, 4])
theta_power_hit_hb = get_gamma_per_window(prm, hit_data_hb[:, 4])
theta_power_miss_ob = get_gamma_per_window(prm, miss_data_ob[:, 4])
theta_power_miss_rz = get_gamma_per_window(prm, miss_data_rz[:, 4])
theta_power_miss_hb = get_gamma_per_window(prm, miss_data_hb[:, 4])
#plot gamma
hit_data_ob = np.vstack((speed_hit_ob, gamma_power_hit_ob))
hit_data_ob = np.transpose(hit_data_ob)
miss_data_ob = np.vstack((speed_miss_ob, gamma_power_miss_ob))
miss_data_ob = np.transpose(miss_data_ob)
hit_data_rz = np.vstack((speed_hit_rz, gamma_power_hit_rz))
hit_data_rz = np.transpose(hit_data_rz)
miss_data_rz = np.vstack((speed_miss_rz, gamma_power_miss_rz))
miss_data_rz = np.transpose(miss_data_rz)
hit_data_hb = np.vstack((speed_hit_hb, gamma_power_hit_hb))
hit_data_hb = np.transpose(hit_data_hb)
miss_data_hb = np.vstack((speed_miss_hb, gamma_power_miss_hb))
miss_data_hb = np.transpose(miss_data_hb)
bins_hit_ob, speed_bins_hit_ob = bin_speed_data(prm, hit_data_ob)
bins_hit_rz, speed_bins_hit_rz = bin_speed_data(prm, hit_data_rz)
bins_hit_hb, speed_bins_hit_hb = bin_speed_data(prm, hit_data_hb)
bins_miss_ob, speed_bins_miss_ob = bin_speed_data(prm, miss_data_ob)
bins_miss_rz, speed_bins_miss_rz = bin_speed_data(prm, miss_data_rz)
bins_miss_hb, speed_bins_miss_hb = bin_speed_data(prm, miss_data_hb)
vr_power_calculations_plots.plot_gamma_power_speed_hit_miss_locations(
prm, bins_hit_ob, speed_bins_hit_ob, bins_miss_ob, speed_bins_miss_ob,
bins_hit_rz, speed_bins_hit_rz, bins_miss_rz, speed_bins_miss_rz,
bins_hit_hb, speed_bins_hit_hb, bins_miss_hb, speed_bins_miss_hb,
channel)