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

View of /FigKernelScripts/make_protein_weights_and_rates.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (as text) (annotate)
Sat Jan 13 03:49:35 2007 UTC (13 years, 4 months ago) by overbeek
Branch: MAIN
weights and rates stuff

use proml;
use gjoseqlib;
use gjonewicklib;

my $usage = "usage: make_protein_weights_and_rates Alignment Tree > weights.and.rates";

(
 ($ali  = shift @ARGV) &&
 ($tree = shift @ARGV)
)
    || die $usage;

my @align = &gjoseqlib::read_fasta($ali);
my $treeP = &gjonewicklib::read_newick_tree($tree);

my($cats,$weights) = proml::estimate_protein_site_rates(\@align,$treeP);
print "$cats->[0]\n",
    join(" ",@{$cats->[1]}),"\n",
    $cats->[2],"\n",
    "$weights\n";

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3