Searching the local neighborhood (again) main content.

Searching the local neighborhood (again)

  1. Although  wipe () seems like a handy command, in interactive tests it would be desirable to have other means to store the state of the program and continue from there later on. The  store () and  use ()commands serve this purpose.
  2. store ("name") will store the current state of the program (that is all the data, trees, and support values) under the name label. Conversely, the command  use ("name") sets the current state of the program to that which was stored under the label name. For example:
                (* These are comments in a POY script *)            
                wipe ()            
                read ("course.fasta", "built_trees.tree")            
                report (data, treestats)            
                (* At this point you can see that we have loaded some molecular data as well as initial trees. We will store it under the label initial_data. *)            
                store ("initial_data")            
                select ()            
                (* We will now modify the data, using a command that will be introduced later. *)            
                transform (static_approximation)            
                report (data, treestats)            
                (* Notice that now we are only analyzing static homology characters. *)            
                store ("static_data")            
                use ("initial_data")            
                report (data, treestats)            
                (* See the characters are the same as the originals? *)            
                use ("static_data")            
                report (data, treestats)            
                (* Is it what you expected? *)
  3. We can now do things more efficiently. Lets compare SPR and TBR again:
                wipe ()            
                read ("course.fasta", "built_trees.tree")           
                store ("initial")            
                swap (tbr)            
                report (treestats)            
                use ("initial")            
                swap (spr)            
                report (treestats)
  4. Finally make sure that the best trees that we found are saved for future use (we will need them later).
                report ("swap_course.tree")            
                exit ()