convert code into matlab code

3 次查看（过去 30 天）
huda nawaf
huda nawaf 2011-12-27
hi, i badly need this code, but it is not matlab code. this is code foe sequence local alignment(smith-waterman), I did not find it as matlab code .I converted most of this code but the last part, I did not know the commands of this programming language can anybody help to convert the following code into matlab?
  1. trace-back
my $align1 = "";
my $align2 = "";
my $j = $max_j;
my $i = $max_i;
while (1) {
last if $matrix[$i][$j]{pointer} eq "none";
if ($matrix[$i][$j]{pointer} eq "diagonal") {
$align1 .= substr($seq1, $j-1, 1);
$align2 .= substr($seq2, $i-1, 1);
$i--; $j--;
}
elsif ($matrix[$i][$j]{pointer} eq "left") {
$align1 .= substr($seq1, $j-1, 1);
$align2 .= "-";
$j--;
}
elsif ($matrix[$i][$j]{pointer} eq "up") {
$align1 .= "-";
$align2 .= substr($seq2, $i-1, 1);
$i--;
}
}
$align1 = reverse $align1;
$align2 = reverse $align2;
print "$align1\n";
print "$align2\n";
thanks in advance

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2011-12-27
Looks like perl.
last if condition
translates to
if condition; break; end
and
A .= B
translates to
A = [A B]
and
$matrix[$i][$j]{pointer}
translates to
matrix(i,j).pointer
and
reverse $align1
translates to
fliplr(align1)
  5 个评论
显示 3更早的评论
huda nawaf
huda nawaf 2011-12-29
hi walter,
this code of perl is very tired:
I totally converted it of course with you help , but I found out something I forgot it so I have not correct result
I looked for in google about what mean 'next' with 'if' in perl , but did not understand it, I did not pick up what is the role of 'next' here:
if ($diagonal_score <= 0 and $up_score <= 0 and $left_score <= 0) {
$matrix[$i][$j]{score} = 0;
$matrix[$i][$j]{pointer} = "none";
next; # terminate this iteration of the loop
}
many thanks in advance
Walter Roberson
Walter Roberson 2011-12-29
"next" in perl is the equivalent of "continue" -- that is, skip any further instructions in the body of the loop and start the next iteration of the loop.
This contrasts with perl's "last", which is the equivalent of "break" -- that is, skip any further instructions in the body of the loop and exit the loop itself, continuing on after the end of the loop.

请先登录，再进行评论。

更多回答（1 个）

Diego
Diego 2011-12-30
If you have the bioinfo toolbox, then you probably have the swalign file that implements the smith-waterman algorithm. On the other hand, I just discovered that Matlab runs perl scripts smoothly! In my case, I have a perl script that uses the bioperl module and Matlab ran it!. I said this because I was facing a similar problem, and these tips may help you to accomplish your work. Regards, Diego
  4 个评论
显示 2更早的评论
huda nawaf
huda nawaf 2011-12-30
hi Diego,
I'm glad to hear that someone use alignment with non biological seq.
in fact, I can not use swalign and nwalign , because these codes defined just for 4 symbols (A,C,G,T) or 20 symols of amino acids. In my situation , need more symbols to represent my data, not just this reason ,but because of substitution matrices, BLOSUM AND PAM , are defined just for Nucleotides and amino acids.
In my situation , I think i have to use unitary matrix instead of BLOSUM and PAM.
currently, i try to write code of smith-waterman . Regarding Needleman I wrote it, it is work correctly.
huda nawaf
huda nawaf 2011-12-31
hi,
is necessary existence perl software to run perl script in matlab?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Data Type Conversion 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by