AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
October 01, 2019 05:33AM
Hi AFNI experts,

I am a very new AFNI user, and I'm trying to estimate my betas through the 3dLSS method. The problem is that in my experimental design (2x2x2x2, 192 total trials distributed in 12 separated runs, 2 trials per condition in each run), the number of stimulus times for a given condition is maximum 2 for run, and can be 1 in some cases (due to participants errors during the task), and also there are some regressor with always 1 stim times (a cue with the block type at the beginning of each run). This creates an error when I'm running the 3dLSS:
"FATAL ERROR: Matrix attribute BasisColumns_000.. shows only 1 column for stim_time_IM. 3dLSS is meant to be used when more than one stimulus time was given, and then it computes the response beta or each stim time separately. If you have only one stim time with -stim_times_IM, you can use the output dataset from 3dDeconvolve (or 3dREMLfit) to get that single beta directly".

Does this mean that I need to run the 3dDeconvolve to estimate the betas of the stim_times with only 1 stimulus time, and after, run the 3dLSS to estimate the betas of the stim times with more than 1 stimulus time? Because this seems complicated, and I wonder if it can be done in some other way.

I attach my script next, in case it is useful or the is any problem somewhere else:

3dDeconvolve -input afniaufIMPRECs01+tlrc.HEAD \
-polort A \
-local_times \
-num_stimts 212 \
-concat '1D: 0 141 282 423 564 705 846 987 1128 1269 1410 1551' \
-stim_times_IM 1 Run1_BlockCue_Impl.txt 'GAM(8.6,0.547)' -stim_label 1 Run1_BlockCue_Impl \
-stim_times_IM 2 Run1_Errors_Grid.txt 'GAM(8.6,0.547)' -stim_label 2 Run1_Errors_Grid \
-stim_times_IM 3 Run1_Errors_Inst.txt 'GAM(8.6,0.547)' -stim_label 3 Run1_Errors_Inst \
-stim_times_IM 4 Run1_Grid_Impl_Faces_1Dim_1Resp.txt 'GAM(8.6,0.547)' -stim_label 4 Run1_Grid_Impl_Faces_1Dim_1Resp \
-stim_times_IM 5 Run1_Grid_Impl_Faces_1Dim_2Resp.txt 'GAM(8.6,0.547)' -stim_label 5 Run1_Grid_Impl_Faces_1Dim_2Resp \
-stim_times_IM 6 Run1_Grid_Impl_Faces_2Dim_1Resp.txt 'GAM(8.6,0.547)' -stim_label 6 Run1_Grid_Impl_Faces_2Dim_1Resp \
-stim_times_IM 7 Run1_Grid_Impl_Faces_2Dim_2Resp.txt 'GAM(8.6,0.547)' -stim_label 7 Run1_Grid_Impl_Faces_2Dim_2Resp \
-stim_times_IM 8 Run1_Grid_Impl_Food_1Dim_1Resp.txt 'GAM(8.6,0.547)' -stim_label 8 Run1_Grid_Impl_Food_1Dim_1Resp \
-stim_times_IM 9 Run1_Grid_Impl_Food_1Dim_2Resp.txt 'GAM(8.6,0.547)' -stim_label 9 Run1_Grid_Impl_Food_1Dim_2Resp \
-stim_times_IM 10 Run1_Grid_Impl_Food_2Dim_1Resp.txt 'GAM(8.6,0.547)' -stim_label 10 Run1_Grid_Impl_Food_2Dim_1Resp \
-stim_times_IM 11 Run1_Grid_Impl_Food_2Dim_2Resp.txt 'GAM(8.6,0.547)' -stim_label 11 Run1_Grid_Impl_Food_2Dim_2Resp \
-stim_times_IM 12 Run1_Inst_Impl_Faces_1Dim_1Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 12 Run1_Inst_Impl_Faces_1Dim_1Resp \
-stim_times_IM 13 Run1_Inst_Impl_Faces_1Dim_2Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 13 Run1_Inst_Impl_Faces_1Dim_2Resp \
-stim_times_IM 14 Run1_Inst_Impl_Faces_2Dim_1Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 14 Run1_Inst_Impl_Faces_2Dim_1Resp \
-stim_times_IM 15 Run1_Inst_Impl_Faces_2Dim_2Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 15 Run1_Inst_Impl_Faces_2Dim_2Resp \
-stim_times_IM 16 Run1_Inst_Impl_Food_1Dim_1Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 16 Run1_Inst_Impl_Food_1Dim_1Resp \
-stim_times_IM 17 Run1_Inst_Impl_Food_1Dim_2Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 17 Run1_Inst_Impl_Food_1Dim_2Resp \
-stim_times_IM 18 Run1_Inst_Impl_Food_2Dim_1Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 18 Run1_Inst_Impl_Food_2Dim_1Resp \
-stim_times_IM 19 Run1_Inst_Impl_Food_2Dim_2Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 19 Run1_Inst_Impl_Food_2Dim_2Resp \
-stim_times_IM 20 Run2_BlockCue_Mem.txt 'GAM(8.6,0.547)' -stim_label 20 Run2_BlockCue_Mem \
-stim_times_IM 21 Run2_Errors_Grid.txt 'GAM(8.6,0.547)' -stim_label 21 Run2_Errors_Grid \
-stim_times_IM 22 Run2_Errors_Inst.txt 'GAM(8.6,0.547)' -stim_label 22 Run2_Errors_Inst \
.
.(the script is very long, so I omitted some for readability)
.
-stim_times_IM 212 Run12_Inst_Mem_Food_2Dim_2Resp.txt 'GAM(8.6,0.547, 2.5)' -stim_label 212 Run12_Inst_Mem_Food_2Dim_2Resp \
-ortvec MovReg.txt'[1..6]' 'MovReg'\
-x1D decon_3ddeconvolve_1XD.x1D \
-x1D_stop \
-xsave \
-bucket post_3ddeconvolve_aufIMPRECs01 \
-allzero_OK \
-GOFORIT 800 \

3dLSS -verb -matrix decon_3ddeconvolve_1XD.x1D \
-input afniaufIMPRECs01+tlrc.HEAD \

Thank you in advance,
Alberto
Subject Author Posted

Problem with 3dLSS after 3dDeconvolve

asobrado October 01, 2019 05:33AM

Re: Problem with 3dLSS after 3dDeconvolve

gang October 01, 2019 04:38PM

Re: Problem with 3dLSS after 3dDeconvolve

asobrado October 04, 2019 11:55AM

Re: Problem with 3dLSS after 3dDeconvolve

gang October 04, 2019 04:30PM

Re: Problem with 3dLSS after 3dDeconvolve

asobrado October 06, 2019 10:26AM

Re: Problem with 3dLSS after 3dDeconvolve

rick reynolds October 07, 2019 05:07PM

Re: Problem with 3dLSS after 3dDeconvolve

asobrado October 08, 2019 04:48AM

Re: Problem with 3dLSS after 3dDeconvolve

rick reynolds October 08, 2019 09:15AM

Re: Problem with 3dLSS after 3dDeconvolve

asobrado October 08, 2019 09:50AM

Re: Problem with 3dLSS after 3dDeconvolve

tommorin November 07, 2019 01:20PM

Re: Problem with 3dLSS after 3dDeconvolve

gang November 07, 2019 05:21PM