Freesurfer

freesurfer@nmr.mgh.harvard.edu

4 participants
30090 discussions

Qdec and GLM ROI analysis
by Jon Wieser 07 May '13

07 May '13

I am interested ing doing a group surfaces analysis like that specified in: ftp.nmr.mgh.harvard.edu/fswiki/FsTutorial/QdecGroupAnalysis, or http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/GroupAnalysis?highlight… but I want to limit the analysis to a specific ROI. how do i do this? -- Jon Wieser Research Specialist UW-Milwaukee Psychology Department, Pearse Hall Rm 366 2441 East Hartford Ave Milwaukee, WI 53211 Phone: 414-229-7145 Fax: 414-229-5219

1 0

segfault during rfx with mri_glmfit
by Alex Kell 07 May '13

07 May '13

hi freesurfers, i am running into a segfault whenever i try to use mri_glmfit to run a random effects analysis. (this is in the volume, not on the surface.) here's my call: mri_glmfit \ --y $ces_fpath \ --X $design_mat_fpath \ --C $c1_fpath \ --C $c2_fpath \ --C $c3_fpath \ --mask $mask_fpath \ --glmdir $glmdir_fpath \ --nii.gz i've tried a number of things to reduce the memory footprint of the call: -- originally i was using WLS, but i was worried about the memory demands of the cesvar, so i switched to OLS -- i introduced a GM mask -- i downsampled the ces and cesvar images (from 1 mm isotropic [cvs space] to 2 mm isotropic) any thoughts on what could be going wrong here? thanks, guys! alex ps. below are my build stamp, my info dump from mri_glmfit, and an example stderr & stdout. my fs build stamp: freesurfer-Linux-centos4_x86_64-dev-20120104 the "all info" dump from mri_glmfit: ProgramName: mri_glmfit ProgramArguments: --all-info ProgramVersion: $Name: $ TimeStamp: 2013/05/07-14:15:28-GMT BuildTimeStamp: Jan 4 2012 05:11:02 CVS: $Id: mri_glmfit.c,v 1.209 2011/12/16 19:02:28 greve Exp $ User: alexkell Machine: ba7 Platform: Linux PlatformVersion: 2.6.38-16-generic CompilerName: GCC CompilerVersion: 30400 an example stderr and stdout: $Id: mri_glmfit.c,v 1.209 2011/12/16 19:02:28 greve Exp $ cwd /mindhive/nklab2/ellison_on_nklab2/adultasd/vss/4_rfx/glm_output/F-O cmdline mri_glmfit --y /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//cvs.all.ces.F-O.runs_1234.2x2x2.nii.gz --wls /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//cvs.all.cesvar.F-O.runs_1234.2x2x2.nii.gz --X /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//design_matrix.dat --C /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//asd.mtx --C /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//tc.mtx --C /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//tc-asd.mtx --mask /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input/cvs.gmMask.2x2x2.nii.gz --glmdir /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_output//F-O --yhat-save --eres-save --save-cond --nii.gz sysname Linux hostname ba7 machine x86_64 user alexkell FixVertexAreaFlag = 1 UseMaskWithSmoothing 1 OneSampleGroupMean 0 y /mindhive/nklab2/ellison_on_nklab2/adultasd/vss/4_rfx/glm_input/cvs.all.ces.F-O.runs_1234.2x2x2.nii.gz logyflag 0 X /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//design_matrix.dat usedti 0 mask /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input/cvs.gmMask.2x2x2.nii.gz maskinv 0 glmdir /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_output//F-O IllCondOK 0 ReScaleX 1 DoFFx 0 wFile /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_input//cvs.all.cesvar.F-O.runs_1234.2x2x2.nii.gz weightinv 1 weightsqrt 1 Creating output directory /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_output//F-O Loading y from /mindhive/nklab2/ellison_on_nklab2/adultasd/vss/4_rfx/glm_input/cvs.all.ces.F-O.runs_1234.2x2x2.nii.gz Saving design matrix to /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_output//F-O/Xg.dat Normalized matrix condition is 1 Matrix condition is 2.21429 Pruning voxels by thr: 0.000000 Found 115353 voxels in mask Saving mask to /mindhive/nklab/projects/ellison/adultasd/vss/4_rfx/glm_output//F-O/mask.nii.gz search space = 922824.000000 Segmentation fault

2 2

Using FSL matrix for Talairach transform
by Johannes Klein 07 May '13

07 May '13

Dear all, I have some tricky subjects that I cannot get to register properly with the standard Talairach process, or even the use-mritotal option. However, with some fiddling, I've produced a reasonable registration matrix with FSL's FLIRT. I'm now looking to convert my FSL registration matrix, and use it as Freesurfer's Talairach transform. From the documentation, I gathered that I need to get it into Freesurfer's xfm format, and I suppose this is a two step process. I converted the FLIRT-provided matrix cor2raw.xfm like so: tkregister2 --mov /usr/local/freesurfer/average/mni305.cor.mgz --fslreg cor2raw.xfm --s mysubject --noedit --reg register.dat I then converted register.dat to the xfm data format: regdat2xfm /usr/local/freesurfer/average/mni305.cor.mgz rawavg.mgz register.dat tal.xfm and finally replaced talairach.xfm in the transforms directory with the output of the last command. Unfortunately, this doesn't work, and the transform is very much off target when I view it in tkregister2. Obviously, I've made a mistake when converting between the different transforms - could anybody point me in the right direction? Thanks all, Johannes

2 1

Probably activation due to vessels
by Jörg Pfannmöller 07 May '13

07 May '13

Dear Freesurfers, I computed the functional activation in a sensory task using freesurfer's fsfast stream. In the cortex surface line like activated regions are present. Their extension is much larger than the extension of the expected regions. Therefore, I suppose that those regions are due to vessels. Is there a possibility in freesurfer to remove them. Ideally there would be something like a volume based mode which does not sxclude the cortex. In this mode the vessels couls be identified and removed. Any help is highly apprexiated. Respectfully yours pfannmoe

2 13

Re: [Freesurfer] ROI analyses
by Douglas N Greve 07 May '13

07 May '13

On 05/07/2013 05:21 AM, Tudor Popescu wrote: > Hi Doug, > I was just wondering whether you received my last email, with > follow-ups on your answers. Sorry for the multitude of questions - as > ever, I really appreciate the help! > Best wishes, > Tudor > > On 4 May 2013 20:31, Tudor Popescu <tudor3(a)gmail.com > <mailto:tudor3@gmail.com>> wrote: > > Sorry Doug, wasn't sure whether you'd seen my reply > Tudor > > > ---------- Forwarded message ---------- > From: *Tudor Popescu* <tudor3(a)gmail.com <mailto:tudor3@gmail.com>> > Date: 2 May 2013 21:54 > Subject: Re: [Freesurfer] ROI analyses > To: freesurfer(a)nmr.mgh.harvard.edu > <mailto:freesurfer@nmr.mgh.harvard.edu> > > > Thanks Doug! If I may follow up: > > > 3) (I can't see whether there's a tutorial for this, as the > main FS > > site seems to be down at the moment) How can a group > comparison of CT > > be done inside a given ROI, and in what ways can the ROI be > specified? > > I guess the command line version of QDEC would produce, in > this case, > > the same contrasts as QDEC, but within the ROI as opposed to > the whole > > brain > You mean you want to do an exploratory analysis within a mask > of an ROI > or that you want to average the CT scan within the ROI and do > a group > analysis of the ROI values? If the former, you can use > mri_glmfit with > the -label or -mask option. If the latter, you can create a > table with > asegstats2table then run mri_glmfit with --table > > 3') Is it the case that mri_glmfit with mask=brainmask.mgz is > equivalent to doing a whole-brain analysis in QDEC? > QDEC is only doing surface-based analysis. brainmask.mgz is a volume. But in principle the idea is similar. > > > 3'') I didn't find, on the mri_glmfit help page, how masks can be > formed out of one of the atlas structures, e.g. > "lh_S_intrapariet_and_P_trans "... The --table argument is also > not present anywhere on that page.. > You can run mri_annotation2label to convert an annotation to label, then pass the label to glmfit. > > > 3''') Would it not be a/parc/stats2table rather than > a/seg/stats2table, if the ROI in question was cortical and not > subcortical? Is the difference between what FS calls > "segmentation" and "parcellation" just terminology, i.e. the > former delineates subcortical structures and the latter delineates > cortical structures? Or does it imply anything else? > Whether you use one or the other depends on whether you generate the individual table with mris_anatomical_stats (aparcstats2table) or mri_segstats (asegstats2table) > > > 4) How do the results of these two different types of ROI group > > analyses differ, and is one of them more "correct" than the > other: > > A) running the command prompt version of QDEC within the > confines > > of a certain atlas-defined ROI, and looking at the resulting > > statistical map (clusters), as per question 3; > > B) extracting the CT values for that ROI for all > subjects using > > aparcstats2table, and doing t-tests to look for a group > difference. > Oops, looks like these are the two I mentioned from #3 above. > The first > is an exploratory analysis in which the groups are compared on a > vertex-by-vertex basis. If there is a subset of vertices that are > different between the groups, it may show up in the exploratory > analysis. However, the effect may be small at each vertex and > averaging > over the vertices may improve your power (unless the effect is > only at a > few vertices). One is not more correct than the other, just > testing > different hypotheses. > > 4') Can the second option (averaging the CT values within the ROI) > not also be done by taking the values from the table into SPSS and > simply doing an independent samples t-test? Is that not equivalent > to what mri_glmfit would do in this case? > Yes, it should be the same. > > > And one more question: > 5) Can the statistical map in QDEC be built using a lower, > uncorrected p value (e.g. 0.0001) instead of an FDR-corrected test > value? > Yes, just set the threshold (FDR is off by default). doug > > > Thanks!! > Tudor > > doug > > > > > _______________________________________________ > > 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 > Fax: 617-726-7422 > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > Outgoing: > ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > 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: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

1 0

Re: [Freesurfer] Issues with Cluster generation
by Douglas N Greve 07 May '13

07 May '13

Hi Joshua, there are several issues. First, you are using a threshold of p<.0001 which is pretty severe. It works for the tutorial, but you need to understand what this threshold is and pick one appropriate for your data. Second, you are only doing 5 iterations. You need to be doing something like 10000. If you are doing a whole cortex analysis, you can use the cached table instead of doing a simulation (with --qache). Running your own simulation will take a long time. The smoothing level is a little high, I usually run at 10 or 15. There's nothing wrong, it is just a little high. doug On 05/05/2013 03:39 AM, Joshua deSouza wrote: > Hi Doug, > > The image has undergone fwhm 20 smoothing. > > The final glmdir file is given in the attachments. RC.txt was the > contrast name and folder in which the cluster correction files are kept. > > The glmdir creation command line : mri_glmfir --y lh.thomasR.20.mgh > --fsgd R.txt dods --C RC.txt --surf fsaverage lh --cortex --glmdir > lh.thomasR.glmdir > > The glmfit-sim command line is : mri_glmfit-sim --glmdir > lh.thomasR.glmdir --sim mc-z 5 4 mc-z.negative --sim-sign neg > --cwpvalthresh 0.999 --overwrite > > The command line to view the image is : tksurfer fsaverage lh inflated > -annot lh.thomasR.glmdir/RC.txt/mz-z.negative.sig.ocn.annot fthresh 1 > -curv -gray > > The version of freesufrer i'm using is 1.379.2.17 > > Thank you helping me out with this. > > > On Sat, May 4, 2013 at 12:12 AM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > Can you send the command line and terminal output and the version > of FS > you are using? > doug > On 05/03/2013 09:01 AM, Joshua deSouza wrote: > > Hello all, > > > > I am facing a problem whereby after running GLMfit and cluserwise > > correction for multiple comparisons using the monte-carlo method, I > > was unable to generate any clusters. I do have significant images > > produced for the glmfit. > > > > All procedures were done according to the tutorial to match the > > effects under FSGD. > > > > There are no clusters available even when I set the cwpvalthresh to > > .999 which is supposed to show all clusters. What is the cluster > > generation dependent upon? As my glmfit does show significant areas. > > > > Please advise > > > > > > _______________________________________________ > > 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: https://gate.nmr.mgh.harvard.edu/filedrop2 > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > Outgoing: > ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > 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: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

1 0

Segmentation fault (core dumped) again in the BETA version
by martina.andellini＠fastwebnet.it 07 May '13

07 May '13

Hi Freesurfer experts, I've got Ubuntu version 13.04 installed on my computer, and after I saw the post : "http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg28267.html" I d ecide to install the BETA version freesurfer-Linux-centos6_x86_64-stable-pub-v5. 3.0-BETA.tar.gz as adviced. I still run into the same error when I run mri_nu_correct.mni and tkmedit: Segmentation fault (core dumped) ERROR: converting to minc and =================================== ERROR: A segfault has occurred. This is not your fault, : but is most likely an unrecoverable error and has : made the program unstable. : : Please send the contents of the file .xdebug_tkmedit : that should be in this directory to freesurfer(a)nmr.mgh.harvard.edu : : Now exiting... : Can anyone please advice on that? thanks a lot Martina

2 1

Neuroimaging Research Technician
by Emily Lindemer 07 May '13

07 May '13

* * * * ** *A Neuroimaging Research Technician position is available in the VA Boston Healthcare System (VABHS) Neuroimaging Research for Veterans (NeRVe) Center. * * The Neuroimaging Research Technician will provide technical support for image preprocessing and analysis (using programs such as FSL and FreeSurfer), maintenance and development of in-house custom processing codes (script writing in unix/linux, matlab, etc), and assist staff members in the conduct of neuroimaging research as well as administrative tasks related to the Center. This includes developing expertise and training other investigators with the operation of current MRI hardware (scanner console, goggles, headphones, BioPac system, etc) and software for data processing and analysis. Additional responsibilities will include management of online documentation and center resources. Technical skills and familiarity with Linix/Unix, Matlab, and imaging software packages such as FSL and FreeSurfer is beneficial. Applicants should have a minimum of a Bachelor’s degree, with a background in Neuroscience, Computer Science, Engineering, Math or Physics. A two year commitment is required. Must be a U.S. citizen. Interested candidates should email a cover letter, resume, and contact information for two professional references to Emily Lindemer at emily.lindemer(a)va.gov . *

1 0

NFS v4 on CentOS 5.4 or Red Hat 6.2
by Garikoitz Lerma-Usabiaga 07 May '13

07 May '13

Hi FS experts, we have changed our file systems to NFS v4 and we found that FS cannot write to it (I mean, FS installed in an app server and subjects in the file system is not working). Last new I have is that FS in CentOS could read but not write and in Red Hat could not write or read. Is it a known problem? Is it something going on with FS or the problem is in our end? (the rest of the programs are working fine: FSL, Matlab...) Many thanks, Gari

1 0

Qdec analysis for only particular regions of interest?
by Klara Mareckova 07 May '13

07 May '13

Hello, I'm new to Freesurfer and I was wondering whether anybody could tell me what steps I need to do in Freesurfer in order to conduct the following analysis? I would like to know which sex differences in cortical thickness (subjects 12-18 years old) are present in the 12 year olds already. I've used qdec to calculate the main effect of sex in a) all participants (12-18) and b) in 12 year olds. Is there a way how to mask the 12-year olds' results with the whole group so that I can calculate which of the whole-group main effects are present before the age of 13? Thanks a lot for any suggestions, Klara

1 0

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