Hello,

I have been trying to set up a working Tracula configuration file and I keep getting an error:

set: Variable name must begin with a letter

I have searched previous message board questions and it seems that others have seen this same error. I have tried implementing some of their solutions, as well as searching my code for a problem, but I cannot find the source. Any advice? I am trying to get this working for a single subject for now. I used the dmrirc.example file as a template. I am calling trac-all from the tcsh command line with the following: trac-all -c dmrirc.example -prep

The contents of my dmrirc file are below, I have deleted most of the options for which I am just using the default, to make things simpler.

Thank you,

Mike

# FreeSurfer SUBJECTS_DIR

# T1 images and FreeSurfer segmentations are expected to be found here

setenv SUBJECTS_DIR /Users/mdatko/Desktop/tracula_tutorial/data/freesurfer

# Output directory where trac-all results will be saved

# Default: Same as SUBJECTS_DIR

set dtroot = /Users/mdatko/Desktop/tracula_tutorial/data/freesurfer

# Subject IDs

set subjlist = (a226A)

# In case you want to analyze only Huey and Louie

# Default: Run analysis on all subjects

set runlist = (1)

# Input diffusion DICOMs (file names relative to dcmroot)

# If original DICOMs don't exist, these can be in other image format

# but then bvecfile and bvalfile must be specified (see below)

set dcmroot = /Users/mdatko/Desktop/tracula_tutorial/data/dti

set dcmlist = (a226A/226A/data.nii.gz)

# Diffusion gradient table

# Must be specified if inputs are not MGH DICOMs

# Three-column format, one row for each volume in the diffusion data set

# Default: Read from DICOM header

set bvecfile = /Users/mdatko/Desktop/tracula_tutorial/data/dti/a226A/226A/bvecs.txt

# Diffusion b-value table

# Must be specified if inputs are not MGH DICOMs

# Single-column format, one value for each volume in the diffusion data set

# Default: Read from DICOM header

set bvalfile = /Users/mdatko/Desktop/tracula_tutorial/data/dti/a226A/226A/bvals.txt

# Perform registration-based eddy-current compensation?

# Default: 1 (yes)

set doeddy = 0

# Perform diffusion-to-T1 registration by bbregister?

# Default: 1 (yes)

set doregbbr = 1

# Perform registration of T1 to MNI template?

# Default: 1 (yes)

set doregmni = 1

# Use brain mask extracted from T1 image instead of low-b diffusion image?

# Has no effect if there is no T1 data

# Default: 1 (yes)

set usemaskanat = 1

# Paths to reconstruct

# Default: All paths in the atlas

set pathlist = ( lh.cst_AS rh.cst_AS \

lh.unc_AS rh.unc_AS \

lh.ilf_AS rh.ilf_AS \

fmajor_PP fminor_PP \

lh.atr_PP rh.atr_PP \

lh.ccg_PP rh.ccg_PP \

lh.cab_PP rh.cab_PP \

lh.slfp_PP rh.slfp_PP \

lh.slft_PP rh.slft_PP )

# Number of path control points

# It can be a single number for all paths or a different number for each of the

# paths specified in pathlist

# Default: 7 for the forceps major, 6 for the corticospinal tract,

# 4 for the angular bundle, and 5 for all other paths

set ncpts = ( 6 6 5 5 5 5 7 5 5 5 5 5 4 4 5 5 5 5 )

# List of training subjects

# This text file lists the locations of training subject directories

# Default: $FREESURFER_HOME/trctrain/trainlist.txt

set trainfile = $FREESURFER_HOME/trctrain/trainlist.txt

# Number of "sticks" (anisotropic diffusion compartments) in the bedpostx

# ball-and-stick model

# Default: 2

set nstick = 2

# Number of MCMC burn-in iterations

# (Path samples drawn initially by MCMC algorithm and discarded)

# Default: 200

set nburnin = 200

# Number of MCMC iterations

# (Path samples drawn by MCMC algorithm and used to estimate path distribution)

# Default: 7500

set nsample = 7500

# Frequency with which MCMC path samples are retained for path distribution

# Default: 5 (keep every 5th sample)

set nkeep = 5