Hi Michael,
I don't really understand what you are trying to do with AlphaSim.
Note that the purpose of AlphaSim is to define a minimal cluster
size to achieve a corrected p-value, given an input grid, mask,
blur widths and uncorrected p-value. Typically, the only dataset
input to AlphaSim is a mask (which would define the voxel grid as
well as the mask itself which denotes the subset of voxels that
you want to consider as your valid world).
So after AlphaSim, one would almost assuredly run 3dmerge, as you
are doing. AlphaSim would not replace 3dmerge.
Anyway, maybe I am just missing the point of what you want...
- rick