Parse a genbank file using regular expressions
...$protein .= $1; $protein =~ s/\s*//g; print "$protein\n"; $trans = 0; } else { print "Problems: end of translation product n...2932 days ago
A multilayer perceptron (MLP) neural network in Perl
...######### my $numInputs = 3; my $numPatterns = 4; ######################################## my $patNum; my $errThisPat; my $outPred; my $RMSerror; my @trainInputs; my...2931 days ago
Generating a random string with Perl
...ngs of a given length sub generate_random_string { my $length_of_randomstring=shift;# the length of...A'..'Z','0'..'9','_'); my $random_string; foreach (1..$length_of_randomstring) { # rand @chars w...2922 days ago
picard tools command to get some insert statistics
...icard tools to get some insert statistics to see whether our reads seem to be in the correct place #module load picard/2.0.1 java -Xmx16g -XX:PermSize=8g -jar $PICARD_HOME/picar...1352 days ago
Collision free write with Perl
...because of multicore usesage sub collision_free_write { my($outFile, $msg) = @_; open my $ofh, ">...flock $ofh, LOCK_EX or die "$0 [$$]: flock: $!"; print $ofh "$msg\n" or die "$0 [$$]: wr...2476 days ago
Genetic Algorithms demonstration with word DNA in Perl
#!/usr/bin/perl -w # GA demonstration with word DNA (512 bits) use strict; use Data::Dumper; # individuals in the population my $popsize = 1024; # a good...2387 days ago
Plot the clock using Lastz -gerenal outfile
...my $R = Statistics::R->new() ; $R->startR ; my $fileN=$ARGV[0]; my $mSize=$ARGV[1]; open (my $fh2,...; if ($arr[7] eq "-") { my $len=$arr[5]-$arr[4]; #next if $len < $mSize; my @chr = split '\_',...2314 days ago
Plot custom gene density with R
library(karyoploteR) pp2264 days ago
Perl script to find palindromic regions in DNA sequences
use strict; use warnings; my $pp = qr/(?: (\w) (?1) \g{-1} | \w? )/ix; my $filename = $ARGV[0]; open(my $fh, '2169 days ago
2159 days ago