How do you get those local maxima? I don't know the appropriate procedure
I think they are reported by default in the feat web output pages but they can also be generated using the cluster command. If I remember correctly, the flag is --olmax. It should be clear from help screen if you type cluster without any arguments.
I don't mean to be critical but if I were a reviewer, I wouldn't be entirely happy with the method you describe. The problem is that the AAL brain regions were determined on the macroanatomy of one brain (Colin27) and therefore are not representative of any other brain, let alone a group sample. Tzourio-Mazoyer and colleagues acknowledge this explicitly in their paper. So the number (or percentage) of voxels in each of the AAL regions is not accurate and the precision of the reporting hides this fact. In addition, if not all AAL regions were imaged (or some where only partially imaged), this too may mislead the reader regarding the the amount of activation in any given region.
Wouldn't it be more appropriate to describe the activations yourself, based on their extent into various macroanatomic regions seen in the group mean structural?
- Joe