Hi all!
Tatiana
I'm running a sensitivity analysis of total evidence and not whether
the weighing of the morphology is the same way that the molecular data
This is my script
read("18s.fas", "rbcl.fas", "rps.fas", "morfo.ss")
set(root:"t1")
transform((names:("18s.fas"), tcm:"111.txt"),(names:("rbcl.fas"),
tcm:"111.txt"),(names:("rps.fas"), tcm:"111.txt"))
store("misdatos")
build(70)
select()
swap()
select()
transform((all, static_approx))
report("test1.mtr", phastwinclad)
echo ("tot Results", output:"Sensitivity_results.txt")
report ("Sensitivity_results.txt", treestats)
select(best:0)
use("misdatos")
select(characters,names:("18s.fas"))
build(70)
select()
swap()
select()
echo("18s.fas Results", output:"Sensitivity_results.txt")
report("Sensitivity_results.txt", treestats)
transform((all, static_approx))
report("test2.mtr", phastwinclad)
use("misdatos")
select(best:0)
select(characters,names:("rbcl.fas"))
build(70)
select()
swap()
select()
echo ("rbcl.fas Results", output:"Sensitivity_results.txt")
report ("Sensitivity_results.txt", treestats)
transform((all, static_approx))
report ("test3.mtr", phastwinclad)
use("misdatos")
select(best:0)
select(characters,names:("rps.fas"))
build(70)
select()
swap()
select()
echo ("rps.fas Results", output:"Sensitivity_results.txt")
report ("Sensitivity_results.txt", treestats)
transform((all, static_approx))
report ("test4.mtr", phastwinclad)
quit()
thanks
Tatiana
ward wheeler
--
You received this message because you are subscribed to the Google Groups "poy4" group.
To post to this group, send email to po...@googlegroups.com.
To unsubscribe from this group, send email to poy4+uns...@googlegroups.com.
For more options, visit this group at http://groups.google.com/group/poy4?hl=en.