Data Quality and Exclusion Criteria in Osprey’s LCModel Output

Hi Admin,

Thank you so much for your detailed response—it’s been incredibly helpful for those of us relatively new to MRS. I also apologize for the delay in following up, as I’ve been testing new parameters for my data and organizing my thoughts. We would be very grateful if you could respond as quickly as you did last time.

Representative Data

All spectra were acquired from the mPFC region, with each subject undergoing two sessions.

Upon review, I noticed that the scan parameters for my MEGA-PRESS data weren’t optimal. Specifically, I found a delta frequency of -2.8 ppm in the water reference scan, a low TR (1500ms), and low averages (128), which differ from the consensus recommendations (e.g., Peek et al., 2023).

Despite this, I’m aiming to make the best of the available data. I tested various parameter combinations and found that changing opts.SpecReg from RobSpecReg to RestrSpecReg was effective. This adjustment, alongside using opts.fit.range [0.5 - 4 ppm] for fitting, generally produced smoother spectra for noisier sessions. It also reduced the overall residual (31.0% vs. 31.7%).

However, there are still cases where RestrSpecReg fails to fully resolve the spectra smoothly.

You can view examples of the improvements here:
RestrSpecReg_Improvement

As requested, I’ve attached some representative data, which includes concatenated output files for both mean and individual session data:
Overview_Data

MM Fitting Strategies

I also have some additional questions regarding signal quality and macromolecule (MM) fitting. I reran the analysis with the 0.55 ppm baseline knot spacing and compared two fitting strategies: 3to2MM and 1to1GABAsoft.

In the 3to2MM model, the macromolecule signal seems to override the expected GABA signal around 3 ppm. Below is a comparison of the fitted difference spectra for the 3to2MM, 1to1GABAsoft, and 1to1GABA models. The red fitted curve for 3to2MM shows only one bump at 3 ppm, whereas the 1to1GABAsoft model shows three bumps, which seems more consistent with the green average. Interestingly, the residuals for 3to2MM (30.0%) are slightly lower than for 1to1GABAsoft (31.0%). The 1to1GABAhard model, while missing the middle bump, has a lower residual than 1to1GABAsoft (31.6%).

You can see the fitting comparison here:
MM_fitting / Fitting Method

I’ve also attached a distribution of individual data points comparing the 3to2MM and 1to1GABAsoft fits (by Max-Min normalization). In about half of the cases under 3to2MM, the GABA concentration is zero, and the GABA CRLB value is set to 999. In contrast, with the 1to1GABAsoft fit, only one data point has a GABA CRLB value of 999, indicating a fit failure. However, two GABA+ CRLB values and one Glx value under 1to1GABAsoft are beyond the 2 SD range.

Distribution for 1to1GABAsoft fitting

Distribution for 3to2MM fitting

Additionally, I’ve uploaded two individual case fitting graphs for your review:
MM_Fitting / Clean_Compare

Each PDF page contains two spectra: 1to1GABAsoft on the left and 3to2MM on the right. These cases have smoother spectra compared to others, and as you can see, the 3to2MM model doesn’t utilize the GABA component for fitting at 3 ppm, being completely dominated by the MM09 peak.

Lastly, referencing a previous forum post, some users have found that 1to1GABAsoft might be more suitable for MEGA-PRESS data. Would this also apply to my case?
Forum Reference

CRLB Screening

We are applying a CRLB ± 2 SD cutoff to filter out participants with extreme data uncertainty. Among the GABA+ CRLB values, two clear outliers stand out. I’ve attached these flagged cases for your review, along with their respective fits under both the 3to2MM and 1to1GABAsoft models. Would you recommend excluding these data points?

Problem_Spectra / CRLB_Flagged

Potentially Noisy Spectra Without CRLB Flags

We’ve also identified several spectra that we are unsure about in terms of quality, even though they weren’t flagged by CRLB values. We’ve categorized these spectra to make it easier for you to assess:

  • Sharp Spikes and Troughs
    Some spectra exhibit rapid, sharp spikes and troughs. Are these spectra acceptable for analysis?
    Problem_Spectra / Noisy GABA
  • Glx Unwanted Oscillations
    For Glx fitting at ~3.7 ppm, some spectra show unwanted oscillations. Are these out-of-voxel artifacts that should be discarded?
    Problem_Spectra / Noisy Glx
  • Baseline Spikes
    In some cases, the baseline between Glx and GABA peaks includes large spikes. Are these spectra suitable for analysis?
    Problem_Spectra / Noisy Baseline

Any insights or recommendations on how to proceed with these cases would be greatly appreciated. Should we exclude certain spectra based on these irregularities, even if they haven’t been flagged by CRLB?

Thanks again for all your help!