X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fschemes%2FHMMERColourScheme.java;h=09092a6e7726225335126802827d1c3ffa8594c9;hb=d47d1eb4b1e5e7a58a220b174bf03b814ecec295;hp=a7fc8b85460df2d4e2617b93c634558e558346bd;hpb=fdda06e4e3b668fb08532ba9192c199bb57fb5bd;p=jalview.git diff --git a/src/jalview/schemes/HMMERColourScheme.java b/src/jalview/schemes/HMMERColourScheme.java index a7fc8b8..09092a6 100644 --- a/src/jalview/schemes/HMMERColourScheme.java +++ b/src/jalview/schemes/HMMERColourScheme.java @@ -51,19 +51,46 @@ public class HMMERColourScheme extends ResidueColourScheme } Double prob; prob = hmm.getMatchEmissionProbability(position, symbol); - double freq = ResidueProperties.aminoBackgroundFrequencies.get(symbol); - Double value = prob - freq; - + double freq = 0; + if ("amino".equals(hmm.getAlphabetType())) + { + if (!ResidueProperties.aminoBackgroundFrequencies.containsKey(symbol)) + { + return Color.white; + } + freq = ResidueProperties.aminoBackgroundFrequencies.get(symbol); + } + else if ("DNA".equals(hmm.getAlphabetType())) + { + if (!ResidueProperties.dnaBackgroundFrequencies.containsKey(symbol)) + { + return Color.white; + } + freq = ResidueProperties.dnaBackgroundFrequencies.get(symbol); + } + else if ("RNA".equals(hmm.getAlphabetType())) + { + if (!ResidueProperties.rnaBackgroundFrequencies.containsKey(symbol)) + { + return Color.white; + } + freq = ResidueProperties.rnaBackgroundFrequencies.get(symbol); + } + if (prob == 0) + { + return Color.red; + } + 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;