[Bioperl-l] phylip - ProtDist Problem
Nathan Haigh
nathanhaigh at ukonline.co.uk
Tue Aug 31 07:26:53 EDT 2004
I am trying to run the following script at the bottom which takes a
SimpleAlign object and runs SeqBoot on it. I would then like to take this
output and run ProtDist, Neighbor and Consense to produce a bootstrap
consensus tree. However, ProtDist just seems to fail when it gets to this
line in ProtDist.pm:
open(PROTDIST,"|".$self->executable); # line 399 in the
latest version of ProtDist.pm from the CVS
I am running this on WinXP from within a cgi script and have phylip v3.61. I
have no idea why this should fail as the similar line in SeqBoot works fine!
Has anyone come across this problem or know why this might be happening!?
Regards
Nathan
# use seqboot to generate bootstrap alignments
my @params_seqboot = ('datatype'=>'SEQUENCE', 'replicates'=>100);
my $seqboot_factory =
Bio::Tools::Run::Phylo::Phylip::SeqBoot->new(@params_seqboot);
my $aln_ref = $seqboot_factory->run($aln);
# create distance matrices and construct trees using neighbor
my $protdist_factory = Bio::Tools::Run::Phylo::Phylip::ProtDist->new();
my $neighbor_factory = Bio::Tools::Run::Phylo::Phylip::Neighbor->new();
my @tree;
foreach my $a (@{$aln_ref}) {
my $matrix = $protdist_factory->run($a);
push @tree, $neighbor_factory->run($matrix);
}
# use consense to get a final tree
my $consense_factory = Bio::Tools::Run::Phylo::Phylip::Consense->new();
my $tree = $consense_factory->run(\@tree);
# output treefile
my $outfilename = "outtrees_ur/homologene.nh"; # unrooted tree
my $outtree = new Bio::TreeIO('-format' => 'newick',
'-file' => ">$outfilename");
$outtree->write_tree($tree);
---
avast! Antivirus: Outbound message clean.
Virus Database (VPS): 0435-2, 28/08/2004
Tested on: 31/08/2004 12:26:50
avast! is copyright (c) 2000-2003 ALWIL Software.
http://www.avast.com
More information about the Bioperl-l
mailing list