I have questions about the strange beta values when I apply 3dDeconvolve to my dataset.
The conditions for the script I tried are as follows:
As I censored the motion in preprocessing step, I used 1d_tool.py with the option '-censor_motion 0.4(“motion_threshold” from here) ’ and then used the output to the 3dDeconvole with the option ‘-censor’.
Used options before setting ‘-num_stimts’ are ‘-censor’, ‘-mask’, ‘-polort A -float’ and ‘-allzero_OK’.
The number of stimulus files is 114 with the option ‘-num_stimts 114’.
3-1. 108 out of 114 are text files in which stimulation onset timing is represented (ex. 1st file has the value 4.03, 2nd file has the value 10.04 in second, and so on. The number of trials of our experiment is 108). Therefore, the statements such as ‘-stim_times 1 1st_file $basis -stim_label 1 Int1’ are repeated from number 1 to number 108. $basis is ‘BLOCK(2,1)’.
3-2. 6 out of 114 are from motion_demean.~.1D file. They are roll, pitch, yaw, dS, dL and dP.
And then, options for output file are followed; ‘-fout’, ‘-tout’, ‘-x1D’, ‘-xjpeg’, ‘-x1D_uncensored’ and ‘-bucket’.
However, beta values for the trials after trying above 3dDeconvolve are strange; beta values for some trials are abnormally large.
For instance, maximum beta values for trials 46-49 and 93-96 are about 10^6 and 10^7, respectively. The order of the minimum beta values is similar.
I thought this might be from subject’s considerable movement and the censored TR, so I tried 1d_tool.py with ‘-censor_motion 0.5’ and then 3dDeconvolve again with newly got censor file. This time, beta values for trials 46-49 are lowered to about 10^5, but those for trials 93-96 is still similar (about 10^7).
I don’t know why those strange beta values are calculated. Could you make suggestions for me to fix this problem? I’d be appreciated to any explanation or suggestions. Thank you.