Beefy Boxes and Bandwidth Generously Provided by pair Networks
XP is just a number

Bioperl _ SeqStats

by MVRS (Acolyte)
on Feb 06, 2013 at 05:35 UTC ( #1017336=perlquestion: print w/replies, xml ) Need Help??
MVRS has asked for the wisdom of the Perl Monks concerning the following question:

i tried running bioperl script to find hydropathicity of a protein sequence , and am unable to fix the error >

Error : Can't use string ("Bio::Tools::SeqStats") as a HASH ref while "strict refs" in use at /usr/share/perl5/Bio/Tools/ line 726.

#!/usr/bin/perl #use SeqStats; use Data::Dumper; use Bio::SeqIO ; use Bio::Seq ; use Bio::Tools::SeqStats; use Bio::Root::Root; # Calculate mean Kyte-Doolittle hydropathicity (aka "gravy" score) my $prot = Bio::PrimarySeq->new(-seq=>'MSFVLVAPDMLATAAADVVQIGSAVSAGS +',-format => 'Fasta',-alphabet=>'protein', -id=>'test'); $gravy = Bio::Tools::SeqStats->hydropathicity($seqobj); print "might be hydropathic" if $gravy > 1;

Replies are listed 'Best First'.
Re: Bioperl _ SeqStats
by vinoth.ree (Monsignor) on Feb 06, 2013 at 06:25 UTC

    You need to give $prot as argument to the hydropathicity() function as follows,

    my $gravy = Bio::Tools::SeqStats->hydropathicity($prot);

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: perlquestion [id://1017336]
Approved by vinoth.ree
and all is quiet...

How do I use this? | Other CB clients
Other Users?
Others having an uproarious good time at the Monastery: (6)
As of 2017-09-23 17:15 GMT
Find Nodes?
    Voting Booth?
    During the recent solar eclipse, I:

    Results (273 votes). Check out past polls.