2 * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
3 * Copyright (C) $$Year-Rel$$ The Jalview Authors
5 * This file is part of Jalview.
7 * Jalview is free software: you can redistribute it and/or
8 * modify it under the terms of the GNU General Public License
9 * as published by the Free Software Foundation, either version 3
10 * of the License, or (at your option) any later version.
12 * Jalview is distributed in the hope that it will be useful, but
13 * WITHOUT ANY WARRANTY; without even the implied warranty
14 * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 * PURPOSE. See the GNU General Public License for more details.
17 * You should have received a copy of the GNU General Public License
18 * along with Jalview. If not, see <http://www.gnu.org/licenses/>.
19 * The Jalview Authors are detailed in the 'AUTHORS' file.
21 package jalview.appletgui;
23 import jalview.analysis.scoremodels.ScoreModels;
24 import jalview.analysis.scoremodels.SimilarityParams;
25 import jalview.api.analysis.ScoreModelI;
26 import jalview.datamodel.Alignment;
27 import jalview.datamodel.AlignmentView;
28 import jalview.datamodel.ColumnSelection;
29 import jalview.datamodel.SeqCigar;
30 import jalview.datamodel.SequenceI;
31 import jalview.util.MessageManager;
32 import jalview.viewmodel.PCAModel;
34 import java.awt.BorderLayout;
35 import java.awt.Button;
36 import java.awt.CheckboxMenuItem;
37 import java.awt.Choice;
38 import java.awt.Color;
39 import java.awt.FlowLayout;
40 import java.awt.Frame;
41 import java.awt.Label;
43 import java.awt.MenuBar;
44 import java.awt.MenuItem;
45 import java.awt.Panel;
46 import java.awt.event.ActionEvent;
47 import java.awt.event.ActionListener;
48 import java.awt.event.ItemEvent;
49 import java.awt.event.ItemListener;
51 public class PCAPanel extends EmbmenuFrame implements Runnable,
52 ActionListener, ItemListener
62 public PCAPanel(AlignViewport viewport)
72 for (int i = 1; i < 8; i++)
74 xCombobox.addItem("dim " + i);
75 yCombobox.addItem("dim " + i);
76 zCombobox.addItem("dim " + i);
80 boolean selected = viewport.getSelectionGroup() != null
81 && viewport.getSelectionGroup().getSize() > 0;
82 AlignmentView seqstrings = viewport.getAlignmentView(selected);
83 boolean nucleotide = viewport.getAlignment().isNucleotide();
87 seqs = viewport.getAlignment().getSequencesArray();
91 seqs = viewport.getSelectionGroup().getSequencesInOrder(viewport.getAlignment());
93 SeqCigar sq[] = seqstrings.getSequences();
94 int length = sq[0].getWidth();
96 for (int i = 0; i < seqs.length; i++)
98 if (sq[i].getWidth() != length)
101 .println("Sequences must be equal length for PCA analysis");
106 ScoreModelI scoreModel = ScoreModels.getInstance().getDefaultModel(
108 pcaModel = new PCAModel(seqstrings, seqs, nucleotide, scoreModel,
109 SimilarityParams.SeqSpace);
111 rc = new RotatableCanvas(viewport);
112 embedMenuIfNeeded(rc);
113 add(rc, BorderLayout.CENTER);
115 jalview.bin.JalviewLite.addFrame(this,
116 MessageManager.getString("label.principal_component_analysis"),
119 Thread worker = new Thread(this);
129 // TODO progress indicator
130 calcSettings.setEnabled(false);
131 rc.setEnabled(false);
134 nuclSetting.setState(pcaModel.isNucleotide());
135 protSetting.setState(!pcaModel.isNucleotide());
142 pcaModel.updateRc(rc);
144 top = pcaModel.getTop();
145 } catch (OutOfMemoryError x)
147 System.err.println("Out of memory when calculating PCA.");
150 calcSettings.setEnabled(true);
152 // TODO revert progress indicator
158 void doDimensionChange()
165 int dim1 = top - xCombobox.getSelectedIndex();
166 int dim2 = top - yCombobox.getSelectedIndex();
167 int dim3 = top - zCombobox.getSelectedIndex();
168 pcaModel.updateRcView(dim1, dim2, dim3);
170 rc.rotmat.setIdentity();
172 rc.paint(rc.getGraphics());
176 public void actionPerformed(ActionEvent evt)
178 if (evt.getSource() == inputData)
182 if (evt.getSource() == resetButton)
189 if (evt.getSource() == values)
191 values_actionPerformed();
196 public void itemStateChanged(ItemEvent evt)
198 if (evt.getSource() == xCombobox)
200 xCombobox_actionPerformed();
202 else if (evt.getSource() == yCombobox)
204 yCombobox_actionPerformed();
206 else if (evt.getSource() == zCombobox)
208 zCombobox_actionPerformed();
210 else if (evt.getSource() == labels)
212 labels_itemStateChanged(evt);
214 else if (evt.getSource() == nuclSetting)
216 if (!pcaModel.isNucleotide())
218 pcaModel.setNucleotide(true);
219 new Thread(this).start();
222 else if (evt.getSource() == protSetting)
224 if (pcaModel.isNucleotide())
226 pcaModel.setNucleotide(false);
227 new Thread(this).start();
232 protected void xCombobox_actionPerformed()
237 protected void yCombobox_actionPerformed()
242 protected void zCombobox_actionPerformed()
247 public void values_actionPerformed()
250 CutAndPasteTransfer cap = new CutAndPasteTransfer(false, null);
251 Frame frame = new Frame();
253 jalview.bin.JalviewLite.addFrame(frame,
254 MessageManager.getString("label.pca_details"), 500, 500);
256 cap.setText(pcaModel.getDetails());
259 void showOriginalData()
261 // decide if av alignment is sufficiently different to original data to
262 // warrant a new window to be created
263 // create new alignmnt window with hidden regions (unhiding hidden regions
264 // yields unaligned seqs)
265 // or create a selection box around columns in alignment view
266 // test Alignment(SeqCigar[])
270 // we try to get the associated view's gap character
271 // but this may fail if the view was closed...
272 gc = av.getGapCharacter();
273 } catch (Exception ex)
277 Object[] alAndColsel = pcaModel.getSeqtrings()
278 .getAlignmentAndColumnSelection(gc);
280 if (alAndColsel != null && alAndColsel[0] != null)
282 Alignment al = new Alignment((SequenceI[]) alAndColsel[0]);
283 AlignFrame af = new AlignFrame(al, av.applet,
284 "Original Data for PCA", false);
286 af.viewport.setHiddenColumns((ColumnSelection) alAndColsel[1]);
290 public void labels_itemStateChanged(ItemEvent itemEvent)
292 rc.showLabels(labels.getState());
295 Panel jPanel2 = new Panel();
297 Label jLabel1 = new Label();
299 Label jLabel2 = new Label();
301 Label jLabel3 = new Label();
303 protected Choice xCombobox = new Choice();
305 protected Choice yCombobox = new Choice();
307 protected Choice zCombobox = new Choice();
309 protected Button resetButton = new Button();
311 FlowLayout flowLayout1 = new FlowLayout();
313 BorderLayout borderLayout1 = new BorderLayout();
315 MenuBar menuBar1 = new MenuBar();
317 Menu menu1 = new Menu();
319 Menu menu2 = new Menu();
321 Menu calcSettings = new Menu();
323 protected CheckboxMenuItem labels = new CheckboxMenuItem();
325 protected CheckboxMenuItem protSetting = new CheckboxMenuItem();
327 protected CheckboxMenuItem nuclSetting = new CheckboxMenuItem();
329 MenuItem values = new MenuItem();
331 MenuItem inputData = new MenuItem();
333 private void jbInit() throws Exception
335 this.setLayout(borderLayout1);
336 jPanel2.setLayout(flowLayout1);
337 jLabel1.setFont(new java.awt.Font("Verdana", 0, 12));
338 jLabel1.setText("x=");
339 jLabel2.setFont(new java.awt.Font("Verdana", 0, 12));
340 jLabel2.setText("y=");
341 jLabel3.setFont(new java.awt.Font("Verdana", 0, 12));
342 jLabel3.setText("z=");
343 jPanel2.setBackground(Color.white);
344 zCombobox.setFont(new java.awt.Font("Verdana", 0, 12));
345 zCombobox.addItemListener(this);
346 yCombobox.setFont(new java.awt.Font("Verdana", 0, 12));
347 yCombobox.addItemListener(this);
348 xCombobox.setFont(new java.awt.Font("Verdana", 0, 12));
349 xCombobox.addItemListener(this);
350 resetButton.setFont(new java.awt.Font("Verdana", 0, 12));
351 resetButton.setLabel(MessageManager.getString("action.reset"));
352 resetButton.addActionListener(this);
353 this.setMenuBar(menuBar1);
354 menu1.setLabel(MessageManager.getString("action.file"));
355 menu2.setLabel(MessageManager.getString("action.view"));
356 calcSettings.setLabel(MessageManager.getString("action.change_params"));
357 labels.setLabel(MessageManager.getString("label.labels"));
358 labels.addItemListener(this);
359 values.setLabel(MessageManager.getString("label.output_values"));
360 values.addActionListener(this);
361 inputData.setLabel(MessageManager.getString("label.input_data"));
362 nuclSetting.setLabel(MessageManager
363 .getString("label.nucleotide_matrix"));
364 nuclSetting.addItemListener(this);
365 protSetting.setLabel(MessageManager.getString("label.protein_matrix"));
366 protSetting.addItemListener(this);
367 this.add(jPanel2, BorderLayout.SOUTH);
368 jPanel2.add(jLabel1, null);
369 jPanel2.add(xCombobox, null);
370 jPanel2.add(jLabel2, null);
371 jPanel2.add(yCombobox, null);
372 jPanel2.add(jLabel3, null);
373 jPanel2.add(zCombobox, null);
374 jPanel2.add(resetButton, null);
377 menuBar1.add(calcSettings);
380 menu1.add(inputData);
381 calcSettings.add(nuclSetting);
382 calcSettings.add(protSetting);
383 inputData.addActionListener(this);