 |
 |
Ok, hopefully third time is a charm. So I've been sorting through the FSL list archives, and I *think* I get what I'm supposed to do to threshold my fdt_paths output. I'm interested in looking at within subject as well as between subject differences. Could someone please tell me if this is correct? For my command "fslmaths fdt_paths -thr 50 -bin -mul dti_FA_str fdt_paths_values" where fdt_paths is the result of the multiple masks command in probtrackx (2 masks as seeding points), I should be setting the -thr for each subject as a set percentage of the total number of voxels in both ROIs, ie. (# of ROI1 voxels + # of ROI2 voxels) x percentage. Thanks! Kristen -----Original Message----- From: FSL - FMRIB's Software Library on behalf of Lindgren, Kristen, Ann Sent: Fri 3/27/2009 8:44 AM To: [log in to unmask] Subject: Re: [FSL] Post probtrackx analyses - thresholds for fdt_paths Whoops, so as followup to my question, I should note that I'm trying to figure out what was meant in the FSL list thread #017802 and its replies. I've looked at the Rilling (2008) paper (trying to get a copy of Matt Glaser's 2008 paper) and read Matt's response: "I think two things are important: 1) That you use a consistent percentage of the total number of samples sent out for each tract across subjects (so that tracts created by larger ROIs, and thus have more total samples sent out, have higher thresholds) and 2) That in setting the percentage that you will use, you try a variety of values and see what seems to produce clean results showing the pathway of interest without many extraneous pathways not clearly connected to the ROI but at the same time does not remove large parts of the pathway of interest (because it is too high)." I used multiple masks for my probtrackx analyses (2 masks only, no waypoints), so would I still find the total number of samples sent out for each tract in the waytotal file? If so, which number is that? There are two numbers in mine. Is the first number the number sent from ROI1 and the second is the number sent from ROI2, in which case I would just add them and make my threshold a consistent percentage of that number? I guess I'm just confused as to how to implement in my data what Matt recommended in the quote above. It seems that Matt found thresholds of 2-6% of total samples sent out worked well for his dataset using structural space seeding points (similar to my study). So if I'm understanding correctly, I would set the threshold for each subject for fdt_paths in fslmaths as a set percentage of the sum of the two values in waypoint.txt. I'm still trying to get a hold of Matt's DTI tractography of language pathways paper to see what he used as a threshold, so any guidance on where to start would be great. Thanks so much! Kristen -----Original Message----- From: FSL - FMRIB's Software Library on behalf of Lindgren, Kristen, Ann Sent: Thu 3/26/2009 10:56 PM To: [log in to unmask] Subject: [FSL] Post probtrackx analyses - thresholds for fdt_paths I was looking through some e-mails that I had saved about thresholding with probtrackx, but I'm still unsure if I'm thresholding my fdt_paths outputs properly so I wanted to double check one last time before committing to my analyses. Any thoughts, concerns, or advice are greatly appreciated. Here's what I've done: 1) I used two freesurfer gray matter parcellation labels as seeding points in probtrackx for each subject. There is a known connection between these two regions. 2) Now I'm pulling out the average FA and volume of the resulting fdt_paths using these formulas (thanks again to Steve for helping me figure out this one!) fslmaths fdt_paths -thr 50 -bin -mul dti_FA_str fdt_paths_values fslstats fdt_paths_values -M -V So here are my concerns: 1) I've been using 50 as my threshold for each subject (used the default value of 5000 samples for probtrackx). Is this correct? Or should I be customizing this per subject since the volumes of the seeding points vary across subjects? If I should be varying it, what should I be basing my thresholds on for each subject? Any references in support of either method? 2) Should I only be including voxels that have an FA value above a certain threshold? If so, anyone know of a paper that has looked at the normal range of FA values for white matter? Sorry for the ton of questions lately, I'm just in the last stretch of my dissertation analyses and I want to make sure I'm doing things correctly before moving forward. Thanks so much for all of your help. I really appreciate it! Kristen