Download Bayesian Analysis Software User Manual

Transcript
THE BAYES ANALYZE ERROR MESSAGES
195
36. TRANS – Program Error Contact Varian
37. ALOGP – Possible Program Error Investigate Further
38. ALOGP2 – Possible Program Error Investigate Further
39. SETMODEL – Program Error Contact Varian
40. SCCETRANS – Program Error Contact Varian
41. SCUETRANS – Program Error Contact Varian
42. GETMODEL – Program Error Contact Varian
43. APPEND TEXT – Program Error Contact Varian
44. READHDR – Program Error Contact Varian
45. WRITEDATA – Program Error Contact Varian
46. WRITEMODEL – Program Error Contact Varian
47. WRITERESID – Program Error Contact Varian
48. COPY PROCPAR – Program Error Contact Varian
49. LIST MODEL – Program Error Contact Varian
The two exceptions that we were discussing are messages numbers 37 and 38. These two messages
are program errors, but they may be caused by the model. For example if the estimated frequencies
become redundant, then the covariance matrix will be singular and these errors are issued. Message
number 37 indicates that this matrix is singular. When this matrix becomes singular it is almost
always because one or more of the resonances have converged to the same resonance frequency. This
behavior may occur under several circumstances: first it is possible that a peak was marked and
this peak was so far from the value indicated by probability theory that Bayes Analyze added one
or more additional resonances near the same location. These resonances are nearly identical and
this can cause the error. The easiest thing to try under these conditions is to remove the offending
resonances and allow Bayes Analyze to run in its automatic model. There is one other condition
where this problem has been known to occur. NMR lines never decay in a perfectly exponential
manner because the magnetic fields are never perfectly uniform. When the signal-to-noise is very
very high Bayes Analyze may try to fit so many resonances at nearly the same resonance frequency
that the matrix it is trying to invert becomes singular. Unfortunately, in this case, there is very
little that can be done. You can try turning on the shimming models and rerunning the analysis,
but if the problem re-occurs then you may not be able to analyze this data.
So far we have been discussing error number 37, but error number 38 is almost the same and it
is caused by essentially the same conditions. The previous error occurs when the matrix inversion
routine detects that the matrix is singular. Here the inversion routine did not detect failure and it
returned a zero status. After the inverse is computed, Bayes Analyze computes the projection of
the data onto the model. What it found was that the projection of the data onto the model was
greater than the projection of the data onto itself. This is simply not possible, and when it happens
it means that the matrix inverse failed.