mailr4686 - /1.3/specific_fns/model_free/model_free.py


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

Header


Content

Posted by edward on January 13, 2008 - 18:50:
Author: bugman
Date: Sun Jan 13 18:39:35 2008
New Revision: 4686

URL: http://svn.gna.org/viewcvs/relax?rev=4686&view=rev
Log:
Converted the overfit_deselect() model-free method to the new relax design.


Modified:
    1.3/specific_fns/model_free/model_free.py

Modified: 1.3/specific_fns/model_free/model_free.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/specific_fns/model_free/model_free.py?rev=4686&r1=4685&r2=4686&view=diff
==============================================================================
--- 1.3/specific_fns/model_free/model_free.py (original)
+++ 1.3/specific_fns/model_free/model_free.py Sun Jan 13 18:39:35 2008
@@ -2583,40 +2583,40 @@
             raise RelaxFault
 
 
-    def overfit_deselect(self, run):
+    def overfit_deselect(self):
         """Function for deselecting residues without sufficient data to 
support minimisation"""
 
         # Test sequence data exists.
-        if not relax_data_store.res.has_key(run):
-            raise RelaxNoSequenceError, run
-
-        # Loop over residue data:
-        for residue in relax_data_store.res[run]:
+        if not exists_mol_res_spin_data():
+            raise RelaxNoSequenceError
+
+        # Loop over the sequence.
+        for spin in spin_loop():
             # Skip unselected data:
-            if not residue.select:
+            if not spin.select:
                 continue
 
             # Check for data structure.
-            if not hasattr(residue, 'relax_data'):
-                residue.select = 0
+            if not hasattr(spin, 'relax_data'):
+                spin.select = 0
                 continue
 
             # Require 3 or more data points
-            if len(residue.relax_data) < 3:
-                residue.select = 0
+            if len(spin.relax_data) < 3:
+                spin.select = 0
                 continue
 
             # Require at least as many data points as params to prevent 
over-fitting
-            if hasattr(residue, 'params'):
-                if len(residue.params) > len(residue.relax_data):
-                    residue.select = 0
+            if hasattr(spin, 'params'):
+                if len(spin.params) > len(spin.relax_data):
+                    spin.select = 0
                     continue
 
             # Test for structural data if required
             if hasattr(relax_data_store, 'diff') and 
relax_data_store.diff.has_key(run):
                 if relax_data_store.diff[run].type == 'spheroid' or 
relax_data_store.diff[run].type == 'ellipsoid':
-                    if not hasattr(residue, 'xh_vect'):
-                        residue.select = 0
+                    if not hasattr(spin, 'xh_vect'):
+                        spin.select = 0
                         continue
 
 




Related Messages


Powered by MHonArc, Updated Sun Jan 13 19:20:18 2008