Rodrigues' formula: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>ClueBot NG
m Reverting possible vandalism by 198.161.203.93 to version by LaguerreLegendre. False positive? Report it. Thanks, ClueBot NG. (1638279) (Bot)
en>Rjwilmsi
m →‎References: Journal cites, Added 1 doi to a journal cite using AWB (10365)
 
Line 1: Line 1:
{{Technical|date=September 2013}}
Today I am sharing a dark experience with you as I am from this darkness today. Whoa! Google Chrome has crashed was the disaster which struck on me early this morning. My blood flow seemed to stop for a limited minutes. Next I regained my senses and decided to fix Chrome crash.<br><br>If you registry gets cluttered up with a lot of junk you don't use, your PC usually run slower. Therefore it's prudent which you frequently receive your registry cleaned.<br><br>It doesn't matter whether you may be not extremely obvious regarding what rundll32.exe is. But remember which it plays an important character in preserving the stability of the computers and the integrity of the system. Whenever certain software or hardware may not respond normally to the system surgery, comes the rundll32 exe error, which might be caused by corrupted files or lost information inside registry. Usually, error content might shows up at booting or the beginning of running a program.<br><br>Windows mistakes is caused by any number of reasons, yet there's virtually constantly one cause. There's a hidden part of the program that is responsible for making 90% of all Windows mistakes, and it's called the 'registry'. This really is the central database for your system plus is where the computer shops all its system files plus settings. It's a fairly important piece of Windows, that is should be able to function. However, it's also among the largest causes of issues on your PC.<br><br>So to fix this, we simply have to be able to make all the registry files non-corrupted again. This might dramatically speed up the loading time of the computer plus allows we to a large number of details on it again. And fixing these files couldn't be simpler - we simply should employ a tool called a [http://bestregistrycleanerfix.com/tune-up-utilities tuneup utilities 2014].<br><br>Another key element whenever we compare registry cleaners is having a facility to manage the start-up jobs. This merely means to select what programs you want to begin whenever you begin a PC. If you have unnecessary programs starting when you boot up a PC this usually cause a slow running computer.<br><br>Reboot PC - Just reboot your PC to see when the error is gone. Frequently, rebooting the PC readjusts the internal settings plus software and hence fixes the problem. If it doesn't then move on to follow the instructions under.<br><br>So inside summary, whenever comparing registry cleaning, make sure the one we choose gives you the following.A backup and restore facility, quickly surgery, automatic deletion center, start-up management, an easy means of contact and a cash back guarantee.
 
The '''Needleman–Wunsch algorithm''' is an [[algorithm]] used in [[bioinformatics]] to [[sequence alignment|align]] [[protein]] or [[nucleotide]] sequences. It was published in 1970 by Saul B. Needleman and Christian D. Wunsch;<ref name=Needleman>{{cite journal | journal=Journal of Molecular Biology | volume=48 | issue=3 | pages=443–53 | year=1970  | author=Needleman, Saul B.; and Wunsch, Christian D. | title=A general method applicable to the search for similarities in the amino acid sequence of two proteins | url=http://linkinghub.elsevier.com/retrieve/pii/0022-2836(70)90057-4 | pmid=5420325 | doi = 10.1016/0022-2836(70)90057-4 }}</ref> it uses [[dynamic programming]], and was the first application of dynamic programming to biological sequence comparison. It is sometimes referred to as the [[optimal matching]] algorithm.
 
==A modern presentation==
Scores for aligned characters are specified by a [[similarity matrix]]. Here, <math>S(a, b)</math> is the similarity of characters ''a'' and ''b''. It uses a linear [[gap penalty]], here called <math>d</math>.
 
For example, if the similarity matrix was
{| style="font-family: monospace; font-size: 150%;" cellpadding="2"
! scope="col" |
! scope="col" | A
! scope="col" | G
! scope="col" | C
! scope="col" | T
|- style="text-align: right;"
! scope="row" | A
| 10 || -1 || -3 || -4
|- style="text-align: right;"
! scope="row" | G
| -1 || 7 || -5 || -3
|- style="text-align: right;"
! scope="row" | C
| -3 || -5 || 9 || 0
|- style="text-align: right;"
! scope="row" | T
| -4 || -3 || 0 || 8
|}
 
then the alignment:
AGACTAGTTAC
CGA---GACGT
with a gap penalty of -5, would have the following score:
:<math>S(A,C) + S(G,G) + S(A,A) + (3\times d) + S(G,G) + S(T,A) + S(T,C) + S(A,G) + S(C,T)</math>
:<math>= -3 + 7 + 10 - (3\times 5) + 7 + -4 + 0 + -1 + 0 = 1</math>
 
To find the alignment with the highest score, a two-dimensional [[Array data structure|array]] (or [[Matrix (mathematics)|matrix]]) ''F'' is allocated. The entry in row ''i'' and column ''j'' is denoted here by
<math>F_{ij}</math>. There is one column for each character in sequence ''A'', and one row for each character in sequence ''B''. Thus, if we are aligning sequences of sizes ''n'' and ''m'', the amount of memory used is in <math>O(nm)</math>. [[Hirschberg's algorithm]] only holds a subset of the array in memory and uses <math>\Theta(\min \{n,m\})</math> space, but is otherwise similar to Needleman-Wunsch (and still requires <math>O(nm)</math> time).
 
As the algorithm progresses, the <math>F_{ij}</math> will be assigned to be the optimal score for the alignment of the first <math>i=0,\dotsc,n</math> characters in ''A'' and the first <math>j=0,\dotsc,m</math> characters in ''B''. The [[Bellman equation#Bellman's Principle of Optimality|principle of optimality]] is then applied as follows:
*Basis:
:<math>F_{0j} = d*j</math>
:<math>F_{i0} = d*i</math>
*Recursion, based on the principle of optimality:
:<math>F_{ij} = \max(F_{i-1,j-1} + S(A_{i}, B_{j}), \; F_{i,j-1} + d, \; F_{i-1,j} + d)</math>
 
The pseudo-code for the algorithm to compute the F matrix therefore looks like this:
'''for''' i=0 '''to''' '''length'''(A)
  F(i,0) ← d*i
'''for''' j=0 '''to''' '''length'''(B)
  F(0,j) ← d*j
'''for''' i=1 '''to''' '''length'''(A)
  '''for''' j=1 '''to''' '''length'''(B)
  {
    Match ← F(i-1,j-1) + S(A<sub>i</sub>, B<sub>j</sub>)
    Delete ← F(i-1, j) + d
    Insert ← F(i, j-1) + d
    F(i,j) ← '''max'''(Match, Insert, Delete)
  }
Once the ''F'' matrix is computed, the entry <math>F_{nm}</math> gives the maximum score among all possible alignments. To compute an alignment that actually gives this score, you start from the bottom right cell, and compare the value with the three possible sources (Match, Insert, and Delete above) to see which it came from. If Match, then <math>A_i</math> and <math>B_j</math> are aligned, if Delete, then <math>A_i</math> is aligned with a gap, and if Insert, then <math>B_j</math> is aligned with a gap. (In general, more than one choices may have the same value, leading to alternative optimal alignments.)
AlignmentA &larr; ""
AlignmentB &larr; ""
i &larr; '''length'''(A)
j &larr; '''length'''(B)
'''while''' (i > 0 '''or''' j > 0)
{
  '''if''' (i > 0 '''and''' j > 0 '''and''' F(i,j) == F(i-1,j-1) + S(A<sub>i</sub>, B<sub>j</sub>))
  {
    AlignmentA &larr; A<sub>i</sub> + AlignmentA
    AlignmentB &larr; B<sub>j</sub> + AlignmentB
    i &larr; i - 1
    j &larr; j - 1
  }
  '''else''' '''if''' (i > 0 '''and''' F(i,j) == F(i-1,j) + d)
  {
    AlignmentA &larr; A<sub>i</sub> + AlignmentA
    AlignmentB &larr; "-" + AlignmentB
    i &larr; i - 1
  }
  '''else''' (j > 0 '''and''' F(i,j) == F(i,j-1) + d)
  {
    AlignmentA &larr; "-" + AlignmentA
    AlignmentB &larr; B<sub>j</sub> + AlignmentB
    j &larr; j - 1
  }
}
 
==Historical notes==
Needleman and Wunsch describe their algorithm explicitly for the case when the alignment is penalized solely by the matches and mismatches, and gaps have no penalty (''d''=0). The original publication<ref name=Needleman /> from 1970 suggests the [[recursion]]
<math>F_{ij} = \max_{h<i,k<j} \{ F_{h,j-1}+S(A_{i},B_{j}), F_{i-1,k}+S(A_i,B_j) \}</math>.
 
The corresponding dynamic programming algorithm takes cubic time. The paper also points out that the recursion can accommodate arbitrary gap penalization formulas:
 
<blockquote>
A penalty factor, a number subtracted for every gap made, may be assessed as a barrier to allowing the gap. The penalty factor could be a function of the size and/or direction of the gap. [page 444]
</blockquote>
 
A better dynamic programming algorithm with quadratic running time for the same problem (no gap penalty) was first introduced<ref name=Sankoff>{{cite journal | doi=10.1073/pnas.69.1.4 | journal=Proceedings of the National Academy of Sciences of the USA | volume=69 | issue=1 | pages=4–6 | year=1972  | author=Sankoff D | title=Matching sequences under deletion/insertion constraints | url=http://www.pnas.org/content/69/1/4.abstract | pmid=4500555 | pmc=427531}}</ref> by David Sankoff in 1972.
Similar quadratic-time algorithms were discovered independently
by T. K. Vintsyuk<ref name=Vintsyuk>{{cite journal | journal=Kibernetika | volume=4 | pages=81–88 | year=1968  | author=Vintsyuk TK | title=Speech discrimination by dynamic programming}}</ref> in 1968 for speech processing
([[Dynamic time warping|"time warping"]]), and by Robert A. Wagner and [[Michael J. Fischer]]<ref name=WagnerFischer>{{cite journal | author=Wagner RA, Fischer MJ | journal = [[Journal of the ACM]] | title=The string-to-string correction problem | volume=21 | issue=1 | year=1974 | pages=168–173 | doi=10.1145/321796.321811}}</ref> in 1974 for string matching.
 
Needleman and Wunsch formulated their problem in terms of maximizing similarity. Another possibility is to minimize the [[Levenshtein distance|edit distance]] between sequences, introduced by [[Vladimir Levenshtein]]. Peter H. Sellers showed<ref name=Sellers>{{cite journal | doi=10.1137/0126070 | title=On the theory and computation of evolutionary distances | author=Sellers PH | journal = SIAM Journal on Applied Mathematics | volume = 26 | issue = 4 | pages = 787–793 | year = 1974}}</ref> in 1974 that the two problems are equivalent.
 
In modern terminology, "Needleman-Wunsch" refers to a [[Sequence alignment#Global and local alignments|global alignment]] algorithm that takes quadratic time for a linear or affine gap penalty.
 
==See also==
*[[Smith-Waterman algorithm]]
*[[BLAST]]
*[[Sequence mining]]
*[[Levenshtein distance]]
*[[Dynamic time warping]]
 
==References==
{{Reflist}}
 
==External links==
* [http://zhanglab.ccmb.med.umich.edu/NW-align NW-align: A protein sequence-to-sequence alignment program by Needleman-Wunsch algorithm (online server & source code)]
* [http://www.bigbold.com/snippets/posts/show/2199 Needleman-Wunsch Algorithm as Ruby Code]
* [http://chneukirchen.org/blog/archive/2006/03/dynamic-programming-in-haskell.html Needleman-Wunsch Algorithm as Haskell Code]
* [http://ds9a.nl/nwunsch A live Javascript-based demo of Needleman-Wunsch]
* [http://baba.sourceforge.net/ B.A.B.A.] &mdash; an applet (with source) which visually explains the algorithm.
* [http://www.ludwig.edu.au/course/lectures2005/Likic.pdf A clear explanation of NW and its applications to sequence alignment {{Dead link|date=July 2013}}]
* [http://technology66.blogspot.com/2008/08/sequence-alignment-techniques.html Sequence Alignment Techniques at Technology Blog]
* [http://opal.przyjaznycms.pl/en/ OPAL] JavaScript implementation of algorithms: Needleman-Wunsch, Needleman-Wunsch-Sellers and Smith-Waterman
* [http://svitsrv25.epfl.ch/R-doc/library/Biostrings/html/00Index.html Biostrings] R package implementing Needleman–Wunsch algorithm among others
* [http://www.cloudbus.org/~alchemi/files/Parallel%20Needleman%20Algo.pdf Parallel Needleman-Wunsch Algorithm for Grid] Implementation by Tahir Naveed, Imitaz Saeed Siddiqui, Shaftab Ahmed - Bahria University
* [https://gist.github.com/jonasmalacofilho/5226596 Needleman-Wunsch Algorithm as Haxe Code]
 
{{DEFAULTSORT:Needleman-Wunsch Algorithm}}
[[Category:Bioinformatics algorithms]]
[[Category:Sequence alignment algorithms]]
[[Category:Computational phylogenetics]]
[[Category:Dynamic programming]]
[[Category:Articles with example pseudocode]]

Latest revision as of 11:06, 16 August 2014

Today I am sharing a dark experience with you as I am from this darkness today. Whoa! Google Chrome has crashed was the disaster which struck on me early this morning. My blood flow seemed to stop for a limited minutes. Next I regained my senses and decided to fix Chrome crash.

If you registry gets cluttered up with a lot of junk you don't use, your PC usually run slower. Therefore it's prudent which you frequently receive your registry cleaned.

It doesn't matter whether you may be not extremely obvious regarding what rundll32.exe is. But remember which it plays an important character in preserving the stability of the computers and the integrity of the system. Whenever certain software or hardware may not respond normally to the system surgery, comes the rundll32 exe error, which might be caused by corrupted files or lost information inside registry. Usually, error content might shows up at booting or the beginning of running a program.

Windows mistakes is caused by any number of reasons, yet there's virtually constantly one cause. There's a hidden part of the program that is responsible for making 90% of all Windows mistakes, and it's called the 'registry'. This really is the central database for your system plus is where the computer shops all its system files plus settings. It's a fairly important piece of Windows, that is should be able to function. However, it's also among the largest causes of issues on your PC.

So to fix this, we simply have to be able to make all the registry files non-corrupted again. This might dramatically speed up the loading time of the computer plus allows we to a large number of details on it again. And fixing these files couldn't be simpler - we simply should employ a tool called a tuneup utilities 2014.

Another key element whenever we compare registry cleaners is having a facility to manage the start-up jobs. This merely means to select what programs you want to begin whenever you begin a PC. If you have unnecessary programs starting when you boot up a PC this usually cause a slow running computer.

Reboot PC - Just reboot your PC to see when the error is gone. Frequently, rebooting the PC readjusts the internal settings plus software and hence fixes the problem. If it doesn't then move on to follow the instructions under.

So inside summary, whenever comparing registry cleaning, make sure the one we choose gives you the following.A backup and restore facility, quickly surgery, automatic deletion center, start-up management, an easy means of contact and a cash back guarantee.