Trouble in running Baseml
Entering edit mode
7.5 years ago
Biocode_user ▴ 30


I'm trying to do a nucleotide substitution analysis using a very simple model like Jukes and Cantor. I have a set of transcripts assembly and based on expression I have chosen the representative component of each and did a reciprocal blastn on the nucleotide sequences which is been jotted down from previous step. I did a pairwise alignment with clustalw for them. I have the alignment file in fasta format and now I'm trying to apply a simple nucleotide substitution model to proceed for further analysis. I'm using the perl wrapper for baseml and have the latest version of PAML and Bioperl installed.

I'm attaching the script that I'm trying to use and the error which I get. It would be greatly appreciable if someone could help on the same.

#!/usr/bin/perl -w

use Bio::Tools::Run::Phylo::PAML::Baseml;
use Bio::AlignIO;
use Data::Dumper;
my $alignio = Bio::AlignIO->new(-format => 'fasta',
             -file   => $ARGV[0]);

my $aln = $alignio->next_aln;
my $baseml = Bio::Tools::Run::Phylo::PAML::Baseml->new();
my ($rc,$parser) = $baseml->run();
my $result = $parser->next_result; 
my $MLmatrix = $result-> get_MLmatrix();
print Dumper($MLmatrix);


Can't call method "get_MLmatrix" on an undefined value at line 27, <GEN1> line 6
PAML RNA-Seq • 1.9k views

Login before adding your answer.

Traffic: 2311 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6