I concatenated two runs each having 147 time points and I was doing true deconvolution with stim_nptr=3. I just figured that the registration wasn't quite satisfactory here. There is obvious difference between the run1 and run2 after
registration. Did I do something wrong for the registration?
Thanks.
george
I guess it's easier to copy the command history to give you more info.
Here is what I did:
3dvolreg -pref SIRP_r1_vr -base SIRP_r1+orig'[90]' SIRP_r1+orig
3dvolreg -pref SIRP_r2_vr -base SIRP_r1+orig'[90]' SIRP_r2+orig
3dTcat -prefix SIRPAll_vr SIRP_r1_vr+orig'[2..148]' SIRP_r2_vr+orig'[2..148]'
3dTshift -prefix SIRPAll_vrTshift SIRPAll_vr+orig
waver -dt 1 -GAM -numout 882 -inline 72@0 37@1 124@0 36@1 172@0 \
187@0 36@1 126@0 36@1 56@0 > SIRP_arrow_1sec_hrf.1D
waver -dt 1 -GAM -numout 882 -inline 186@0 37@1 63@0 36@1 119@0 \
85@0 37@1 117@0 37@1 165@0 > SIRP_prob1_1sec_hrf.1D
waver -dt 1 -GAM -numout 882 -inline 26@0 37@1 280@0 37@1 61@0 \
144@0 37@1 223@0 37@1 > SIRP_prob3_1sec_hrf.1D
waver -dt 1 -GAM -numout 882 -inline 137@0 36@1 232@0 36@1 \
32@0 37@1 235@0 37@1 100@0 > SIRP_prob5_1sec_hrf.1D
waver -dt 1 -GAM -numout 882 -inline 1@0 12@1 50@0 6@1 40@0 8@1 57@0 3@1 46@0 \
6@1 40@0 8@1 45@0 6@1 52@0 4@1 57@0 \
1@0 12@1 56@0 6@1 47@0 8@1 51@0 4@1 38@0 \
6@1 47@0 8@1 57@0 6@1 38@0 4@1 52@0 > SIRP_cross_1sec_hrf.1D
waver -dt 1 -GAM -numout 882 -inline 16@0 8@1 97@0 13@1 47@0 3@1 76@0 3@1 50@0 \
8@1 47@0 14@1 39@0 \
16@0 14@1 49@0 3@1 52@0 8@1 92@0 3@1 50@0 \
15@1 92@0 8@1 39@0 > SIRP_encod_1sec_hrf.1D
echo 0 > SIRPconcat.1D
echo 147 >> SIRPconcat.1D
3dAutomask -prefix SIRPAll_vrTshift_mask SIRPAll_vrTshift+orig
3dDeconvolve -input SIRPAll_vrTshift+orig -xout -nfirst 0 -nlast 293 \
-concat SIRPconcat.1D -mask SIRPAll_vrTshift_mask+orig -num_stimts 6 \
-stim_file 1 SIRP_arrow_1sec_hrf.1D -stim_label 1 arrow -stim_maxlag 1 3 -stim_nptr 1 3 \
-stim_file 2 SIRP_prob1_1sec_hrf.1D -stim_label 2 prob1 -stim_maxlag 2 3 -stim_nptr 2 3 \
-stim_file 3 SIRP_prob3_1sec_hrf.1D -stim_label 3 prob3 -stim_maxlag 3 3 -stim_nptr 3 3 \
-stim_file 4 SIRP_prob5_1sec_hrf.1D -stim_label 4 prob5 -stim_maxlag 4 3 -stim_nptr 4 3 \
-stim_file 5 SIRP_cross_1sec_hrf.1D -stim_label 5 cross -stim_maxlag 5 3 -stim_nptr 5 3 \
-stim_file 6 SIRP_encod_1sec_hrf.1D -stim_label 6 encod -stim_maxlag 6 3 -stim_nptr 6 3 \
-full_first -fout -tout -rout -fitts SIRPAll_vrTshift_3ddeconv_fitts -bucket SIRPAll_vrTshift_3ddeconv