X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fanalysis%2Fscoremodels%2FFeatureDistanceModel.java;h=c20d4f1a76353ee54adcb41221ce5175c607057a;hb=3d0101179759ef157b088ea135423cd909512d9f;hp=04a7b143330693a2e922e89bed83105fd2d9da57;hpb=fcb39fa3bc47777bf4e0eb209f765dd254dc3cb9;p=jalview.git

diff --git a/src/jalview/analysis/scoremodels/FeatureDistanceModel.java b/src/jalview/analysis/scoremodels/FeatureDistanceModel.java
index 04a7b14..c20d4f1 100644
--- a/src/jalview/analysis/scoremodels/FeatureDistanceModel.java
+++ b/src/jalview/analysis/scoremodels/FeatureDistanceModel.java
@@ -20,11 +20,15 @@
  */
 package jalview.analysis.scoremodels;
 
-import jalview.api.analysis.DistanceModelI;
-import jalview.api.analysis.ViewBasedAnalysisI;
+import jalview.api.AlignmentViewPanel;
+import jalview.api.FeatureRenderer;
+import jalview.api.analysis.ScoreModelI;
+import jalview.api.analysis.SimilarityParamsI;
 import jalview.datamodel.AlignmentView;
 import jalview.datamodel.SeqCigar;
 import jalview.datamodel.SequenceFeature;
+import jalview.math.Matrix;
+import jalview.math.MatrixI;
 import jalview.util.SetUtils;
 
 import java.util.HashMap;
@@ -33,13 +37,40 @@ import java.util.List;
 import java.util.Map;
 import java.util.Set;
 
-public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
+public class FeatureDistanceModel extends DistanceScoreModel
 {
-  jalview.api.FeatureRenderer fr;
+  private static final String NAME = "Sequence Feature Similarity";
+
+  private String description;
+
+  FeatureRenderer fr;
+
+  /**
+   * Constructor
+   */
+  public FeatureDistanceModel()
+  {
+  }
 
   @Override
-  public boolean configureFromAlignmentView(
-          jalview.api.AlignmentViewPanel view)
+  public ScoreModelI getInstance(AlignmentViewPanel view)
+  {
+    FeatureDistanceModel instance;
+    try
+    {
+      instance = this.getClass().newInstance();
+      instance.configureFromAlignmentView(view);
+      return instance;
+    } catch (InstantiationException | IllegalAccessException e)
+    {
+      System.err.println("Error in " + getClass().getName()
+              + ".getInstance(): " + e.getMessage());
+      return null;
+    }
+  }
+
+  boolean configureFromAlignmentView(AlignmentViewPanel view)
+
   {
     fr = view.cloneFeatureRenderer();
     return true;
@@ -51,18 +82,30 @@ public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
    * features each sequence pair has at each column, ignore feature types they
    * have in common, and count the rest. The totals are normalised by the number
    * of columns processed.
+   * <p>
+   * The parameters argument provides settings for treatment of gap-residue
+   * aligned positions, and whether the score is over the longer or shorter of
+   * each pair of sequences
+   * 
+   * @param seqData
+   * @param params
    */
   @Override
-  public float[][] findDistances(AlignmentView seqData)
+  public MatrixI findDistances(AlignmentView seqData,
+          SimilarityParamsI params)
   {
-    List<String> dft = fr.getDisplayedFeatureTypes();
     SeqCigar[] seqs = seqData.getSequences();
     int noseqs = seqs.length;
     int cpwidth = 0;// = seqData.getWidth();
-    float[][] distance = new float[noseqs][noseqs];
-    if (dft.isEmpty())
+    double[][] distances = new double[noseqs][noseqs];
+    List<String> dft = null;
+    if (fr != null)
+    {
+      dft = fr.getDisplayedFeatureTypes();
+    }
+    if (dft == null || dft.isEmpty())
     {
-      return distance;
+      return new Matrix(distances);
     }
 
     // need to get real position for view position
@@ -79,10 +122,10 @@ public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
         cpwidth++;
 
         /*
-         * first pass: record features types in column for each sequence
+         * first record feature types in this column for each sequence
          */
-        Map<SeqCigar, Set<String>> sfap = findFeatureTypesAtColumn(
-                seqs, cpos);
+        Map<SeqCigar, Set<String>> sfap = findFeatureTypesAtColumn(seqs,
+                cpos);
 
         /*
          * count feature types on either i'th or j'th sequence but not both
@@ -92,9 +135,23 @@ public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
         {
           for (int j = i + 1; j < noseqs; j++)
           {
-            int seqDistance = SetUtils.countDisjunction(sfap.get(seqs[i]),
-                    sfap.get(seqs[j]));
-            distance[i][j] += seqDistance;
+            SeqCigar sc1 = seqs[i];
+            SeqCigar sc2 = seqs[j];
+            Set<String> set1 = sfap.get(sc1);
+            Set<String> set2 = sfap.get(sc2);
+            boolean gap1 = set1 == null;
+            boolean gap2 = set2 == null;
+
+            /*
+             * gap-gap always scores zero
+             * residue-residue is always scored
+             * include gap-residue score if params say to do so
+             */
+            if ((!gap1 && !gap2) || params.includeGaps())
+            {
+              int seqDistance = SetUtils.countDisjunction(set1, set2);
+              distances[i][j] += seqDistance;
+            }
           }
         }
       }
@@ -103,21 +160,24 @@ public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
     /*
      * normalise the distance scores (summed over columns) by the
      * number of visible columns used in the calculation
+     * and fill in the bottom half of the matrix
      */
+    // TODO JAL-2424 cpwidth may be out by 1 - affects scores but not tree shape
     for (int i = 0; i < noseqs; i++)
     {
       for (int j = i + 1; j < noseqs; j++)
       {
-        distance[i][j] /= cpwidth;
-        distance[j][i] = distance[i][j];
+        distances[i][j] /= cpwidth;
+        distances[j][i] = distances[i][j];
       }
     }
-    return distance;
+    return new Matrix(distances);
   }
 
   /**
-   * Builds and returns a list (one per SeqCigar) of visible feature types at
-   * the given column position
+   * Builds and returns a map containing a (possibly empty) list (one per
+   * SeqCigar) of visible feature types at the given column position. The map
+   * has no entry for sequences which are gapped at the column position.
    * 
    * @param seqs
    * @param columnPosition
@@ -129,18 +189,18 @@ public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
     Map<SeqCigar, Set<String>> sfap = new HashMap<SeqCigar, Set<String>>();
     for (SeqCigar seq : seqs)
     {
-      Set<String> types = new HashSet<String>();
       int spos = seq.findPosition(columnPosition);
       if (spos != -1)
       {
+        Set<String> types = new HashSet<String>();
         List<SequenceFeature> sfs = fr.findFeaturesAtRes(seq.getRefSeq(),
                 spos);
         for (SequenceFeature sf : sfs)
         {
           types.add(sf.getType());
         }
+        sfap.put(seq, types);
       }
-      sfap.put(seq, types);
     }
     return sfap;
   }
@@ -148,7 +208,13 @@ public class FeatureDistanceModel implements DistanceModelI, ViewBasedAnalysisI
   @Override
   public String getName()
   {
-    return "Sequence Feature Similarity";
+    return NAME;
+  }
+
+  @Override
+  public String getDescription()
+  {
+    return description;
   }
 
   @Override