Hi,
I’m wondering if there is any way to get the Cramer-Rao lower bounds for a fit in Osprey. Alternatively, what type of error estimations in the fit parameters are available in Osprey?.
Thanks,
/Greger
Hi Greger,
We’ve been working hard on getting CRLBs returned correctly, but there are still glitches when we calculate the Jacobian matrix… it’s something that keeps haunting us. We’re aware that’s quite the shortcoming for now, but hope to have it fixed at some point. You’ll read it here 
Best,
Georg
Not at this point. We’re planning some larger changes to the fit module this year, and we’re certainly eager to include CRLBs in those changes, too.
Very excited to see what you all have planned!
1 Like
Hello,
I was looking for info on CRLB and came across this thread. Is there any news regarding this topic?
BW,
Andreia
Not yet, and I suspect that @admin will support this with articles (below) explaining that CRLB is not an adequate approximation to standard deviation in MRS. 
Best,
Frederic
Roland Kreis.The trouble with quality filtering based on relative Cramér-Rao lower bounds.2016
Karl Landheer, Christoph Juchem.Are Cramér-Rao lower bounds an accurate estimate for standard deviations in in vivo magnetic resonance spectroscopy?2021
Hi Frederic,
Thanks for your reply, it makes sense indeed!
@admins Are there any plans to include some type of error estimation for metabolite quantification? Or how do you recommend that one report on Osprey’s fitting quality?
Thanks,
Andreia
@admin Can I please check how you would recommend reporting quality in Osprey? I understand that CRLB% is not appropriate, but should absolute CRLB still be reported (and if so will Osprey provide it)? Is reporting Cr SNR, Cr FWHM, water FWHM, and freqShift adequate? Many thanks.
1 Like
Hi, I also wanted to follow up on this topic and I am interested in reporting data quality. Is there also an option for linebroading?
Thank you. Best,
Stephanie
-
The current version of Osprey does not provide CRLB, unfortunately. We are testing a new release currently which will come with a completely new modeling module including 2D fitting and, finally, properly validated CRLB. @jbarnacle Yes, I would report SNR and FWHM; Osprey also gives you a relative fit quality metric (calculated by comparing the maximum residual to the standard deviation of the noise).
-
There is no option for linebroadening. That is on purpose. Linebroadening (and, btw, zero-filling beyond a factor of 2) destroys the independence of data points and therefore violates assumptions made during linear-combination modeling (more precisely, the least-squares solution is identical to the maximum-likelihood estimate only if the noise is independent and identically distributed - to be fair, I’m not aware how bad it is to violate these assumptions and how far off the solutions would be - someone should write a paper about this).
2 Likes