The Richard Gilder Graduate School embraces graduate training, post- doctoral fellowships, and undergraduate training programs at the Museum, through both independent activities and partnerships with universities.

Searching the local neighborhood (again)

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.

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? *)

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)

Finally make sure that the best trees that we found are saved for future use (we will need them later).