**--in <file>**	input data file ( NIFTI for volumetric, **GIFTI for surface** )
Optional arguments (You may optionally specify one or more of):
	--ac		perform autocorrelation estimation only
	--thr <num>	initial threshold to apply to input data
	--ar		fits autoregressive model - default is to use tukey with M=sqrt(numvols)
	--help		prints this message
	--noest		do not estimate auto corrs
	--outputPWdata		output prewhitened data and average design matrix
	--pava		estimates autocorr using PAVA - default is to use tukey with M=sqrt(numvols)
	--sa		smooths auto corr estimates
	-v		outputs full data
	--rn <file>	directory name to store results in, default is results
	--mode <mode>	analysis mode, options are volumetric ( default ) or surface. Caution: surface-based functionality is still BETA
	**--in2 <file>**	input surface for autocorr smoothing in surface-based analyses
	--mf <file>	re-estimate mean_func baseline - for use with perfusion subtraction
	--mft <file>	minimum timepoint file
	--pd <file>	paradigm file
	--con <file>	t-contrasts file
	--fcon <file>	f-contrasts file
	**--epith <num>**	susan brightness threshold for volumetric analysis/smoothing sigma for surface analysis
	**--ms <num>**	susan mask size for volumetric analysis/smoothing extent for surface analysis
	--tukey <num>	uses tukey window to estimate autocorr with window size num - default is to use tukey with M=sqrt(numvols)
	--mt <num>	uses multitapering with slepian tapers and num is the time-bandwidth product - default is to use tukey with M=sqrt(numvols)
	--ven		list of numbers indicating voxelwise EVs position in the design matrix (list order corresponds to files in vxf option). Caution BETA option, only use with volumetric analysis.
	--vef		list of 4D images containing voxelwise EVs (list order corresponds to numbers in vxl option). Caution BETA option, only use with volumetric analysis.