Use NULL (default) to just return the template
objects directly. Otherwise, use a character vector of length 3 or 4 of file
path(s) to save the output to:
the mean template, the variance template, the variance decomposition, and
the FC template if present, in that order. If one file name is provided,
it will be appended with
"_mean.[file_ext]" for the template mean map,
"_var.[file_ext]" for the template variance map,
"_varDecomp.rds" for the variance decomposition, and
"_FC.rds" where [file_ext]
will be "dscalar.nii" for CIFTI input, "nii" for NIFTI input,
and "rds" for data input.