...
10 80
12 100
14 99
Make clusters
Once you have found a set of stable habitats you want to go with, then make the clusters (this is another set of commands, could be an old format, I might also be using habitats with 2, 3, 4 and 5 predicted habitats):
foreach f (All_simple_u_4 All_simple_u_3 All_simple_u_5 All_simple_u_2)
foreach? /home/software/python/python-2.5.1/bin/python ~/AdaptML_dir/latest_AdaptML/habitats/trunk/AdaptML_Analyze.py tree=./${f}_noroot.phy_phyml_tree.txt3 hab0=16 outgroup=ECK_1 write=./${f}_dir thresh=0.05
foreach? perl ~/bin/migration2color_codes.pl ./${f}_dir/migration.matrix ~/bin/color_template.txt > ./${f}_dir/color.file
foreach? mkdir ./${f}_dir/rand_trials
foreach? /home/software/python/python-2.5.1/bin/python ~/AdaptML_dir/latest_AdaptML/clusters/getstats/rand_JointML.py habitats=./${f}_dir/migration.matrix mu=./${f}_dir/mu.val tree=./${f}_noroot.phy_phyml_tree.txt3 outgroup=ECK_1 write=./${f}_dir/rand_trials/
foreach? /home/software/python/python-2.5.1/bin/python ~/AdaptML_dir/latest_AdaptML/clusters/getstats/GetLikelihoods.py ./${f}_dir/rand_trials/ ./${f}_dir/
foreach? /home/software/python/python-2.5.1/bin/python ~/AdaptML_dir/latest_AdaptML/clusters/trunk/JointML.py habitats=./${f}_dir/migration.matrix mu=./${f}_dir/mu.val color=./${f}_dir/color.file tree=./${f}_noroot.phy_phyml_tree.txt3 write=./${f}_dir/ outgroup=ECK_1 thresh=./${f}_dir/9900.file
foreach? end