Configure mkcontrast-sess

Contrasts are embodiments of hypotheses that you want to test. They are linear combinations of conditions. In order to construct a contrast, you need to know the numeric ID associated with each condition; this is specified in your paradigm file.

Suppose, for example, your par-files for the LDT experiment were encoded with the following conditions:

1 nonword
2 low-frequency word
3 high-frequency word

Contrasts between conditions are vector weightings of each condition.

The key thing is that the activated and contrast weightings must sum to zero. Continuing with this example, to configure the words > nonwords contrast for the analysis specified by the previous call to mkanalysis-sess, you would run the following command within your $SUBJECTS_DIR:

mkcontrast-sess -analysis LDT.sm6.lh -contrast words-v-nonwords -a 2 -a 3 -c 1

Pertinent to some of our borrowed data, we have the following conditions for the children's dyslexia study:

1 O+P+
2 O+P-
3 O-P+
4 O-P-
5 Null (Fixation)
6 Perceptual

Contrasts between all first four lexical conditions and the fixation condition (i.e., LEX-FIX) would be coded:

mkcontrast-sess -analysis booth500.sm6.lh -contrast LEX_v_FIX -a 1 -a 2 -a 3 -a 4 -c 5

Because this needs to be repeated for each surface/volume for which you wish to run an analysis, you might find it helpful to get a little practice automating this command with a shell script. Below is a little script I wrote that takes a couple parameters that correspond to how I generated the analysis directories when I ran mkanalysis-sess, and uses them to run mkcontrast-sess for the lh, rh and mni305 data:

!#/bin/bash
PARADIGM=$1 #the first parameter is the PAR-file name
#e.g., when using FAM.par, I would pass "FAM" as the first parameter
shift 
# the shift command moves all the parameters down the line
# like a conveyor belt. Now the first parameter "falls off", and the former
# 2nd parameter is now the first parameter
SMOOTHING=$1
# my analysis folders also include the smoothing kernel in the name (e.g., FAM.sm4.lh)
# so I will need to also use this information when running the command
shift

# finally, I will list all the hemispheres/volumes to loop through
HEMIS=( lh rh mni )

ANROOT="${PARADIGM}" #base name for the analysis directories
for hemi in "${HEMIS[@]}"
do
  mkcontrast-sess -analysis ${ANROOT}.${SMOOTHING}.${hemi} -contrast task -a 1$
  mkcontrast-sess -analysis ${ANROOT}.${SMOOTHING}.${hemi} -contrast fam-v-unf$
  mkcontrast-sess -analysis ${ANROOT}.${SMOOTHING}.${hemi} -contrast fam-v-bas$
  mkcontrast-sess -analysis ${ANROOT}.${SMOOTHING}.${hemi} -contrast unfam-v-b$
done