Trusted, Bureaucrats
4,223
edits
Changes
m
Switch to the {{gna task link}} template to remove dead Gna! links.
{{lowercase title}}
{{infobox relax release| version = 3.3.9| prev = 3.3.8| next = 4.0.0| keywords = N-state model, relaxation dispersion| type = Minor feature| date = 30 September 2015| manual = yes}} == Description ==
<section begin=description/>
This is a minor feature release with improvements to the automatic relaxation dispersion protocol for repeated CPMG data, support for Monte Carlo or Bootstrap simulating RDC and PCS errorsQ factors, a huge speedup of Monte Carlo simulations in the N-state model analysis, and geometric mean and standard deviation functions added to the relax library.
<section end=description/>
== Download ==
<section begin=download/>
<section end=download/>
== CHANGES file ==
<section begin=metadata/>
<section end=metadata/>
=== Features ===
<section begin=features/>
* Improvements to the automatic relaxation dispersion protocol for repeated CPMG data.
* Support for Monte Carlo or Bootstrap simulating the RDC and PCS Q factors.
* Huge speedup of Monte Carlo simulations in the N-state model analysis.
* Geometric mean and standard deviation functions added to the relax library.
<section end=features/>
=== Changes ===
<section begin=changes/>
* Wrote a method to store parameter data and dispersion curves, for the protocol of repeated CPMG analysis. This is to prepare for analysis in other programs. The method loops through the data pipes, and writes the data out. It then writes a bash script that will concatenate the data in an matrix array style, for reading and processing in other programs. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
* Added to write out a collection script for χ<sup>2</sup> and rate parameters. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
* In the collection bash script, removes spins which have not been fitted. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
* Fix for use of " instead of ' in bash script. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
* 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. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
* Printout when minimising Monte Carlo simulations. {{gna task link|7826|text=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].
* Small fix for the shell script to collect data files, and not use the program "column" in the end. The line width becomes to large to handle for column. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
* 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.
* Silenced a warning in the N-state model optimisation if the verbosity is set to zero. This removes a repetitive warning from the Monte Carlo or Bootstrap simulations.
* 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].
* Made a summarize function to compare results. {{gna task link|7826|text=Task #7826: Write an Python class for the repeated analysis of dispersion data}}.
<section end=changes/>
=== Bugfixes ===
<section begin=bugfixes/>
* Fix committed, where an empty spin placeholder has the select flag set to False. [https://gna.org/bugs/?23642 Bug #23642: When deleting all spins for a residue, an empty placeholder is where select=True].
<section end=bugfixes/>
== Links ==
<section begin=links/>
For reference, the announcement for this release can also be found at following links:
* [http://wiki.nmr-relax.com/Relax_3.3.8 9 Official release notes on the relax wiki].* [https://{{gna link|url=gna.org/forum/forum.php?forum_id=2496 2513|text=Gna! news item]}}.
* [http://article.gmane.org/gmane.science.nmr.relax.announce/67 Gmane mailing list archive].
* [http://www.mail-archive.com/relax-announce%40gna.org/msg00057msg00058.html The Mail Archive].* [https://mail.gna.org/public/relax-announce/2015-0410/msg00000.html Local archives].* [https://marc.info/?l=relax-announce&m=142628131518605144376869310288&w=2 Mailing list ARChives (MARC)].
{{:Softpedia}}
<section end=links/>
== Announcements ==
{{:relax release announcements}}
== See also ==
* [http://www.nmr-relax.com/api/3.3/ The relax 3.3 API documentation]
{{:relax release see also}}
[[Category:Relaxation dispersion analysis]]