Author: tlinnet Date: Tue Sep 2 14:03:37 2014 New Revision: 25543 URL: http://svn.gna.org/viewcvs/relax?rev=25543&view=rev Log: Fix for dimension extractions, when estimating the parameter errors.
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=25543&r1=25542&r2=25543&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 14:03:37 2014 @@ -250,7 +250,7 @@ # NM: Number of spectrometer frequencies. # NO: Maximum number of offsets. # ND: Number of dispersion(data) points. - NJ, NE, NS, NM, NO, ND = func_jacobian.shape + ND, NE, NS, NM, NO, NJ = func_jacobian.shape if NS != 1: raise RelaxError("The number of spins does not fit.") @@ -268,7 +268,7 @@ cur_weights = weights[ei, si, mi, oi] # Extract every column/row from the first to last columns. Is this correct? - cur_jacobian = func_jacobian[0:NJ:1, ei, si, mi, oi] + cur_jacobian = func_jacobian[0:ND:1, ei, si, mi, oi] # Get the co-variance pcov = multifit_covar(J=cur_jacobian, weights=cur_weights) _______________________________________________ relax (http://www.nmr-relax.com) This is the relax-commits 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-commits
