Hi Kaitlin,

You could always load these into Matlab using something like:

 filename=sprintf('pFDR_corr_Subject1_Condition001vsCondition002_Source%03d.nii',source);
 a = spm_vol(filename);
 pvalues = spm_read_vols(a);

but if you are mostly interested in displaying these results perhaps the simplest way would be to first threshold these maps, for example by adding in the original code something like the following (within the for loop):

 filename=sprintf('pFDR_thresholded_corr_Subject1_Condition001vsCondition002_Source%03d.nii',source);
 V=struct('mat',a(1).mat,'dim',a(1).dim,'fname',filename,'pinfo',[1;0;0],'n',[1,1],'dt',[spm_type('uint8') spm_platform('bigend')]);
 spm_write_vol(V,double(p<.05));

and then explore them using the following (or your preferred nifti visualizer)

 conn_mesh_display(filename,'');

or 

 conn_mesh_display('',filename);

to view those suprathreshold clusters using a surface-projected view or a standard 3d-volume view, respectively.

Hope this helps
Alfonso
Originally posted by Kaitlin Cassady: