Listed below are the data quality metrics that Gannet computes during data processing, signal fitting, and tissue segmentation.

Linewidth

Linewidth is calculated as the full-width half-maximum (FWHM) (in Hz) of fitted model signals. When reporting linewidths of datasets, it is recommended to use the FWHM of Cr or NAA from the OFF spectrum, or the FWHM of the water reference (if a water reference is provided).

Signal-to-noise ratio (SNR)

The SNR of fitted signals is calculated as the amplitude of the given signal model divided by twice the standard deviation of the noise signal. To estimate the noise signal, Gannet takes two independent segments of the OFF or DIFF spectrum (as appropriate to the modeled signal of interest) between 8–9 ppm and 9–10 ppm, and detrends them using a second-order polynomial function. Detrending is performed to remove baseline artifacts (often related to the residual water signal). The standard deviation of each detrended noise segment is then calculated. The smaller of the two standard deviations is then used as the estimate of noise, which is then multiplied by 2.

Formulaically, this is defined as:

\[ SNR_{signal} = \frac{A_{signal}}{2\cdot\mathrm{std}(noise)} \]

where:

Parameter Definition
\(A_{signal}\) Signal model amplitude
\(noise\) Detrended noise signal between either 8–9 or 9–10 ppm (whichever produces a lower standard deviation) in the appropriate spectrum (i.e., either the OFF or DIFF spectrum)

Frequency offsets (frequency drift/motion)

To estimate the degree of frequency offsets that result from scanner-related frequency drift1 and participant motion2, Gannet calculates an average frequency offset3 \(\overline{\Delta\delta_{0}}\). This is calculated as the mean (over the course of the acquisition) difference between the observed frequency of the residual water signal in the pre-frequency-corrected subspectra and the nominal water frequency \(\delta_{0}\) at 4.68 ppm (4.8 ppm for room-temperature phantoms), or the nominal Cr frequency at 3.02 ppm for HERMES acquisitions. It should be noted that using the mean of offset differences does not fully characterize frequency offsets but is a useful heuristic.

Formulaically, this is defined as:

\[ \overline{\Delta\delta_{0}} = \frac{1}{m}\sum{\widehat{\delta_{0,m}} - \delta_{0}} \] where:

Parameter Definition
\(m\) Each individual subspectrum index number
\(\widehat{\delta_{0,m}}\) Observed water or Cr frequency in each individual subspectrum
\(\delta_{0}\) Nominal water or Cr frequency

Fit error

When fitting signal functions to metabolite peaks, Gannet will also estimate the error of the fit. This is defined as the standard deviation of the residuals of the signal model fit divided by the signal model amplitude and multiplied by 100 to give a percentage.

For metabolite peak fits, this is:

\[ \epsilon_{metab} = 100\cdot\frac{\mathrm{std}(resid_{metab})}{A_{metab}} \]

where:

Parameter Definition
\(resid_{metab}\) Signal model fit residuals
\(A_{metab}\) Signal model amplitude

Similarly, for reference signal fits, this is:

\[ \epsilon_{ref} = 100\cdot\frac{\mathrm{std}(resid_{ref})}{A_{ref}} \]

Since all metabolites are normalized to a reference signal (either Cr or unsuppressed water) and reported as such, the fit error that really should be considered (and reported) is the combined error of the metabolite and reference signal model fits, which add up in quadrature.

Formulaically, this is defined as:

\[ \epsilon_{metab,ref} = \sqrt{\epsilon_{metab}^2 + \epsilon_{ref}^2} \]

References

1.
Hui SCN, Mikkelsen M, Zöllner HJ, et al. Frequency drift in MR spectroscopy at 3T. NeuroImage. 2021;241(21):118430. doi:10.1016/j.neuroimage.2021.118430
2.
Evans CJ, Puts NAJ, Robson SE, et al. Subtraction artifacts and frequency (Mis-)alignment in J-difference GABA editing. Journal of Magnetic Resonance Imaging. 2013;38(4):970-975. doi:10.1002/jmri.23923
3.
Mikkelsen M, Barker PB, Bhattacharyya PK, et al. Big GABA: Edited MR spectroscopy at 24 research sites. NeuroImage. 2017;159:32-45. doi:10.1016/j.neuroimage.2017.07.021
LS0tCnRpdGxlOiAiRGF0YSBxdWFsaXR5IG1ldHJpY3MiCmRhdGU6ICJMYXN0IHVwZGF0ZWQ6IGByIGZvcm1hdChTeXMudGltZSgpLCAnJUIgJWQsICVZJylgIgpiaWJsaW9ncmFwaHk6IGJpYmxpb2dyYXBoeS5iaWIKY3NsOiBhbWVyaWNhbi1tZWRpY2FsLWFzc29jaWF0aW9uLmNzbApsaW5rLWNpdGF0aW9uczogeWVzCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiBUUlVFCiAgICB0b2NfZGVwdGg6IDIKICAgIHRvY19mbG9hdDoKICAgICAgY29sbGFwc2VkOiBGQUxTRQotLS0KCmBgYHtyIHNldHVwLCBpbmNsdWRlID0gRkFMU0V9CmtuaXRyOjpvcHRzX2NodW5rJHNldChlY2hvID0gVFJVRSkKYGBgCgpgYGB7ciwgY2hpbGQgPSAianMvYmFjay10by10b3AuanMifQpgYGAKCmBgYHtjc3MsIGVjaG8gPSBGQUxTRX0KLndhcm5pbmcgewogIG1hcmdpbi1ib3R0b206IDIwcHg7Cn0KYGBgCgo8YnI+CgpMaXN0ZWQgYmVsb3cgYXJlIHRoZSBkYXRhIHF1YWxpdHkgbWV0cmljcyB0aGF0IEdhbm5ldCBjb21wdXRlcyBkdXJpbmcgZGF0YSBwcm9jZXNzaW5nLCBzaWduYWwgZml0dGluZywgYW5kIHRpc3N1ZSBzZWdtZW50YXRpb24uCgojIyBMaW5ld2lkdGgKCkxpbmV3aWR0aCBpcyBjYWxjdWxhdGVkIGFzIHRoZSBmdWxsLXdpZHRoIGhhbGYtbWF4aW11bSAoRldITSkgKGluIEh6KSBvZiBmaXR0ZWQgbW9kZWwgc2lnbmFscy4gV2hlbiByZXBvcnRpbmcgbGluZXdpZHRocyBvZiBkYXRhc2V0cywgaXQgaXMgcmVjb21tZW5kZWQgdG8gdXNlIHRoZSBGV0hNIG9mIENyIG9yIE5BQSBmcm9tIHRoZSBPRkYgc3BlY3RydW0sIG9yIHRoZSBGV0hNIG9mIHRoZSB3YXRlciByZWZlcmVuY2UgKGlmIGEgd2F0ZXIgcmVmZXJlbmNlIGlzIHByb3ZpZGVkKS4KCiMjIFNpZ25hbC10by1ub2lzZSByYXRpbyAoU05SKQoKVGhlIFNOUiBvZiBmaXR0ZWQgc2lnbmFscyBpcyBjYWxjdWxhdGVkIGFzIHRoZSBhbXBsaXR1ZGUgb2YgdGhlIGdpdmVuIHNpZ25hbCBtb2RlbCBkaXZpZGVkIGJ5IHR3aWNlIHRoZSBzdGFuZGFyZCBkZXZpYXRpb24gb2YgdGhlIG5vaXNlIHNpZ25hbC4gVG8gZXN0aW1hdGUgdGhlIG5vaXNlIHNpZ25hbCwgR2FubmV0IHRha2VzIHR3byBpbmRlcGVuZGVudCBzZWdtZW50cyBvZiB0aGUgT0ZGIG9yIERJRkYgc3BlY3RydW0gKGFzIGFwcHJvcHJpYXRlIHRvIHRoZSBtb2RlbGVkIHNpZ25hbCBvZiBpbnRlcmVzdCkgYmV0d2VlbiA44oCTOSBwcG0gYW5kIDnigJMxMCBwcG0sIGFuZCBkZXRyZW5kcyB0aGVtIHVzaW5nIGEgc2Vjb25kLW9yZGVyIHBvbHlub21pYWwgZnVuY3Rpb24uIERldHJlbmRpbmcgaXMgcGVyZm9ybWVkIHRvIHJlbW92ZSBiYXNlbGluZSBhcnRpZmFjdHMgKG9mdGVuIHJlbGF0ZWQgdG8gdGhlIHJlc2lkdWFsIHdhdGVyIHNpZ25hbCkuIFRoZSBzdGFuZGFyZCBkZXZpYXRpb24gb2YgZWFjaCBkZXRyZW5kZWQgbm9pc2Ugc2VnbWVudCBpcyB0aGVuIGNhbGN1bGF0ZWQuIFRoZSBzbWFsbGVyIG9mIHRoZSB0d28gc3RhbmRhcmQgZGV2aWF0aW9ucyBpcyB0aGVuIHVzZWQgYXMgdGhlIGVzdGltYXRlIG9mIG5vaXNlLCB3aGljaCBpcyB0aGVuIG11bHRpcGxpZWQgYnkgMi4KCkZvcm11bGFpY2FsbHksIHRoaXMgaXMgZGVmaW5lZCBhczoKCiQkClNOUl97c2lnbmFsfSA9IFxmcmFje0Ffe3NpZ25hbH19ezJcY2RvdFxtYXRocm17c3RkfShub2lzZSl9CiQkCgp3aGVyZToKCnwgPHU+UGFyYW1ldGVyPC91PiB8IDx1PkRlZmluaXRpb248L3U+IHwKfCA6LSB8IDotLS0tLS0tLS0gfAp8ICRBX3tzaWduYWx9JCB8IFNpZ25hbCBtb2RlbCBhbXBsaXR1ZGUgfAp8ICRub2lzZSQgfCBEZXRyZW5kZWQgbm9pc2Ugc2lnbmFsIGJldHdlZW4gZWl0aGVyIDjigJM5IG9yIDnigJMxMCBwcG0gKHdoaWNoZXZlciBwcm9kdWNlcyBhIGxvd2VyIHN0YW5kYXJkIGRldmlhdGlvbikgaW4gdGhlIGFwcHJvcHJpYXRlIHNwZWN0cnVtIChpLmUuLCBlaXRoZXIgdGhlIE9GRiBvciBESUZGIHNwZWN0cnVtKSB8CgojIyBGcmVxdWVuY3kgb2Zmc2V0cyAoZnJlcXVlbmN5IGRyaWZ0L21vdGlvbikKClRvIGVzdGltYXRlIHRoZSBkZWdyZWUgb2YgZnJlcXVlbmN5IG9mZnNldHMgdGhhdCByZXN1bHQgZnJvbSBzY2FubmVyLXJlbGF0ZWQgZnJlcXVlbmN5IGRyaWZ0IFtASHVpMjAyMWFdIGFuZCBwYXJ0aWNpcGFudCBtb3Rpb24gW0BFdmFuczIwMTNdLCBHYW5uZXQgY2FsY3VsYXRlcyBhbiBhdmVyYWdlIGZyZXF1ZW5jeSBvZmZzZXQgW0BNaWtrZWxzZW4yMDE3XSAkXG92ZXJsaW5le1xEZWx0YVxkZWx0YV97MH19JC4gVGhpcyBpcyBjYWxjdWxhdGVkIGFzIHRoZSBtZWFuIChvdmVyIHRoZSBjb3Vyc2Ugb2YgdGhlIGFjcXVpc2l0aW9uKSBkaWZmZXJlbmNlIGJldHdlZW4gdGhlIG9ic2VydmVkIGZyZXF1ZW5jeSBvZiB0aGUgcmVzaWR1YWwgd2F0ZXIgc2lnbmFsIGluIHRoZSBwcmUtZnJlcXVlbmN5LWNvcnJlY3RlZCBzdWJzcGVjdHJhIGFuZCB0aGUgbm9taW5hbCB3YXRlciBmcmVxdWVuY3kgJFxkZWx0YV97MH0kIGF0IDQuNjggcHBtICg0LjggcHBtIGZvciByb29tLXRlbXBlcmF0dXJlIHBoYW50b21zKSwgb3IgdGhlIG5vbWluYWwgQ3IgZnJlcXVlbmN5IGF0IDMuMDIgcHBtIGZvciBIRVJNRVMgYWNxdWlzaXRpb25zLiBJdCBzaG91bGQgYmUgbm90ZWQgdGhhdCB1c2luZyB0aGUgbWVhbiBvZiBvZmZzZXQgZGlmZmVyZW5jZXMgZG9lcyBub3QgZnVsbHkgY2hhcmFjdGVyaXplIGZyZXF1ZW5jeSBvZmZzZXRzIGJ1dCBpcyBhIHVzZWZ1bCBoZXVyaXN0aWMuCgpGb3JtdWxhaWNhbGx5LCB0aGlzIGlzIGRlZmluZWQgYXM6CgokJApcb3ZlcmxpbmV7XERlbHRhXGRlbHRhX3swfX0gPSBcZnJhY3sxfXttfVxzdW17XHdpZGVoYXR7XGRlbHRhX3swLG19fSAtIFxkZWx0YV97MH19CiQkCndoZXJlOgoKfCA8dT5QYXJhbWV0ZXI8L3U+IHwgPHU+RGVmaW5pdGlvbjwvdT4gfAp8IDotIHwgOi0tLS0tLS0tLSB8CnwgJG0kIHwgRWFjaCBpbmRpdmlkdWFsIHN1YnNwZWN0cnVtIGluZGV4IG51bWJlciB8CnwgJFx3aWRlaGF0e1xkZWx0YV97MCxtfX0kIHwgT2JzZXJ2ZWQgd2F0ZXIgb3IgQ3IgZnJlcXVlbmN5IGluIGVhY2ggaW5kaXZpZHVhbCBzdWJzcGVjdHJ1bSB8CnwgJFxkZWx0YV97MH0kIHwgTm9taW5hbCB3YXRlciBvciBDciBmcmVxdWVuY3kgfAoKIyMgRml0IGVycm9yCgpXaGVuIGZpdHRpbmcgc2lnbmFsIGZ1bmN0aW9ucyB0byBtZXRhYm9saXRlIHBlYWtzLCBHYW5uZXQgd2lsbCBhbHNvIGVzdGltYXRlIHRoZSBlcnJvciBvZiB0aGUgZml0LiBUaGlzIGlzIGRlZmluZWQgYXMgdGhlIHN0YW5kYXJkIGRldmlhdGlvbiBvZiB0aGUgcmVzaWR1YWxzIG9mIHRoZSBzaWduYWwgbW9kZWwgZml0IGRpdmlkZWQgYnkgdGhlIHNpZ25hbCBtb2RlbCBhbXBsaXR1ZGUgYW5kIG11bHRpcGxpZWQgYnkgMTAwIHRvIGdpdmUgYSBwZXJjZW50YWdlLgoKRm9yIG1ldGFib2xpdGUgcGVhayBmaXRzLCB0aGlzIGlzOgoKJCQKXGVwc2lsb25fe21ldGFifSA9IDEwMFxjZG90XGZyYWN7XG1hdGhybXtzdGR9KHJlc2lkX3ttZXRhYn0pfXtBX3ttZXRhYn19CiQkCgp3aGVyZToKCnwgPHU+UGFyYW1ldGVyPC91PiB8IDx1PkRlZmluaXRpb248L3U+IHwKfCA6LSB8IDotLS0tLS0tLS0gfAp8ICRyZXNpZF97bWV0YWJ9JCB8IFNpZ25hbCBtb2RlbCBmaXQgcmVzaWR1YWxzIHwKfCAkQV97bWV0YWJ9JCB8IFNpZ25hbCBtb2RlbCBhbXBsaXR1ZGUgfAoKU2ltaWxhcmx5LCBmb3IgcmVmZXJlbmNlIHNpZ25hbCBmaXRzLCB0aGlzIGlzOgoKJCQKXGVwc2lsb25fe3JlZn0gPSAxMDBcY2RvdFxmcmFje1xtYXRocm17c3RkfShyZXNpZF97cmVmfSl9e0Ffe3JlZn19CiQkCgpTaW5jZSBhbGwgbWV0YWJvbGl0ZXMgYXJlIG5vcm1hbGl6ZWQgdG8gYSByZWZlcmVuY2Ugc2lnbmFsIChlaXRoZXIgQ3Igb3IgdW5zdXBwcmVzc2VkIHdhdGVyKSBhbmQgcmVwb3J0ZWQgYXMgc3VjaCwgdGhlIGZpdCBlcnJvciB0aGF0IHJlYWxseSBzaG91bGQgYmUgY29uc2lkZXJlZCAoYW5kIHJlcG9ydGVkKSBpcyB0aGUgY29tYmluZWQgZXJyb3Igb2YgdGhlIG1ldGFib2xpdGUgYW5kIHJlZmVyZW5jZSBzaWduYWwgbW9kZWwgZml0cywgd2hpY2ggYWRkIHVwIGluIHF1YWRyYXR1cmUuCgpGb3JtdWxhaWNhbGx5LCB0aGlzIGlzIGRlZmluZWQgYXM6CgokJApcZXBzaWxvbl97bWV0YWIscmVmfSA9IFxzcXJ0e1xlcHNpbG9uX3ttZXRhYn1eMiArIFxlcHNpbG9uX3tyZWZ9XjJ9CiQkCgojIyMgUmVmZXJlbmNlcwoKCgo=


Built with R Markdown in RStudio

Copyright © 2020–2024, Mark Mikkelsen