{
double iep = hmm.getInsertEmissionProbability(2, 'A');
assertEquals(iep, Math.pow(Math.E, -2.68618), 0.0001d);
+ // symbol is not case-sensitive
+ assertEquals(iep, hmm.getInsertEmissionProbability(2, 'a'));
iep = hmm.getInsertEmissionProbability(5, 'T');
assertEquals(iep, Math.pow(Math.E, -2.77519), 0.0001d);
}
@Test(groups = "Functional")
- public void testGetConsensusAtAlignColumn()
- {
- assertEquals(hmm.getConsensusAtAlignColumn(10), 's');
- assertEquals(hmm.getConsensusAtAlignColumn(50), 'k');
- hmm.setConsensusResidueStatus(false);
- assertEquals(hmm.getConsensusAtAlignColumn(100), 'l');
- assertEquals(hmm.getConsensusAtAlignColumn(400), 'k');
- }
-
- @Test(groups = "Functional")
public void testGetConsensusSequence()
{
SequenceI seq = hmm.getConsensusSequence();
String subStr = seq.getSequenceAsString().substring(0, 10);
- assertEquals(subStr, "YELLEKLGSG");
+ assertEquals(subStr, "yelleklGsG");
subStr = seq.getSequenceAsString().substring(150, 161);
- assertEquals(subStr, "-DLLK------");
+ assertEquals(subStr, "-dllk------");
}
/**
.get("amino");
int col = 4;
float expected = 0f;
- for (char aa : hmm.getSymbols())
+ for (char aa : hmm.getSymbols().toCharArray())
{
double mep = hmm.getMatchEmissionProbability(col, aa);
float background = uniprotFreqs.get(aa);