Menu

Commit [r25238]  Maximize  Restore  History

Modified systemtest Relax_disp.test_estimate_r2eff.

This is to compare against errors simulated with 2000 MC.

The paramaters are comparable, but not equal.
Mostly, it seems that the errors from scipy.optimize.leastsq, are twice as high than the Monte Carlo simulations.
This affect model fitting, and the calculated chi2 value.

Left column is 2000 Monte Carlo, right column is scipy.optimize.leastsq.

Optimised parameters for spin: 52V @N
Model: No Rex
Parameter: r1 Value: 1.46138806 - 1.46328102
Parameter: r2 Value: 11.48392438 - 11.48040934
Parameter: chi2 Value: 848.42015672 - 3363.95829122

Model: DPL94
Parameter: r1 Value: 1.44845743 - 1.45019848
Parameter: r2 Value: 10.15688373 - 10.16304892
Parameter: phi_ex Value: 0.07599563 - 0.07561937
Parameter: kex Value: 4460.43707304 - 4419.03906628
Parameter: chi2 Value: 179.47041255 - 710.24767560

Model: TP02
Parameter: r1 Value: 1.54354392 - 1.54352369
Parameter: r2 Value: 9.72654895 - 9.72772727
Parameter: pA Value: 0.88827039 - 0.88807488
Parameter: dw Value: 1.08875836 - 1.08765645
Parameter: kex Value: 4921.28597928 - 4904.70134941
Parameter: chi2 Value: 29.33882481 - 114.47142772

Model: TAP03
Parameter: r1 Value: 1.54356410 - 1.54354368
Parameter: r2 Value: 9.72641885 - 9.72759371
Parameter: pA Value: 0.88828925 - 0.88809317
Parameter: dw Value: 1.08837248 - 1.08726695
Parameter: kex Value: 4926.42974479 - 4909.86896567
Parameter: chi2 Value: 29.29050624 - 114.27987534

Model: MP05
Parameter: r1 Value: 1.54356415 - 1.54354372
Parameter: r2 Value: 9.72641730 - 9.72759220
Parameter: pA Value: 0.88828927 - 0.88809322
Parameter: dw Value: 1.08837250 - 1.08726707
Parameter: kex Value: 4926.44228958 - 4909.88128236
Parameter: chi2 Value: 29.29054252 - 114.28002272

Model: NS R1rho 2-site
Parameter: r1 Value: 1.41359226 - 1.41321968
Parameter: r2 Value: 9.34531364 - 9.34602793
Parameter: pA Value: 0.94504369 - 0.94496541
Parameter: dw Value: 1.56001843 - 1.55833321
Parameter: kex Value: 5628.66529504 - 5610.20221435
Parameter: chi2 Value: 34.44010458 - 134.14368365

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.

tlinnet 2014-08-24

changed /trunk/test_suite/system_tests/relax_disp.py
/trunk/test_suite/system_tests/relax_disp.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.