Changes

* Adding option to [http://www.nmr-relax.com/api/3.3/auto_analyses.relax_disp_repeat_cpmg.Relax_disp_rep-class.html#minimise_execute minimise class function], to perform Monte Carlo error analysis. [https://gna.org/task/?7826 Task #7826: Write an Python class for the repeated analysis of dispersion data].

* Added additional test to [http://www.nmr-relax.com/api/3.3/test_suite.system_tests.relax_disp.Relax_disp-class.html#test_bug_23186_cluster_error_calc_dw system test Relax_disp.test_bug_23186_cluster_error_calc_dw() ] to prove that Bug #23619 is invalid. [https://gna.org/bugs/?23619 Bug #23619: Stored chi2 sim values from Monte Carlo simulations does not equal normal chi2 values].

* Added a unit test that triggers the bug. Test added in [http://www.nmr-relax.com/api/3.3/test_suite.unit_tests.spin_testing_base.Spin_base_class-class.html#test_delete_spin_all test_delete_spin_all], and can be accessed with: relax -u _pipe_control.test_spin. [https://gna.org/bugs/?23642 Bug #23642]: When deleting all spins for a residue, an empty placeholder is where select=True.* Added sample data and analysis script, that will eventually show that there is not much difference in the sample statistics used for comparing the output of two very similar datasets. This is a multiple comparison test with many T-tests at once, where the familywise error is controlled by the Holm method. Even if the values are close to equal, and within the standard deviation, this procedure will reject up to 20% of the null hypothesis. This is not deemed as a suitable method. [https://gna.org/bugs/?23644 Bug #23644: monte_carlo.error_analysis() does not update the mean value/expectation value from simulations].* Added Monte Carlo simulations to the [http://www.nmr-relax.com/api/3.3/test_suite.system_tests.n_state_model.N_state_model-class.html#test_absolute_T N_state_model.test_absolute_T system test]. This is to demonstrate a failure of the simulations in certain N-state model setups.* Added a missing call to [http://www.nmr-relax.com/manual/monte_carlo_initial_values.html monte_carlo.initial_values ] in the [http://www.nmr-relax.com/api/3.3/test_suite.system_tests.n_state_model.N_state_model-class.html#test_absolute_T N_state_model.test_absolute_T system test]. This fixes the N_state_model.test_absolute_T system test, showing that there is not a problem with the Monte Carlo simulations.* Added Monte Carlo and Bootstrap simulation support for the RDC and PCS Q factor calculations. The [http://www.nmr-relax.com/api/3.3/pipe_control.rdc-module.html#q_factors pipe_control.rdc.q_factors() ] and [http://www.nmr-relax.com/api/3.3/pipe_control.pcs-module.html#q_factors pipe_control.pcs.q_factors() ] functions have been modified to support Monte Carlo and Bootstrap simulations. The sim_index argument has been added to allow the Q factor for the given simulation number to be calculated. All of the Q factor data structures in the base data pipe now have *_sim equivalents for permanently storing the simulation values. For the simulation values, all the warnings have been silenced.* Added simulation support for the RDC and PCS Q factors in the N-state model analysis. This is for both Monte Carlo and Bootstrap simulation. The simulation RDC and PCS values, as well as the simulation back calculated values are now stored via the [http://www.nmr-relax.com/api/3.3/specific_analyses.n_state_model.optimisation-module.html#minimise_bc_data minimise_bc_data() function ] of [http://www.nmr-relax.com/api/3.3/specific_analyses.n_state_model.optimisation -module.html specific_analyses.n_state_model.optimisation] in the respective spin or interatomic data containers. The analysis specific API methods now send the sim_index value into minimise_bc_data(), as well as the [http://www.nmr-relax.com/api/3.3/pipe_control.rdc-module.html#q_factors pipe_control.rdc.q_factors() ] and [http://www.nmr-relax.com/api/3.3/pipe_control.pcs-module.html#q_factors pipe_control.pcs.q_factors() ] functions.

* Huge speed up for the Monte Carlo simulations in the N-state model analyses. This speed up is also for Bootstrap simulations and the frame order analysis. The change affects the [http://www.nmr-relax.com/manual/monte_carlo_initial_values.html monte_carlo.initial_values] user function. The alignment tensor [http://www.nmr-relax.com/api/3.3/data_store.align_tensor.AlignTensorData-class.html#_update_object _update_object() method ] was very inefficient when updating the Monte Carlo simulation data structures. For each simulation, each of the alignment tensor data structures were being updated for all simulations. Now only the current simulations is being updated. This speeds up the user function by many orders of magnitude.* Added functions for calculating the geometric mean and standard deviation to the relax library. These are the [http://www.nmr-relax.com/api/3.3/lib.statistics-module.html#geometric_mean geometric_mean() ] and [http://www.nmr-relax.com/api/3.3/lib.statistics-module.html#geometric_std geometric_std() ] functions of the [http://www.nmr-relax.com/api/3.3/lib.statistics-module.html lib.statistics module]. The implementation is designed to be fast, using numpy array arithmetic rather than Python loops.* Created a simple unit test for the new [http://www.nmr-relax.com/api/3.3/lib.statistics-module.html#geometric_mean lib.statistics.geometric_mean() function].* Added a unit test for the new [http://www.nmr-relax.com/api/3.3/lib.statistics-module.html#geometric_std lib.statistics.geometric_std() function].