14.2.3. Wardle et al. (2025). Brief encounters with real objects modulate medial parietal …¶
Introduction¶
Here we present commands used in the following paper:
- Wardle SG, Rispoli B, Roopchansingh V, Baker CI (2025). Brief encounters with real objects modulate medial parietal but not occipitotemporal cortex. (submitted) bioRxiv doi: https://doi.org/10.1101/2024.08.05.606667
Abstract: Humans are skilled at recognizing everyday objects from pictures, even if we have never encountered the depicted object in real life. But if we have encountered an object, how does that real-world experience affect the representation of its photographic image in the human brain? We developed a paradigm that involved brief real-world manual exploration of everyday objects prior to the measurement of brain activity with fMRI while viewing pictures of those objects. We found that while object-responsive regions in lateral occipital and ventral temporal cortex were visually driven and contained highly invariant representations of specific objects, those representations were not modulated by this brief real-world exploration. However, there was an effect of visual experience in object-responsive regions in the form of repetition suppression of the BOLD response over repeated presentations of the object images. Real-world experience with an object did, however, produce foci of increased activation in medial parietal and posterior cingulate cortex, regions that have previously been associated with the encoding and retrieval of remembered items in explicit memory paradigms. Our discovery that these regions are engaged during spontaneous recognition of real-world objects from their 2D image demonstrates that modulation of activity in medial regions by familiarity is neither stimulus nor task-specific. Overall, our results support separable coding in the human brain of the visual appearance of an object from the associations gained via real-world experience. The richness of object representations beyond their photographic image has important implications for understanding object recognition in both the human brain and in computational models.
Main programs:
afni_proc.py, 3dDeconvolve
Download scripts¶
To download, either:
... click the link(s) in the following table (perhaps Rightclick -> “Save Link As…”):
afni_proc.pycommand for multiecho FMRI3dDeconvolvecommand for the functional localizer data3dDeconvolvecommand for the experimental runs3dDeconvolvecommand 3dDeconvolve for the experimental runs, with individual modulation (1 beta per condition PER RUN)... or copy+paste into a terminal:
curl -O https://afni.nimh.nih.gov/pub/dist/doc/htmldoc/codex/fmri/media/2025_WardleEtal/s1.afni_proc.tcsh curl -O https://afni.nimh.nih.gov/pub/dist/doc/htmldoc/codex/fmri/media/2025_WardleEtal/s2a.3dDeconvolve_localizer.tcsh curl -O https://afni.nimh.nih.gov/pub/dist/doc/htmldoc/codex/fmri/media/2025_WardleEtal/s2b.3dDeconvolve_experimental.tcsh curl -O https://afni.nimh.nih.gov/pub/dist/doc/htmldoc/codex/fmri/media/2025_WardleEtal/s2c.3dDeconvolve_experimental_IM.tcsh
View scripts¶
s1.afni_proc.tcsh¶
This script contains an afni_proc.py command for ME-FMRI processing, with detailed regression modeling done in a separate follow-up script.
1#!/bin/tcsh
2
3# For this study, all analyses are conducted in each participant's
4# native brain space (i.e. aligned to in-session anatomical MPRAGE)
5# and displayed on the common mesh in SUMA for any group analyses
6
7# This script contains the afni_proc.py command used to process
8# task-based multi-echo multiband fMRI data in the below study. A
9# dummy 'regress' block is included to generate AFNI QC files, however
10# 3dDecolvolve is run separately in later scripts as the betas are
11# analyzed differently depending on the analysis and run type
12# (i.e. block design of localizer runs vs. event related design of
13# experimental runs, single beta estimate per condition for univariate
14# analyses vs. one beta estimate per condition per run for
15# multivariate analyses). Here the raw data is preprocessed in the
16# order collected for motion correction and alignment to the
17# in-session anatomical.
18
19# Used for processing in:
20#
21# Wardle, S. G., Rispoli, B., Roopchansingh, V. & Baker, C. (2024)
22# Brief encounters with real objects modulate medial parietal but
23# not occipitotemporal cortex. bioRxiv. 2024.08.05.606667
24# https://doi.org/10.1101/2024.08.05.606667
25
26# To run for a single participant, type (while providing an actual
27# value for P_ID):
28#
29# tcsh s1.afni_proc.tcsh P_ID
30
31# =============================================================================
32
33# collect user input
34if ( $#argv > 0 ) then
35 set pname = $argv[1]
36else
37 echo 'WARNING: NO PARTICIPANT IS ENTERED'
38endif
39
40# specify project directory
41set myroot = <INSERT PROJECT HOME DIRECTORY>
42
43# define directories
44set input_root = ${myroot}/rawdata/rawMRI
45set subdir = ${input_root}/${pname}
46set output_root = ${myroot}/preprocessed
47set suboutdir = ${output_root}/${pname}
48
49# print into terminal
50echo ${subdir}
51echo ${pname}
52echo ${input_root}
53
54# create the output root directory if it doesn't exist
55# nb: don't create the subject directory here because afni_proc.py will
56\mkdir -pv $output_root
57echo $output_root
58
59cd $subdir
60
61# generate the afni_proc.py script for this participant, for multiecho
62# FMRI processing
63
64# Notes
65#
66# + This command includes reverse phase encoding for EPI distortion
67# correction (`-blip_* ..` options)
68#
69# + AFNI's formulation for the optimal combination (OC; Posse et al.,
70# 1999) of multiple echos is used (`-combine_method OC`)
71#
72# + The 'regress' block is included to enable the QC HTML to be made,
73# but the actual regression commands are run later/separately, so no
74# `-regress_* ..` options were used here
75#
76
77afni_proc.py \
78 -subj_id ${pname} \
79 -script ${subdir}proc.${pname} \
80 -out_dir ${suboutdir} \
81 -dsets_me_echo ${pname}*tSeries*e01*orig.HEAD \
82 -dsets_me_echo ${pname}*tSeries*e02*orig.HEAD \
83 -dsets_me_echo ${pname}*tSeries*e03*orig.HEAD \
84 -echo_times 12.9 32.228 51.556 \
85 -reg_echo 2 \
86 -copy_anat ${subdir}/*anat.nii \
87 -blocks tshift align volreg blur mask combine regress \
88 -tcat_remove_first_trs 8 8 8 8 8 8 8 8 8 8 \
89 -blip_forward_dset ${pname}blip_forward-e02*orig.HEAD \
90 -blip_reverse_dset ${pname}blip_reverse-e02*orig.HEAD \
91 -volreg_align_to MIN_OUTLIER \
92 -volreg_align_e2a \
93 -blur_size 4 \
94 -combine_method OC \
95
96# run the afni_proc.py script for this participant
97tcsh -xef ${subdir}proc.${pname} |& tee ${subdir}output.proc.${pname}
s2a.3dDeconvolve_localizer.tcsh¶
1#!/bin/tcsh
2
3# In this script, 3dDeconvolve is run for the functional localizer
4# data (objects vs. scrambled objects, in a block design).
5#
6# The contrast [objects - scrambled objects] is later used to
7# functionally define the object-responsive regions of interest.
8#
9# LO and VTC in each participant's native brain space.
10
11# Used in:
12#
13# Wardle, S. G., Rispoli, B., Roopchansingh, V. & Baker, C. (2024)
14# Brief encounters with real objects modulate medial parietal but
15# not occipitotemporal cortex. bioRxiv. 2024.08.05.606667
16# https://doi.org/10.1101/2024.08.05.606667
17
18# To run, type:
19#
20# tcsh s2a.3dDeconvolve_localizer.tcsh
21
22# =============================================================================
23
24# specify participant number
25set pname = <ENTER SUBJECT NUMBER>
26
27# specify project directory
28set myroot = <INSERT PROJECT HOME DIRECTORY>
29
30# define directories
31set predir = ${myroot}/preprocessed/${pname} # preprocessed data
32set timedir = ${myroot}/timestamps/afni/${pname} # timestamps
33
34# create output directory
35set betadir = ${myroot}/betas_LOC/${pname}
36\mkdir -pv ${betadir}
37
38echo "++ Beta dir: ${betadir}"
39
40# ========================================
41# Run 3dDeconvolve for the functional localizer runs (block design)
42
43set motion_file = dfileLOC_rall.1D
44
45set epi_files = ( ${predir}/pb05.${pname}.r*.combineLOC+orig.HEAD )
46
47echo "++ EPI files: ${epi_files}"
48
493dDeconvolve \
50 -overwrite \
51 -input ${epi_files} \
52 -polort A \
53 -local_times \
54 -mask ${predir}/full_mask.${pname}+orig.HEAD \
55 -jobs 16 \
56 -xjpeg designLOC \
57 -x1D ${pname}.LOC \
58 -ortvec ${motion_file} motion \
59 -num_stimts 2 \
60 -stim_times 1 ${timedir}/${pname}_LOCtimestamps_objects.txt \
61 'BLOCK(16,1)' \
62 -stim_label 1 Objects \
63 -stim_times 2 ${timedir}/${pname}_LOCtimestamps_scrambled.txt \
64 'BLOCK(16,1)' \
65 -stim_label 2 Scrambled \
66 -num_glt 1 \
67 -gltsym 'SYM: +Objects -Scrambled' \
68 -glt_label 1 Objects-Scrambled \
69 -full_first \
70 -fout \
71 -tout \
72 -cbucket ${betadir}/${pname}.LOCstats.cbucket \
73 -bucket ${betadir}/${pname}.LOCstats.bucket
74
75tcsh ${betadir}/${pname}.REML_cmd
s2b.3dDeconvolve_experimental.tcsh¶
1#!/bin/tcsh
2
3# Run 3dDeconvolve for the experimental runs (event-related design, 48
4# conditions). This version produces 1 beta per condition (for use in
5# univariate analyses).
6
7# Used in:
8
9# Wardle, S. G., Rispoli, B., Roopchansingh, V. & Baker, C. (2024)
10# Brief encounters with real objects modulate medial parietal but
11# not occipitotemporal cortex. bioRxiv. 2024.08.05.606667
12# https://doi.org/10.1101/2024.08.05.606667
13
14# To run, type:
15#
16# tcsh s2b.3dDeconvolve_experimental.tcsh
17
18# ===============================================================================
19
20# specify participant number
21set pname = <ENTER SUBJECT NUMBER>
22
23# specify project directory
24set myroot = <INSERT PROJECT HOME DIRECTORY>
25
26# define directories
27set predir = ${myroot}/preprocessed/${pname} # preprocessed data
28set timedir = ${myroot}/timestamps/afni/${pname} # timestamps
29
30# create output directory
31set betadir=${myroot}/betas_EXP/$pname/
32\mkdir -pv ${betadir}
33
34echo "++ Beta dir: ${betadir}"
35
36set motion_file = dfileEXP_rall.1D
37
38set epi_files = ( ${predir}/pb05.$pname.r*.combine+orig.HEAD )
39
40echo "++ EPI files: ${epi_files}"
41
423dDeconvolve \
43 -overwrite \
44 -input ${EXP_epi_files} \
45 -polort A \
46 -local_times \
47 -mask ${predir}/full_mask.${pname}+orig.HEAD \
48 -jobs 16 \
49 -xjpeg designEXP \
50 -x1D ${pname}.EXP \
51 -ortvec ${EXP_motion_file} motion \
52 -num_stimts 48 \
53 -stim_times 1 ${timedir}/${pname}_Aind01.txt 'GAM' -stim_label 1 Aind01 \
54 -stim_times 2 ${timedir}/${pname}_Aind02.txt 'GAM' -stim_label 2 Aind02 \
55 -stim_times 3 ${timedir}/${pname}_Aind03.txt 'GAM' -stim_label 3 Aind03 \
56 -stim_times 4 ${timedir}/${pname}_Aind04.txt 'GAM' -stim_label 4 Aind04 \
57 -stim_times 5 ${timedir}/${pname}_Aind05.txt 'GAM' -stim_label 5 Aind05 \
58 -stim_times 6 ${timedir}/${pname}_Aind06.txt 'GAM' -stim_label 6 Aind06 \
59 -stim_times 7 ${timedir}/${pname}_Aind07.txt 'GAM' -stim_label 7 Aind07 \
60 -stim_times 8 ${timedir}/${pname}_Aind08.txt 'GAM' -stim_label 8 Aind08 \
61 -stim_times 9 ${timedir}/${pname}_Aind09.txt 'GAM' -stim_label 9 Aind09 \
62 -stim_times 10 ${timedir}/${pname}_Aind10.txt 'GAM' -stim_label 10 Aind10 \
63 -stim_times 11 ${timedir}/${pname}_Aind11.txt 'GAM' -stim_label 11 Aind11 \
64 -stim_times 12 ${timedir}/${pname}_Aind12.txt 'GAM' -stim_label 12 Aind12 \
65 -stim_times 13 ${timedir}/${pname}_Aout01.txt 'GAM' -stim_label 13 Aout01 \
66 -stim_times 14 ${timedir}/${pname}_Aout02.txt 'GAM' -stim_label 14 Aout02 \
67 -stim_times 15 ${timedir}/${pname}_Aout03.txt 'GAM' -stim_label 15 Aout03 \
68 -stim_times 16 ${timedir}/${pname}_Aout04.txt 'GAM' -stim_label 16 Aout04 \
69 -stim_times 17 ${timedir}/${pname}_Aout05.txt 'GAM' -stim_label 17 Aout05 \
70 -stim_times 18 ${timedir}/${pname}_Aout06.txt 'GAM' -stim_label 18 Aout06 \
71 -stim_times 19 ${timedir}/${pname}_Aout07.txt 'GAM' -stim_label 19 Aout07 \
72 -stim_times 20 ${timedir}/${pname}_Aout08.txt 'GAM' -stim_label 20 Aout08 \
73 -stim_times 21 ${timedir}/${pname}_Aout09.txt 'GAM' -stim_label 21 Aout09 \
74 -stim_times 22 ${timedir}/${pname}_Aout10.txt 'GAM' -stim_label 22 Aout10 \
75 -stim_times 23 ${timedir}/${pname}_Aout11.txt 'GAM' -stim_label 23 Aout11 \
76 -stim_times 24 ${timedir}/${pname}_Aout12.txt 'GAM' -stim_label 24 Aout12 \
77 -stim_times 25 ${timedir}/${pname}_Bind13.txt 'GAM' -stim_label 25 Bind13 \
78 -stim_times 26 ${timedir}/${pname}_Bind14.txt 'GAM' -stim_label 26 Bind14 \
79 -stim_times 27 ${timedir}/${pname}_Bind15.txt 'GAM' -stim_label 27 Bind15 \
80 -stim_times 28 ${timedir}/${pname}_Bind16.txt 'GAM' -stim_label 28 Bind16 \
81 -stim_times 29 ${timedir}/${pname}_Bind17.txt 'GAM' -stim_label 29 Bind17 \
82 -stim_times 30 ${timedir}/${pname}_Bind18.txt 'GAM' -stim_label 30 Bind18 \
83 -stim_times 31 ${timedir}/${pname}_Bind19.txt 'GAM' -stim_label 31 Bind19 \
84 -stim_times 32 ${timedir}/${pname}_Bind20.txt 'GAM' -stim_label 32 Bind20 \
85 -stim_times 33 ${timedir}/${pname}_Bind21.txt 'GAM' -stim_label 33 Bind21 \
86 -stim_times 34 ${timedir}/${pname}_Bind22.txt 'GAM' -stim_label 34 Bind22 \
87 -stim_times 35 ${timedir}/${pname}_Bind23.txt 'GAM' -stim_label 35 Bind23 \
88 -stim_times 36 ${timedir}/${pname}_Bind24.txt 'GAM' -stim_label 36 Bind24 \
89 -stim_times 37 ${timedir}/${pname}_Bout13.txt 'GAM' -stim_label 37 Bout13 \
90 -stim_times 38 ${timedir}/${pname}_Bout14.txt 'GAM' -stim_label 38 Bout14 \
91 -stim_times 39 ${timedir}/${pname}_Bout15.txt 'GAM' -stim_label 39 Bout15 \
92 -stim_times 40 ${timedir}/${pname}_Bout16.txt 'GAM' -stim_label 40 Bout16 \
93 -stim_times 41 ${timedir}/${pname}_Bout17.txt 'GAM' -stim_label 41 Bout17 \
94 -stim_times 42 ${timedir}/${pname}_Bout18.txt 'GAM' -stim_label 42 Bout18 \
95 -stim_times 43 ${timedir}/${pname}_Bout19.txt 'GAM' -stim_label 43 Bout19 \
96 -stim_times 44 ${timedir}/${pname}_Bout20.txt 'GAM' -stim_label 44 Bout20 \
97 -stim_times 45 ${timedir}/${pname}_Bout21.txt 'GAM' -stim_label 45 Bout21 \
98 -stim_times 46 ${timedir}/${pname}_Bout22.txt 'GAM' -stim_label 46 Bout22 \
99 -stim_times 47 ${timedir}/${pname}_Bout23.txt 'GAM' -stim_label 47 Bout23 \
100 -stim_times 48 ${timedir}/${pname}_Bout24.txt 'GAM' -stim_label 48 Bout24 \
101 -num_glt 1 \
102 -gltsym 'SYM: Aind01 Aind02 Aind03 Aind04 Aind05 Aind06 \
103 Aind07 Aind08 Aind09 Aind10 Aind11 Aind12 \
104 Aout01 Aout02 Aout03 Aout04 Aout05 Aout06 \
105 Aout07 Aout08 Aout09 Aout10 Aout11 Aout12 \
106 -Bind13 -Bind14 -Bind15 -Bind16 -Bind17 -Bind18 \
107 -Bind19 -Bind20 -Bind21 -Bind22 -Bind23 -Bind24 \
108 -Bout13 -Bout14 -Bout15 -Bout16 -Bout17 -Bout18 \
109 -Bout19 -Bout20 -Bout21 -Bout22 -Bout23 -Bout24' \
110 -glt_label 1 SetA-SetB \
111 -full_first \
112 -fout \
113 -tout \
114 -cbucket ${betadir}/${pname}.EXPstats.cbucket \
115 -bucket ${betadir}/${pname}.EXPstats.bucket
116
117
118tcsh ${betadir}/${pname}.REML_cmd
s2c.3dDeconvolve_experimental_IM.tcsh¶
1#!/bin/tcsh
2
3# Run 3dDeconvolve for the experimental runs (event-related design, 48
4# conditions). This version is individually modulated (for use in
5# multivariate analyses, e.g. decoding) = 1 beta per condition PER RUN.
6
7# Used in:
8#
9# Wardle, S. G., Rispoli, B., Roopchansingh, V. & Baker, C. (2024)
10# Brief encounters with real objects modulate medial parietal but
11# not occipitotemporal cortex. bioRxiv. 2024.08.05.606667
12# https://doi.org/10.1101/2024.08.05.606667
13
14# To run, type:
15#
16# tcsh s2c.3dDeconvolve_experimental_IM.tcsh
17
18# ==============================================================================
19
20# specify participant number
21set pname = <ENTER SUBJECT NUMBER>
22
23# specify project directory
24set myroot = <INSERT PROJECT HOME DIRECTORY>
25
26# define directories
27set predir = ${myroot}/preprocessed/${pname} # preprocessed data
28set timedir = ${myroot}/timestamps/afni/${pname} # timestamps
29
30# create output directory
31set betadir = ${myroot}/betas_EXP_IM/${pname}
32\mkdir -pv ${betadir}
33
34echo "++ Beta dir: ${betadir}"
35
36set motion_file = dfileEXP_rall.1D
37
38# for 3T (MULTIband Multiecho) data, combine is the last step == take
39# only the EXP files here (no LOC in filename).
40set epi_files = ( ${predir}/pb05.${pname}.r*.combine+orig.HEAD )
41
42echo "++ EPI files: ${epi_files}"
43
44
453dDeconvolve \
46 -overwrite \
47 -input $epi_files \
48 -polort A \
49 -local_times \
50 -mask ${predir}/full_mask.${pname}+orig.HEAD \
51 -jobs 16 \
52 -xjpeg designEXP \
53 -x1D ${pname}.EXP \
54 -ortvec ${motion_file} motion \
55 -num_stimts 48 \
56 -stim_times_IM 1 ${timedir}/${pname}_Aind01.txt 'GAM' -stim_label 1 Aind01 \
57 -stim_times_IM 2 ${timedir}/${pname}_Aind02.txt 'GAM' -stim_label 2 Aind02 \
58 -stim_times_IM 3 ${timedir}/${pname}_Aind03.txt 'GAM' -stim_label 3 Aind03 \
59 -stim_times_IM 4 ${timedir}/${pname}_Aind04.txt 'GAM' -stim_label 4 Aind04 \
60 -stim_times_IM 5 ${timedir}/${pname}_Aind05.txt 'GAM' -stim_label 5 Aind05 \
61 -stim_times_IM 6 ${timedir}/${pname}_Aind06.txt 'GAM' -stim_label 6 Aind06 \
62 -stim_times_IM 7 ${timedir}/${pname}_Aind07.txt 'GAM' -stim_label 7 Aind07 \
63 -stim_times_IM 8 ${timedir}/${pname}_Aind08.txt 'GAM' -stim_label 8 Aind08 \
64 -stim_times_IM 9 ${timedir}/${pname}_Aind09.txt 'GAM' -stim_label 9 Aind09 \
65 -stim_times_IM 10 ${timedir}/${pname}_Aind10.txt 'GAM' -stim_label 10 Aind10 \
66 -stim_times_IM 11 ${timedir}/${pname}_Aind11.txt 'GAM' -stim_label 11 Aind11 \
67 -stim_times_IM 12 ${timedir}/${pname}_Aind12.txt 'GAM' -stim_label 12 Aind12 \
68 -stim_times_IM 13 ${timedir}/${pname}_Aout01.txt 'GAM' -stim_label 13 Aout01 \
69 -stim_times_IM 14 ${timedir}/${pname}_Aout02.txt 'GAM' -stim_label 14 Aout02 \
70 -stim_times_IM 15 ${timedir}/${pname}_Aout03.txt 'GAM' -stim_label 15 Aout03 \
71 -stim_times_IM 16 ${timedir}/${pname}_Aout04.txt 'GAM' -stim_label 16 Aout04 \
72 -stim_times_IM 17 ${timedir}/${pname}_Aout05.txt 'GAM' -stim_label 17 Aout05 \
73 -stim_times_IM 18 ${timedir}/${pname}_Aout06.txt 'GAM' -stim_label 18 Aout06 \
74 -stim_times_IM 19 ${timedir}/${pname}_Aout07.txt 'GAM' -stim_label 19 Aout07 \
75 -stim_times_IM 20 ${timedir}/${pname}_Aout08.txt 'GAM' -stim_label 20 Aout08 \
76 -stim_times_IM 21 ${timedir}/${pname}_Aout09.txt 'GAM' -stim_label 21 Aout09 \
77 -stim_times_IM 22 ${timedir}/${pname}_Aout10.txt 'GAM' -stim_label 22 Aout10 \
78 -stim_times_IM 23 ${timedir}/${pname}_Aout11.txt 'GAM' -stim_label 23 Aout11 \
79 -stim_times_IM 24 ${timedir}/${pname}_Aout12.txt 'GAM' -stim_label 24 Aout12 \
80 -stim_times_IM 25 ${timedir}/${pname}_Bind13.txt 'GAM' -stim_label 25 Bind13 \
81 -stim_times_IM 26 ${timedir}/${pname}_Bind14.txt 'GAM' -stim_label 26 Bind14 \
82 -stim_times_IM 27 ${timedir}/${pname}_Bind15.txt 'GAM' -stim_label 27 Bind15 \
83 -stim_times_IM 28 ${timedir}/${pname}_Bind16.txt 'GAM' -stim_label 28 Bind16 \
84 -stim_times_IM 29 ${timedir}/${pname}_Bind17.txt 'GAM' -stim_label 29 Bind17 \
85 -stim_times_IM 30 ${timedir}/${pname}_Bind18.txt 'GAM' -stim_label 30 Bind18 \
86 -stim_times_IM 31 ${timedir}/${pname}_Bind19.txt 'GAM' -stim_label 31 Bind19 \
87 -stim_times_IM 32 ${timedir}/${pname}_Bind20.txt 'GAM' -stim_label 32 Bind20 \
88 -stim_times_IM 33 ${timedir}/${pname}_Bind21.txt 'GAM' -stim_label 33 Bind21 \
89 -stim_times_IM 34 ${timedir}/${pname}_Bind22.txt 'GAM' -stim_label 34 Bind22 \
90 -stim_times_IM 35 ${timedir}/${pname}_Bind23.txt 'GAM' -stim_label 35 Bind23 \
91 -stim_times_IM 36 ${timedir}/${pname}_Bind24.txt 'GAM' -stim_label 36 Bind24 \
92 -stim_times_IM 37 ${timedir}/${pname}_Bout13.txt 'GAM' -stim_label 37 Bout13 \
93 -stim_times_IM 38 ${timedir}/${pname}_Bout14.txt 'GAM' -stim_label 38 Bout14 \
94 -stim_times_IM 39 ${timedir}/${pname}_Bout15.txt 'GAM' -stim_label 39 Bout15 \
95 -stim_times_IM 40 ${timedir}/${pname}_Bout16.txt 'GAM' -stim_label 40 Bout16 \
96 -stim_times_IM 41 ${timedir}/${pname}_Bout17.txt 'GAM' -stim_label 41 Bout17 \
97 -stim_times_IM 42 ${timedir}/${pname}_Bout18.txt 'GAM' -stim_label 42 Bout18 \
98 -stim_times_IM 43 ${timedir}/${pname}_Bout19.txt 'GAM' -stim_label 43 Bout19 \
99 -stim_times_IM 44 ${timedir}/${pname}_Bout20.txt 'GAM' -stim_label 44 Bout20 \
100 -stim_times_IM 45 ${timedir}/${pname}_Bout21.txt 'GAM' -stim_label 45 Bout21 \
101 -stim_times_IM 46 ${timedir}/${pname}_Bout22.txt 'GAM' -stim_label 46 Bout22 \
102 -stim_times_IM 47 ${timedir}/${pname}_Bout23.txt 'GAM' -stim_label 47 Bout23 \
103 -stim_times_IM 48 ${timedir}/${pname}_Bout24.txt 'GAM' -stim_label 48 Bout24 \
104 -full_first \
105 -fout \
106 -tout \
107 -cbucket ${betadir}/${pname}.EXPstatsIM.cbucket \
108 -bucket ${betadir}/${pname}.EXPstatsIM.bucket
109
110
111tcsh ${betadir}/${pname}.REML_cmd
