Author: bugman Date: Thu May 1 17:05:05 2008 New Revision: 6046 URL: http://svn.gna.org/viewcvs/relax?rev=6046&view=rev Log: Updated the modsel.py script to the new relax design. Modified: 1.3/sample_scripts/modsel.py Modified: 1.3/sample_scripts/modsel.py URL: http://svn.gna.org/viewcvs/relax/1.3/sample_scripts/modsel.py?rev=6046&r1=6045&r2=6046&view=diff ============================================================================== --- 1.3/sample_scripts/modsel.py (original) +++ 1.3/sample_scripts/modsel.py Thu May 1 17:05:05 2008 @@ -24,30 +24,25 @@ ######################################## -# Nuclei type -nuclei('N') +# Set the data pipe names. +pipes = ['m0', 'm1', 'm2', 'm3', 'm4', 'm5', 'm6', 'm7', 'm8', 'm9'] -# Set the run names. -runs = ['m0', 'm1', 'm2', 'm3', 'm4', 'm5', 'm6', 'm7', 'm8', 'm9'] - -# Loop over the run names. -for name in runs: +# Loop over the data pipe names. +for name in pipes: print "\n\n# " + name + " #" - # Create the run. + # Create the data pipe. pipe.create(name, 'mf') - # Reload precalculated results from the file 'm1/results', etc. - results.read(run=name, file='results', dir=name) + # Reload precalculated results from the files 'm1/results', etc. + results.read(file='results', dir=name) # Model elimination. eliminate() # Model selection. -pipe.create('aic', 'mf') -model_selection('AIC', 'aic') +model_selection(method='AIC', modsel_pipe='aic') # Write the results. -state.save('save', force=1) -results.write(run='aic', file='results', force=1) - +state.save('save', force=True) +results.write(file='results', force=True)