<div dir="ltr">Hi David,<div><br></div><div>So I guess this is where the issue of quantification becomes a matter of debate. As was shown quite comprehensively by <a href="http://www.ncbi.nlm.nih.gov/pubmed/22846632">Derek Jones et al</a>, streamlines counts are not in themselves particularly meaningful - and hence can't be viewed as quantitative. There is also evidence that they're not particularly reproducible (at least when using probabilistic streamlines without an enormous number of streamlines) - see e.g. <a href="http://www.ncbi.nlm.nih.gov/pubmed/24246491">our TDI paper on the issue</a>. This is what I meant when I said there isn't much quantification you can do in MRtrix 0.2. You have a lot more options in MRtrix3, including <a href="http://www.ncbi.nlm.nih.gov/pubmed/22705374">ACT</a>, <a href="http://www.ncbi.nlm.nih.gov/pubmed/23238430">SIFT</a>, <a href="http://www.ncbi.nlm.nih.gov/pubmed/26163802">SIFT2</a>, and <a href="http://www.ncbi.nlm.nih.gov/pubmed/22036682">AFD</a> - SIFT being particularly <a href="http://www.ncbi.nlm.nih.gov/pubmed/25312774">beneficial for connectomic analyses</a>. </div><div><br></div><div>As to your question: I don't think these minor inaccuracies in the responses are going to make much difference at all to the tractography itself. These are minor noise ripples in the low signal region of the response - what makes a big difference is the overall size of the response (this influences the scaling of the FOD), and how broad it is. But even if it does have a minor effect on the FOD, this is unlikely to affect the tractography a great deal, and is even less likely to influence any comparisons you might make between subjects - I don't see how this could introduce a bias between your groups, at worst you might lose a negligible amount of statistical power. This of course all assumes that the procedure you're using is the same for all subjects.</div><div><br></div><div>So I would run at lmax=8 without worrying too much about it... If you were to switch to use some of the methods in MRtrix3, then this would indeed be a different matter.</div><div><br></div><div>Cheers,</div><div>Donald.</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 1 April 2016 at 18:02, David Grayson <span dir="ltr"><<a href="mailto:dgrayson@ucdavis.edu" target="_blank">dgrayson@ucdavis.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
Hey Donald,<br>
<br>
I am in fact doing quantitative analyses of number of pobabilistic streamlines connecting a given pair of regions, across subjects. So I would prefer to minimize variation in the rf across subjects.<br>
<br>
>From your experience, would you expect one of these rf's to be more consistent across subjects?<br>
<hr style="display:inline-block;width:98%">
<div dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> J-Donald Tournier <<a href="mailto:jdtournier@gmail.com" target="_blank">jdtournier@gmail.com</a>><br>
<b>Sent:</b> Friday, April 1, 2016 12:15:30 AM<br>
<b>To:</b> David Grayson<br>
<b>Cc:</b> <a href="mailto:mrtrix-discussion@www.nitrc.org" target="_blank">mrtrix-discussion@www.nitrc.org</a><div><div class="h5"><br>
<b>Subject:</b> Re: [Mrtrix-discussion] Response functions at lmax=6 vs lmax=8</div></div></font>
<div> </div>
</div><div><div class="h5">
<div>
<p dir="ltr">Hi David,</p>
<p dir="ltr">Yes, these look fine. That's often a little protuberance like this in the response, it'll be due to some extent to the Rician bias driving the signal away from zero in this region. It's a minor inaccuracy that shouldn't make a great deal of difference
to the output of CSD. In fact, some people reckon it provides a better fit to the actual signal (since it contains Rician bias anyway), and yields sharper FODs... Not entirely sold on that one, but just to say that it really shouldn't impact on any downstream
inferences. </p>
<p dir="ltr">On that note, these subtle differences and inaccuracies in the response may conceivably make a difference if you're doing something quantitative with the tracks themselves across subjects. There's really very little you can do on that front with
MRtrix 0.2. If you'd been using MRtrix3, there would be many more options available for quantification, at which point you'd need to worry about this. But even then, the main concern is to make sure the same response is used across subjects, rather than whether
the response is perfect...</p>
<p dir="ltr">Hope this helps,<br>
Donald</p>
<p dir="ltr">--<br>
Dr J-Donald Tournier (PhD)</p>
<p dir="ltr">Senior Lecturer, Biomedical Engineering<br>
Division of Imaging Sciences & Biomedical Engineering<br>
King's College London</p>
<p dir="ltr">A: Department of Perinatal Imaging & Health, 1st Floor South Wing, St Thomas' Hospital, London. SE1 7EH<br>
T: +44 (0)20 7188 7118 ext 53613<br>
W: <a href="http://www.kcl.ac.uk/medicine/research/divisions/imaging/departments/biomedengineering" target="_blank">http://www.kcl.ac.uk/medicine/research/divisions/imaging/departments/biomedengineering</a></p>
<div class="gmail_quote">On 1 Apr 2016 4:39 am, "David Grayson" <<a href="mailto:dgrayson@ucdavis.edu" target="_blank">dgrayson@ucdavis.edu</a>> wrote:<br type="attribution">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Thanks for the responses here.</p>
<p><br>
</p>
<p>I should have mentioned, I have been using the freesurfer WM mask (co-registered to the B0) to exclude noisy voxels. Looks like there was an issue with co-registration. Now that this is fixed, indeed the lmax=8 response looks (very slightly) sharper than
lmax=6.</p>
<p><br>
</p>
<p>However, I'm wondering if you can have a look at the following comparisons of lmax=6 vs lmax=8, using an eroded WM mask and using the full WM mask. The lmax=8 still produces the protruberance on the z-axis, which is exaggerated using the eroded WM mask.
Does this suggest that I should perhaps stick with lmax=6, and avoid eroding any further the WM mask? Would you consider these RF's acceptable now?</p>
<p><br>
</p>
<p><a href="http://imgur.com/42NCktt" title="http://imgur.com/42NCktt
Cmd+Click or tap to follow the link" target="_blank">http://imgur.com/42NCktt</a><br>
<br>
</p>
<p><span style="font-size:12pt">Thank you very much,</span><br>
</p>
<p>David</p>
<br>
<br>
<div style="color:rgb(0,0,0)">
<hr style="display:inline-block;width:98%">
<div dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b>
<a href="mailto:albnet@gmail.com" target="_blank">albnet@gmail.com</a> <<a href="mailto:albnet@gmail.com" target="_blank">albnet@gmail.com</a>> on behalf of Dorian P. <<a href="mailto:alb.net@gmail.com" target="_blank">alb.net@gmail.com</a>><br>
<b>Sent:</b> Tuesday, March 22, 2016 3:10 PM<br>
<b>To:</b> J-Donald Tournier<br>
<b>Cc:</b> David Grayson; <a href="mailto:mrtrix-discussion@www.nitrc.org" target="_blank">
mrtrix-discussion@www.nitrc.org</a><br>
<b>Subject:</b> Re: [Mrtrix-discussion] Response functions at lmax=6 vs lmax=8</font>
<div> </div>
</div>
<div>
<div dir="ltr">I always used to correct manually the single fiber mask. The automatic thresholding includes out of brain voxels in inferior slices. A simple threshold is not to be trusted, I think.
<div><br>
</div>
<div>Dorian</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Tue, Mar 22, 2016 at 6:01 PM, J-Donald Tournier <span dir="ltr">
<<a href="mailto:jdtournier@gmail.com" target="_blank">jdtournier@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Hi David,
<div><br>
</div>
<div>I've never heard of this being referred to as a 'nubbin', but it's clear enough... ;)</div>
<div><br>
</div>
<div>Looking at your responses though, the 'nubbin' isn't all that worrying. What's more worrying is the fact that the lmax=6 response is sharper than the lmax=8 response, which seems wrong. Your lmax=8 response is certainly broader than I'd expect at b=3000. </div>
<div><br>
</div>
<div>So all this suggests that the response function estimation isn't working all that well, which is typically a symptom of a poor single-fibre mask. You used an FA threshold of 0.7 here, but bear in mind that this is by no means a 'default' - this is more
of a guideline, <a href="http://jdtournier.github.io/mrtrix-0.2/tractography/preprocess.html#csd" target="_blank">
as stated in the documentation</a>. Here's the relevant excerpt:</div>
<div><br>
</div>
<div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Note that this value is a guide only - feel free to use a different value if this does not produce satisfactory results. Ideally, you should now have a mask containing a few hundred voxels, all located within high FA white matter regions. <b>It is very important
to check that the single-fibre mask is suitable, as otherwise the response function produced in the following step may be totally inappropriate, which would seriously affect the quality of the CSD output</b>. If needed, you can edit this mask image to remove
unwanted voxels using the <a href="http://jdtournier.github.io/mrtrix-0.2/general/mrview.html#roi" style="text-decoration:none;color:teal" target="_blank">ROI analysis</a> sidebar tool within <a href="http://jdtournier.github.io/mrtrix-0.2/general/mrview.html" style="text-decoration:none;color:teal" target="_blank">MRview</a>.</blockquote>
</div>
<div><br>
</div>
<div>I would have a good look at the single-fibre mask used in this step, and check against the above. The chances are you might find a lot of noisy high-FA edge voxels got included or something, despite the erosion step. This all depends on how good the initial
brain mask was (I often found it difficult to exclude the nasal sinuses, for example). Unfortunately, this step can be a bit fiddly. </div>
<div><br>
</div>
<div>Alternatively, you could upgrade to <a href="http://www.mrtrix.org" target="_blank">
MRtrix3</a> - we've put in quite a bit of work on <a href="http://mrtrix.readthedocs.org/en/latest/concepts/response_function_estimation.html" target="_blank">estimating the response function</a>...</div>
<div><br>
</div>
<div>All the best,</div>
<div>Donald.</div>
<div><br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">
<div>
<div>On 22 March 2016 at 19:47, David Grayson <span dir="ltr"><<a href="mailto:dgrayson@ucdavis.edu" target="_blank">dgrayson@ucdavis.edu</a>></span> wrote:<br>
</div>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div>
<div lang="EN-US">
<div>
<p class="MsoNormal">Hi MRtrixers,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have HARDI datasets on a large sample of young adults (7-16yrs) taken with 72-dir and b0=3000. I am wondering whether I should proceed with lmax=6 or 8. I am using MRtrix 0.2.12.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I am concerned about using lmax=8 because in the response function there is often a ‘nubbin’ that appears right in the center of the Z-axis line that doesn’t appear at lmax=6. The below link shows an example of the RF’s generated for a
single subject taken at lmax=6 (on the left) and lmax=8 (on the right). These are generated using default parameters (FA threshold of 0.7 for the single-fiber-orientation mask).<u></u><u></u></p>
<p class="MsoNormal"><a href="http://imgur.com/a/ODrcf" target="_blank">http://imgur.com/a/ODrcf</a><u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Hopefully it’s obvious what I’m talking about. Sometimes this nubbin is more pronounced than what I’m showing here. I’m not sure how concerned I should be about this, but based on my intuition about what the RF means, the Z-axis itself
should represent the global minimum of the DW signal, which is obviously not the case with lmax=8. Should I avoid any appearance of this effect at all costs, or is there some rule of thumb for how big it can be before I should be concerned about it?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks,<span><font color="#888888"><u></u><u></u></font></span></p>
<span><font color="#888888">
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">David<u></u><u></u></p>
</font></span></div>
</div>
<br>
</div>
</div>
_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">Mrtrix-discussion@www.nitrc.org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/mrtrix-discussion" rel="noreferrer" target="_blank">http://www.nitrc.org/mailman/listinfo/mrtrix-discussion</a><br>
<br>
</blockquote>
</div>
<span><font color="#888888"><br>
<br clear="all">
<div><br>
</div>
-- <br>
<div>
<div dir="ltr"><b><font color="#990000">Dr J-Donald Tournier (PhD)</font></b><br>
<div><font color="#990000"><br>
</font></div>
<i><font color="#990000">Senior Lecturer, </font></i><i><font color="#990000">Biomedical Engineering</font></i>
<div><i><font color="#990000">Division of Imaging Sciences & Biomedical Engineering<br>
King's College London</font></i>
<div><i><font color="#990000"><br>
</font></i></div>
<div><i><font color="#990000"><b style="font-family:Calibri,sans-serif;font-size:15px"><span style="font-size:10pt">A:</span></b><span style="font-family:Calibri,sans-serif;font-size:10pt"> Department of Perinatal Imaging & Health, 1<sup>st</sup> Floor South
Wing, St Thomas' Hospital, London. SE1 7EH</span><br>
</font></i></div>
<div><i><font color="#990000"><b>T:</b> <a href="tel:%2B44%20%280%2920%207188%207118%20ext%2053613" value="+442071887118" target="_blank">
+44 (0)20 7188 7118 ext 53613</a></font></i></div>
</div>
<div><i><font color="#990000"><b>W:</b> <a href="http://www.kcl.ac.uk/medicine/research/divisions/imaging/departments/biomedengineering" target="_blank">http://www.kcl.ac.uk/medicine/research/divisions/imaging/departments/biomedengineering</a></font></i><br>
</div>
</div>
</div>
</font></span></div>
<br>
_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">Mrtrix-discussion@www.nitrc.org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/mrtrix-discussion" rel="noreferrer" target="_blank">http://www.nitrc.org/mailman/listinfo/mrtrix-discussion</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
<br>
_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">Mrtrix-discussion@www.nitrc.org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/mrtrix-discussion" rel="noreferrer" target="_blank">http://www.nitrc.org/mailman/listinfo/mrtrix-discussion</a><br>
<br>
</blockquote>
</div>
</div>
</div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><b><font color="#990000">Dr J-Donald Tournier (PhD)</font></b><br><div><font color="#990000"><br></font></div><i><font color="#990000">Senior Lecturer, </font></i><i><font color="#990000">Biomedical Engineering</font></i><div><i><font color="#990000">Division of Imaging Sciences & Biomedical Engineering<br>King's College London</font></i><div><i><font color="#990000"><br></font></i></div><div><i><font color="#990000"><b style="font-family:Calibri,sans-serif;font-size:15px"><span style="font-size:10pt">A:</span></b><span style="font-family:Calibri,sans-serif;font-size:10pt"> Department of Perinatal Imaging & Health, 1<sup>st</sup> Floor South Wing, St Thomas' Hospital, London. SE1 7EH</span><br></font></i></div><div><i><font color="#990000"><b>T:</b> +44 (0)20 7188 7118 ext 53613</font></i></div></div><div><i><font color="#990000"><b>W:</b> <a href="http://www.kcl.ac.uk/medicine/research/divisions/imaging/departments/biomedengineering" target="_blank">http://www.kcl.ac.uk/medicine/research/divisions/imaging/departments/biomedengineering</a></font></i><br></div></div></div>
</div>