mailr6046 - /1.3/sample_scripts/modsel.py


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

Header


Content

Posted by edward on May 01, 2008 - 17:16:
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)




Related Messages


Powered by MHonArc, Updated Thu May 01 17:20:16 2008