Closed
Description
The docs seem to suggest the function could do gamma convolution.
I set the attribute at hand accrodingly:
level1design = pe.Node(interface=Level1Design(), name="level1design")
level1design.inputs.bases = {'gamma': {'derivs':False}}
And sure enough I got the following in my node's _report/report.rst file:
Original Inputs
---------------
* bases : {'gamma': {'derivs': False}}
* contrasts : [('allStim', 'T', ['s1', 's2', 's3', 's4', 's5', 's6'], [1.0, 1.0, 1.0, 1.0, 1.0, 1.0])]
* ignore_exception : False
* interscan_interval : 1.0
* model_serial_correlations : True
* session_info : [{'regress': [], 'cond': [{'duration': [30.0], 'onset': [172.77000000000001], 'name': 's1'}, {'duration': [30.0], 'onset': [322.76999999999998], 'name': 's2'}, {'duration': [30.0], 'onset': [472.76999999999998], 'name': 's3'}, {'duration': [30.0], 'onset': [622.76999999999998], 'name': 's4'}, {'duration': [30.0], 'onset': [772.76999999999998], 'name': 's5'}, {'duration': [30.0], 'onset': [922.76999999999998], 'name': 's6'}], 'hpf': 180.0, 'scans': '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/specify_model/corr_10_trans_filt.nii.gz'}]
Execution Inputs
----------------
* bases : {'gamma': {'derivs': False}}
* contrasts : [('allStim', 'T', ['s1', 's2', 's3', 's4', 's5', 's6'], [1.0, 1.0, 1.0, 1.0, 1.0, 1.0])]
* ignore_exception : False
* interscan_interval : 1.0
* model_serial_correlations : True
* session_info : [{'regress': [], 'cond': [{'duration': [30.0], 'onset': [172.77000000000001], 'name': 's1'}, {'duration': [30.0], 'onset': [322.76999999999998], 'name': 's2'}, {'duration': [30.0], 'onset': [472.76999999999998], 'name': 's3'}, {'duration': [30.0], 'onset': [622.76999999999998], 'name': 's4'}, {'duration': [30.0], 'onset': [772.76999999999998], 'name': 's5'}, {'duration': [30.0], 'onset': [922.76999999999998], 'name': 's6'}], 'hpf': 180.0, 'scans': '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/specify_model/corr_10_trans_filt.nii.gz'}]
Execution Outputs
-----------------
* ev_files : ['/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/ev_s1_0_1.txt', '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/ev_s2_0_2.txt', '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/ev_s3_0_3.txt', '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/ev_s4_0_4.txt', '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/ev_s5_0_5.txt', '/mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/ev_s6_0_6.txt']
* fsf_files : /mnt/data/NIdata/ofM.erc/GLM/level1_ext_gamma/first_level/_condition_ERC_ofM_subject_5503/_scan_type_EPI_CBV_alej/level1design/run0.fsf
So far so good.
Now, if I try to browse to the produced .fsf file, I see the EVs still use "dgamma" convolution, e.g.:
# EV title
set fmri(evtitle6) "s6"
# Basic waveform shape (EV 6)
#0 : Square
#1 : Sinusoid
#2 : Custom (1 entry per volume)
#3 : Custom (3 column format)
#4 : Interaction
#10 : Empty (all zeros)
set fmri(shape6) 3
# Convolution (EV 6)
#0 : None
#1 : Gaussian
#2 : Gamma
#3 : Double-Gamma HRF
#4 : Gamma basis functions
#5 : Sine basis functions
#6 : FIR basis functions
set fmri(convolve6) 3
I have tried to browse the sources to figure out what's wrong, and I found that the file needed for the gamma option exists, but is not used. Apprently the dgamma file is added by default unless the convolution is set to "none".