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  

|
May 31, 2017 03:48PM
Hello,

I'd like to create an impulse response model. My functional data are not continuous (e.g. sparse) and thus it is not appropriate to insert stimulus timing files. For example, my timing files begin at the beginning of the experiment, but the functional data are not continuous. So if I use TENT or CSPLIN then it would take the timing file and generate my impulses accordingly, except that it wouldn't be in the same temporal domain as the functional data (because the fMRI data isn't continuous).

So I have my own impulse response matrix that I would like to provide where columns refer to an event type and rows refer to a specific volume.

For example, here is a sample of 6 rows:

1 0 0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0 0 0 0 1

Here is how I am getting the matrix into 3dDeconvolve.
3dDeconvolve 												                \
			-input $runs -jobs 24										\
			-polort A 	 										                \
			-ortvec MoPar_demean.1D 	demean 						        \
			-ortvec MoPar_derv.1D 	        derv 						        \
			-ortvec FIR.1D 		        ir	# impulse response matrix	        \
			-censor censor.1D 						                                        \
			-num_glt 1 											        \
			-gltsym 'SYM: ir[9..11] -ir[0..2]'  -glt_label 1 sound-nosound 			\

Obviously there are more deconvolve options for statistics and such, I just wanted to provide how I'm setting up the relevant parts.
I want to be able to submit some of the impulse response columns to a glt test. In the above example, lets say that I want to do a t test between columns 9 through 11 minus columns 0 through 2. How would I specify this?

If I specify it in the above way I receive an error:

** ERROR: row #1 of matrix 'SYM: ir{9]}' is all zero!
** ERROR: -gltsym errors immediately above from file 'SYM: ir{9]}'
** FATAL ERROR: Can't continue after the above -gltsym problems!

My rationale for this specification is because this is how I'd do it in TENT or CSPLIN. Thank you for your time.

Dustin
Subject Author Posted

On ortvec and glts

dmoracze May 31, 2017 03:48PM

Re: On ortvec and glts

rick reynolds June 02, 2017 09:18AM

Re: On ortvec and glts

dmoracze June 02, 2017 12:37PM

Re: On ortvec and glts

dmoracze June 02, 2017 02:17PM

Re: On ortvec and glts

rick reynolds June 02, 2017 03:14PM

Re: On ortvec and glts

rick reynolds June 06, 2017 11:07AM

Re: On ortvec and glts

dmoracze June 06, 2017 11:17AM

Re: On ortvec and glts

dmoracze June 14, 2017 01:49PM

Re: On ortvec and glts

dmoracze June 15, 2017 04:04PM

Re: On ortvec and glts

rick reynolds June 20, 2017 05:12PM

Re: On ortvec and glts

rick reynolds June 21, 2017 11:51AM

Re: On ortvec and glts Attachments

dmoracze June 21, 2017 04:44PM

Re: On ortvec and glts

rick reynolds June 28, 2017 08:57AM

Re: On ortvec and glts Attachments

dmoracze July 14, 2017 11:50AM