Re: Question about NiedlemanWunsch class

Jose Manuel Duarte <[email protected]> Wed, 05 Nov 2014 12:10:13 +0100
Newsgroups gmane.comp.java.bio.general
Message-ID <[email protected]>
Hi Alexander

Maybe it would be helpful if you can post your code and also even the 
sequences you are looking at, just copy/paste them into the email text. 
Then it would be possible to have a look and see whether it is a bug or 
a feature...

By the way, what version of Biojava are you using? is it 3.1? or are you 
directly getting the latest snapshots from github?

Cheers

Jose



On 04/11/14 23:30, Andreas Prlic wrote:
> Hi Alexander,
>
> Not sure, perhaps the javadoc is misleading (it comes from a high 
> level interface and the lower level implementations might do something 
> slightly different). Does the alignment image look correct? You might 
> want to test this with two shorter and perhaps quite similar 
> sequences. Also, take a look at the SmithWaterman approach, which does 
> local alignments, rather than global. Perhaps your two sequences are 
> too far apart and the global alignment does not yield good results...
>
> Andreas
>
>
>
> On Mon Nov 03 2014 at 11:01:59 AM Александр Бескровный 
> <[email protected] 
> <mailto:[email protected]>> wrote:
>
>     Hello!
>     I have a problem with aligning the sequences with the help of
>     NiedlemanWunsch class. The problem is following: I found the
>     sequences, where method "getSimilarity" got negative numbers, but
>     the description said: "It returns score as a similarity between
>     0.0 and 1.0"
>     I looked at the output of methods "getMaxScore", "getMinScore" and
>     "getScore", which are involved in calculations of method
>     "getSimilarity", and I have noticed, that "getMinScore" and
>     "getMaxScore" got strange numbers (maxScore was less than minScore
>     or both were negative), but output of "getScore" seemed adequate.
>     I used SimpleGapPenalty and SubstitutionMatrixHelper.getNuc4_4()
>     as a parameters of constructor of NiedlemanWunsch class.
>     According to rules, i can’t attach file with sequences to my
>     letter. But in case it is possible, I would kindly ask you to
>     explain me, how to do it correctly,and then i'll send you the
>     sequences. They are too long to be send like a text in a letter.
>
>     Thanks for your attention!
>
>     -- 
>     Regards,
>     Alexander Beskrovnyy
>     _______________________________________________
>     Biojava-l mailing list  - [email protected]
>     <mailto:[email protected]>
>     http://mailman.open-bio.org/mailman/listinfo/biojava-l
>
>
>
> _______________________________________________
> Biojava-l mailing list  -  [email protected]
> http://mailman.open-bio.org/mailman/listinfo/biojava-l

_______________________________________________
Biojava-l mailing list  -  [email protected]
http://mailman.open-bio.org/mailman/listinfo/biojava-l