}
Double prob;
prob = hmm.getMatchEmissionProbability(position, symbol);
+ if (prob == 0)
+ {
+ return Color.red;
+ }
double freq = ResidueProperties.aminoBackgroundFrequencies.get(symbol);
- Double value = prob - freq;
-
+ Double value = Math.log(prob / freq);
Color colour = null;
- if (value >= 0)
+ if (value > 0)
{
colour = ColorUtils.getGraduatedColour(value.floatValue(), 0,
- Color.WHITE, 1f, Color.green);
+ Color.WHITE, 3f, Color.blue);
}
else if (value < 0)
{
- return Color.YELLOW;
+ return Color.ORANGE;
}
return colour;