Freesurfer

freesurfer@nmr.mgh.harvard.edu

16 participants
30084 discussions

Hippocampus subfields outputs missing
by Muhammad Naveed Iqbal Qureshi 06 Dec '19

06 Dec '19

External Email - Use Caution Hi I ran the following command not the subject on HPC module load matlab/r2013a module load freesurfer/6.0.0 source $FREESURFER_HOME/SetUpFreeSurfer.sh fixup_mni_paths recon-all -all -no-isrunning -s /hpf/projects/glongoni/ADHD-200/$file -subjid ADHD-${PBS_ARRAYID} -hippocampal-subfields-T1 -sd /hpf/projects/glongoni/ADHD-200/ -openmp 8 in the output directory I get only following files for both left and right hemisphere lh.hippoSfLabels-T1.v10.FSvoxelSpace.mgz lh.hippoSfLabels-T1.v10.mgz lh.hippoSfVolumes-T1.v10.txt How can I get following files? lh.hippoAmygLabels-T1.v21.CA.FSvoxelSpace.mgz lh.hippoAmygLabels-T1.v21.CA.mgz lh.hippoAmygLabels-T1.v21.FS60.FSvoxelSpace.mgz lh.hippoAmygLabels-T1.v21.HBT.FSvoxelSpace.mgz I try to find this command but it is not available in my installation segmentHA_T1.sh Regards, Naveed

1 0

4-year PhD studentship ENIGMA Antisocial Behavior
by Stephane De Brito 06 Dec '19

06 Dec '19

External Email - Use Caution Dear all, We have a 4-year fully funded PhD studentship available to work on our ENIGMA working group focusing on different antisocial phenotypes. https://www.findaphd.com/phds/project/enhancing-neuroimaging-genetics-throu… More information about ENIGMA available here http://enigma.ini.usc.edu/ongoing/enigma-antisocial-behavior/ and about the ENIGMA Antisocial Behavior working group here http://enigma.ini.usc.edu/ongoing/enigma-antisocial-behavior/ Best wishes, Stephane --------------------------------------------------- Stephane De Brito, PhD Senior Lecturer Centre for Human Brain Health School of Psychology University of Birmingham, B15 2TT, UK Tel: +44(0)121 414 2641 (internal: 42641) Robert Aitken Building, Room 337a Web: Dr Stephane De Brito - School of Psychology - University of Birmingham<http://www.birmingham.ac.uk/staff/profiles/psychology/de-brito-stephane.aspx> Current projects: http://www.femnat-cd.eu/ and ENIGMA Antisocial Behavior<http://enigma.ini.usc.edu/ongoing/enigma-antisocial-behavior/> Twitter: @SCAN_Lab<https://twitter.com/scan_lab>

1 0

Postdoctoral Researcher in Neuroimaging at Karolinska Institute
by Joana Braga Pereira 06 Dec '19

06 Dec '19

External Email - Use Caution Dear all, We are looking for a Postdoc to work in an interdisciplinary and internationally competitive environment at Karolinska Institute, one of the top medical universities in the world and the home for the Nobel Prize in Medicine. We combine imaging analyses, statistical tools, machine learning algorithms, biofluid assays and genetics to study brain aging and the earliest signs of Alzheimer’s and Parkinson’s diseases. This project will be carried out under the supervision of Assistant Professor Joana B. Pereira and it will involve analyzing cross-sectional and longitudinal imaging data, combining clinical measures and cerebrospinal fluid biomarkers, as well as writing manuscripts and presenting results in international conferences. The fellow is expected to be independent, develop his or her own projects, supervise students and show an interest in approaches based on graph theory. In addition, he/she will have the opportunity to develop new skills, analyze large datasets, apply for grants and develop a wide network of collaborations. The ideal candidate should have the following qualities: - A PhD in neuroimaging, computational science, neuropsychology, physics, statistics or a related field; - Experience in statistical techniques, programming, or scripting; - Experience with structural MRI, DTI, fMRI or PET imaging analyses; - Clinical or research experience with aging and/or neurodegenerative disorders. - Excellent record of scientific publications. The Postdoctoral studies will be performed for 1 year, with possible extension for an additional year. There are also possibilities for a longer appointment after this 2-year period. To apply please go to: https://ki.varbi.com/what:job/jobID:305265/ For further information or general inquires, please contact Joana B. Pereira (joana.pereira(a)ki.se) Best wishes, Joana

1 0

topological defect error
by Laurel Quinlan 05 Dec '19

05 Dec '19

External Email - Use Caution Hi Freesurfer Support, One of my colleagues receives the following error after running a recon: exited with errors: "excessive topologic defect encountered: could not allocate 289309485 edges for retessellation" I checked the inflated view and there does appear to be significant portions of the brain missing (however I cannot see this in 2D). I've attached a screenshot of the inflated view. How would I go about fixing the error? Thank you! Laurel Quinlan Associate Research Specialist Center for Healthy Minds University of Wisconsin-Madison

3 6

Reporting Results in a three group analysis
by cody samth 05 Dec '19

05 Dec '19

External Email - Use Caution Hi, I have a statistical question about how to approach reporting results from FreeSurfer analyses containing three groups. I ran a group effect (F-test) and then post-hoc tests looking at pair-wise comparisons between the three groups. My question is why is that we run separate vertex-wise analyses for the post-hocs rather than extracting the values from significant clusters and running post-hocs in a statistical software for just the regions where a significance difference was found (ie the clusters)? As a post-hoc vertex wide analysis can lead to different results. For example in my group effect contrasts I found a cluster in the parietal lobe. Whereas in my post hoc-vertex wide analyses one of them found two clusters 1) within the parietal 2) within the frontal lobe. If I choose to run these post-hoc analyses via freesurfer (ie vertex-wise) rather than extracting the results to analyze in a statistical program, is it standard to report the second cluster? Even though it didn't come up in the group model? If so is there a paper that people reference that uses this approach? Thanks, Cody

2 4

surfaceholes in aseg.stats
by Laurel Quinlan 05 Dec '19

05 Dec '19

External Email - Use Caution Hi Freesurfer Support, I'm wondering how ?h.surfaceholes prior to fixing and "Total number of defect holes in surfaces prior to fixing" are calculated in the aseg.stats file for each brain. I compared aseg.stats files for the same brain edited by two different people, and I noticed that there were different numbers of holes listed for each person, despite the brain being the same. For example: lhSurfaceHoles 130 vs 61 rhSurfaceHoles 102 vs 80 Total number of defect holes 232 vs 141 Why would those numbers be so different? Thanks, Laurel Quinlan Associate Research Specialist Center for Healthy Minds University of Wisconsin-Madison

3 3

Multiple comparison using Monte Carlo
by Ting Li 05 Dec '19

05 Dec '19

External Email - Use Caution Dear Freesurfer Expert: I have used the simulation code as show below: mri_glmfit-sim \ --glmdir lh.gender_age.glmdir \ --cache 4 neg \ --cwp 0.05\ --2spaces My question is how many iterations are done during this simulation? I didn’t specify the iteration number. Why I didn’t have the csd file? Thanks a lot for your response! Best regards, Ting

2 11

negative SUVr using PETsurfer
by miracle ozzoude 05 Dec '19

05 Dec '19

External Email - Use Caution Hello Experts, I ran the PETsurfer pipeline using AV45 pet and performed pvc using pons as reference point. When i looked at the gtm.stats.dat file for one of my subjects, the SUVr for right and left nucleus-accumbens areas were in the negative (left = -0.361; right = -0.046). This is strange because SUVR shouldn't be negative. How do i go about fixing this error? Below, is my pvc command. mri_gtmpvc --i pet.nii --reg register.lta --psf 6 --seg gtmseg.mgz --default--seg-merge --auto-mask PSF .01 --mgx .01 --o subject_gtmpvc --rescale 174 Thanks. Best, Paul

2 3

T2 in recon-all: command to correct pial surfaces errors
by Steve Petersen 05 Dec '19

05 Dec '19

External Email - Use Caution Dear Freesurfer experts, I used t2 image as input in subjects that have already been processed without them following the instructions of the website. Although pial surface segmentation has improved there is still some error. After edit the brainmask.mgz file in Freeview (recon edit option), what command should I use to correct the pial surfaces obtained after the -T2pial flag? I tried to use the -autorecon-pial command to correct the edits but I think that this command only takes into account the T1, is that right? Best regards,

2 1

Re: [Freesurfer] Question on GLM - correlation
by Greve, Douglas N.,Ph.D. 05 Dec '19

05 Dec '19

That does not have the pvr in it. You can extract the pvr for that vertex and add it as a column to the X matrix with pvr = fast_vol2mat(MRIread('pvr.mgz')); pvrvtx = pvr(:,vertexno+1); ps. Please remember to post to the list On 12/3/2019 1:29 PM, Barletta, Valeria wrote: Ok done. This is the Xg.mat 1 29 1 44 1 38 1 47 1 34 1 25 1 30 1 22 1 52 1 31 1 37 1 42 1 28 1 39 1 44 1 44 1 32 1 32 1 36 1 47 1 52 And the cond(Xg'*Xg) is again: 3.0555e+04 Looks like there is no extra column for the pvr in the Xg.mat ________________________________ From: Greve, Douglas N.,Ph.D. <DGREVE(a)mgh.harvard.edu><mailto:DGREVE@mgh.harvard.edu> Sent: Tuesday, December 3, 2019 11:03 AM To: Barletta, Valeria <VBARLETTA(a)mgh.harvard.edu><mailto:VBARLETTA@mgh.harvard.edu>; Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu><mailto:freesurfer@nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] Question on GLM - correlation I'm not sure those condition numbers include the pvr. Try loading the Xg.mat into matlab and computing cond(Xg'*Xg) Make sure Xg has the extra column for the pvr On 12/3/2019 10:42 AM, Barletta, Valeria wrote: Ok it worked. So I tried with a non-zero vertex and I got: Normalized matrix condition is 81.1382 Matrix condition is 30558.5 Found 149926 points in label. Found 149926 voxels in mask Saving mask to /autofs/space/lipari_001/users/valeria/NODDI/GLM//GLM_CME_ficvf_doug/rh_contrast-MS_CME_all_FIS_ICA1/mask.mgh Reshaping mriglm->mask... search space = 74490.928733 DOF = 18 Dumping voxel 1024 0 0 to /autofs/space/lipari_001/users/valeria/NODDI/GLM//GLM_CME_ficvf_doug/rh_contrast-MS_CME_all_FIS_ICA1/voxdump-1024-0-0 ...and with a zero vertex and got: Normalized matrix condition is 81.1382 Matrix condition is 30558.5 Found 149926 points in label. Found 149926 voxels in mask Saving mask to /autofs/space/lipari_001/users/valeria/NODDI/GLM//GLM_CME_ficvf_doug2/rh_contrast-MS_CME_all_FIS_ICA1/mask.mgh Reshaping mriglm->mask... search space = 74490.928733 DOF = 18 Dumping voxel 150252 0 0 to /autofs/space/lipari_001/users/valeria/NODDI/GLM//GLM_CME_ficvf_doug2/rh_contrast-MS_CME_all_FIS_ICA1/voxdump-150252-0-0 ________________________________ From: freesurfer-bounces(a)nmr.mgh.harvard.edu<mailto:freesurfer-bounces@nmr.mgh.harvard.edu> <freesurfer-bounces(a)nmr.mgh.harvard.edu><mailto:freesurfer-bounces@nmr.mgh.harvard.edu> on behalf of Greve, Douglas N.,Ph.D. <DGREVE(a)mgh.harvard.edu><mailto:DGREVE@mgh.harvard.edu> Sent: Monday, December 2, 2019 2:56 PM To: freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu> <freesurfer(a)nmr.mgh.harvard.edu><mailto:freesurfer@nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] Question on GLM - correlation oh, sorry, did not see the pvr. It could very well be the pvr that is causing the patchiness. Unfortunately, you have to debug it vertex by vertex. I would choose a few vertices (one with signal, on without) and run mri_glmfit with --voxdump vertexno 0 0 This will create a subfolder with the GLM for that vertex. I'm not sure what to tell you to look for, but I would at least check the condition numbers in each vertex. What is the nature of ficvf? On 12/2/2019 1:22 PM, Barletta, Valeria wrote: My other continuous variable is --pvr ${hemi}.ficvf.mgh Not sure about --no-prune, I found it in my script as a default... Shall I remove it? ________________________________ From: freesurfer-bounces(a)nmr.mgh.harvard.edu<mailto:freesurfer-bounces@nmr.mgh.harvard.edu> <freesurfer-bounces(a)nmr.mgh.harvard.edu><mailto:freesurfer-bounces@nmr.mgh.harvard.edu> on behalf of Barletta, Valeria <VBARLETTA(a)mgh.harvard.edu><mailto:VBARLETTA@mgh.harvard.edu> Sent: Friday, November 29, 2019 3:08 PM To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu><mailto:freesurfer@nmr.mgh.harvard.edu> Subject: [Freesurfer] Question on GLM - correlation Dear Freesurfers, I am running glmfit to assess for correlation between y and another continue variable across the cortex (ficvf), correcting for the effect of age, using the command: mri_glmfit \ --y {surfname_cme}.mgh \ --fsgd {contrast}.fsgd \ --C {contrast}.mat \ --pvr ${hemi}.ficvf.mgh" \ --surf fsaverage $hemi \ --glmdir ${glm_dir}/${hemi}_contrast-${contrast}_${surfname_cme} \ --no-prune This is my matrix: 1 0 0 And my fsgd file: Title CME Class MS Variable age Input MS069_R01_01 MS 29 Input MS070_R01_02_R1year MS 44 Input MS072_R01_04 MS 38 Input MS073_R01_05 MS 47 Input MS075_R01_07_R1year MS 34 Input MS076_R01_08 MS 25 Input MS077_R01_09 MS 30 Input MS078_R01_10 MS 22 Input MS080_R01_12 MS 52 Input MS081_R01_13 MS 31 Input MS082_R01_14 MS 37 Input MS084_R01_16 MS 42 Input MS086_R01_18 MS 28 Input MS087_R01_19_R1year MS 39 Input MS088_R01_20 MS 44 Input MS089_R01_29 MS 44 Input MS090_R01_21 MS 32 Input MS091_R01_22 MS 32 Input MS092_R01_23 MS 36 Input MS094_R01_25 MS 47 Input MS104_R01_36 MS 52 This is the result: [cid:part10.AAFB4919.2AE170D2@mgh.harvard.edu][cid:part11.F1A0EFEC.57C56985@mgh.harvard.edu] Am I proceeding correctly? Thank you, Valeria _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

2 4

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