From c521a798765bb3672e89ea8cda3241a97d6c6e73 Mon Sep 17 00:00:00 2001
From: jprocter <Jim Procter>
Date: Mon, 30 May 2005 11:01:57 +0000
Subject: [PATCH] Modified output on stdout for similarity matrix.

---
 src/jalview/gui/PairwiseAlignPanel.java |   61 +++++++++++++++++--------------
 1 file changed, 34 insertions(+), 27 deletions(-)

diff --git a/src/jalview/gui/PairwiseAlignPanel.java b/src/jalview/gui/PairwiseAlignPanel.java
index 972d4f5..76f4eb5 100755
--- a/src/jalview/gui/PairwiseAlignPanel.java
+++ b/src/jalview/gui/PairwiseAlignPanel.java
@@ -17,40 +17,47 @@ public class PairwiseAlignPanel extends GPairwiseAlignPanel
     {
       super();
       this.av = av;
-      float scores[][] = new float[av.getAlignment().getHeight()][av.getAlignment().getHeight()];
+      Vector selsubset = new Vector();
+      for (int i=0, j=av.getSelectionGroup().getSize(); i<j; i++)
+        if (av.getAlignment().getSequences().contains(av.getSelectionGroup().getSequenceAt(i)))
+          selsubset.add(av.getSelectionGroup().getSequenceAt(i));
+      float scores[][] = new float[selsubset.size()][selsubset.size()];
       double totscore = 0;
-      int count = av.getSelectionGroup().getSize();
+      int count = selsubset.size();
 
       int acount = 0;
-        for (int i = 1; i < count; i++)
+      for (int i = 1; i < count; i++)
+      {
+        for (int j = 0; j < i; j++)
         {
-          for (int j = 0; j < i; j++)
-          {
-            acount++;
-            AlignSeq as = new AlignSeq(av.getSelectionGroup().getSequenceAt(i),
-                                       av.getSelectionGroup().getSequenceAt(j),"pep");
-
-            as.calcScoreMatrix();
-            as.traceAlignment();
-            as.printAlignment();
-            scores[i][j] = (float)as.getMaxScore()/(float)as.getASeq1().length;
-            totscore = totscore + scores[i][j];
-
-            textarea.append(as.getOutput());
-            sequences.add( new Sequence( as.getS1().getName(), as.getAStr1()) );
-            sequences.add( new Sequence( as.getS2().getName(), as.getAStr2()) );
-
-
-          }
+          acount++;
+          AlignSeq as = new AlignSeq( (SequenceI) selsubset.elementAt(i),
+                                     (SequenceI) selsubset.elementAt(j), "pep");
+          as.calcScoreMatrix();
+          as.traceAlignment();
+          as.printAlignment();
+          scores[i][j] = (float) as.getMaxScore() /
+              (float) as.getASeq1().length;
+          totscore = totscore + scores[i][j];
+
+          textarea.append(as.getOutput());
+          sequences.add(new Sequence(as.getS1().getName(), as.getAStr1()));
+          sequences.add(new Sequence(as.getS2().getName(), as.getAStr2()));
         }
+      }
+      if (count > 2)
+      {
+        System.out.println("Pairwise alignment scaled similarity score matrix\n");
+        for (int i = 0; i < count;i++)
+          jalview.util.Format.print(System.out, "%s \n", (""+i)+" "+((SequenceI) selsubset.elementAt(i)).getName());
+        System.out.println("\n");
+        for (int i = 0; i < count;i++)
+          for (int j = 0; j < i; j++)
+            jalview.util.Format.print(System.out,"%7.3f",scores[i][j]/totscore);
+        System.out.println("\n");
 
-        if (count > 2)
-        {
-          for (int i = 0; i < count;i++)
-            for (int j = 0; j < i; j++)
-              jalview.util.Format.print(System.out,"%7.3f",scores[i][j]/totscore);
-        }
       }
+    }
 
 
   protected void viewInEditorButton_actionPerformed(ActionEvent e)
-- 
1.7.10.2