mailr25527 - /branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py


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

Header


Content

Posted by tlinnet on September 02, 2014 - 03:27:
Author: tlinnet
Date: Tue Sep  2 03:27:50 2014
New Revision: 25527

URL: http://svn.gna.org/viewcvs/relax?rev=25527&view=rev
Log:
Possible fix for the weight, but as the error is 1, then it does not matter.

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/specific_analyses/relax_disp/estimate_r2eff.py

Modified: 
branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py?rev=25527&r1=25526&r2=25527&view=diff
==============================================================================
--- branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py     
  (original)
+++ branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py     
  Tue Sep  2 03:27:50 2014
@@ -258,7 +258,7 @@
             param_key = generate_r20_key(exp_type=exp_type, frq=frq)
 
             # Extract weights.
-            cur_weights = weights[ei, si, mi, oi]
+            cur_weights = 1/weights[ei, si, mi, oi]**2
 
             # Extract every column/row from the first to last columns. Is 
this correct?
             cur_jacobian = jacobian[0:NJ:1, ei, si, mi, oi]




Related Messages


Powered by MHonArc, Updated Tue Sep 02 09:40:07 2014