def simulated_test(units, stimulus_list):
# assert len(next(iter(units.values())).spike_train)==2200
test_pipeline = glia.compose(glia.f_create_experiments(stimulus_list),
glia.f_has_stimulus_type(["GRATING"]),
glia.f_group_by_stimulus(),
glia.f_calculate_firing_rate_by_stimulus())
firing_rates = glia.apply_pipeline(test_pipeline, units, progress=True)
for stimulus, rates in next(iter(firing_rates.values())).items():
for rate in rates:
assert np.isclose(rate, 60, 1)
def simulated_test(units, stimulus_list): assert len(next(iter(units.values())).spike_train)==2200 test_pipeline = glia.compose( glia.f_create_experiments(stimulus_list), glia.f_has_stimulus_type(["GRATING"]), glia.f_group_by_stimulus(), glia.f_calculate_firing_rate_by_stimulus() ) firing_rates = glia.apply_pipeline(test_pipeline, units) for stimulus,rates in next(iter(firing_rates.values())).items(): for rate in rates: assert rate==1
def get_fr_dsi_osi(units, stimulus_list):
get_bar_firing_rate = glia.compose(
glia.f_create_experiments(stimulus_list),
glia.f_has_stimulus_type(["BAR"]),
glia.f_group_by_stimulus(),
glia.f_calculate_firing_rate_by_stimulus(),
)
bar_firing_rate = glia.apply_pipeline(get_bar_firing_rate,
units,
progress=True)
get_bar_dsi = glia.compose(glia.by_speed_width_then_angle,
glia.calculate_dsi_by_speed_width)
bar_dsi = glia.apply_pipeline(get_bar_dsi, bar_firing_rate, progress=True)
get_bar_osi = glia.compose(glia.by_speed_width_then_angle,
glia.calculate_osi_by_speed_width)
bar_osi = glia.apply_pipeline(get_bar_osi, bar_firing_rate, progress=True)
return (bar_firing_rate, bar_dsi, bar_osi)
def get_fr_dsi_osi(units, stimulus_list):
get_bar_firing_rate = glia.compose(
glia.f_create_experiments(stimulus_list),
glia.f_has_stimulus_type(["BAR"]),
glia.f_group_by_stimulus(),
glia.f_calculate_firing_rate_by_stimulus(),
)
bar_firing_rate = glia.apply_pipeline(get_bar_firing_rate,units)
get_bar_dsi = glia.compose(
glia.by_speed_width_then_angle,
glia.calculate_dsi_by_speed_width
)
bar_dsi = glia.apply_pipeline(get_bar_dsi,bar_firing_rate)
get_bar_osi = glia.compose(
glia.by_speed_width_then_angle,
glia.calculate_osi_by_speed_width
)
bar_osi = glia.apply_pipeline(get_bar_osi,bar_firing_rate)
return (bar_firing_rate, bar_dsi, bar_osi)
