r25524 - /branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/cpmg_synthetic.py -- September 02, 2014 - 03:05

Author: tlinnet
Date: Tue Sep  2 03:05:57 2014
New Revision: 25524

URL: http://svn.gna.org/viewcvs/relax?rev=25524&view=rev
Log:
Modified data generator script, to have 5% error, and add the error to the 
data.

task #7824(https://gna.org/task/index.php?7824): Model parameter ERROR 
estimation from Jacobian and Co-variance matrix of dispersion models.

Modified:
    
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/cpmg_synthetic.py

Modified: 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/cpmg_synthetic.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/cpmg_synthetic.py?rev=25524&r1=25523&r2=25524&view=diff
==============================================================================
--- 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/cpmg_synthetic.py
 (original)
+++ 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/cpmg_synthetic.py
 Tue Sep  2 03:05:57 2014
@@ -23,6 +23,7 @@
 
 # Python module imports.
 from math import sqrt
+from random import gauss
 
 # relax module imports.
 from lib.io import open_write_file
@@ -168,7 +169,7 @@
 # The set r2eff err, for defining the error on the graphs and in the fitting 
weight.
 # We set the error to be the same for all ncyc eksperiments.
 if not hasattr(ds, 'r2eff_err'):
-    ds.r2eff_err = 0.1
+    ds.r2eff_err = 1.
 
 # The print result info.
 if not hasattr(ds, 'print_res'):
@@ -438,7 +439,9 @@
             set_r2eff_err = r2eff_errs[di]
 
             # Add the defined error to the calculated error.
-            r2eff_w_err = r2eff + set_r2eff_err
+            #r2eff_w_err = r2eff + set_r2eff_err
+            r2eff_w_err = gauss(r2eff, ds.r2eff_err)
+            print r2eff_w_err, r2eff
 
             string = "%.15f %.15f %.3f %.15f\n"%(point, r2eff_w_err, 
ds.r2eff_err, r2eff)
             file.write(string)
@@ -602,7 +605,7 @@
 
 min_algor = 'simplex'
 constraints = False
-number = 10000
+number = 2000
 
 monte_carlo.setup(number=number)
 monte_carlo.create_data()