Hi all,
I have defined some array variables a
***************************************** -
@L=('CTT','CTC','CTA','CTG','TTA','TTG');
-
-
@S=('TCT','TCC','TCA','TCG','AGT','AGC');
-
-
@R=('CGT','CGC','CGA','CGG','AGA','AGG');
-
-
@V=('GTT','GTC','GTA','GTG');
-
-
@A=('GCT','GCC','GCA','GCG');
-
-
@G=('GGT','GGC','GGA','GGG');
-
-
@P=('CCT','CCC','CCA','CCG');
-
-
@T=('ACT','ACC','ACA','ACG');
-
-
@I=('ATT','ATC','ATA');
-
-
@F=('TTT','TTC');
-
-
@C=('TGT','TGC');
-
-
@Y=('TAT','TAC');
-
-
@Q=('CAA','CAG');
-
-
@N=('AAT','AAC');
-
-
@H=('CAT','CAC');
-
-
@E=('GAA','GAG');
-
-
@D=('GAT','GAC');
-
-
@K=('AAA','AAG');
-
-
@M=('ATG');
-
-
@W=('TGG');
-
************************************************** *****
and also created a hash to store the arrays in it and finally I wanted to do a foreach loop over the desired input say for example input=VEFCDAMP and to create all possible permutation of the word taken from the definition from dictionary as defined in the array.
The following is the perl script that I conceived to do the job but I am confused here how to use the forloop here to iterate through this and make permutation,
************************************************** ********** -
print "\n\n\t\#################### AA 2 PROTEIN #################### \n\n";
-
-
print "This script will convert your amino acid sequence to DNA Sequence\n\n";
-
-
print "ENTER THE FILENAME OF THE Amino acid SEQUENCE:= ";
-
-
$aafilename = <STDIN>;
-
-
chomp $aafilename;
-
-
unless ( open(aaFILE, $aafilename) ) {
-
-
print "Cannot open file \"$aafilename\"\n\n";
-
-
}
-
-
@aa = <aaFILE>;
-
-
close aaFILE;
-
-
$aa = join( '', @aa);
-
-
print " \nThe original AA file is:\n$aa \n";
-
-
$aa =~ s/\s//g;
-
-
@aa=$aa;
-
-
my $protein='';
-
-
my $codon;
-
-
@L=('CTT','CTC','CTA','CTG','TTA','TTG');
-
-
@S=('TCT','TCC','TCA','TCG','AGT','AGC');
-
-
@R=('CGT','CGC','CGA','CGG','AGA','AGG');
-
-
@V=('GTT','GTC','GTA','GTG');
-
-
@A=('GCT','GCC','GCA','GCG');
-
-
@G=('GGT','GGC','GGA','GGG');
-
-
@P=('CCT','CCC','CCA','CCG');
-
-
@T=('ACT','ACC','ACA','ACG');
-
-
@I=('ATT','ATC','ATA');
-
-
@F=('TTT','TTC');
-
-
@C=('TGT','TGC');
-
-
@Y=('TAT','TAC');
-
-
@Q=('CAA','CAG');
-
-
@N=('AAT','AAC');
-
-
@H=('CAT','CAC');
-
-
@E=('GAA','GAG');
-
-
@D=('GAT','GAC');
-
-
@K=('AAA','AAG');
-
-
@M=('ATG');
-
-
@W=('TGG');
-
-
-
-
$L=\@L;
-
-
$S=\@S;
-
-
$R=\@R;
-
-
$V=\@V;
-
-
$A=\@A;
-
-
$G=\@G;
-
-
$P=\@P;
-
-
$T=\@T;
-
-
$I=\@I;
-
-
$F=\@F;
-
-
$C=\@C;
-
-
$Y=\@Y;
-
-
$Q=\@Q;
-
-
$N=\@N;
-
-
$H=\@H;
-
-
$E=\@E;
-
-
$D=\@D;
-
-
$K=\@K;
-
-
$M=\@M;
-
-
$W=\@W;
-
-
-
-
%hash = (L=>$L,S=>$S,R=>$R,V=>$V,A=>$A,G=>$G,P=>$P,T=>$T,I=>$I,F=>$F,C=>$C,Y=>$Y,Q=>$Q,N=>$N,H=>$H,E=>$E,D=>$D,K=>$K,M=>$M,W=>$W);
-
-
-
-
$hash_ref=\%hash;
-
-
-
my @aa_split = split //, $aa;
-
-
-
-
foreach my $val (@{hash{@aa_split}})
-
-
{
-
print "@{$val} ";
-
-
}
-
print " \n";
-
************************************************** ***************
to run this program please create and save a txt file with the content of word VEFCDAMP in it.
I want my program to perform simillar task what the following program does it for
************************************************** ************** -
use strict;
-
use warnings;
-
my @V=('GTT','GTC','GTA','GTG');
-
my @A=('GCT','GCC','GCA','GCG');
-
my @E=('GAA','GAG');
-
my @F=('TTT','TTC');
-
my @G=('GGT','GGC','GGA','GGG');
-
my @H=('CAT','CAC');
-
my $i = 1;
-
foreach my $v (@V) {
-
foreach my $a (@A) {
-
foreach my $e (@E) {
-
foreach my $f (@F) {
-
foreach my $g (@G) {
-
foreach my $h (@H) {
-
print "$i $v$a$e$f$g$h\n";
-
$i++;
-
}
-
}
-
}
-
}
-
}
-
}
-
************************************************** ************
Thanks in advance,
B.Nataraj
28 3140
And PLEASE use code tags around any and all code you post in the forums. The link that my colleague provided show's their use as well.
Regards,
Moderator
Assuming everything else is good with your code: -
open (OUT, ">>" , 'path/to/outfile.txt') or die "$!";
-
my @aa_split = split //, $aa;
-
foreach my $val (@{hash{@aa_split}}){
-
print OUT "@{$val} ";
-
}
-
print " \n";
-
close OUT;
-
Assuming everything else is good with your code: -
open (OUT, ">>" , 'path/to/outfile.txt') or die "$!";
-
my @aa_split = split //, $aa;
-
foreach my $val (@{hash{@aa_split}}){
-
print OUT "@{$val} ";
-
}
-
print " \n";
-
close OUT;
-
Hi Kevin,
The code that you added, does not give the desired output of all permutation of the word "VEFCDAMP". This particular code is the one that I wrote as an extention of your previous code given in the thread ( http://bytes.com/forum/thread847179.html ) for fixed input like "VEFCDAMP" (your code is also pasted in this topic) . I tried here to extend the code for all possible input word (generalization) and the number of possible letter is only 20 (for my case) and the maximum size of any word never going to exceed 10 letters, thats why I thought of iterating only the desired letter in foreach loop to save computational time, otherwise I tried for all possible foreach loop for 20 letters and it goes on and on in my system.
Hope you can understand my problem by reading the previous thread and also this thread.
Thanks in advance,
B.Nataraj
I'm sorry but your current requirements are out of my range of experience. Trying to make a list of all possible permutations from an arbitrarily deep list of arrays is something I have no experience with.
Nataraj,
When you are taking input dynamically, it would be difficult to get all permutations if the length of input string varies. But if it is fixed at 8 as in VEFCDAMP, then you should be able to do it using 8 foreach loops as you have done in the second script.
Nataraj,
When you are taking input dynamically, it would be difficult to get all permutations if the length of input string varies. But if it is fixed at 8 as in VEFCDAMP, then you should be able to do it using 8 foreach loops as you have done in the second script.
Hi nithin & Kevin,
Ok then, Is there any other way to address this problem ? I mean otherthan using foreach loop.
Thanks ,
B.Nataraj
This is what I came up with although I am not sure its what you want and it means you will have to install the List::Permutor module and any of its dependencies.
This is a script with all extraneous code removed just so we can get the permutations given a list of letters. -
use List::Permutor;
-
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @loops = split '', $seq;
-
my @AoA;
-
foreach my $letter (@loops) {
-
push @AoA, @{$hash{$letter}};
-
}
-
my $p = List::Permutor->new(@AoA);
-
{
-
local $";
-
while(my @set = $p->next) {
-
print "@set\n";
-
}
-
}
-
I used a very small list of letters (VF) because for a long string of letters all the possible premutations is a large number and takes a little while to run. Basically 'VF' has 6 unique codons, so there is a possible 720 unique permutations (UP) just for those:
UP = 1*2*3*4*5*6 = 720
if you add one more unique codon the number obviously grows exponentially to 5,040. You will probably be OK until you get into the range of 10 unique codons in which case the UP is 3,628,800. A fast computer could figure out the permutations quickly, but writing them all to a file is what will take time.
This is what I came up with although I am not sure its what you want and it means you will have to install the List::Permutor module and any of its dependencies.
This is a script with all extraneous code removed just so we can get the permutations given a list of letters. -
use List::Permutor;
-
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @loops = split '', $seq;
-
my @AoA;
-
foreach my $letter (@loops) {
-
push @AoA, @{$hash{$letter}};
-
}
-
my $p = List::Permutor->new(@AoA);
-
{
-
local $";
-
while(my @set = $p->next) {
-
print "@set\n";
-
}
-
}
-
I used a very small list of letters (VF) because for a long string of letters all the possible premutations is a large number and takes a little while to run. Basically 'VF' has 6 unique codons, so there is a possible 720 unique permutations (UP) just for those:
UP = 1*2*3*4*5*6 = 720
if you add one more unique codon the number obviously grows exponentially to 5,040. You will probably be OK until you get into the range of 10 unique codons in which case the UP is 3,628,800. A fast computer could figure out the permutations quickly, but writing them all to a file is what will take time.
Hi Kevin,
Thanks once again for your help. It took some time for me to install the List::Permutor module into my system and so a delay in replying you, now your code is working fine but as you said it goes on exponential time for a larger size input. I have not much gone in detail of your code, Hope I soon understand it and will come back to you.
With Thanks,
B.Nataraj
Hi Kevin,
Thanks once again for your help. It took some time for me to install the List::Permutor module into my system and so a delay in replying you, now your code is working fine but as you said it goes on exponential time for a larger size input. I have not much gone in detail of your code, Hope I soon understand it and will come back to you.
With Thanks,
B.Nataraj
Hi Kevin,
Sorry about my last reply, I just carried away by seeing new code and implementing new module and running the code. Infact your new code is not doing the right job, To say it short to you , I wanted the output for the word "VF" (you took this word as sample in your code) like the one come out of the following code (This is also your code) , The total permutation would be 512. - use strict;
-
use warnings;
-
my @V=('GTT','GTC','GTA','GTG');
-
my @A=('GCT','GCC','GCA','GCG');
-
my @E=('GAA','GAG');
-
my @F=('TTT','TTC');
-
my @G=('GGT','GGC','GGA','GGG');
-
my @H=('CAT','CAC');
-
my $i = 1;
-
foreach my $v (@V) {
-
foreach my $a (@A) {
-
foreach my $e (@E) {
-
foreach my $f (@F) {
-
foreach my $g (@G) {
-
foreach my $h (@H) {
-
print "$i $v$f\n";
-
$i++;
-
}
-
}
-
}
-
}
-
}
-
}
-
-
Hope you can now better understand the problem. In mean time I am trying to change your code the way I wanted it is to be.
Thanks for your all help,
B.Nataraj
Hi Kevin,
Sorry about my last reply, I just carried away by seeing new code and implementing new module and running the code. Infact your new code is not doing the right job, To say it short to you , I wanted the output for the word "VF" (you took this word as sample in your code) like the one come out of the following code (This is also your code) , The total permutation would be 512.
Thanks for your all help,
B.Nataraj
The total permutation for 'VF' should not be 512. It will be 8. The script you have shown will just create multiple duplicates to produce 512 results, because of usage of unnecessary foreach loops. -
use strict;
-
use warnings;
-
my @V=('GTT','GTC','GTA','GTG');
-
my @A=('GCT','GCC','GCA','GCG');
-
my @E=('GAA','GAG');
-
my @F=('TTT','TTC');
-
my @G=('GGT','GGC','GGA','GGG');
-
my @H=('CAT','CAC');
-
my $i = 1;
-
-
foreach my $v (@V ){
-
foreach my $f (@F) {
-
print "$i $v$f\n";
-
$i++;
-
}
-
}
-
-
However, I haven't tested the output using List::Permutor. I will check it in my free-time and compare the output.
-Nithin
The total permutation for 'VF' should not be 512. It will be 8. The script you have shown will just create multiple duplicates to produce 512 results, because of usage of unnecessary foreach loops. -
use strict;
-
use warnings;
-
my @V=('GTT','GTC','GTA','GTG');
-
my @A=('GCT','GCC','GCA','GCG');
-
my @E=('GAA','GAG');
-
my @F=('TTT','TTC');
-
my @G=('GGT','GGC','GGA','GGG');
-
my @H=('CAT','CAC');
-
my $i = 1;
-
-
foreach my $v (@V ){
-
foreach my $f (@F) {
-
print "$i $v$f\n";
-
$i++;
-
}
-
}
-
-
However, I haven't tested the output using List::Permutor. I will check it in my free-time and compare the output.
-Nithin
Hi Nithin,
Sorry about my careless mistake and thanks for correcting it. Looking forward your next reply.
B.Nataraj
The total permutation for 'VF' should not be 512. It will be 8. The script you have shown will just create multiple duplicates to produce 512 results, because of usage of unnecessary foreach loops. -
use strict;
-
use warnings;
-
my @V=('GTT','GTC','GTA','GTG');
-
my @A=('GCT','GCC','GCA','GCG');
-
my @E=('GAA','GAG');
-
my @F=('TTT','TTC');
-
my @G=('GGT','GGC','GGA','GGG');
-
my @H=('CAT','CAC');
-
my $i = 1;
-
-
foreach my $v (@V ){
-
foreach my $f (@F) {
-
print "$i $v$f\n";
-
$i++;
-
}
-
}
-
-
However, I haven't tested the output using List::Permutor. I will check it in my free-time and compare the output.
-Nithin
8? Then you are not talking about permutations. VF will produce 512 unique permutations with no duplicates assuming all 6 codons can be in all six positions of each permutation.
Using VF as an example, I think your requirements are:
V (first position)
F (second position)
So all the codons in @V can only ever occupy the first position of the permutations you desire. All codons in @F can only occupy the second position. There are only two positions because that is the number of letters/arrays. Using VFas the example, there can not be a permutation with TTT or TTC in the first position, for example:
TTTGTT
TTCGTT
Is that correct?
8? Then you are not talking about permutations. VF will produce 512 unique permutations with no duplicates assuming all 6 codons can be in all six positions of each permutation.
Using VF as an example, I think your requirements are:
V (first position)
F (second position)
So all the codons in @V can only ever occupy the first position of the permutations you desire. All codons in @F can only occupy the second position. There are only two positions because that is the number of letters/arrays. Using VFas the example, there can not be a permutation with TTT or TTC in the first position, for example:
TTTGTT
TTCGTT
Is that correct?
Hi Kevin,
Exactly... yea the position should not be interchanged.
Thanks,
B.Nataraj
Hi Kevin,
Exactly... yea the position should not be interchanged.
Thanks,
B.Nataraj
List::Permutor is not the module of choice then. It will produce all possible combinations/ orders of list elements. E.g for a list with elements - 'A', 'B', 'C'.,
the possible combinations are:
A, B, C
A, C, B
B, A, C
B, C, A
C, A, B
C, B, A
But, what you are looking for is something like A can be constituted by three possible sets of codons, B by two etc.
I believe there are Bio-Perl modules which do the job that you are looking for.
List::Permutor is not the module of choice then. It will produce all possible combinations/ orders of list elements. E.g for a list with elements - 'A', 'B', 'C'.,
the possible combinations are:
A, B, C
A, C, B
B, A, C
B, C, A
C, A, B
C, B, A
But, what you are looking for is something like A can be constituted by three possible sets of codons, B by two etc.
I believe there are Bio-Perl modules which do the job that you are looking for.
Hi Nithin,
Before coming to this forum, I tried possible code for the same in bioperl module and all possible code depository for any snippet to do my objective but my search ended in vein. Technically the approach is called "Back-translation"(bioinformatics term). More the Bioinformatics problem that I try to address here is so unique and rare and it may not be a regular requirement for many , that’s why I suppose that this is not there to find in bioperl. that’s how I started to develop my own code and sought the help from this forum. Anyhow I would try again there then.
Thanks ,
B.Nataraj
Hi Nithin,
Before coming to this forum, I tried possible code for the same in bioperl module and all possible code depository for any snippet to do my objective but my search ended in vein. Technically the approach is called "Back-translation"(bioinformatics term). More the Bioinformatics problem that I try to address here is so unique and rare and it may not be a regular requirement for many , that’s why I suppose that this is not there to find in bioperl. that’s how I started to develop my own code and sought the help from this forum. Anyhow I would try again there then.
Thanks ,
B.Nataraj
Here is one way of doing it. This code dynamically generates required foreach loops and executes the command. -
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @loops = split '', $seq;
-
my $cmd;
-
local $i=0;
-
my $vars;
-
-
foreach my $letter (@loops) {
-
my ($str,$v)= createloop($letter);
-
$vars.=$v;
-
$cmd.=$str;
-
}
-
$cmd=$cmd.'print'. " \"$vars\\n\";"; ## adding print statement
-
$cmd.= '}'x$i; ### ending statement blocks
-
-
eval {eval $cmd; }; ## first eval returns statement block
-
-
##creating foreach loops
-
sub createloop {
-
$i++;
-
my $var ='$a'.$i;
-
return ("foreach $var (qw(@{$hash{$_[0]}})) {",$var) ;
-
-
}
-
Here is one way of doing it. This code dynamically generates required foreach loops and executes the command. -
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @loops = split '', $seq;
-
my $cmd;
-
local $i=0;
-
my $vars;
-
-
foreach my $letter (@loops) {
-
my ($str,$v)= createloop($letter);
-
$vars.=$v;
-
$cmd.=$str;
-
}
-
$cmd=$cmd.'print'. " \"$vars\\n\";"; ## adding print statement
-
$cmd.= '}'x$i; ### ending statement blocks
-
-
eval {eval $cmd; }; ## first eval returns statement block
-
-
##creating foreach loops
-
sub createloop {
-
$i++;
-
my $var ='$a'.$i;
-
return ("foreach $var (qw(@{$hash{$_[0]}})) {",$var) ;
-
-
}
-
Hi Nithin,
Thank you very much for the code, it does great and meet my objective.
I really thankful for you and Kevin for extending your kind help to me.
You guys are doing really a fantastic job.
Thanks once again,
B.Nataraj
nithinpes has posted a very interesting solution, but I urge caution using it since eval will run any code it evaluates. Just make sure you don't accidently introduce something unwise into the eval block. Its fine as its written, but read up on eval before making any changes that might have serious side effects.
Heres a solution that does not use eval, but I would benchmark my code and nithinpes code when you get a chance and see if one is better than the other. -
use warnings;
-
use strict;
-
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @letters = split '', $seq;
-
my $loops = 1;
-
foreach my $L (@letters) {
-
$loops *= @{$hash{$L}};
-
}
-
my @perms;
-
foreach my $S (@letters) {
-
my @t = @{$hash{$S}};
-
for my $i (0..$loops-1) {
-
my $c = shift @t;
-
$perms[$i] .= $c;
-
push @t, $c;
-
}
-
}
-
my $i = 1;
-
foreach my $s (@perms) {
-
print $i++, " $s\n";
-
}
-
I'm not sure how efficient my use of "shift" and "push" is to loop through the arrays over and over the number of times necessary.
Heres a solution that does not use eval, but I would benchmark my code and nithinpes code when you get a chance and see if one is better than the other. -
use warnings;
-
use strict;
-
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @letters = split '', $seq;
-
my $loops = 1;
-
foreach my $L (@letters) {
-
$loops *= @{$hash{$L}};
-
}
-
my @perms;
-
foreach my $S (@letters) {
-
my @t = @{$hash{$S}};
-
for my $i (0..$loops-1) {
-
my $c = shift @t;
-
$perms[$i] .= $c;
-
push @t, $c;
-
}
-
}
-
my $i = 1;
-
foreach my $s (@perms) {
-
print $i++, " $s\n";
-
}
-
I'm not sure how efficient my use of "shift" and "push" is to loop through the arrays over and over the number of times necessary.
Hi Kevin,
Thanks for your new code and pointing out possible pitfall. In fact I am half way down my real work and another half way to pass on. I am making the remaining code and hopefully that wont be so difficult logic to make it out and all about string manipulation and once entire code is ready and I will be in a position to check yours and nithin for a real time problem and I hope can give you my feedback very soon.
B.Nataraj
Hi Kevin,
The last code you have posted, is in fact repeating the first four sequence as output but it make up to the desired 512 (for a input seqnence of VAEFGH) sequence as per rule. But the alternative code by Nithin is working just fine for creating all possible combination and there is no repetition.
I just happened to compare this bug in your code at the last stage of implementing in my entire code using your code as a starting point. I could not fix in it and possibly you can do it quickly so I am posting it to you once again.
Thanks,
B.Nataraj
(I tried to delete this message, but I could not do it)
Hi Kevin,
Sorry I did not mention my input sequence in my previous message, here it is "VAEFGH".
Thanks,
B.Nataraj
(I tried to delete this message, but I could not do it)
Hi Kevin,
Sorry I did not mention my input sequence in my previous message, here it is "VAEFGH".
Thanks,
B.Nataraj
hmmm...guess I made a logical error in the code I wrote. I'll look at it later today and see if I can find the error.
well, I can't get it to work. I think you are will need to use nithinpes code or wait and see if someone else can help or post on www.perlmonks.com and see if one of the members there has a suggestion.
Hi Kevin,
Ok, I will use nithin's code then.
Thanks,
B.Nataraj
Here is one way of doing it. This code dynamically generates required foreach loops and executes the command. -
my $seq = 'VF';
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
-
my @loops = split '', $seq;
-
my $cmd;
-
local $i=0;
-
my $vars;
-
-
foreach my $letter (@loops) {
-
my ($str,$v)= createloop($letter);
-
$vars.=$v;
-
$cmd.=$str;
-
}
-
$cmd=$cmd.'print'. " \"$vars\\n\";"; ## adding print statement
-
$cmd.= '}'x$i; ### ending statement blocks
-
-
eval {eval $cmd; }; ## first eval returns statement block
-
-
##creating foreach loops
-
sub createloop {
-
$i++;
-
my $var ='$a'.$i;
-
return ("foreach $var (qw(@{$hash{$_[0]}})) {",$var) ;
-
-
}
-
Hi Nithinpes,
I am using your code posted by you in this thread for my work. Now I had a little problem in using your code to direct the output to a file. I could not do it since I could not understand the code properly, especially the eval statement and its contribution. Kindly help me to add an extra code to direct the generated sequence to save in a text file.
Thanks,
B.Nataraj
Hi Nithinpes,
I am using your code posted by you in this thread for my work. Now I had a little problem in using your code to direct the output to a file. I could not do it since I could not understand the code properly, especially the eval statement and its contribution. Kindly help me to add an extra code to direct the generated sequence to save in a text file.
Thanks,
B.Nataraj
Hi Nithinpes,
Yes, I done it myself by following way. Please igonre my post. - print "\n\n\t\#################### AA 2 PROTEIN #################### \n\n";
-
print "This script will convert your amino acid sequence to DNA Sequence\n\n";
-
print "ENTER THE FILENAME OF THE Amino acid SEQUENCE:= ";
-
$aafilename = <STDIN>;
-
chomp $aafilename;
-
-
unless ( open(aaFILE, $aafilename) ) {
-
print "Cannot open file \"$aafilename\"\n\n";
-
}
-
@aa = <aaFILE>;
-
close aaFILE;
-
$aa = join( '', @aa);
-
-
my $seq = "$aa";
-
my %hash = (
-
L => ['CTT','CTC','CTA','CTG','TTA','TTG'],
-
S => ['TCT','TCC','TCA','TCG','AGT','AGC'],
-
R => ['CGT','CGC','CGA','CGG','AGA','AGG'],
-
V => ['GTT','GTC','GTA','GTG'],
-
A => ['GCT','GCC','GCA','GCG'],
-
G => ['GGT','GGC','GGA','GGG'],
-
P => ['CCT','CCC','CCA','CCG'],
-
T => ['ACT','ACC','ACA','ACG'],
-
I => ['ATT','ATC','ATA'],
-
F => ['TTT','TTC'],
-
C => ['TGT','TGC'],
-
Y => ['TAT','TAC'],
-
Q => ['CAA','CAG'],
-
N => ['AAT','AAC'],
-
H => ['CAT','CAC'],
-
E => ['GAA','GAG'],
-
D => ['GAT','GAC'],
-
K => ['AAA','AAG'],
-
M => ['ATG'],
-
W => ['TGG'],
-
);
-
open(OUT, ">>dna_outfile.txt");
-
-
my @loops = split '', $seq;
-
my $cmd;
-
local $i=0;
-
my $vars;
-
foreach my $letter (@loops) {
-
my ($str,$v)= createloop($letter);
-
$vars.=$v;
-
$cmd.=$str;
-
}
-
$cmd=$cmd.'print OUT'. " \"$vars\\n\";"; ## adding print statement
-
#$cmd=$cmd.'print'. " \"$vars\\n\";"; ## adding print statement
-
$cmd.= '}'x$i; ; ### ending statement blocks
-
-
eval {eval $cmd; }; ## first eval returns statement block
-
-
##creating foreach loops
-
sub createloop {
-
$i++;
-
my $var ='$a'.$i;
-
return ("foreach $var (qw(@{$hash{$_[0]}})) {",$var) ;
-
-
}
-
-
-
-
close OUT;
Sign in to post your reply or Sign up for a free account.
Similar topics
by: Bart Nessux |
last post by:
How is this done in php? I've tried several variations of the following:
contact_array =
($_POST,
$_POST,
$_POST,
$_POST,
$_POST);
The ultimate goal is to insert the array into a table in...
|
by: Gary |
last post by:
I haver a table of students - Say 100 students that I need to be able to
update/delete and amend. I know I can do this one student at a time which is
simple but lets say I want to see all the...
|
by: Sjoerd |
last post by:
Summary: Use foreach(..) instead of while(list(..)=each(..)).
--==--
Foreach is a language construct, meant for looping through arrays.
There are two syntaxes; the second is a minor but useful...
|
by: khalidbaloch |
last post by:
hi every one, how are you folf , hope fine
dear Friends i want to get values of multi-dimensional arrays using foreach loop and after that print out the html using an other while loop ,
i tried...
|
by: Akira |
last post by:
I noticed that using foreach is much slower than using for-loop, so I
want to change our current code from foreach to for-loop.
But I can't figure out how.
Could someone help me please?
Current...
|
by: cheesecaker |
last post by:
Okay well, let's say you have a script that evaluates exam questions. Arrays are used for both the answer key and the test taker's provided answers. Repeated comparison is required, so obviously, a...
|
by: nirmal1349 |
last post by:
Hi,
I have a hash in hash file in perl which looks some thing like this:
$FIELDS = {
'abc' => {
'Description' => {
'Purpose' => 'some data is present',...
|
by: fjm |
last post by:
For some reason, I have always had a hard time understanding arrays as they pertain to php and databases. I understand associative arrays just fine but when there are multidimensional arrays, I kinda...
|
by: Dheeraj Joshi |
last post by:
Hi.. I wrote following code.
iindex = 0;
foreach (XmlNode Labelnodes1 in Labelnode)
{
countlab = Labelnodes1.Attributes.Count;
for (int i = 0; i < countlab; i++)
{
|
by: emmanuelkatto |
last post by:
Hi All, I am Emmanuel katto from Uganda. I want to ask what challenges you've faced while migrating a website to cloud.
Please let me know.
Thanks!
Emmanuel
|
by: BarryA |
last post by:
What are the essential steps and strategies outlined in the Data Structures and Algorithms (DSA) roadmap for aspiring data scientists? How can individuals effectively utilize this roadmap to progress...
|
by: Sonnysonu |
last post by:
This is the data of csv file
1 2 3
1 2 3
1 2 3
1 2 3
2 3
2 3
3
the lengths should be different i have to store the data by column-wise with in the specific length.
suppose the i have to...
|
by: Hystou |
last post by:
There are some requirements for setting up RAID:
1. The motherboard and BIOS support RAID configuration.
2. The motherboard has 2 or more available SATA protocol SSD/HDD slots (including MSATA, M.2...
|
by: Hystou |
last post by:
Most computers default to English, but sometimes we require a different language, especially when relocating. Forgot to request a specific language before your computer shipped? No problem! You can...
|
by: Oralloy |
last post by:
Hello folks,
I am unable to find appropriate documentation on the type promotion of bit-fields when using the generalised comparison operator "<=>".
The problem is that using the GNU compilers,...
|
by: tracyyun |
last post by:
Dear forum friends,
With the development of smart home technology, a variety of wireless communication protocols have appeared on the market, such as Zigbee, Z-Wave, Wi-Fi, Bluetooth, etc. Each...
|
by: agi2029 |
last post by:
Let's talk about the concept of autonomous AI software engineers and no-code agents. These AIs are designed to manage the entire lifecycle of a software development project—planning, coding, testing,...
|
by: isladogs |
last post by:
The next Access Europe User Group meeting will be on Wednesday 1 May 2024 starting at 18:00 UK time (6PM UTC+1) and finishing by 19:30 (7.30PM).
In this session, we are pleased to welcome a new...
| |