Hi,
I'll answer below:
I performed 'cpmg fixed' a relaxation dispersion analysis, for dataset with
68 residues.
Having 22 intensity files, with 4 triple replications.
I took from 17 pm to 13 pm following day, app 20 hours.
The analysed models were: R2eff', 'No Rex','LM63','CR72'
Is 20 hours for an analysis expected?
, or should I look for some errors somewhere?
This depends. By default relax uses much higher precision
optimisation than most other softwares. This is based on the
philosophy that more accurate results are worth the wait, especially
considering that this time is relatively small in comparison to the
measurement and processing time (and adding the inevitable
re-measurements). In addition, lots of Monte Carlo simulations are
used for really accurately determining your parameter errors. For
comparison, note that a full a Lipari and Szabo model-free analysis
can take between 1 day and 2 weeks to complete.
For initial analyses where errors are not so important, the number of
simulations can be dropped massively to speed things up. This can
also be done for students. The result is that the error estimates of
the parameters are horrible but, in some cases, excluding publication,
that is not such a problem. This will not affect model selection
either. Therefore if errors are not important for specific cases,
then set the number of MC sims to 3. Then watch how much quicker
things run.
Oh, if this is for students using the GUI, you could hack a special
version of relax to have the auto-analysis perform much lower
precision optimisation. It should be possible to make things run very
quick for them. If they are using scripting, but still the relaxation
dispersion auto-analysis, then no hacking is necessary. The function
tolerance and maximum number of iterations can be set using special
tricks ;)
Relax made a model selection.
model_selection(method='AIC', modsel_pipe='final', bundle=None, pipes=['No
Rex', 'LM63', 'CR72'])
How can I inspect which model is then chosen?
Well, the hard way would be to open the relax saved state in the
'final' directory. Or to look at the logs. That information is
unfortunately not presented in a text file. How would you suggest
that such information is presented?
Then I would like to make a global fit / cluster analysis.
Is this implemented yet?
Yes. You need to use the relax_disp.cluster user function.
And should I use the:
User functions (n-z) -> relax_disp -> cluster ?
I am still thinking about how to bring this into the relaxation
dispersion auto-analysis GUI element without requiring me to write a
lot of code!
By the way.
The sherekhan input function is super great!
:)
The CPMGFit and NESSY input user functions will hopefully also be of
use to some people. However if you start adding some new models to
relax, then relax will soon have all of the capabilities that these
these programs possess. And if Paul Schanda's numerical integration
Python code is merged into relax, then relax will soon surpass all of
these softwares.
Thanks a lot :-)
This is super great
You're welcome. I'm hoping that you can help to make this even greater ;)
Regards,
Edward
Re: Time of running, Model selection and global/cluster analysis for relaxation dispersion analysis