Beefy Boxes and Bandwidth Generously Provided by pair Networks
Keep It Simple, Stupid
 
PerlMonks  

Re: alphabet counting

by Cristoforo (Deacon)
on Jun 03, 2012 at 00:07 UTC ( #974073=note: print w/replies, xml ) Need Help??


in reply to alphabet counting

Here is a solution to counting the letters taken from the BioPerl webpage, the 'HOWTO:Beginners' link on the opening page, then link 17 - (Obtaining basic sequence statistics). That led to the docs on Bio::Tools::SeqStats.

#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; use Bio::Tools::SeqStats; my @prot = qw/ A C D E F G H I K L M N P Q R S T V W Y /; my $outputfile = "countaa"; open my $OUT, ">", $outputfile or die "Can't open file \"$outputfile\" to write to $!\n\n"; my $proteinio=Bio::SeqIO->new (-file=>"ec 1.1.1.fasta",-format=>'fasta +'); while(my $seq = $proteinio->next_seq() ) { my $seq_stats = Bio::Tools::SeqStats->new(-seq => $seq); my $count = $seq_stats->count_monomers(); print $OUT join(' ', map {$_ || 0} @$count{ @prot }), "\n"; } close $OUT or die $!;

Hope this is of some help,

Chris

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://974073]
help
Chatterbox?
[shmem]: Discipulus: while (readdir DIR) works fine
[Discipulus]: Eataly vs World(coffe) ;=)
[Discipulus]: shemem i wrote it directly.. but when i tested it i got a lot of 'undefined value..'
[Corion]: Oh, poor Trump can't get a break. Even the Vatican only sends the second in command to greet him...
[shmem]: Discipulus: readdir has the same semantics as readline (a.k.a <>). Have another coffee ☕ ;-)
[Discipulus]: ah! just with my old system perl (5.8.8) was needed
[Corion]: Has anybody toyed around with Mastodon? There is a Perl module for it, but I'm not sure that it's a yak worth chasing
[Corion]: On the upside, maybe running that for my family would be interesting, or at least for myself
[Corion]: I'm not fond of using a server hosted by other people, even if they claim to be GNU ;)

How do I use this? | Other CB clients
Other Users?
Others wandering the Monastery: (8)
As of 2017-05-26 07:31 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?