We have been working on timecourse data that has had the prestim removed
(at the mkcontrast-sess level). The timecourses obtained with
roisummary-sess do not seem to correspond to the timecourses plotted with
tksurfer. The discrepancy is two-fold: the peak latency and the overall
timecourse shape. What could be the source of this discrepancy? This is
particularly obvious for the ALC (alcohol) and PLAC (placebo) conditions -
their peak activity is at the same timepoint in tksurfer but at different
points in the text output (as shown in the Excel file)
Included are:
- The original text ouput file
- Excel file with the text output along with the graph,
- screen shots of the tksurfer timecourses and activity - on the left
side is the timecourse of the activity with prestim removed and on the
right is the timecourse with all the points
We have also included some of the relevant script
commands we used throughout the analysis.
The entire working directory can be found here...
/space/veles/3/users/sheeva/TEST-S51
This is the analysis I used...
analysis.info
analysis SNfir-test
TR 2
designtype event-related
nconditions 4
parname para.para
fsd bold_SN
funcstem fmcsm8
tpexclude
-timewindow 18
-prestim 4
-polyfit 2
-TER 2
-nskip 0
-fwhm 0
-extreg mcextreg
-nextreg 3
example of mkcont command...removing the prestim
entire script can be found
here.../space/veles/3/users/sheeva/TEST-S51/mkcont-sess.csh
setenv aa 1
setenv cc 0
echo '0 0 0 1 1 1 1 1 1' | mkcontrast-sess -analysis SNfir-test -contrast
CvFix-RMPS -a $aa -c $cc -rmprestim
stxgrinder-sess -analysis SNfir-test -contrast CvFix-RMPS -sf
51-SN-intragp
Finally, here is the csh that I used for making the timecourses....
setenv SUBJECTS_DIR /space/perun/3/users/cylix/Subjects
foreach roi (lh-AC-med)
foreach num (51)
foreach ltr (a j)
mri_label2label --srclabel
/space/veles/3/users/sheeva/fMRIanal/fsavg-labels/RMPS/surface/LH/${roi}.label
--srcsubject fsaverage --trgsubject sub${num}_recon --trglabel
sub${num}_recon/label/${roi}.label --hemi lh --regmethod surface
func2roi-sess -roidef ${roi}.label -labelfile
/space/veles/3/users/sheeva/fMRIanal/fsavg-labels/RMPS/surface/LH/${roi}.label
-d . -s sub${num}${ltr} -analysis SNfir-test
roisummary-sess -analysis SNfir-test -roidef ${roi}.label -d . -sf
SN51-a-j-RMPS -sumfile TEST-S51aj/Unmasked/LH/${roi}-unmasked.txt
end
end
end
Thank you
Elizabeth
sub51a sub51j
27 27
27 27
459.1759 342.7553
7.4028 5.0723
384.000000 384.000000
2.000000 2.000000
4.000000 4.000000
4 4
9 9
-0.8460 -0.0432
-0.0129 -0.2949
0.8421 0.1958
0.1448 1.5211
0.3556 0.5827
0.5880 -0.4853
-0.2520 -0.3832
-0.5757 0.0549
-1.0957 0.1103
-0.4473 0.0758
0.0342 -0.3010
0.9660 -0.1193
-0.1295 0.9076
0.0990 -0.3110
0.2888 -1.1811
0.8339 -0.1293
0.5635 -0.5169
-0.0138 -1.2245
-1.2507 0.6735
-0.6719 0.2709
-0.5685 0.2360
-0.2412 0.4716
-0.5896 0.3330
-0.8171 -0.5989
-0.7288 -0.3390
-0.1770 -0.5283
-0.7550 -0.7002
-0.7418 0.5066
-0.2738 -0.1271
0.1742 0.2551
-0.5554 1.1041
-0.9806 0.7834
-0.4049 -0.0559
-0.9098 -0.2644
-0.3250 -0.2043
-0.6926 0.3521
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer