create_afni_timing_faces_over_time.py

# Purpose: take teh first and last trials of the runs
# there's 3 trials/run - get only first 2 trials, last 2 trials

import os
import glob
from shutil import copyfile
import pandas as pd
import json
import numpy as np

save_path = '/data/jux/cnds/amennen/rtAttenPenn/fmridata/Nifti/derivatives/afni/first_level/timing_files';
all_categories = ['fearful','happy', 'neutral']
subjects = np.array([1,2,3,4,5,6,7,8,9,10,11,101, 102,103,104,105,106, 107,108,109,110,111,112,113,114])
days=[1,3]
for subjectNum in subjects:
    bids_id = 'sub-{0:03d}'.format(subjectNum)
    for d in days:
        subjectDay=d
        ses_id = 'ses-{0:02d}'.format(subjectDay)
        # then for each category, we want the first, last, and 1 TR after
        for c in np.arange(len(all_categories)):
            category = all_categories[c]
            category_str = category + '.txt'
            file_name = os.path.join(save_path,bids_id,ses_id, category_str)
            t = pd.read_fwf(file_name,header=None)
            timing = t.values # now 2 x 18 array
            # all first values are every 6
            ntrials = 6
            nhalves = 2
            # first half: colums 0 & 1
            # last half: columns 4 & 5
            #timing[0,trial] timing[0,trial+ntrials] timing[0,trial+(ntrials*2)]
            nreps = 3 # each stimulus gets shown in 3 blocks
            counter = 0
            for r in np.arange(nreps):
                for h in np.arange(nhalves):

                    trial_1 = (h*4) + (ntrials*r)
                    trial_2 = trial_1 + 1
                    print(trial_1)
                    print(trial_2)
                    category_str_out = category + '_' + 'rep_' + str(r) + '_' + 'half' + '_' + str(h) + '.txt'
                    print(category_str_out)
                    output_file = os.path.join(save_path,bids_id,ses_id,category_str_out)
                    outF = open(output_file,"w")
                    outF.write('%8.4f %8.4f\n' % (timing[0,trial_1],timing[0,trial_2]))
                    outF.write('%8.4f %8.4f' % (timing[1,trial_1],timing[1,trial_2]))
                    outF.close()
                    counter+=1