Hi Gang,
Thanks for your reply.
Unfortunately the '3drefit ' command couldn't append the fico parameters in the header of the dataset. But I have another question with this approach. I used the following 3dLocalBistat to output the dataset with Pearson correlation,
3dLocalBistat -nbhd 'RHDD(6)' -stat pearson -mask MASK+orig -prefix OUT A+orig B+orig
Since the volume of the rhombic dodecahedron with radius of 6mm is 432mm3 and the voxel size of A+orig and B+orig is 1x1x1mm3, the number of voxel pairs used in correlation analysis at the edge of the mask will be <432. Using a universal number of voxel pairs to derive the p-values may have a problem.
Any suggestions?
Phoebe.