Freesurfer

freesurfer@nmr.mgh.harvard.edu

30132 discussions

[Major Upgrade] - iSurf Brainview 4.0.0
by Pedro Paulo de Magalhães Oliveira Junior 20 Oct '11

20 Oct '11

Today it was released the version 4.0.0 of iSurf Brainview, it's a majorversion release for IPad, IPhone and IPod Touch. It's also a landmark release. Major New Features - 3D Visualization of Cortical structure - 3D Model is fully "touchable" and displays the name of touched structure using aparc.annot Minor Changes - iOS5 compliance - Memory usage improvement - Versioning As usual suggestions and feedbacks are welcome. The link for download is: http://itunes.apple.com/WebObjects/MZStore.woa/wa/viewSoftware?id=381072423… A personal note: This App is being used in the Radiology course. All med students in São Paulo University now have basic MRI neuroanatomy training using iSurf --------------------------------------------------------------------- Pedro Paulo de Magalhães Oliveira Junior Netfilter & SpeedComm Telecom -- www.netfilter.com.br -- For mobile: http://itunes.apple.com/br/artist/netfilter/id365306441

1 0

Subcortical segmentation atlas
by cloud.ctrl＠me.com 19 Oct '11

19 Oct '11

Hello, I am working on comparing hand-drawn subcortical volumes to those output automatically by FreeSurfer. While reading about the segmentation process on the wiki, I came across this line: "The final segmentation is based on both a subject-independent probabilistic atlas and subject-specific measured values. The atlas is built from a training set, i.e., a set of subjects whose brains (surfaces or volumes) have been labeled by hand. These labels are then mapped into a common space (Talairach space for volumes..." (http://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferAnalysisPipelineOverview) Because of the way we decided to delineate our ROI, we need to compare the shapes we defined in AFNI to those defined by this training set and output by FreeSurfer. I've been looking for more information about this probabilistic atlas, but I cannot find any more information about it. The only thing I found was the last line of this thread from the mail archives: Am 23.08.2011 um 04:03 schrieb Anthony Dick: > The subcortical segmentation identifies subcortical structures such as > thalamus, basal ganglia structures, hippocampus, amygdala, and > cerebellum. The cortical parcellation delineates different areas of the > cortex, such as the superior temporal sulcus, the central sulcus, the > pars triangularis, opercularis, and orbitalis of the inferior frontal > gyrus. A recommendation is to get a copy of the Duvernoy atlas, which > covers definitions of different parts of the cortex. The Mai atlas is > good for learning subcortical structures. Was the Mai atlas used to create the define the structures in the training set? Could you provide me with any more information regarding the training set?

2 1

Re: [Freesurfer] Group Average of Retinotopy Results
by Michelle Umali 19 Oct '11

19 Oct '11

Hi Doug and the Freesurfers, Thanks for your help. I have 2 questions.: 1) I did the retinotopy analysis on the fsaverage surface as you suggested below and then ran isxconcat-sess for my eccen and polar data: e.g. isxconcat-sess -analysis rtopy.fsaverage.lh -contrast polar -sf sessid -o group_polar_lh but how do you do this for the fieldsign analysis, since the cess.nii and other files don't exist? 2) So when trying to run a group mri_glmfit on the polar and eccen maps with: mri_glmfit --y ces.nii.gz --wls cesvar.nii.gz --osgm --surface fsaverage lh --glmdir group_eccen_lh.wls --nii.gz I got an error, because I have two ces.nii.gz and cesvar.nii.gz files each. I am not sure how this happened. I only have one run each of polar and wedge data. What may have generated these two files and which one do I use? Thanks. Michelle On Thu, Oct 13, 2011 at 11:36 AM, Douglas N Greve <greve(a)nmr.mgh.harvard.edu > wrote: > When you run preproc-sess use -surface fsaverage lhrh instead of -surface > self lhrh. Then create a new analysis with the same parameters, except use > -surface fsaverage lh instead of -surface self lh. > doug > > Michelle Umali wrote: > >> Hi Doug, >> I'm a little confused as to how to assign functional data to the fsaverage >> space. >> 1) Do you mean just change each subject's subjectname file to fsaverage >> and then do the same thing as before for each person? >> 2) Also, when I specify the analysis do I change -rtopy.self.lh to >> rtopy.fsaverage.lh, i.e. will isconcat later be sensitive to analysis name? >> >> e.g. mkanalysis-sess -a rtopy.self.lh -surface self lh -TR 2 -retinotopy >> 48 -paradigm rtopy.par -nskip 4 -fwhm 0 -fsd bold -force >> >> then the same for selxavg3-sess and fieldsign-sess >> >> and then isxconcat-sess would be: >> isxconcat-sess -sf sessid -analysis rtopy.fsaverage.lh -o group_lh >> >> >> Thanks. >> Michelle >> >> On Wed, Oct 12, 2011 at 3:54 PM, Douglas N Greve < >> greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.**edu<greve(a)nmr.mgh.harvard.edu>>> >> wrote: >> >> If you want to combine across subjects, then you need to run the >> analysis on fsaverage (ie, rerun preproc-sess, create a new >> analysis, run selxavg3-sess, and then isxconcat-sess). >> doug >> >> Michelle Umali wrote: >> >> Dear Freesurfers, >> I've generated individual polar, eccentricity, and fieldsign >> maps for each of my subjects. I would like to generate group >> average maps of these that I can visualize on the fsaverage brain. >> >> When I ran: >> isxconcat-sess -sf sessid -analysis rtopy.self.lh -o group_lh >> >> I got: >> ERROR: analysis space is self surface, not supported >> >> How do I do a group average for the 3 different maps? >> >> Thanks. >> Michelle >> >> For each person I ran recon-all and flattened occipital patches. >> Then: >> mkanalysis-sess -a rtopy.self.lh -surface self lh -TR 2 >> -retinotopy 48 -paradigm rtopy.par -nskip 4 -fwhm 0 -fsd bold >> -force >> mkanalysis-sess -a rtopy.self.rh -surface self rh -TR 2 >> -retinotopy 48 -paradigm rtopy.par -nskip 4 -fwhm 0 -fsd bold >> -force >> >> preproc-sess -surface self lhrh -fwhm 5 -per-run -s sjXX -fsd >> bold -force >> >> selxavg3-sess -a rtopy.self.lh -s sjXX -force >> selxavg3-sess -a rtopy.self.rh -s sjXX -force >> >> fieldsign-sess -a rtopy.self.lh -occip -fwhm 20 -s sjXX >> fieldsign-sess -a rtopy.self.rh -occip -fwhm 20 -s sjXX >> >> ------------------------------**------------------------------** >> ------------ >> >> ______________________________**_________________ >> Freesurfer mailing list >> Freesurfer(a)nmr.mgh.harvard.edu >> <mailto:Freesurfer@nmr.mgh.**harvard.edu<Freesurfer(a)nmr.mgh.harvard.edu> >> > >> >> https://mail.nmr.mgh.harvard.**edu/mailman/listinfo/**freesurfer<https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer> >> >> >> -- Douglas N. Greve, Ph.D. >> MGH-NMR Center >> greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.**edu<greve(a)nmr.mgh.harvard.edu> >> > >> Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422 >> <tel:617-726-7422> >> >> Bugs: surfer.nmr.mgh.harvard.edu/**fswiki/BugReporting<http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> >> <http://surfer.nmr.mgh.**harvard.edu/fswiki/**BugReporting<http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> >> > >> FileDrop: www.nmr.mgh.harvard.edu/**facility/filedrop/index.html<http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> >> <http://www.nmr.mgh.harvard.**edu/facility/filedrop/index.**html<http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> >> > >> >> >> >> >> The information in this e-mail is intended only for the person to >> whom it is >> addressed. If you believe this e-mail was sent to you in error and >> the e-mail >> contains patient information, please contact the Partners >> Compliance HelpLine at >> http://www.partners.org/**complianceline<http://www.partners.org/complianceline>. If the e-mail was sent to >> you in error >> but does not contain patient information, please contact the >> sender and properly >> dispose of the e-mail. >> >> >> > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu > Phone Number: 617-724-2358 Fax: 617-726-7422 > > Bugs: surfer.nmr.mgh.harvard.edu/**fswiki/BugReporting<http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: www.nmr.mgh.harvard.edu/**facility/filedrop/index.html<http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > >

2 1

FreeSurfer Questions
by Christopher McCarthy 19 Oct '11

19 Oct '11

Dear FreeSurfer Gurus, We are looking for a way to normalize the intensity histogram for our subjects. In the past we have used Brain Image Java, using FSL-Fast inhomogeneity corrections. Unfortunately we are unable to open the output files (.img, .hdr) in FreeSurfer. When trying to import the files (after converting them to nii) we get a talaraich error: ERROR: talairach_afd: Talairach Transform: transforms/talairach.xfm ***FAILED*** (p=0.0000, pval=0.0000 < threshold=0.0050) here is the mri_info for a subject we get this error with: Volume information for S029_nii type: nii dimensions: 256 x 256 x 100 voxel sizes: 0.9375, 0.9375, 1.5000 type: SHORT (4) fov: 240.000 dof: 0 xstart: -120.0, xend: 120.0 ystart: -120.0, yend: 120.0 zstart: -75.0, zend: 75.0 TR: 0.00 msec, TE: 0.00 msec, TI: 0.00 msec, flip angle: 0.00 degrees nframes: 1 PhEncDir: UNKNOWN ras xform present xform info: x_r = -1.0000, y_r = 0.0000, z_r = 0.0000, c_r = 0.0000 : x_a = 0.0000, y_a = 1.0000, z_a = 0.0000, c_a = 0.0000 : x_s = 0.0000, y_s = 0.0000, z_s = 1.0000, c_s = 0.0000 Orientation : LAS Primary Slice Direction: axial voxel to ras transform: -0.9375 0.0000 0.0000 120.0000 0.0000 0.9375 0.0000 -120.0000 0.0000 0.0000 1.5000 -75.0000 0.0000 0.0000 0.0000 1.0000 voxel-to-ras determinant -1.31836 ras to voxel transform: -1.0667 0.0000 0.0000 128.0000 -0.0000 1.0667 -0.0000 128.0000 -0.0000 -0.0000 0.6667 50.0000 0.0000 0.0000 0.0000 1.0000 We are either looking to get the files we currently have to import correctly in FreeSurfer, or import our DICOM files into a program and normalize the intensity histogram and output a file freesurfer can read. We are looking to get rid of hyper intensities in the temporal and frontal lobes that cause main and pial surface errors. Thanks for whatever help you can give, now and in the future. Chris McCarthy

2 1

problem with generating cesvar.nii
by SHAHIN NASR 19 Oct '11

19 Oct '11

Hi, I have generated two different functional connectivity maps for two different ROI-based seeds. These maps are based on group-average of the same subjects (generated by using separate isxconcat-sess commands ), and now I want to see the difference map. To do so, I used mris_calc as below. > mris_calc -o diff_map/rh.ces.nii map1/rh.ces.nii sub map2/rh.ces.nii Also, to generate the rh.cesvar.nii file I used this command: > mris_calc -o diff_map1/rh.cesvar.nii map1/rh.cesvar.nii sub map2/rh.cesvar.nii But then when I used the glmfit with WLS option as below: > mri_glmfit --y rh.ces.nii --osgm --glmdir rh.ffx.osgm.wls --nii --mask ../rh.mask.nii --yffxvar rh.cesvar.nii --ffxdofdat ../ffxdof.dat --surf fsaverage rh --wls rh.cesvar.nii I faced an error saying: "Error: MRInormWeights: values less than or eq to 0" which I assumed is related to an error in rh.cesvar.nii file (am I right?). Would you tell me how can I fix this problem? P.S.: please note that sequence of subjects and all MRI parameters are exactly the same between the two maps. -- Shahin Nasr PhD in Cognitive Neuroscience Martinos Imaging Center, MGH Harvard Medical School

2 5

Re: [Freesurfer] reg-feat2anat failure
by Douglas N Greve 19 Oct '11

19 Oct '11

This file is created when you run reg-feat2anat and should be the one that gets changed with --manxfm func2anat. You can handle this differently if you 1. Manually edit the anat2exf.register.dat so that it is relatively close (no need to be perfect) cd feat/reg/freesurfer tkregister2 --mov ../../example_func.nii.gz --reg anat2exf.register.dat --surfs 2. Run bbregister: cp anat2exf.register.dat backup.anat2exf.register.dat # Make a backup bbregister --mov ../../example_func.nii.gz --init-reg anat2exf.register.dat --reg anat2exf.register.dat --bold this second step will fine-tune your hand-edited version. Then go back and run aparc2feat, etc. doug Michelle Umali wrote: > Hi Doug, > The file that changed was exf2std.register.dat. How do I specify > that I wanted that anat2exf.register.dat is supposed to be edited? At > what point in the processing stream does this occur? > > Did editing exf2std.register.dat do anything? > > Thanks. > Michelle > > > > On Tue, Oct 18, 2011 at 2:17 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > The registration that should be used for all of these is > anat2exf.register.dat found in feat/reg/freesurfer. Can you verify > that that was the one that you manually changed? Look at the time > stamp as compared to the other files in that folder. > doug > > Michelle Umali wrote: > > Hi Doug, > I am using reg-feat2anat on the feat directory of each run. > reg-feat2anat --feat feadir.feat --subject XX > Then I checked and adjusted the registration with: > reg-feat2anat --feat featdir.feat --manual > Next I checked and adjusted the registration to the standard with: > reg-feat2anat --feat featdir.feat --manxfm func2std > > To generate cortical binary masks I used: > aparc2feat --feat featdir.feat --annot BA > fslmaths featdir.feat/reg/freesurfer/lh.BA.nii.gz -thr 10 > -uthr 10 featdir.feat/reg/freesurfer/V1_l.nii.gz > > To generate subcortical binary masks I did: > aseg2feat --feat featdir.feat --aseg aparc+aseg > fslmaths featdir.feat/reg/freesurfer/aparc+aseg.nii.gz > -thr 54 -uthr 54 featdir.feat/reg/freesurfer/amygdala_r.nii.gz > > When I looked at the masks in FSLview and checked the masks > against the other atlases, everything was off even after > painstakingly editing the registrations. > Is there something I am doing wrong? > > Thanks. > Michelle > > On Mon, Oct 17, 2011 at 1:36 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > <mailto:greve@nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>>> wrote: > > How are you generating the masks? > doug > > Michelle Umali wrote: > > Hi Freesurfers, > I've run reg-feat2anat on several subjects in order to > generate binary masks to use in an FSL ppi. > For 3 subjects, this completely failed. > > When I looked at the fsl's own registration of fsl feat 2 > standard (via tkregister2), the registration is fine, but > after running reg-feat2anat, the registration of the > feat to > the freesurfer anatomical and the feat to standard are > completely off. > > I manually fixed and saved the registration edits via > reg-feat2anat --manual, but when I generated the masks > again, > they were totally wrong still. > When looking at freesurfer's feat to standard > registration, it > looks like all the brains weren't tall enough in the > saggital > view. Am I supposed to stretch the feat brain to fit? > > How can I fix this so that the masks are accurate? > > Thanks. > Michelle > > ------------------------------------------------------------------------ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > <mailto:Freesurfer@nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu>> > > > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > -- Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu> > <mailto:greve@nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>> > Phone Number: 617-724-2358 <tel:617-724-2358> > <tel:617-724-2358 <tel:617-724-2358>> Fax: 617-726-7422 > <tel:617-726-7422> > <tel:617-726-7422 <tel:617-726-7422>> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > > > > The information in this e-mail is intended only for the > person to > whom it is > addressed. If you believe this e-mail was sent to you in > error and > the e-mail > contains patient information, please contact the Partners > Compliance HelpLine at > http://www.partners.org/complianceline . If the e-mail was > sent to > you in error > but does not contain patient information, please contact the > sender and properly > dispose of the e-mail. > > > ------------------------------------------------------------------------ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422 > <tel:617-726-7422> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html

1 0

Re: [Freesurfer] QDec questions
by Douglas N Greve 18 Oct '11

18 Oct '11

Yes, see the answer to your first question below. doug krista kelly wrote: > Am I able to do vertex-by-vertex analyses using freesurfer if I have > more than 2 discrete variables? If so, how? > > Is it also ok to use the unsmoothed data when doing surface area label > analyses? > > Thanks! > > On Tue, Oct 18, 2011 at 3:14 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > Sorry, I meant to say that you cannot use QDEC with more than 2 > discrete variables (you have 3). You do not need to extract data > from the smoothed thickness files, but the results you get from > your ROI analysis may then differ a bit from that you get from the > QDEC analysis. > doug > > krista kelly wrote: > > Thanks Doug! Just a few more things: > > I only have two groups, so my contrast was wrong. Should it > only be 1 -1, even if I have other variables I want to control > for? Or should I have a contrast file for each discrete variable? > > How would I extract the smoothed data using mri_anatomical_stats? > > I realize the importance of smoothing data while doing a > vertex-by-vertex analysis in QDec, but if I'm doing an ROI > analysis by extracting the average thickness in a certain > label (i.e. v1), is it necessary to do the stats on the > smoothed data from each participant? > > Thanks! > > > On Tue, Oct 18, 2011 at 2:23 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > <mailto:greve@nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>>> wrote: > > > > krista kelly wrote: > > Hi, just wondering if anyone has any advice on the > following > questions: > > I'm using QDEC to analyze cortical thickness differences > between groups and I have a few questions: > > 1) There are a few variables that I need to control > for. First > variable is scanner - I used two different scanners of > different field strengths (1.5T GE and 3T Siemens) so I > wanted > to control for scanner. I also wanted to control for > age and > gender. Here's my ideal design setup: > > discrete variables - groups, scanner, gender > nuisance variable - age > > contrast would be 1 -1 0 0 0 0 > I realize that scanner, gender, and group are all discrete > variables, and that age would be a continuous one. > However, in > the design tab, I am not able to select all 3 discrete > variables (groups, scanner, and gender). Should I designate > scanner or gender as a nuisance variable instead? Would it > matter which one I use as the nuisance variable? > Possible setup: > > discrete vars - groups, scanner > nuisance vars - age, gender > > contrast - 1 -1 0 0 (Correct?) > > Unfortunately, you cannot use QDEC if you have three groups. > You'll have to create an FSGD file and contrasts by hand, > then run > mri_glmfit. Do a search for FSGD on our wiki and go to the > examples page. I think there is an example there for what > you want > to do. > > > 2) If I use DODS, should I demean the ages if age is a > nuisance variable? > > You should actually test whether there is an interaction > between > age and your group contrast of interest (using the DODS > model). If > there is not interaction, then switch to a DOSS where demeaning > does not matter. > > > 3) If I use a smoothing kernel during the analysis, then > create an ROI and map that ROI to each individual subject, > will the resulting .stats file in each subjects stats > folder > reflect the smoothed data or does mapping onto each ROI > just > extract the raw unsmoothed data from each vertex? > > It depends on how you do the extraction. If you are using > mris_anatomical_stats on each individual's unsmoothed data, > then > it will reflect the unsmoothed. > > doug > > > Thanks! > Krista > > ------------------------------------------------------------------------ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > <mailto:Freesurfer@nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu>> > > > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > -- Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu> > <mailto:greve@nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>> > Phone Number: 617-724-2358 <tel:617-724-2358> > <tel:617-724-2358 <tel:617-724-2358>> Fax: 617-726-7422 > <tel:617-726-7422> > <tel:617-726-7422 <tel:617-726-7422>> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > > > > The information in this e-mail is intended only for the > person to > whom it is > addressed. If you believe this e-mail was sent to you in > error and > the e-mail > contains patient information, please contact the Partners > Compliance HelpLine at > http://www.partners.org/complianceline . If the e-mail was > sent to > you in error > but does not contain patient information, please contact the > sender and properly > dispose of the e-mail. > > > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422 > <tel:617-726-7422> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html

1 0

Re: [Freesurfer] Surface area measure
by Douglas N Greve 18 Oct '11

18 Oct '11

I don't know if that has ever been explored. doug krista kelly wrote: > Great! Do you know if it is better to use one surface area (i.e. WM or > pial) over the other when investigating differences between groups? > > On Tue, Oct 18, 2011 at 2:24 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > The white surface. > doug > > krista kelly wrote: > > When I extract surface area data using the following > commands, which measure of surface area am I getting, the pial > surface, the WM surface, or the total of both combined? > > aparcstats2table --subjects 1 2 3 4 --parc aparc.a2009 --hemi > rh --meas area --tablefile aparc_stats_rh.txt > > Thanks! > Krista > || > ------------------------------------------------------------------------ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422 > <tel:617-726-7422> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > > > The information in this e-mail is intended only for the person to > whom it is > addressed. If you believe this e-mail was sent to you in error and > the e-mail > contains patient information, please contact the Partners > Compliance HelpLine at > http://www.partners.org/complianceline . If the e-mail was sent to > you in error > but does not contain patient information, please contact the > sender and properly > dispose of the e-mail. > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html

1 0

Re: [Freesurfer] QDec questions
by Douglas N Greve 18 Oct '11

18 Oct '11

Sorry, I meant to say that you cannot use QDEC with more than 2 discrete variables (you have 3). You do not need to extract data from the smoothed thickness files, but the results you get from your ROI analysis may then differ a bit from that you get from the QDEC analysis. doug krista kelly wrote: > Thanks Doug! Just a few more things: > > I only have two groups, so my contrast was wrong. Should it only be 1 > -1, even if I have other variables I want to control for? Or should I > have a contrast file for each discrete variable? > > How would I extract the smoothed data using mri_anatomical_stats? > > I realize the importance of smoothing data while doing a > vertex-by-vertex analysis in QDec, but if I'm doing an ROI analysis by > extracting the average thickness in a certain label (i.e. v1), is it > necessary to do the stats on the smoothed data from each participant? > > Thanks! > > > On Tue, Oct 18, 2011 at 2:23 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > > > krista kelly wrote: > > Hi, just wondering if anyone has any advice on the following > questions: > > I'm using QDEC to analyze cortical thickness differences > between groups and I have a few questions: > > 1) There are a few variables that I need to control for. First > variable is scanner - I used two different scanners of > different field strengths (1.5T GE and 3T Siemens) so I wanted > to control for scanner. I also wanted to control for age and > gender. Here's my ideal design setup: > > discrete variables - groups, scanner, gender > nuisance variable - age > > contrast would be 1 -1 0 0 0 0 > I realize that scanner, gender, and group are all discrete > variables, and that age would be a continuous one. However, in > the design tab, I am not able to select all 3 discrete > variables (groups, scanner, and gender). Should I designate > scanner or gender as a nuisance variable instead? Would it > matter which one I use as the nuisance variable? Possible setup: > > discrete vars - groups, scanner > nuisance vars - age, gender > > contrast - 1 -1 0 0 (Correct?) > > Unfortunately, you cannot use QDEC if you have three groups. > You'll have to create an FSGD file and contrasts by hand, then run > mri_glmfit. Do a search for FSGD on our wiki and go to the > examples page. I think there is an example there for what you want > to do. > > > 2) If I use DODS, should I demean the ages if age is a > nuisance variable? > > You should actually test whether there is an interaction between > age and your group contrast of interest (using the DODS model). If > there is not interaction, then switch to a DOSS where demeaning > does not matter. > > > 3) If I use a smoothing kernel during the analysis, then > create an ROI and map that ROI to each individual subject, > will the resulting .stats file in each subjects stats folder > reflect the smoothed data or does mapping onto each ROI just > extract the raw unsmoothed data from each vertex? > > It depends on how you do the extraction. If you are using > mris_anatomical_stats on each individual's unsmoothed data, then > it will reflect the unsmoothed. > > doug > > > Thanks! > Krista > ------------------------------------------------------------------------ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422 > <tel:617-726-7422> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > > > The information in this e-mail is intended only for the person to > whom it is > addressed. If you believe this e-mail was sent to you in error and > the e-mail > contains patient information, please contact the Partners > Compliance HelpLine at > http://www.partners.org/complianceline . If the e-mail was sent to > you in error > but does not contain patient information, please contact the > sender and properly > dispose of the e-mail. > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html

1 0

eTIV in qdec
by dolphinede＠aol.com 18 Oct '11

18 Oct '11

Dear FreeSurfer experts, my idea was to put the eTIV as a variable in qdec. I have put this in my qdec.table.dat and when I run qdec and choose eTIV as variable I get the following error and qdec closes some seconds later automatically: Creating output directory /Users/dolphinede/Desktop/Probanden/qdec/Untitled Loading y from /Users/dolphinede/Desktop/Probanden/qdec/Untitled/y.mgh INFO: gd2mtx_method is dods Saving design matrix to /Users/dolphinede/Desktop/Probanden/qdec/Untitled/Xg.dat Normalized matrix condition is 389.695 Matrix condition is 1e+08 Found 148151 points in label. Pruning voxels by thr: 0.000000 Found 148142 voxels in mask Saving mask to /Users/dolphinede/Desktop/Probanden/qdec/Untitled/mask.mgh Reshaping mriglm->mask... search space = 73700.662359 DOF = 77 Starting fit and test Fit completed in 0.0607667 minutes Computing spatial AR1 on surface Residual: ar1mn=0.995148, ar1std=0.005328, gstd=5.802773, fwhm=13.664487 Writing results lh-Avg-Intercept-thickness maxvox sig=35.9893 F=546.618 at index 1515 0 0 seed=1319335878 lh-Avg-thickness-eTIV-Cor maxvox sig=6.48644 F=31.2833 at index 51144 0 0 seed=1319335878 mri_glmfit done ninputs = 2 Checking inputs nframestot = 2 Allocing output Done allocing nframes = 2 Writing to /Users/dolphinede/Desktop/Probanden/qdec/Untitled/contrasts.sig.mgh reading colortable from annotation file... colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt) reading colortable from annotation file... colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt) gdfReadHeader: reading /Users/dolphinede/Desktop/Probanden/qdec/Untitled/y.fsgd INFO: ignoring tag Creator INFO: ignoring tag SUBJECTS_DIR INFO: ignoring tag SynthSeed INFO: NOT demeaning continuous variables Continuous Variable Means (all subjects) 0 eTIV 1.45233e+06 147517 Class Means of each Continuous Variable 1 Main 1452332.4051 MatrixMultiply: m1 is null! No such file or directory I have no idea what is wrong. Best regards Daniel Klein

2 1

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer