Re^8: how to output the final Tree into a file?

#!/usr/bin/perl

$ENV{CLUSTALDIR} = 'C:\\\\\\\\\\\\\\\\\\;
use Bio::Tools::Run::Alignment::Clustalw;

my $file = <>; # Get file name (fasta file) from command prompt.

#  Build a clustalw alignment factory
my $factory = Bio::Tools::Run::Alignment::Clustalw->new(-matrix => 'BL
+OSUM',-outfile=>'out', -output =>'clustalw'); 
my $ktuple = 2;
$factory->ktuple($ktuple);

my $inseq = Bio::SeqIO->new(
                    -file=> "<$file",
                    -format=> 'fasta'
                   );
my $seq;
my @seq_array;
while ($seq = $inseq->next_seq) {
    push(@seq_array, $seq);
}
# Now we do the actual alignment.
my $seq_array_ref = \@seq_array;
my $aln = $factory->align($seq_array_ref);

#constarcting the Tree
my $alnio = Bio::AlignIO->new(-file => 'out', -format=>'clustalw');

my $dfactory = Bio::Tree::DistanceFactory->new(-method => 'NJ');

my $stats = Bio::Align::ProteinStatistics->new;

my $treeout = Bio::TreeIO->new(-format => 'newick');

while( my $aln = $alnio->next_aln ) 
{
 my $mat = $stats->distance(-method => 'Kimura',
                            -align  => $aln);
 my $tree = $dfactory->make_tree($mat);
 $treeout->write_tree($tree);
}
my $outfile = 'C:\Users\Ioulia\Desktop\outree.tre';
my $treeout = Bio::TreeIO->new(-format => 'newick', -file => ">$outfil
+e");
exit;
[download]

Comment on Re^8: how to output the final Tree into a file? Download Code

Replies are listed 'Best First'.
Re^9: how to output the final Tree into a file? by poj (Abbot) on Mar 16, 2016 at 07:04 UTC
I think you need to create the `$treeout` object with `-file` first so that it uses it as the stream for the `write_tree` method `my $outfile = 'C:\Users\AnyUser\Desktop\outree.tre'; my $treeout = Bio::TreeIO->new( -format => 'newick', -file => ">$outfile" ); while( my $aln = $alnio->next_aln ) { my $mat = $stats->distance(-method => 'Kimura', -align => $aln); my $tree = $dfactory->make_tree($mat); $treeout->write_tree($tree); } exit;` [download] poj	[reply] [d/l] [select]
Re^10: how to output the final Tree into a file? by dimitris852 (Acolyte) on Mar 16, 2016 at 08:40 UTC
You are absolutely right! Bless you my friend!	[reply]