[Bio] / FigKernelScripts / make_protein_weights_and_rates.pl Repository:
ViewVC logotype

Annotation of /FigKernelScripts/make_protein_weights_and_rates.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (view) (download) (as text)

1 : overbeek 1.1
2 :     use proml;
3 :     use gjoseqlib;
4 :     use gjonewicklib;
5 :    
6 :     my $usage = "usage: make_protein_weights_and_rates Alignment Tree > weights.and.rates";
7 :    
8 :     (
9 :     ($ali = shift @ARGV) &&
10 :     ($tree = shift @ARGV)
11 :     )
12 :     || die $usage;
13 :    
14 :     my @align = &gjoseqlib::read_fasta($ali);
15 :     my $treeP = &gjonewicklib::read_newick_tree($tree);
16 :    
17 :     my($cats,$weights) = proml::estimate_protein_site_rates(\@align,$treeP);
18 :     print "$cats->[0]\n",
19 :     join(" ",@{$cats->[1]}),"\n",
20 :     $cats->[2],"\n",
21 :     "$weights\n";

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3