) args = parser.parse_args() # --------------------------------------------------------------------------- # Process argv # --------------------------------------------------------------------------- # Replace this with good arg processing.... # basename, dataname, rois, cond, tr, filtfile = process_exp_argv(sys.argv) data = get_data(args.data) _, rois = load_roifile(args.roifile) ## roifile # --------------------------------------------------------------------------- # Setup exp # --------------------------------------------------------------------------- spacetime = SelectSpace(PCA(6, whiten=True), fir, mode="decompose") exp = DecomposeExp( spacetime, data, window=args.window, nsig=3, tr=args.tr ) # --------------------------------------------------------------------------- # And run each roi # --------------------------------------------------------------------------- for n, roi in enumerate(rois): print("{3}: {0} ({1}/{2})".format(roi, n+1, len(rois), args.data)) exp.run( args.name, roi, args.cond, smooth=False, filtfile=args.filtfile, event=True )
from fmrilearn.analysis import fir from sklearn.decomposition import PCA from wheelerexp.base import Space from wheelerexp.base import DecomposeExp from wheelerexp.common import process_exp_argv from wheelerdata.load.meta import get_data # --------------------------------------------------------------------------- # Process argv # --------------------------------------------------------------------------- basename, dataname, rois, cond, tr, filtfile = process_exp_argv(sys.argv) data = get_data(dataname) # --------------------------------------------------------------------------- # Setup exp # --------------------------------------------------------------------------- spacetime = Space(PCA(5, whiten=True), fir, mode="decompose") exp = DecomposeExp(spacetime, data, window=15, nsig=3, tr=tr) # --------------------------------------------------------------------------- # And run each roi # --------------------------------------------------------------------------- for n, roi in enumerate(rois): print("{3}: {0} ({1}/{2})".format(roi, n+1, len(rois), dataname)) exp.run(basename, roi, cond, smooth=False, filtfile=filtfile, event=True)