CGI/Perl Guide | Learning Center | Forums | Advertise | Login
Site Search: in

  Main Index MAIN
INDEX
Search Posts SEARCH
POSTS
Who's Online WHO'S
ONLINE
Log in LOG
IN

Home: Perl Programming Help: Intermediate: Re: [Kevone07] DNA Sequence Count Perl Program HELP!: Edit Log



Chris Charley
User

Oct 26, 2012, 8:17 AM


Views: 7813
Re: [Kevone07] DNA Sequence Count Perl Program HELP!

Maybe the attached files will give you an idea. t33.pl is the program file, fasta_dat.txt is the data file and junk.txt is the output from the program.

I used a smaller window size because my sample fasta file sequences weren't as long as the ones you have and the skip is smaller too. You would just need to make those adjustments for your program.

Look it over and ask any questions.

Hope this helps.

Update: Made a change to t33.pl

if ($begin_pos + $window_size > $length)

to

if ($begin_pos + $window_size >= $length)

Update2: Just to show a possible output ordered by skew, (uses SQLite database, pretty simple). File is bio_stat_db.pl.


Update3: Made adjustments to t33.pl and bio_stat_db.pl. Changed from $window_size to $size_read for skew calculations. Laurent caught this error (later in this thread). Output is now at the end of the 2 program files (after the __END__ token).


(This post was edited by Chris Charley on Oct 29, 2012, 8:20 AM)
Attachments: fasta_dat.txt (19.4 KB)
  bio_stat_db.pl (4.90 KB)
  t33.pl (9.98 KB)


Edit Log:
Post edited by Chris Charley (User) on Oct 26, 2012, 8:18 AM
Post edited by Chris Charley (User) on Oct 26, 2012, 8:20 AM
Post edited by Chris Charley (User) on Oct 26, 2012, 8:23 AM
Post edited by Chris Charley (User) on Oct 26, 2012, 8:26 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 6:26 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 6:29 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 7:48 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 7:50 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 10:48 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 11:13 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 11:16 AM
Post edited by Chris Charley (User) on Oct 27, 2012, 1:17 PM
Post edited by Chris Charley (User) on Oct 27, 2012, 7:00 PM
Post edited by Chris Charley (User) on Oct 28, 2012, 8:31 AM
Post edited by Chris Charley (User) on Oct 29, 2012, 8:20 AM


Search for (options) Powered by Gossamer Forum v.1.2.0

Web Applications & Managed Hosting Powered by Gossamer Threads
Visit our Mailing List Archives