removing the prestim will not affect any of the graphs. Is the discrepancy between group time course (you show on fsaverage) and the individual (from roisummary-sess)?

doug

erickenb@nmr.mgh.harvard.edu wrote:
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

-- 
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
Phone Number: 617-724-2358 
Fax: 617-726-7422

In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting