def test_make_regressor_1():
""" test the generated regressor
"""
condition = ([1, 20, 36.5], [2, 2, 2], [1, 1, 1])
frame_times = np.linspace(0, 69, 70)
hrf_model = 'spm'
reg, reg_names = compute_perfusion_regressors([condition], ['cond'],
hrf_model,
frame_times)
assert_equal(len(reg_names), 2)
assert_equal(reg_names[0], 'perfusion_baseline')
assert_equal(reg_names[1], 'perfusion_cond')
assert_almost_equal(np.sum(reg), -3, 1)
def test_make_regressor_3():
""" test the generated regressor
"""
condition = ([1, 20, 36.5], [0, 0, 0], [1, 1, 1])
frame_times = np.linspace(0, 138, 70)
hrf_model = 'fir'
reg, reg_names = compute_perfusion_regressors([condition], ['cond'],
hrf_model,
frame_times,
fir_delays=np.arange(4))
assert_equal(len(reg_names), 5)
assert_array_equal(np.unique(reg), np.array([-.5, 0, .5]))
assert_array_equal(np.sum(reg, 1), np.array([ 0., -1.5, 1.5, -1.5, 1.5]))
