I am trying to perform the population level analysis and need help understanding how to do the trial level modeling (without RT). I went through your code. I am not familiar with R and I am unable to understand the input data structure.

Can you help me understand the input data structure of the task including the how the header detail need to be provided?

And is it also possible for you to upload the input files used in the population analysis for the above mentioned paper (file name : 'RT-task-data.txt' or ‘task-data.txt’)

The paper you referred to aims to explore the extent and impact of cross-trial variability in a typical fMRI dataset. To better understand your study, could you clarify the following:

Are you interested in cross-trial variability?

What is your research hypothesis for the 2 × 3 experimental design?

How many trials are there per condition?

Are you planning to perform the population-level analysis at the voxel level or the region level?

I am working with Srikanth Padmala to understand the cross-trial variability in our dataset.

To answer your question:

Yes, we feel cross trial variability is something good to look at in our dataset.

So, we are interested in understanding the interaction in our 2 × 3 experimental design. We have 2 Reward (High, Low) x 3 Emotion (Positive, Neutral, Negative).

I have 24 trials per condition. Is that too low? But I have 60 subject data.

For now, I am looking at performing the population-level analysis at the ROI level

The header for the ROI data should look like this:

subjID ROI Reward Emotion trial beta SE

Please note that the estimates for the BOLD response might be significantly biased (e.g., underestimated) because the adopted HRF is likely not well accommodative for some regions.

You can add -tout in your 3dREMLfit script to obtain the t-statistics that are associated with the estimated regression coefficients from -Rbuck. Then the standard errors can be obtained through the following formula SE=\frac{\widehat\beta}{t}.

Can I perform the same REMLfit at the voxel level for the wholeBrain ?

Sure. However, it would be computationally challenging if you plan to implement the Bayesian multilevel modeling approach at the voxel level.

Gang Chen

The
National Institute of Mental Health (NIMH) is part of the National Institutes of
Health (NIH), a component of the U.S. Department of Health and Human
Services.