mailHybrid selection problem ...


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by Alexandar Hansen on November 09, 2006 - 19:01:
Hi,

I'm trying to run a hybrid model selection on my data and have been getting a Traceback error.  This is my script and the output.

script = 'my_hybrid.py'
----------------------------------------------------------------------------------------------------

run.create('Dom1','mf')
run.create('Dom2','mf')
run.create('hybrid_select','mf')
results.read('Dom1','results.bz2','15N_final_domain1/final')
results.read('Dom2','results.bz2','15N_final_domain2/final')

run.hybridise('hybrid',['Dom1','Dom2'])
model_selection('AIC','hybrid_select',['Dom1','Dom2','hybrid'])
----------------------------------------------------------------------------------------------------

relax> run.create(run='Dom1', run_type='mf')
relax> run.create(run='Dom2', run_type='mf')
relax> run.create(run='hybrid_select', run_type='mf')
relax> results.read(run='Dom1', file='results.bz2', dir='15N_final_domain1/final', format='columnar')
Opening the file '15N_final_domain1/final/results.bz2' for reading.
Loading all structures from the PDB file.
Structure('etar_au22_phage_df_pt.pdb', model=1):
  Nucleotide chain of length 68

relax> results.read(run='Dom2', file='results.bz2', dir='15N_final_domain2/final', format='columnar')
Opening the file '15N_final_domain2/final/results.bz2' for reading.
Using the structures from the run 'Dom1'.
Structure('etar_au22_phage_df_pt.pdb', model=1):
  Nucleotide chain of length 68

relax> run.hybridise(hybrid='hybrid', runs=['Dom1', 'Dom2'])
relax> model_selection(method='AIC', modsel_run='hybrid_select', runs=['Dom1', 'Dom2', 'hybrid'])
AIC model selection.

Instance 0.

Run                  Num_params_(k)       Num_data_sets_(n)    Chi2                 Criterion          
Dom1                 9                    12                   3.70303              21.70303            
Dom2                 10                   12                   1.89550              21.89550           
Traceback (most recent call last):
  File "relax", line 454, in ?
    Relax()
  File "relax", line 166, in __init__
    self.interpreter.run()
  File "/local/home/viochem/relax/prompt/interpreter.py", line 213, in run
    run_script(intro=self.relax.intro_string, local= self.local, script_file=self.relax.script_file, quit=1)
  File "/local/home/viochem/relax/prompt/interpreter.py", line 388, in run_script
    console.interact(intro, local, script_file, quit)
  File "/local/home/viochem/relax/prompt/interpreter.py", line 340, in interact_script
    execfile(script_file, local)
  File "my_hybrid.py", line 9, in ?
    model_selection('AIC','hybrid_select',['Dom1','Dom2','hybrid'])
  File "/local/home/viochem/relax/prompt/model_selection.py", line 121, in model_selection
    self.relax.generic.model_selection.select(method=method, modsel_run=modsel_run, runs=runs)
  File "/local/home/viochem/relax/generic_fns/model_selection.py", line 217, in select
    k, n, chi2 = self.model_statistics [run](run=run, instance=i, global_stats=global_stats)
  File "/local/home/viochem/relax/specific_fns/hybrid.py", line 117, in model_statistics
    for i in xrange(num):
NameError: global name 'num' is not defined


Is this a bug, or am I doing something totally wrong?



Alex



Related Messages


Powered by MHonArc, Updated Fri Nov 10 10:00:29 2006