Hi Gang,

Many thanks for the offer. I'm asking since we're starting a new project and I need to decide whether the group level analysis will be done in FSL, pure R or an R/Afni tool (of the sort you're developing). I can certainly wait till it's out.

A few questions: will you be using T values estimated by 3dREMLfit? Also, I am curious to know if you are concerned about heteroskedasticity between beta and residuals on

(a) the single subject level (i.e., in subject Joe's brain, higher betas are associated with larger var(residuals) -- I've seen this in a few cases, and it's interesting [sometimes]) and/or

(b) heteroskedasticity on the group level -- voxels with higher mean beta on the group level being associated with higher var(beta) or higher var(residuals) in that voxel

In any case, if you have a reference to the approach you're implementing I'd appreciate a pointer.

All the best, and thank you for your assistance.
Oori