mailRe: [bug #22554] The distribution of intensity with errors in Monte-Carlo simulations are markedly more narrow than expected.


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by Edward d'Auvergne on August 31, 2014 - 23:32:
Hi Troels,

Before this drives you crazy, the reason is dead simple.  The Monte Carlo
simulations use a Gaussian distribution.  This is correct.  But you are
comparing it to a uniform distribution, which is not correct and is much
broader.  There's no need to worry.

Regards,

Edward




On Sunday, 31 August 2014, Troels E. Linnet <
NO-REPLY.INVALID-ADDRESS@xxxxxxx> wrote:

URL:
  <http://gna.org/bugs/?22554>

                 Summary: The distribution of intensity with errors in
Monte-Carlo simulations are markedly more narrow than expected.
                 Project: relax
            Submitted by: tlinnet
            Submitted on: Sun 31 Aug 2014 07:03:07 PM UTC
                Category: relax's source code
Specific analysis category: Relaxation dispersion
                Priority: 9 - Immediate
                Severity: 4 - Important
                  Status: None
             Assigned to: None
         Originator Name:
        Originator Email:
             Open/Closed: Open
                 Release: Repository: trunk
         Discussion Lock: Any
        Operating System: All systems

    _______________________________________________________

Details:

It seems that relax draw a lot more narrow distribution of Intensity with
errors, than python module "random".
This has an influence on estimated parameter error.

This is a potential huge error in relax.
A possible example of a catastrophic implementation of Monte-Carlo
simulations.

See:
http://thread.gmane.org/gmane.science.nmr.relax.scm/23239
http://thread.gmane.org/gmane.science.nmr.relax.scm/23238

Attached image that show distribution narrowing for Monte-Carlo
simulations in
relax.

Also attached 2 images, that show there is perfect correlation for
estimation
of errors between boot strapping and from Co-variance matrix.




    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: Sun 31 Aug 2014 07:03:07 PM UTC  Name:
estimate_errors_analyse_relax.py.png  Size: 47kB   By: tlinnet

<http://gna.org/bugs/download.php?file_id=21975>
-------------------------------------------------------
Date: Sun 31 Aug 2014 07:03:07 PM UTC  Name: estimate_errors_data_nt5.png
Size: 73kB   By: tlinnet

<http://gna.org/bugs/download.php?file_id=21976>
-------------------------------------------------------
Date: Sun 31 Aug 2014 07:03:07 PM UTC  Name: estimate_errors_data_nt10.png
Size: 74kB   By: tlinnet

<http://gna.org/bugs/download.php?file_id=21977>

    _______________________________________________________

Reply to this item at:

  <http://gna.org/bugs/?22554>

_______________________________________________
  Message sent via/by Gna!
  http://gna.org/


_______________________________________________
relax (http://www.nmr-relax.com)

This is the relax-devel mailing list
relax-devel@xxxxxxx <javascript:;>

To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-devel



Related Messages


Powered by MHonArc, Updated Mon Sep 01 00:00:38 2014