Beefy Boxes and Bandwidth Generously Provided by pair Networks
more useful options
 
PerlMonks  

comment on
( [id://3333]=superdoc: print w/replies, xml ) Need Help??
The code that i wrote for mutation of a hardcoded string based on random number isnt working properly. Heres what the code does. The string to be mutated has already been provided inside the code(line 30). The string is made of 20 alphabets(line 5-24). %nuc (line 4)is a 24x24 matrix made in hash. The program generates a random number between 0 and 1 for length times the input string (length $string). It then checks this value in the corresponding hash %nuc. In our case since the $string has been hardcoded consisting of only "A", therefore each time a random number is generated the program goes to line 17 only. It then traverses the line 17 horizontally until a number greater than the random number is encountered. The value of the counter is noted when a number greater than the random number is encountered. For example in case our our $string which has only "A", it will mutate to another alphabet only if the random number is greater than .9799 otherwise it will remain as "A" only. The biggest problem with this code is that hashes dont get printed in the same order you store them,this shuffling leads to erratic results. I tried rearranging the rows of hash %nuc but it still didnt help. Another error is that Im getting an alphabet more in the output string than the hardcoded string. 
#!/usr/bin/perl -w

use Time::HiRes qw(usleep nanosleep);



%nuc = ( 

F  => [qw( .0001 .0001 .0001 .0000 .0000 .0000 .0000 .0001 .0002 .0008
+ .0006 .0000 .0004 .9944 .0000 .0000 .0001 .0003 .0028 .0000)],

T  => [qw( .0022 .0002 .0013 .0004 .0001 .0003 .0002 .0002 .0001 .0011
+ .0002 .0008 .0006 .0001 .0005 .0032 .9874 .0000 .0002 .0009)],

N  => [qw( .0004 .0001 .9867 .0036 .0000 .0004 .0006 .0006 .0021 .0003
+ .0001 .0013 .0000 .0001 .0002 .0020 .0009 .0001 .0004 .0001)],

V  => [qw( .0013 .0002 .0001 .0001 .0003 .0002 .0002 .0003 .0003 .0057
+ .0011 .0001 .0017 .0001 .0003 .0002 .0010 .0000 .0002 .9866)], 

K  => [qw( .0002 .0037 .0025 .0006 .0000 .0012 .0007 .0002 .0002 .0004
+ .0001 .9858 .0020 .0000 .0003 .0008 .0011 .0000 .0001 .0001)],

E  => [qw( .0010 .0000 .0007 .0056 .0000 .0035 .9865 .0004 .0002 .0003
+ .0001 .0004 .0001 .0000 .0003 .0004 .0002 .0000 .0001 .0002)],

Y  => [qw( .0001 .0000 .0003 .0000 .0003 .0000 .0001 .0000 .0004 .0001
+ .0001 .0000 .0000 .0021 .0000 .0001 .0001 .0002 .9960 .0001)],

Q  => [qw( .0003 .0010 .0004 .0005 .0000 .9901 .0027 .0001 .0024 .0001
+ .0003 .0006 .0004 .0000 .0006 .0002 .0002 .0000 .0000 .0001)],

I  => [qw( .0002 .0002 .0003 .0001 .0002 .0001 .0002 .0000 .0000 .9915
+ .0009 .0002 .0012 .0007 .0000 .0001 .0007 .0000 .0001 .0033)],

C  => [qw( .0001 .0001 .0000 .0000 .9987 .0000 .0000 .0000 .0001 .0001
+ .0000 .0000 .0000 .0000 .0001 .0005 .0001 .0000 .0000 .0002)],

L  => [qw( .0003 .0001 .0003 .0000 .0000 .0006 .0001 .0001 .0004 .0022
+ .9871 .0002 .0045 .0013 .0003 .0001 .0003 .0004 .0002 .0015)],

M  => [qw( .0001 .0001 .0000 .0000 .0000 .0002 .0000 .0000 .0000 .0005
+ .0009 .0005 .9968 .0001 .0000 .0001 .0002 .0000 .0000 .0005)],

A  => [qw( .9799 .0002 .0009 .0010 .0003 .0008 .0017 .0021 .0002 .0006
+ .0004 .0002 .0006 .0002 .0022 .0035 .0032 .0000 .0002 .0018)],

W  => [qw( .0000 .0002 .0000 .0000 .0000 .0000 .0000 .0000 .0000 .0000
+ .0000 .0000 .0000 .0001 .0000 .0001 .0000 .9995 .0001 .0000)],

S  => [qw( .0028 .0011 .0034 .0007 .0011 .0004 .0006 .0016 .0002 .0002
+ .0001 .0007 .0004 .0003 .0017 .9800 .0038 .0005 .0002 .0002)],

P  => [qw( .0013 .0005 .0002 .0001 .0001 .0008 .0003 .0002 .0005 .0001
+ .0002 .0002 .0001 .0001 .9935 .0012 .0004 .0000 .0000 .0002)],

H  => [qw( .0001 .0008 .0018 .0003 .0001 .0020 .0001 .0000 .9933 .0000
+ .0001 .0001 .0000 .0002 .0003 .0001 .0001 .0001 .0004 .0001)],

D  => [qw( .0006 .0000 .0042 .9869 .0000 .0006 .0053 .0006 .0004 .0001
+ .0000 .0003 .0000 .0000 .0001 .0005 .0003 .0000 .0000 .0001)],

R  => [qw( .0001 .9929 .0001 .0000 .0001 .0010 .0000 .0000 .0010 .0003
+ .0001 .0019 .0004 .0001 .0004 .0006 .0001 .0008 .0000 .0001)],

G  => [qw( .0021 .0001 .0012 .0011 .0001 .0003 .0007 .9906 .0001 .0000
+ .0001 .0002 .0001 .0001 .0003 .0021 .0003 .0000 .0000 .0005)]);









#SVEQRISTDIGQAYQLQGLGSNLRSIRSKTGAGEVNYIDAAKSVNDNQLLAEIG

$string='AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA';

@seq=split(//,$string);

$size=scalar(@seq);

print "$size\n";



foreach $seq (@seq)

{ 

&mutate($seq) ;

print "##########SEQ=$seq\tRANDOM=$ran\tK=$k################";

usleep (1000000);

 }



$k=0;

$p=0;

sub mutate # Start Of Sub

{

$ran=`awk 'BEGIN {{srand()} {print rand()}}'`;

for ($j = 1; $j <= 20;$j++)

   { $mut[$j]=$nuc{$seq}[$j]; }

for ($i=1; $i<=20; $i++)

{

        print "$mut[$i]\t$ran\t";   

    if ($ran < $mut[$i])

    {#print "\n$mut[$i] \t hello";

        push (@names, $i);

        print "i=$i";

         #push (@names, '-');

    }

    last if $ran < $mut[$i];

        print"\tran=$ran\t muti=$mut[$i]";

        $k =$mut[$i];

        print "K = $k";

} $p++; #counter for serial number

} # End Of Sub





# $nam=join ' ',@names;

 # for ($j = 0; $j <=$size;$j++)

 # {print $names[$j];}

 print"\n";

for ($j = 0; $j < $size;$j++)

{

 if($names[$j]==1)

  {$names[$j]="A";}

  elsif($names[$j]==2)

  {$names[$j]="R";}

  elsif($names[$j]==3)

  {$names[$j]="N";}

  elsif($names[$j]==4)

  {$names[$j]="D";}

  elsif($names[$j]==5)

  {$names[$j]="C";}

  elsif($names[$j]==6)

  {$names[$j]="Q";}

  elsif($names[$j]==7)

  {$names[$j]="E";}

  elsif($names[$j]==8)

  {$names[$j]="G";}

  elsif($names[$j]==9)

  {$names[$j]="H";}

  elsif($names[$j]==10)

  {$names[$j]="I";}

  elsif($names[$j]==11)

  {$names[$j]="L";}

  elsif($names[$j]==12)

  {$names[$j]="K";}

  elsif($names[$j]==13)

  {$names[$j]="M";}

  elsif($names[$j]==14)

  {$names[$j]="F";}

  elsif($names[$j]==15)

  {$names[$j]="P";}

  elsif($names[$j]==16)

  {$names[$j]="S";}

  elsif($names[$j]==17)

  {$names[$j]="T";}

  elsif($names[$j]==18)

  {$names[$j]="W";}

  elsif($names[$j]==19)

  {$names[$j]="Y";}

  else {$names[$j]="V";}

  }

   



print "orig = $string   \nmuta = ";

for ($j = 0; $j <=$size;$j++)

 {print $names[$j];}

 @arr = %nuc;

 print "@arr";

[download]

In reply to string mutation script behaving erraticaly by zing

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post; it's "PerlMonks-approved HTML":


  • Are you posting in the right place? Check out Where do I post X? to know for sure.
  • Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
    <code> <a> <b> <big> <blockquote> <br /> <dd> <dl> <dt> <em> <font> <h1> <h2> <h3> <h4> <h5> <h6> <hr /> <i> <li> <nbsp> <ol> <p> <small> <strike> <strong> <sub> <sup> <table> <td> <th> <tr> <tt> <u> <ul>
  • Snippets of code should be wrapped in <code> tags not <pre> tags. In fact, <pre> tags should generally be avoided. If they must be used, extreme care should be taken to ensure that their contents do not have long lines (<70 chars), in order to prevent horizontal scrolling (and possible janitor intervention).
  • Want more info? How to link or How to display code and escape characters are good places to start.
Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?

www.com | www.net | www.org
Chatterbox?
and the web crawler heard nothing...

How do I use this? • Last hour • Other CB clients
Other Users?
Others rifling through the Monastery: (6)
As of 2024-04-25 18:06 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found