X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fgui%2FPairwiseAlignPanel.java;h=2c65a49595cda0170925365493744472a7dad8d9;hb=e1e98ff6e24e5f26f07ab94cdd94799891662319;hp=248a54c4ba013bb37b28942ab6bd2b654341e975;hpb=f3bc67915e61855d87866b4a614fe7a953813408;p=jalview.git diff --git a/src/jalview/gui/PairwiseAlignPanel.java b/src/jalview/gui/PairwiseAlignPanel.java index 248a54c..2c65a49 100755 --- a/src/jalview/gui/PairwiseAlignPanel.java +++ b/src/jalview/gui/PairwiseAlignPanel.java @@ -1,160 +1,156 @@ -/* - * Jalview - A Sequence Alignment Editor and Viewer - * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle - * - * This program is free software; you can redistribute it and/or - * modify it under the terms of the GNU General Public License - * as published by the Free Software Foundation; either version 2 - * of the License, or (at your option) any later version. - * - * This program is distributed in the hope that it will be useful, - * but WITHOUT ANY WARRANTY; without even the implied warranty of - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the - * GNU General Public License for more details. - * - * You should have received a copy of the GNU General Public License - * along with this program; if not, write to the Free Software - * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA - */ -package jalview.gui; - -import jalview.analysis.*; - -import jalview.datamodel.*; - -import jalview.jbgui.*; - -import java.awt.event.*; - -import java.util.*; - - -/** - * DOCUMENT ME! - * - * @author $author$ - * @version $Revision$ - */ -public class PairwiseAlignPanel extends GPairwiseAlignPanel -{ - - AlignViewport av; - Vector sequences; - - /** - * Creates a new PairwiseAlignPanel object. - * - * @param av DOCUMENT ME! - */ - public PairwiseAlignPanel(AlignViewport av) - { - super(); - this.av = av; - - sequences = new Vector(); - - SequenceI [] seqs; - String [] seqStrings = av.getViewAsString(true); - - if(av.getSelectionGroup()==null) - { - seqs = av.alignment.getSequencesArray(); - } - else - { - seqs = av.getSelectionGroup().getSequencesInOrder(av.alignment); - } - - - float[][] scores = new float[seqs.length][seqs.length]; - double totscore = 0; - int count = seqs.length; - - Sequence seq; - - for (int i = 1; i < count; i++) - { - for (int j = 0; j < i; j++) - { - - AlignSeq as = new AlignSeq(seqs[i], seqStrings[i], - seqs[j], seqStrings[j], "pep"); - - if(as.s1str.length()==0 || as.s2str.length()==0) - { - continue; - } - - as.calcScoreMatrix(); - as.traceAlignment(); - - - as.printAlignment(System.out); - scores[i][j] = (float) as.getMaxScore() / (float) as.getASeq1().length; - totscore = totscore + scores[i][j]; - - textarea.append(as.getOutput()); - seq = new Sequence(as.getS1().getName(), - as.getAStr1(), - as.getS1().getStart(), - as.getS1().getEnd() - ); - sequences.add(seq); - - seq = new Sequence(as.getS2().getName(), - as.getAStr2(), - as.getS2().getStart(), - as.getS2().getEnd() ); - sequences.add(seq); - } - } - - 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) + " " + - seqs[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"); - } - } - - /** - * DOCUMENT ME! - * - * @param e DOCUMENT ME! - */ - protected void viewInEditorButton_actionPerformed(ActionEvent e) - { - Sequence[] seq = new Sequence[sequences.size()]; - - for (int i = 0; i < sequences.size(); i++) - { - seq[i] = (Sequence) sequences.elementAt(i); - } - - AlignFrame af = new AlignFrame(new Alignment(seq), - AlignFrame.DEFAULT_WIDTH, - AlignFrame.DEFAULT_HEIGHT); - - Desktop.addInternalFrame(af, "Pairwise Aligned Sequences", - AlignFrame.DEFAULT_WIDTH, AlignFrame.DEFAULT_HEIGHT); - } -} +/* + * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$) + * Copyright (C) $$Year-Rel$$ The Jalview Authors + * + * This file is part of Jalview. + * + * Jalview is free software: you can redistribute it and/or + * modify it under the terms of the GNU General Public License + * as published by the Free Software Foundation, either version 3 + * of the License, or (at your option) any later version. + * + * Jalview is distributed in the hope that it will be useful, but + * WITHOUT ANY WARRANTY; without even the implied warranty + * of MERCHANTABILITY or FITNESS FOR A PARTICULAR + * PURPOSE. See the GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with Jalview. If not, see . + * The Jalview Authors are detailed in the 'AUTHORS' file. + */ +package jalview.gui; + +import jalview.analysis.AlignSeq; +import jalview.datamodel.Alignment; +import jalview.datamodel.Sequence; +import jalview.datamodel.SequenceI; +import jalview.jbgui.GPairwiseAlignPanel; +import jalview.util.MessageManager; +import jalview.viewmodel.AlignmentViewport; + +import java.awt.event.ActionEvent; +import java.util.Vector; + +/** + * DOCUMENT ME! + * + * @author $author$ + * @version $Revision$ + */ +public class PairwiseAlignPanel extends GPairwiseAlignPanel +{ + + AlignmentViewport av; + + Vector sequences; + + /** + * Creates a new PairwiseAlignPanel object. + * + * @param av + * DOCUMENT ME! + */ + public PairwiseAlignPanel(AlignmentViewport av) + { + super(); + this.av = av; + + sequences = new Vector(); + + SequenceI[] seqs; + String[] seqStrings = av.getViewAsString(true); + + if (av.getSelectionGroup() == null) + { + seqs = av.getAlignment().getSequencesArray(); + } + else + { + seqs = av.getSelectionGroup().getSequencesInOrder(av.getAlignment()); + } + + String type = (av.getAlignment().isNucleotide()) ? AlignSeq.DNA + : AlignSeq.PEP; + + float[][] scores = new float[seqs.length][seqs.length]; + double totscore = 0; + int count = seqs.length; + + Sequence seq; + + for (int i = 1; i < count; i++) + { + for (int j = 0; j < i; j++) + { + + AlignSeq as = new AlignSeq(seqs[i], seqStrings[i], seqs[j], + seqStrings[j], type); + + if (as.s1str.length() == 0 || as.s2str.length() == 0) + { + continue; + } + + as.calcScoreMatrix(); + as.traceAlignment(); + + as.printAlignment(System.out); + scores[i][j] = (float) as.getMaxScore() + / (float) as.getASeq1().length; + totscore = totscore + scores[i][j]; + + textarea.append(as.getOutput()); + sequences.add(as.getAlignedSeq1()); + sequences.add(as.getAlignedSeq2()); + } + } + + 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) + " " + + seqs[i].getName()); + } + + System.out.println("\n"); + + for (int i = 0; i < count; i++) + { + for (int j = 0; j < i; j++) + { + jalview.util.Format.printDouble(System.out, "%7.3f", scores[i][j] + / totscore); + } + } + + System.out.println("\n"); + } + } + + /** + * DOCUMENT ME! + * + * @param e + * DOCUMENT ME! + */ + protected void viewInEditorButton_actionPerformed(ActionEvent e) + { + Sequence[] seq = new Sequence[sequences.size()]; + + for (int i = 0; i < sequences.size(); i++) + { + seq[i] = (Sequence) sequences.elementAt(i); + } + + AlignFrame af = new AlignFrame(new Alignment(seq), + AlignFrame.DEFAULT_WIDTH, AlignFrame.DEFAULT_HEIGHT); + + Desktop.addInternalFrame(af, + MessageManager.getString("label.pairwise_aligned_sequences"), + AlignFrame.DEFAULT_WIDTH, AlignFrame.DEFAULT_HEIGHT); + } +}