meta-retention index

6 views
Skip to first unread message

Caroline Greve

unread,
Nov 25, 2009, 8:07:17 AM11/25/09
to po...@googlegroups.com
Hello All
I´m performing a sensitivity analysis using the script sugested in the "commands" file of POY 4.0:
(* sensitivity analysis *)
read("9.fas")
set(root:"t1")
store("original_data")
transform(tcm:"111.txt")
build(100)
swap(timeout:3600)
select()
report("111.tre",trees:(total) ,"111con.tre",consensus,
"111con.pdf",graphconsensus)
use("original_data")
transform(tcm:"112.txt")
build(100)
swap(timeout:3600)
select()
report("112.tre",trees:(total),"112con.tre",consensus,
"112con.pdf",graphconsensus)
exit()
How can I calculate the MRI  with these results?  Specifically it´s not clear to me how to obtain the maximum cost of fragments.

Ward Wheeler

unread,
Nov 25, 2009, 11:57:11 AM11/25/09
to po...@googlegroups.com
Caroline,
    You can set the weight to (-1) for the fragments and run them separately.  That way you are 
searching for the 'worst' tree.
W


--

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.



--
Ward Wheeler
Division of Invertebrate Zoology
American Museum of Natural History
Central Park West at 79th Street
New York, NY 10024-5192
USA
whe...@amnh.org
http://research.amnh.org/scicomp/ward_wheeler.html

Caroline Greve

unread,
Nov 26, 2009, 7:40:02 AM11/26/09
to po...@googlegroups.com
Thank you very much for your help!

2009/11/25 Ward Wheeler <ward.w...@gmail.com>
Reply all
Reply to author
Forward
0 new messages