/* * 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.appletgui; import jalview.bin.JalviewLite; import jalview.datamodel.Alignment; import jalview.datamodel.AlignmentView; import jalview.datamodel.ColumnSelection; import jalview.datamodel.SeqCigar; import jalview.datamodel.SequenceI; import jalview.util.MessageManager; import jalview.viewmodel.PCAModel; import java.awt.BorderLayout; import javax.swing.JButton; import javax.swing.JCheckBoxMenuItem; import java.awt.Color; import java.awt.FlowLayout; import javax.swing.JFrame; import javax.swing.JLabel; import javax.swing.JMenu; import javax.swing.JMenuBar; import javax.swing.JMenuItem; import javax.swing.JPanel; import java.awt.event.ActionEvent; import java.awt.event.ActionListener; import java.awt.event.ItemEvent; import java.awt.event.ItemListener; public class PCAPanel extends EmbmenuFrame implements Runnable, ActionListener, ItemListener { RotatableCanvas rc; AlignViewport av; PCAModel pcaModel; int top = 0; public PCAPanel(AlignViewport av) { try { jbInit(); } catch (Exception e) { e.printStackTrace(); } for (int i = 1; i < 8; i++) { xCombobox.addItem("dim " + i); yCombobox.addItem("dim " + i); zCombobox.addItem("dim " + i); } this.av = av; boolean selected = av.getSelectionGroup() != null && av.getSelectionGroup().getSize() > 0; AlignmentView seqstrings = av.getAlignmentView(selected); boolean nucleotide = av.getAlignment().isNucleotide(); SequenceI[] seqs; if (!selected) { seqs = av.getAlignment().getSequencesArray(); } else { seqs = av.getSelectionGroup().getSequencesInOrder(av.getAlignment()); } SeqCigar sq[] = seqstrings.getSequences(); int length = sq[0].getWidth(); for (int i = 0; i < seqs.length; i++) { if (sq[i].getWidth() != length) { System.out .println("Sequences must be equal length for PCA analysis"); return; } } pcaModel = new PCAModel(seqstrings, seqs, nucleotide); rc = new RotatableCanvas(av); embedMenuIfNeeded(rc); add(rc, BorderLayout.CENTER); JalviewLite.addFrame(this, MessageManager.getString("label.principal_component_analysis"), 475, 400); Thread worker = new Thread(this); worker.start(); } /** * DOCUMENT ME! */ public void run() { // TODO progress indicator calcSettings.setEnabled(false); rc.setEnabled(false); try { nuclSetting.setState(pcaModel.isNucleotide()); protSetting.setState(!pcaModel.isNucleotide()); pcaModel.run(); // //////////////// xCombobox.select(0); yCombobox.select(1); zCombobox.select(2); pcaModel.updateRc(rc); // rc.invalidate(); top = pcaModel.getTop(); } catch (OutOfMemoryError x) { System.err.println("Out of memory when calculating PCA."); return; } calcSettings.setEnabled(true); // TODO revert progress indicator rc.setEnabled(true); rc.repaint(); this.repaint(); } void doDimensionChange() { if (top == 0) { return; } int dim1 = top - xCombobox.getSelectedIndex(); int dim2 = top - yCombobox.getSelectedIndex(); int dim3 = top - zCombobox.getSelectedIndex(); pcaModel.updateRcView(dim1, dim2, dim3); rc.img = null; rc.rotmat.setIdentity(); rc.initAxes(); rc.paint(rc.getGraphics()); } public void actionPerformed(ActionEvent evt) { if (evt.getSource() == inputData) { showOriginalData(); } if (evt.getSource() == resetButton) { xCombobox.select(0); yCombobox.select(1); zCombobox.select(2); doDimensionChange(); } if (evt.getSource() == values) { values_actionPerformed(); } } public void itemStateChanged(ItemEvent evt) { if (evt.getSource() == xCombobox) { xCombobox_actionPerformed(); } else if (evt.getSource() == yCombobox) { yCombobox_actionPerformed(); } else if (evt.getSource() == zCombobox) { zCombobox_actionPerformed(); } else if (evt.getSource() == labels) { labels_itemStateChanged(evt); } else if (evt.getSource() == nuclSetting) { if (!pcaModel.isNucleotide()) { pcaModel.setNucleotide(true); new Thread(this).start(); } } else if (evt.getSource() == protSetting) { if (pcaModel.isNucleotide()) { pcaModel.setNucleotide(false); new Thread(this).start(); } } } protected void xCombobox_actionPerformed() { doDimensionChange(); } protected void yCombobox_actionPerformed() { doDimensionChange(); } protected void zCombobox_actionPerformed() { doDimensionChange(); } public void values_actionPerformed() { CutAndPasteTransfer cap = new CutAndPasteTransfer(false, null); JFrame frame = new JFrame(); frame.add(cap); JalviewLite.addFrame(frame, MessageManager.getString("label.pca_details"), 500, 500); cap.setText(pcaModel.getDetails()); } void showOriginalData() { // decide if av alignment is sufficiently different to original data to // warrant a new window to be created // create new alignmnt window with hidden regions (unhiding hidden regions // yields unaligned seqs) // or create a selection box around columns in alignment view // test Alignment(SeqCigar[]) char gc = '-'; try { // we try to get the associated view's gap character // but this may fail if the view was closed... gc = av.getGapCharacter(); } catch (Exception ex) { } ; Object[] alAndColsel = pcaModel.getSeqtrings() .getAlignmentAndColumnSelection(gc); if (alAndColsel != null && alAndColsel[0] != null) { Alignment al = new Alignment((SequenceI[]) alAndColsel[0]); AlignFrame af = new AlignFrame(al, av.applet, "Original Data for PCA", false); af.viewport.setHiddenColumns((ColumnSelection) alAndColsel[1]); } } public void labels_itemStateChanged(ItemEvent itemEvent) { rc.showLabels(labels.isSelected()); } JPanel jPanel2 = new JPanel(); JLabel jLabel1 = new JLabel(); JLabel jLabel2 = new JLabel(); JLabel jLabel3 = new JLabel(); protected Choice xCombobox = new Choice(); protected Choice yCombobox = new Choice(); protected Choice zCombobox = new Choice(); protected JButton resetButton = new JButton(); FlowLayout flowLayout1 = new FlowLayout(); BorderLayout borderLayout1 = new BorderLayout(); JMenuBar menuBar1 = new JMenuBar(); JMenu menu1 = new JMenu(); JMenu menu2 = new JMenu(); JMenu calcSettings = new JMenu(); protected JCheckBoxMenuItem labels = new JCheckBoxMenuItem(); protected JCheckBoxMenuItem protSetting = new JCheckBoxMenuItem(); protected JCheckBoxMenuItem nuclSetting = new JCheckBoxMenuItem(); JMenuItem values = new JMenuItem(); JMenuItem inputData = new JMenuItem(); private void jbInit() throws Exception { this.setLayout(borderLayout1); jPanel2.setLayout(flowLayout1); jLabel1.setFont(new java.awt.Font("Verdana", 0, 12)); jLabel1.setText("x="); jLabel2.setFont(new java.awt.Font("Verdana", 0, 12)); jLabel2.setText("y="); jLabel3.setFont(new java.awt.Font("Verdana", 0, 12)); jLabel3.setText("z="); jPanel2.setBackground(Color.white); zCombobox.setFont(new java.awt.Font("Verdana", 0, 12)); zCombobox.addItemListener(this); yCombobox.setFont(new java.awt.Font("Verdana", 0, 12)); yCombobox.addItemListener(this); xCombobox.setFont(new java.awt.Font("Verdana", 0, 12)); xCombobox.addItemListener(this); resetButton.setFont(new java.awt.Font("Verdana", 0, 12)); resetButton.setLabel(MessageManager.getString("action.reset")); resetButton.addActionListener(this); this.setJMenuBar(menuBar1); menu1.setLabel(MessageManager.getString("action.file")); menu2.setLabel(MessageManager.getString("action.view")); calcSettings.setLabel(MessageManager.getString("action.change_params")); labels.setLabel(MessageManager.getString("label.labels")); labels.addItemListener(this); values.setLabel(MessageManager.getString("label.output_values")); values.addActionListener(this); inputData.setLabel(MessageManager.getString("label.input_data")); nuclSetting.setLabel(MessageManager .getString("label.nucleotide_matrix")); nuclSetting.addItemListener(this); protSetting.setLabel(MessageManager.getString("label.protein_matrix")); protSetting.addItemListener(this); this.add(jPanel2, BorderLayout.SOUTH); jPanel2.add(jLabel1, null); jPanel2.add(xCombobox, null); jPanel2.add(jLabel2, null); jPanel2.add(yCombobox, null); jPanel2.add(jLabel3, null); jPanel2.add(zCombobox, null); jPanel2.add(resetButton, null); menuBar1.add(menu1); menuBar1.add(menu2); menuBar1.add(calcSettings); menu2.add(labels); menu1.add(values); menu1.add(inputData); calcSettings.add(nuclSetting); calcSettings.add(protSetting); inputData.addActionListener(this); } }