[Biojava-l] How to read in an aligment and then start viterbi

Fabian Schreiber fab.schreiber at gmx.de
Wed May 18 12:17:12 EDT 2005


Hello!
I created a Hmm with an Alphabet, that
contains 9 different Proteins such as ALA, ARG, etc.
After i created this model, i want to read in an
aligment from a file and start Viterbi with the aligment as input.
When i do so, i always get the following error:
[code]
org.biojava.bio.symbol.IllegalSymbolException: Symbol ARG not found in 
alphabet DnaAlphabet
	at 
org.biojava.bio.symbol.AbstractAlphabet.validate(AbstractAlphabet.java:278)
	at 
org.biojava.bio.symbol.LinearAlphabetIndex.indexForSymbol(LinearAlphabetIndex.java:117)
	at 
org.biojava.bio.dist.SimpleDistribution.getWeightImpl(SimpleDistribution.java:131)
	at 
org.biojava.bio.dist.AbstractDistribution.getWeight(AbstractDistribution.java:197)
	at 
org.biojava.bio.dp.ScoreType$Probability.calculateScore(ScoreType.java:54)
	at org.biojava.bio.dp.onehead.SingleDP.getEmission(SingleDP.java:100)
	at org.biojava.bio.dp.onehead.SingleDP.viterbi(SingleDP.java:553)
	at org.biojava.bio.dp.onehead.SingleDP.viterbi(SingleDP.java:488)
[/code]

It looks as if the alphabets (one from the hmm and the other from the 
alignment) differ, but i debugged it carefully and also checked
the indices of both alphabets, which are equal, but could not find
a problem.
Is it a bug?
Did anyone of you experienced a similiar problem?

Can someone help me please?

Its really important.

Thanks a lot

Greetings

Fabian



More information about the Biojava-l mailing list