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 08, 2018 12:04PM
When you say "many runs", how many is "many"? The reason I ask is that it sounds like each run has one movie-event only, and averaging across events is part of the usual FMRI analysis strategy.

The answer to your specific question depends on the question you are asking about the time series shape during the movie interval, versus the pre- and post-movie rest intervals. If all you are asking is "the size of the average BOLD effect during the movie", then 'BLOCK(300,1)" is a reasonable choice. But that does not mean you are expecting the BLOCK shape to fit the data well -- it is just that the beta coefficient for the BLOCK shape for the 300 s movie interval will provide an estimate of the average signal change (up or down) during that interval. If you want a better fit, you could use 'TENTzero(0,310,32)' to provide an undulating fit -- in this case, you would get 30 betas, and then you would have to average them to get the average response -- which should be close to the result you would get from 'BLOCK(300,1)' -- this averaging could be done with the proper 3dTstat command, if you parse out the details of the format of the output statistics dataset using 3dinfo. Negative average responses would correspond to "deactivation".

Other details that occur to me, if you are using afni_proc.py:
  • Use option -regress_motion_per_run if you are analyzing all runs for one subject in a single afni_proc.py command.
  • Alternatively, analyze each run in a separate afni_proc.py run, and then combine the results later somehow.
I have never analyzed datasets like this, so to some extent you'll have to feel your way through the analysis pathways.
Subject Author Posted

basis functions for movie data

insularcortex May 07, 2018 05:02PM

Re: basis functions for movie data

RWCox May 08, 2018 12:04PM