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 <
[email protected]> 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
> [email protected] <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
>
_______________________________________________
relax (http://www.nmr-relax.com)

This is the relax-devel mailing list
[email protected]

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

Reply via email to