Menu

Commit [r25532]  Maximize  Restore  History

Modified module which estimate model parameter errors, to scale to units in relax.

The R2eff errors come in rad/s. The Jacobian come in rad/s.
The co-variance matrix come in (rad/s)^2. The standard deviation of the parameters the comes in rad/s.

But dw is stored as ppm, therefore a scaling is necessary, and is performed with the scaling matrix from the target function.

Now there is excellent agreement.

Systemtest:
relax -s Relax_disp.test_task_model_par_est_tsmfk01 -d

param: r2a_err, with err: 0.48131663, compared to MC: 0.49354342, with 2 boot 0.00000000
param: dw_err, with err: 0.30874093, compared to MC: 0.32813703, with 2 boot 0.00000000
param: k_AB_err, with err: 0.64606994, compared to MC: 0.65384783, with 2 boot 0.00000000

Juhuuu!

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

tlinnet 2014-09-02

changed /branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py
/branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py Diff Switch to side-by-side view
Loading...
Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.