TY - JOUR
T1 - Maximum likelihood and the EM algorithm for 2D NMR spectroscopy
AU - Miller, Michael I.
AU - Chen, Shou C.
AU - Kuefler, David A.
AU - D’avignon, D. André
PY - 1993/10/1
Y1 - 1993/10/1
N2 - By extending the maximum-likelihood method of Miller and Greene for one-dimensional nuclear magnetic resonance spectroscopy, maximum-likelihood (ML) methods for two-dimensional spectroscopy are developed. The time-domain free-induction decay is modeled as a 2D signal f(t1, t2), the superposition of exponentially decaying complex exponentials in both dimensions of unknown amplitudes, frequencies, and decays. The measurement is modeled as a 2D Gaussian process with mean f(t1, t2) and fixed variance. The ML estimation problem is to choose the set of amplitudes, frequencies, and decays which minimizes the squared error between the measured FID and the parameterized signal model. An iterative expectation-maximization algorithm is proposed for solving the nonlinear least-squares problem, providing a method for estimating the signal parameters directly from the 2D FID. The ML method and implementation algorithm are tested by analyzing both simulated and experimental N-butanol data collected via the hypercomplex data-collection protocol.
AB - By extending the maximum-likelihood method of Miller and Greene for one-dimensional nuclear magnetic resonance spectroscopy, maximum-likelihood (ML) methods for two-dimensional spectroscopy are developed. The time-domain free-induction decay is modeled as a 2D signal f(t1, t2), the superposition of exponentially decaying complex exponentials in both dimensions of unknown amplitudes, frequencies, and decays. The measurement is modeled as a 2D Gaussian process with mean f(t1, t2) and fixed variance. The ML estimation problem is to choose the set of amplitudes, frequencies, and decays which minimizes the squared error between the measured FID and the parameterized signal model. An iterative expectation-maximization algorithm is proposed for solving the nonlinear least-squares problem, providing a method for estimating the signal parameters directly from the 2D FID. The ML method and implementation algorithm are tested by analyzing both simulated and experimental N-butanol data collected via the hypercomplex data-collection protocol.
UR - http://www.scopus.com/inward/record.url?scp=38248999845&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=38248999845&partnerID=8YFLogxK
U2 - 10.1006/jmra.1993.1221
DO - 10.1006/jmra.1993.1221
M3 - Article
AN - SCOPUS:38248999845
SN - 1064-1858
VL - 104
SP - 247
EP - 257
JO - Journal of Magnetic Resonance - Series A
JF - Journal of Magnetic Resonance - Series A
IS - 3
ER -