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...2921 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...2920 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...2911 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...1341 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...2465 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...2376 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 '\_',...2303 days ago
Plot custom gene density with R
library(karyoploteR) pp2254 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, '2158 days ago
2148 days ago