def plot_thresh_zstat(self): """Plot the thresholded zstat.""" m = Mosaic(stat=self.inputs.zstat_thresh_file, mask=self.inputs.mask_file) m.plot_activation(pos_cmap="OrRd_r", vfloor=3.3, alpha=.9) out_fname = nii_to_png(self.inputs.zstat_thresh_file) self.out_files.append(out_fname) m.savefig(out_fname) m.close()
def plot_zstats(self): """Plot the positive and negative z stats with a low threshold.""" for z_file in self.inputs.zstat_files: m = Mosaic(self.mean, z_file, self.mask, step=1) m.plot_activation(pos_cmap="Reds_r", neg_cmap="Blues", thresh=1.7, alpha=.85) png_name = nii_to_png(z_file) m.savefig(png_name) m.close() self.out_files.append(png_name)
def write_zstat_images(self): """Show the fixed effects inference maps.""" for fname in self.inputs.zstat_files: m = Mosaic(self.inputs.anatomy, fname, self.every_mask, step=self.step) m.plot_activation(pos_cmap="Reds_r", neg_cmap="Blues", thresh=2.3, alpha=.85) contrast = fname.split("/")[-2] os.mkdir(contrast) out_fname = op.join(contrast, "zstat1.png") self.zstat_files.append(op.abspath(contrast)) m.savefig(out_fname) m.close()