Hi Hongyan,
It looks to me like you don't have enough data to support optimisation
of m6-m8. These models have 4 or 5 parameters, so you need >3 data
points per residue in order to use them (eg. R1, R2 and NOE at more than
one field). Because of this, relax won't optimise these models if you
only provide 3 data sets and so you see no fit parameters.
To limit the models considered in each relax script, you simply need to
edit the line:
runs = ['m0', 'm1', 'm2', 'm3', 'm4', 'm5', 'm6', 'm7', 'm8', 'm9']
to include only the models you want to analyse.
I'm fairly sure the error you get on model selection is a bug trigered
by the fact that some models have not been optimised - if you run the
modsel.py script with:
runs = ['m0', 'm1', 'm2', 'm3', 'm4', 'm5', 'm9']
does that solve the problem?
Hope that helps,
Chris
On Tue, 2006-12-19 at 12:44 +0800, Hongyan Li wrote:
Dear relax-users,
I am a new user of relax, I have tried to run modelfree analysis and get
m0, m1,
to m9 models when I cheked the results it seems m6,m7,m8 with no parameters
fitted and m9 with Rex fitted. I don't think there is m9 in modelfree. Can
we
just fit m1,m2,m3,m4,and m5 as modelfree does? and when I run modsel.py
using
sample scripts I got error and asked me to re-run with debug flag. By typing
relax --debug modsel.py and I got error again and it looks like that:
relax> state.save(file='relax_state_20061912_124038', dir=None, force=0,
compress_type=1)
Opening the file 'relax_state_20061912_124038.bz2' for writing.
Traceback (most recent call last):
File "modsel.py", line 20, in ?
eliminate()
File "/data2/usr-local/relax-1.2.9/prompt/eliminate.py", line 97, in
eliminate
self.relax.generic.eliminate.eliminate(run=run, function=function,
args=args)
File "/data2/usr-local/relax-1.2.9/generic_fns/eliminate.py", line 56, in
eliminate
for i in xrange(num_instances(self.run)):
File "/data2/usr-local/relax-1.2.9/specific_fns/model_free.py", line
2605, in
num_instan
ces
raise RelaxFault
RelaxFault: RelaxError: Impossible to be here, please re-run relax with the
'--debug' flag
and summit a bug report at https://gna.org/projects/relax/.
Does anybody know how to do modelselection?
Very much appreciate any suggestion!
Cheers!
Hongyan
_______________________________________________
relax (http://nmr-relax.com)
This is the relax-devel mailing list
relax-devel@xxxxxxx
To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-devel
Re: error in modsel.py