2 // FORESTER -- software libraries and applications
\r
3 // for evolutionary biology research and applications.
\r
5 // Copyright (C) 2008-2009 Christian M. Zmasek
\r
6 // Copyright (C) 2008-2009 Burnham Institute for Medical Research
\r
7 // Copyright (C) 2003-2007 Ethalinda K.S. Cannon
\r
8 // All rights reserved
\r
10 // This library is free software; you can redistribute it and/or
\r
11 // modify it under the terms of the GNU Lesser General Public
\r
12 // License as published by the Free Software Foundation; either
\r
13 // version 2.1 of the License, or (at your option) any later version.
\r
15 // This library is distributed in the hope that it will be useful,
\r
16 // but WITHOUT ANY WARRANTY; without even the implied warranty of
\r
17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
\r
18 // Lesser General Public License for more details.
\r
20 // You should have received a copy of the GNU Lesser General Public
\r
21 // License along with this library; if not, write to the Free Software
\r
22 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
\r
24 // Contact: phylosoft @ gmail . com
\r
25 // WWW: https://sites.google.com/site/cmzmasek/home/software/forester
\r
27 package org.forester.archaeopteryx;
\r
29 import java.awt.BorderLayout;
\r
30 import java.awt.Font;
\r
31 import java.awt.event.ActionEvent;
\r
32 import java.awt.event.ComponentAdapter;
\r
33 import java.awt.event.ComponentEvent;
\r
34 import java.awt.event.WindowAdapter;
\r
35 import java.awt.event.WindowEvent;
\r
36 import java.io.File;
\r
37 import java.io.FileInputStream;
\r
38 import java.io.IOException;
\r
39 import java.io.InputStream;
\r
40 import java.net.MalformedURLException;
\r
41 import java.net.URL;
\r
42 import java.util.ArrayList;
\r
43 import java.util.HashSet;
\r
44 import java.util.List;
\r
45 import java.util.Set;
\r
47 import javax.swing.ButtonGroup;
\r
48 import javax.swing.JCheckBoxMenuItem;
\r
49 import javax.swing.JFileChooser;
\r
50 import javax.swing.JMenu;
\r
51 import javax.swing.JMenuBar;
\r
52 import javax.swing.JMenuItem;
\r
53 import javax.swing.JOptionPane;
\r
54 import javax.swing.JRadioButtonMenuItem;
\r
55 import javax.swing.UIManager;
\r
56 import javax.swing.UnsupportedLookAndFeelException;
\r
57 import javax.swing.WindowConstants;
\r
58 import javax.swing.event.ChangeEvent;
\r
59 import javax.swing.event.ChangeListener;
\r
60 import javax.swing.plaf.synth.SynthLookAndFeel;
\r
62 import org.forester.analysis.TaxonomyDataManager;
\r
63 import org.forester.archaeopteryx.AptxUtil.GraphicsExportType;
\r
64 import org.forester.archaeopteryx.Options.CLADOGRAM_TYPE;
\r
65 import org.forester.archaeopteryx.Options.NODE_LABEL_DIRECTION;
\r
66 import org.forester.archaeopteryx.Options.PHYLOGENY_GRAPHICS_TYPE;
\r
67 import org.forester.archaeopteryx.tools.AncestralTaxonomyInferrer;
\r
68 import org.forester.archaeopteryx.tools.InferenceManager;
\r
69 import org.forester.archaeopteryx.tools.PhyloInferenceDialog;
\r
70 import org.forester.archaeopteryx.tools.PhylogeneticInferenceOptions;
\r
71 import org.forester.archaeopteryx.tools.PhylogeneticInferrer;
\r
72 import org.forester.archaeopteryx.tools.SequenceDataRetriver;
\r
73 import org.forester.archaeopteryx.webservices.PhylogeniesWebserviceClient;
\r
74 import org.forester.archaeopteryx.webservices.WebservicesManager;
\r
75 import org.forester.io.parsers.FastaParser;
\r
76 import org.forester.io.parsers.GeneralMsaParser;
\r
77 import org.forester.io.parsers.PhylogenyParser;
\r
78 import org.forester.io.parsers.nexus.NexusPhylogeniesParser;
\r
79 import org.forester.io.parsers.nhx.NHXParser;
\r
80 import org.forester.io.parsers.nhx.NHXParser.TAXONOMY_EXTRACTION;
\r
81 import org.forester.io.parsers.phyloxml.PhyloXmlDataFormatException;
\r
82 import org.forester.io.parsers.phyloxml.PhyloXmlParser;
\r
83 import org.forester.io.parsers.phyloxml.PhyloXmlUtil;
\r
84 import org.forester.io.parsers.tol.TolParser;
\r
85 import org.forester.io.parsers.util.ParserUtils;
\r
86 import org.forester.io.writers.PhylogenyWriter;
\r
87 import org.forester.io.writers.SequenceWriter;
\r
88 import org.forester.msa.Msa;
\r
89 import org.forester.msa.MsaFormatException;
\r
90 import org.forester.phylogeny.Phylogeny;
\r
91 import org.forester.phylogeny.PhylogenyMethods;
\r
92 import org.forester.phylogeny.PhylogenyNode;
\r
93 import org.forester.phylogeny.PhylogenyNode.NH_CONVERSION_SUPPORT_VALUE_STYLE;
\r
94 import org.forester.phylogeny.data.Confidence;
\r
95 import org.forester.phylogeny.data.PhylogenyDataUtil;
\r
96 import org.forester.phylogeny.data.Sequence;
\r
97 import org.forester.phylogeny.data.Taxonomy;
\r
98 import org.forester.phylogeny.factories.ParserBasedPhylogenyFactory;
\r
99 import org.forester.phylogeny.factories.PhylogenyFactory;
\r
100 import org.forester.phylogeny.iterators.PhylogenyNodeIterator;
\r
101 import org.forester.sequence.MolecularSequence;
\r
102 import org.forester.util.BasicDescriptiveStatistics;
\r
103 import org.forester.util.BasicTable;
\r
104 import org.forester.util.BasicTableParser;
\r
105 import org.forester.util.DescriptiveStatistics;
\r
106 import org.forester.util.ForesterUtil;
\r
108 public final class MainFrameApplication extends MainFrame {
\r
110 public static MainFrameApplication createInstance( final Phylogeny[] phys, final Configuration config ) {
\r
111 return new MainFrameApplication( phys, config );
\r
114 static MainFrame createInstance( final Phylogeny[] phys, final Configuration config, final String title ) {
\r
115 return new MainFrameApplication( phys, config, title );
\r
118 public static MainFrame createInstance( final Phylogeny[] phys,
\r
119 final Configuration config,
\r
120 final String title,
\r
121 final File current_dir ) {
\r
122 return new MainFrameApplication( phys, config, title, current_dir );
\r
125 static MainFrame createInstance( final Phylogeny[] phys, final String config_file_name, final String title ) {
\r
126 return new MainFrameApplication( phys, config_file_name, title );
\r
129 static void warnIfNotPhyloXmlValidation( final Configuration c ) {
\r
130 if ( !c.isValidatePhyloXmlAgainstSchema() ) {
\r
132 .showMessageDialog( null,
\r
134 .wordWrap( "phyloXML XSD-based validation is turned off [enable with line 'validate_against_phyloxml_xsd_schem: true' in configuration file]",
\r
137 JOptionPane.WARNING_MESSAGE );
\r
140 static final String INFER_ANCESTOR_TAXONOMIES = "Infer Ancestor Taxonomies";
\r
141 static final String OBTAIN_DETAILED_TAXONOMIC_INFORMATION = "Obtain Detailed Taxonomic Information";
\r
142 private final static int FRAME_X_SIZE = 800;
\r
143 private final static int FRAME_Y_SIZE = 800;
\r
144 // Filters for the file-open dialog (classes defined in this file)
\r
145 private static final long serialVersionUID = -799735726778865234L;
\r
146 private static final boolean PREPROCESS_TREES = false;
\r
147 private final JFileChooser _values_filechooser;
\r
148 private final JFileChooser _sequences_filechooser;
\r
149 private final JFileChooser _open_filechooser;
\r
150 private final JFileChooser _msa_filechooser;
\r
151 private final JFileChooser _seqs_pi_filechooser;
\r
152 private final JFileChooser _open_filechooser_for_species_tree;
\r
153 // Application-only print menu items
\r
154 private JMenuItem _collapse_below_threshold;
\r
155 private JMenuItem _collapse_below_branch_length;
\r
156 private ButtonGroup _radio_group_1;
\r
157 private ButtonGroup _radio_group_2;
\r
159 double _min_not_collapse = Constants.MIN_NOT_COLLAPSE_DEFAULT;
\r
160 double _min_not_collapse_bl = 0.001;
\r
161 // Phylogeny Inference menu
\r
162 private JMenu _inference_menu;
\r
163 private JMenuItem _inference_from_msa_item;
\r
164 private JMenuItem _inference_from_seqs_item;
\r
165 // Phylogeny Inference
\r
166 private PhylogeneticInferenceOptions _phylogenetic_inference_options = null;
\r
167 private Msa _msa = null;
\r
168 private File _msa_file = null;
\r
169 private List<MolecularSequence> _seqs = null;
\r
170 private File _seqs_file = null;
\r
171 JMenuItem _read_values_jmi;
\r
172 JMenuItem _read_seqs_jmi;
\r
174 private MainFrameApplication( final Phylogeny[] phys, final Configuration config ) {
\r
175 _configuration = config;
\r
176 if ( _configuration == null ) {
\r
177 throw new IllegalArgumentException( "configuration is null" );
\r
179 setVisible( false );
\r
180 setOptions( Options.createInstance( _configuration ) );
\r
181 _mainpanel = new MainPanel( _configuration, this );
\r
182 _open_filechooser = null;
\r
183 _open_filechooser_for_species_tree = null;
\r
184 _save_filechooser = null;
\r
185 _writetopdf_filechooser = null;
\r
186 _writetographics_filechooser = null;
\r
187 _msa_filechooser = null;
\r
188 _seqs_pi_filechooser = null;
\r
189 _values_filechooser = null;
\r
190 _sequences_filechooser = null;
\r
191 _jmenubar = new JMenuBar();
\r
194 _contentpane = getContentPane();
\r
195 _contentpane.setLayout( new BorderLayout() );
\r
196 _contentpane.add( _mainpanel, BorderLayout.CENTER );
\r
198 setSize( MainFrameApplication.FRAME_X_SIZE, MainFrameApplication.FRAME_Y_SIZE );
\r
199 // The window listener
\r
200 setDefaultCloseOperation( WindowConstants.DO_NOTHING_ON_CLOSE );
\r
201 addWindowListener( new WindowAdapter() {
\r
204 public void windowClosing( final WindowEvent e ) {
\r
208 // setVisible( true );
\r
209 if ( ( phys != null ) && ( phys.length > 0 ) ) {
\r
210 AptxUtil.addPhylogeniesToTabs( phys, "", null, _configuration, _mainpanel );
\r
212 getMainPanel().getControlPanel().showWholeAll();
\r
213 getMainPanel().getControlPanel().showWhole();
\r
215 //activateSaveAllIfNeeded();
\r
216 // ...and its children
\r
217 _contentpane.repaint();
\r
220 private MainFrameApplication( final Phylogeny[] phys, final Configuration config, final String title ) {
\r
221 this( phys, config, title, null );
\r
224 private MainFrameApplication( final Phylogeny[] phys,
\r
225 final Configuration config,
\r
226 final String title,
\r
227 final File current_dir ) {
\r
229 _configuration = config;
\r
230 if ( _configuration == null ) {
\r
231 throw new IllegalArgumentException( "configuration is null" );
\r
234 boolean synth_exception = false;
\r
235 if ( Constants.__SYNTH_LF ) {
\r
237 final SynthLookAndFeel synth = new SynthLookAndFeel();
\r
238 synth.load( MainFrameApplication.class.getResourceAsStream( "/resources/synth_look_and_feel_1.xml" ),
\r
239 MainFrameApplication.class );
\r
240 UIManager.setLookAndFeel( synth );
\r
242 catch ( final Exception ex ) {
\r
243 synth_exception = true;
\r
244 ForesterUtil.printWarningMessage( Constants.PRG_NAME,
\r
245 "could not create synth look and feel: "
\r
246 + ex.getLocalizedMessage() );
\r
249 if ( !Constants.__SYNTH_LF || synth_exception ) {
\r
250 if ( _configuration.isUseNativeUI() ) {
\r
251 UIManager.setLookAndFeel( UIManager.getSystemLookAndFeelClassName() );
\r
254 UIManager.setLookAndFeel( UIManager.getCrossPlatformLookAndFeelClassName() );
\r
257 //UIManager.setLookAndFeel( "com.sun.java.swing.plaf.nimbus.NimbusLookAndFeel" );
\r
259 catch ( final UnsupportedLookAndFeelException e ) {
\r
260 AptxUtil.dieWithSystemError( "unsupported look and feel: " + e.toString() );
\r
262 catch ( final ClassNotFoundException e ) {
\r
263 AptxUtil.dieWithSystemError( "class not found exception: " + e.toString() );
\r
265 catch ( final InstantiationException e ) {
\r
266 AptxUtil.dieWithSystemError( "instantiation exception: " + e.toString() );
\r
268 catch ( final IllegalAccessException e ) {
\r
269 AptxUtil.dieWithSystemError( "illegal access exception: " + e.toString() );
\r
271 if ( ( current_dir != null ) && current_dir.canRead() && current_dir.isDirectory() ) {
\r
272 setCurrentDir( current_dir );
\r
274 // hide until everything is ready
\r
275 setVisible( false );
\r
276 setOptions( Options.createInstance( _configuration ) );
\r
277 setInferenceManager( InferenceManager.createInstance( _configuration ) );
\r
278 setPhylogeneticInferenceOptions( PhylogeneticInferenceOptions.createInstance( _configuration ) );
\r
279 // _textframe = null; #~~~~
\r
281 setTitle( Constants.PRG_NAME + " " + Constants.VERSION + " (" + Constants.PRG_DATE + ")" );
\r
282 _mainpanel = new MainPanel( _configuration, this );
\r
283 // The file dialogs
\r
284 _open_filechooser = new JFileChooser();
\r
285 _open_filechooser.setCurrentDirectory( new File( "." ) );
\r
286 _open_filechooser.setMultiSelectionEnabled( false );
\r
287 _open_filechooser.addChoosableFileFilter( MainFrame.xmlfilter );
\r
288 _open_filechooser.addChoosableFileFilter( MainFrame.nhxfilter );
\r
289 _open_filechooser.addChoosableFileFilter( MainFrame.nhfilter );
\r
290 _open_filechooser.addChoosableFileFilter( MainFrame.nexusfilter );
\r
291 _open_filechooser.addChoosableFileFilter( MainFrame.tolfilter );
\r
292 _open_filechooser.addChoosableFileFilter( _open_filechooser.getAcceptAllFileFilter() );
\r
293 _open_filechooser.setFileFilter( MainFrame.defaultfilter );
\r
294 _open_filechooser_for_species_tree = new JFileChooser();
\r
295 _open_filechooser_for_species_tree.setCurrentDirectory( new File( "." ) );
\r
296 _open_filechooser_for_species_tree.setMultiSelectionEnabled( false );
\r
297 _open_filechooser_for_species_tree.addChoosableFileFilter( MainFrame.xmlfilter );
\r
298 _open_filechooser_for_species_tree.addChoosableFileFilter( MainFrame.tolfilter );
\r
299 _open_filechooser_for_species_tree.setFileFilter( MainFrame.xmlfilter );
\r
300 _save_filechooser = new JFileChooser();
\r
301 _save_filechooser.setCurrentDirectory( new File( "." ) );
\r
302 _save_filechooser.setMultiSelectionEnabled( false );
\r
303 _save_filechooser.setFileFilter( MainFrame.xmlfilter );
\r
304 _save_filechooser.addChoosableFileFilter( MainFrame.nhfilter );
\r
305 _save_filechooser.addChoosableFileFilter( MainFrame.nexusfilter );
\r
306 _save_filechooser.addChoosableFileFilter( _save_filechooser.getAcceptAllFileFilter() );
\r
307 _writetopdf_filechooser = new JFileChooser();
\r
308 _writetopdf_filechooser.addChoosableFileFilter( MainFrame.pdffilter );
\r
309 _writetographics_filechooser = new JFileChooser();
\r
310 _writetographics_filechooser.addChoosableFileFilter( MainFrame.graphicsfilefilter );
\r
312 _msa_filechooser = new JFileChooser();
\r
313 _msa_filechooser.setName( "Read Multiple Sequence Alignment File" );
\r
314 _msa_filechooser.setCurrentDirectory( new File( "." ) );
\r
315 _msa_filechooser.setMultiSelectionEnabled( false );
\r
316 _msa_filechooser.addChoosableFileFilter( _msa_filechooser.getAcceptAllFileFilter() );
\r
317 _msa_filechooser.addChoosableFileFilter( MainFrame.msafilter );
\r
319 _seqs_pi_filechooser = new JFileChooser();
\r
320 _seqs_pi_filechooser.setName( "Read Sequences File" );
\r
321 _seqs_pi_filechooser.setCurrentDirectory( new File( "." ) );
\r
322 _seqs_pi_filechooser.setMultiSelectionEnabled( false );
\r
323 _seqs_pi_filechooser.addChoosableFileFilter( _seqs_pi_filechooser.getAcceptAllFileFilter() );
\r
324 _seqs_pi_filechooser.addChoosableFileFilter( MainFrame.seqsfilter );
\r
326 _values_filechooser = new JFileChooser();
\r
327 _values_filechooser.setCurrentDirectory( new File( "." ) );
\r
328 _values_filechooser.setMultiSelectionEnabled( false );
\r
330 _sequences_filechooser = new JFileChooser();
\r
331 _sequences_filechooser.setCurrentDirectory( new File( "." ) );
\r
332 _sequences_filechooser.setMultiSelectionEnabled( false );
\r
333 // build the menu bar
\r
334 _jmenubar = new JMenuBar();
\r
335 if ( !_configuration.isUseNativeUI() ) {
\r
336 _jmenubar.setBackground( getConfiguration().getGuiMenuBackgroundColor() );
\r
339 if ( Constants.__ALLOW_PHYLOGENETIC_INFERENCE ) {
\r
340 buildPhylogeneticInferenceMenu();
\r
342 buildAnalysisMenu();
\r
345 buildFontSizeMenu();
\r
346 buildOptionsMenu();
\r
349 setJMenuBar( _jmenubar );
\r
350 _jmenubar.add( _help_jmenu );
\r
351 _contentpane = getContentPane();
\r
352 _contentpane.setLayout( new BorderLayout() );
\r
353 _contentpane.add( _mainpanel, BorderLayout.CENTER );
\r
355 setSize( MainFrameApplication.FRAME_X_SIZE, MainFrameApplication.FRAME_Y_SIZE );
\r
356 // addWindowFocusListener( new WindowAdapter() {
\r
359 // public void windowGainedFocus( WindowEvent e ) {
\r
360 // requestFocusInWindow();
\r
363 // The window listener
\r
364 setDefaultCloseOperation( WindowConstants.DO_NOTHING_ON_CLOSE );
\r
365 addWindowListener( new WindowAdapter() {
\r
368 public void windowClosing( final WindowEvent e ) {
\r
369 if ( isUnsavedDataPresent() ) {
\r
370 final int r = JOptionPane.showConfirmDialog( null,
\r
371 "Exit despite potentially unsaved changes?",
\r
373 JOptionPane.YES_NO_OPTION );
\r
374 if ( r != JOptionPane.YES_OPTION ) {
\r
379 final int r = JOptionPane.showConfirmDialog( null,
\r
380 "Exit Archaeopteryx?",
\r
382 JOptionPane.YES_NO_OPTION );
\r
383 if ( r != JOptionPane.YES_OPTION ) {
\r
390 // The component listener
\r
391 addComponentListener( new ComponentAdapter() {
\r
394 public void componentResized( final ComponentEvent e ) {
\r
395 if ( _mainpanel.getCurrentTreePanel() != null ) {
\r
396 _mainpanel.getCurrentTreePanel().calcParametersForPainting( _mainpanel.getCurrentTreePanel()
\r
398 _mainpanel.getCurrentTreePanel()
\r
403 requestFocusInWindow();
\r
404 // addKeyListener( this );
\r
405 setVisible( true );
\r
406 if ( ( phys != null ) && ( phys.length > 0 ) ) {
\r
407 AptxUtil.addPhylogeniesToTabs( phys, title, null, _configuration, _mainpanel );
\r
409 getMainPanel().getControlPanel().showWholeAll();
\r
410 getMainPanel().getControlPanel().showWhole();
\r
412 activateSaveAllIfNeeded();
\r
413 // ...and its children
\r
414 _contentpane.repaint();
\r
418 private MainFrameApplication( final Phylogeny[] phys, final String config_file, final String title ) {
\r
419 // Reads the config file (false, false => not url, not applet):
\r
420 this( phys, new Configuration( config_file, false, false, true ), title );
\r
424 public void actionPerformed( final ActionEvent e ) {
\r
426 super.actionPerformed( e );
\r
427 final Object o = e.getSource();
\r
428 // Handle app-specific actions here:
\r
429 if ( o == _open_item ) {
\r
430 readPhylogeniesFromFile();
\r
432 if ( o == _open_url_item ) {
\r
433 readPhylogeniesFromURL();
\r
435 else if ( o == _new_item ) {
\r
438 else if ( o == _save_all_item ) {
\r
441 else if ( o == _close_item ) {
\r
442 closeCurrentPane();
\r
444 else if ( o == _write_to_jpg_item ) {
\r
445 writeToGraphicsFile( _mainpanel.getCurrentPhylogeny(), GraphicsExportType.JPG );
\r
447 else if ( o == _write_to_gif_item ) {
\r
448 writeToGraphicsFile( _mainpanel.getCurrentPhylogeny(), GraphicsExportType.GIF );
\r
450 else if ( o == _write_to_tif_item ) {
\r
451 writeToGraphicsFile( _mainpanel.getCurrentPhylogeny(), GraphicsExportType.TIFF );
\r
453 else if ( o == _write_to_bmp_item ) {
\r
454 writeToGraphicsFile( _mainpanel.getCurrentPhylogeny(), GraphicsExportType.BMP );
\r
456 else if ( o == _print_item ) {
\r
459 else if ( o == _load_species_tree_item ) {
\r
460 readSpeciesTreeFromFile();
\r
462 else if ( o == _lineage_inference ) {
\r
463 if ( isSubtreeDisplayed() ) {
\r
464 JOptionPane.showMessageDialog( this,
\r
465 "Subtree is shown.",
\r
466 "Cannot infer ancestral taxonomies",
\r
467 JOptionPane.ERROR_MESSAGE );
\r
470 executeLineageInference();
\r
472 else if ( o == _obtain_detailed_taxonomic_information_jmi ) {
\r
473 if ( isSubtreeDisplayed() ) {
\r
476 obtainDetailedTaxonomicInformation();
\r
478 else if ( o == _obtain_detailed_taxonomic_information_deleting_jmi ) {
\r
479 if ( isSubtreeDisplayed() ) {
\r
482 obtainDetailedTaxonomicInformationDelete();
\r
484 else if ( o == _obtain_seq_information_jmi ) {
\r
485 obtainSequenceInformation();
\r
487 else if ( o == _read_values_jmi ) {
\r
488 if ( isSubtreeDisplayed() ) {
\r
491 addExpressionValuesFromFile();
\r
493 else if ( o == _read_seqs_jmi ) {
\r
494 if ( isSubtreeDisplayed() ) {
\r
497 addSequencesFromFile();
\r
499 else if ( o == _move_node_names_to_tax_sn_jmi ) {
\r
500 moveNodeNamesToTaxSn();
\r
502 else if ( o == _move_node_names_to_seq_names_jmi ) {
\r
503 moveNodeNamesToSeqNames();
\r
505 else if ( o == _extract_tax_code_from_node_names_jmi ) {
\r
506 extractTaxDataFromNodeNames();
\r
508 else if ( o == _internal_number_are_confidence_for_nh_parsing_cbmi ) {
\r
509 updateOptions( getOptions() );
\r
511 else if ( o == _replace_underscores_cbmi ) {
\r
512 if ( ( _extract_taxonomy_no_rbmi != null ) && !_extract_taxonomy_no_rbmi.isSelected() ) {
\r
513 _extract_taxonomy_no_rbmi.setSelected( true );
\r
515 updateOptions( getOptions() );
\r
517 else if ( o == _allow_errors_in_distance_to_parent_cbmi ) {
\r
518 updateOptions( getOptions() );
\r
520 else if ( o == _collapse_below_threshold ) {
\r
521 if ( isSubtreeDisplayed() ) {
\r
524 collapseBelowThreshold();
\r
526 else if ( o == _collapse_below_branch_length ) {
\r
527 if ( isSubtreeDisplayed() ) {
\r
530 collapseBelowBranchLengthThreshold();
\r
532 else if ( ( o == _extract_taxonomy_pfam_strict_rbmi ) || ( o == _extract_taxonomy_pfam_relaxed_rbmi )
\r
533 || ( o == _extract_taxonomy_agressive_rbmi ) ) {
\r
534 if ( _replace_underscores_cbmi != null ) {
\r
535 _replace_underscores_cbmi.setSelected( false );
\r
537 updateOptions( getOptions() );
\r
539 else if ( o == _extract_taxonomy_no_rbmi ) {
\r
540 updateOptions( getOptions() );
\r
542 else if ( o == _inference_from_msa_item ) {
\r
543 executePhyleneticInference( false );
\r
545 else if ( o == _inference_from_seqs_item ) {
\r
546 executePhyleneticInference( true );
\r
548 _contentpane.repaint();
\r
550 catch ( final Exception ex ) {
\r
551 AptxUtil.unexpectedException( ex );
\r
553 catch ( final Error err ) {
\r
554 AptxUtil.unexpectedError( err );
\r
558 private void addExpressionValuesFromFile() {
\r
559 if ( ( getCurrentTreePanel() == null ) || ( getCurrentTreePanel().getPhylogeny() == null ) ) {
\r
560 JOptionPane.showMessageDialog( this,
\r
561 "Need to load evolutionary tree first",
\r
562 "Can Not Read Expression Values",
\r
563 JOptionPane.WARNING_MESSAGE );
\r
566 final File my_dir = getCurrentDir();
\r
567 if ( my_dir != null ) {
\r
568 _values_filechooser.setCurrentDirectory( my_dir );
\r
570 final int result = _values_filechooser.showOpenDialog( _contentpane );
\r
571 final File file = _values_filechooser.getSelectedFile();
\r
572 if ( ( file != null ) && ( file.length() > 0 ) && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
573 BasicTable<String> t = null;
\r
575 t = BasicTableParser.parse( file, '\t' );
\r
576 if ( t.getNumberOfColumns() < 2 ) {
\r
577 t = BasicTableParser.parse( file, ',' );
\r
579 if ( t.getNumberOfColumns() < 2 ) {
\r
580 t = BasicTableParser.parse( file, ' ' );
\r
583 catch ( final IOException e ) {
\r
584 JOptionPane.showMessageDialog( this,
\r
586 "Could Not Read Expression Value Table",
\r
587 JOptionPane.ERROR_MESSAGE );
\r
590 if ( t.getNumberOfColumns() < 2 ) {
\r
591 JOptionPane.showMessageDialog( this,
\r
592 "Table contains " + t.getNumberOfColumns() + " column(s)",
\r
593 "Problem with Expression Value Table",
\r
594 JOptionPane.ERROR_MESSAGE );
\r
597 if ( t.getNumberOfRows() < 1 ) {
\r
598 JOptionPane.showMessageDialog( this,
\r
599 "Table contains zero rows",
\r
600 "Problem with Expression Value Table",
\r
601 JOptionPane.ERROR_MESSAGE );
\r
604 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
605 if ( t.getNumberOfRows() != phy.getNumberOfExternalNodes() ) {
\r
606 JOptionPane.showMessageDialog( this,
\r
607 "Table contains " + t.getNumberOfRows() + " rows, but tree contains "
\r
608 + phy.getNumberOfExternalNodes() + " external nodes",
\r
610 JOptionPane.WARNING_MESSAGE );
\r
612 final DescriptiveStatistics stats = new BasicDescriptiveStatistics();
\r
614 for( final PhylogenyNodeIterator iter = phy.iteratorPreorder(); iter.hasNext(); ) {
\r
615 final PhylogenyNode node = iter.next();
\r
616 final String node_name = node.getName();
\r
617 if ( !ForesterUtil.isEmpty( node_name ) ) {
\r
620 row = t.findRow( node_name );
\r
622 catch ( final IllegalArgumentException e ) {
\r
624 .showMessageDialog( this,
\r
626 "Error Mapping Node Identifiers to Expression Value Identifiers",
\r
627 JOptionPane.ERROR_MESSAGE );
\r
631 if ( node.isExternal() ) {
\r
636 final List<Double> l = new ArrayList<Double>();
\r
637 for( int col = 1; col < t.getNumberOfColumns(); ++col ) {
\r
640 d = Double.parseDouble( t.getValueAsString( col, row ) );
\r
642 catch ( final NumberFormatException e ) {
\r
643 JOptionPane.showMessageDialog( this,
\r
644 "Could not parse \"" + t.getValueAsString( col, row )
\r
645 + "\" into a decimal value",
\r
646 "Issue with Expression Value Table",
\r
647 JOptionPane.ERROR_MESSAGE );
\r
650 stats.addValue( d );
\r
653 if ( !l.isEmpty() ) {
\r
654 if ( node.getNodeData().getProperties() != null ) {
\r
655 node.getNodeData().getProperties()
\r
656 .removePropertiesWithGivenReferencePrefix( PhyloXmlUtil.VECTOR_PROPERTY_REF );
\r
658 node.getNodeData().setVector( l );
\r
662 if ( not_found > 0 ) {
\r
663 JOptionPane.showMessageDialog( this, "Could not fine expression values for " + not_found
\r
664 + " external node(s)", "Warning", JOptionPane.WARNING_MESSAGE );
\r
666 getCurrentTreePanel().setStatisticsForExpressionValues( stats );
\r
670 private void addSequencesFromFile() {
\r
671 if ( ( getCurrentTreePanel() == null ) || ( getCurrentTreePanel().getPhylogeny() == null ) ) {
\r
672 JOptionPane.showMessageDialog( this,
\r
673 "Need to load evolutionary tree first",
\r
674 "Can Not Read Sequences",
\r
675 JOptionPane.WARNING_MESSAGE );
\r
678 final File my_dir = getCurrentDir();
\r
679 if ( my_dir != null ) {
\r
680 _sequences_filechooser.setCurrentDirectory( my_dir );
\r
682 final int result = _sequences_filechooser.showOpenDialog( _contentpane );
\r
683 final File file = _sequences_filechooser.getSelectedFile();
\r
684 List<MolecularSequence> seqs = null;
\r
685 if ( ( file != null ) && !file.isDirectory() && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
687 if ( FastaParser.isLikelyFasta( new FileInputStream( file ) ) ) {
\r
688 seqs = FastaParser.parse( new FileInputStream( file ) );
\r
691 JOptionPane.showMessageDialog( this,
\r
692 "Format does not appear to be Fasta",
\r
693 "Multiple sequence file format error",
\r
694 JOptionPane.ERROR_MESSAGE );
\r
698 catch ( final MsaFormatException e ) {
\r
700 JOptionPane.showMessageDialog( this,
\r
701 e.getLocalizedMessage(),
\r
702 "Multiple sequence file format error",
\r
703 JOptionPane.ERROR_MESSAGE );
\r
706 catch ( final IOException e ) {
\r
708 JOptionPane.showMessageDialog( this,
\r
709 e.getLocalizedMessage(),
\r
710 "Failed to read multiple sequence file",
\r
711 JOptionPane.ERROR_MESSAGE );
\r
714 catch ( final Exception e ) {
\r
716 e.printStackTrace();
\r
717 JOptionPane.showMessageDialog( this,
\r
718 e.getLocalizedMessage(),
\r
719 "Unexpected error during reading of multiple sequence file",
\r
720 JOptionPane.ERROR_MESSAGE );
\r
723 if ( ( seqs == null ) || ( seqs.size() < 1 ) ) {
\r
724 JOptionPane.showMessageDialog( this,
\r
725 "Multiple sequence file is empty",
\r
726 "Empty multiple sequence file",
\r
727 JOptionPane.ERROR_MESSAGE );
\r
732 if ( seqs != null ) {
\r
733 for( final MolecularSequence seq : seqs ) {
\r
734 System.out.println( seq.getIdentifier() );
\r
736 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
737 int total_counter = 0;
\r
738 int attached_counter = 0;
\r
739 for( final MolecularSequence seq : seqs ) {
\r
741 final String seq_name = seq.getIdentifier();
\r
742 if ( !ForesterUtil.isEmpty( seq_name ) ) {
\r
743 List<PhylogenyNode> nodes = phy.getNodesViaSequenceName( seq_name );
\r
744 if ( nodes.isEmpty() ) {
\r
745 nodes = phy.getNodesViaSequenceSymbol( seq_name );
\r
747 if ( nodes.isEmpty() ) {
\r
748 nodes = phy.getNodesViaGeneName( seq_name );
\r
750 if ( nodes.isEmpty() ) {
\r
751 nodes = phy.getNodes( seq_name );
\r
753 if ( nodes.size() > 1 ) {
\r
754 JOptionPane.showMessageDialog( this,
\r
755 "Sequence name \"" + seq_name + "\" is not unique",
\r
756 "Sequence name not unique",
\r
757 JOptionPane.ERROR_MESSAGE );
\r
761 final String[] a = seq_name.split( "\\s" );
\r
762 if ( nodes.isEmpty() && ( a.length > 1 ) ) {
\r
763 final String seq_name_split = a[ 0 ];
\r
764 nodes = phy.getNodesViaSequenceName( seq_name_split );
\r
765 if ( nodes.isEmpty() ) {
\r
766 nodes = phy.getNodesViaSequenceSymbol( seq_name_split );
\r
768 if ( nodes.isEmpty() ) {
\r
769 nodes = phy.getNodes( seq_name_split );
\r
771 if ( nodes.size() > 1 ) {
\r
772 JOptionPane.showMessageDialog( this, "Split sequence name \"" + seq_name_split
\r
773 + "\" is not unique", "Sequence name not unique", JOptionPane.ERROR_MESSAGE );
\r
778 if ( nodes.size() == 1 ) {
\r
779 ++attached_counter;
\r
780 final PhylogenyNode n = nodes.get( 0 );
\r
781 if ( !n.getNodeData().isHasSequence() ) {
\r
782 n.getNodeData().addSequence( new org.forester.phylogeny.data.Sequence() );
\r
784 n.getNodeData().getSequence().setMolecularSequence( seq.getMolecularSequenceAsString() );
\r
785 if ( ForesterUtil.isEmpty( n.getNodeData().getSequence().getName() ) ) {
\r
786 n.getNodeData().getSequence().setName( seq_name );
\r
791 if ( attached_counter > 0 ) {
\r
793 int ext_nodes_with_seq = 0;
\r
794 for( final PhylogenyNodeIterator iter = phy.iteratorExternalForward(); iter.hasNext(); ) {
\r
796 final PhylogenyNode n = iter.next();
\r
797 if ( n.getNodeData().isHasSequence()
\r
798 && !ForesterUtil.isEmpty( n.getNodeData().getSequence().getMolecularSequence() ) ) {
\r
799 ++ext_nodes_with_seq;
\r
803 if ( ext_nodes == ext_nodes_with_seq ) {
\r
804 s = "All " + ext_nodes_with_seq + " external nodes now have a molecular sequence attached to them.";
\r
807 s = ext_nodes_with_seq + " out of " + ext_nodes
\r
808 + " external nodes now have a molecular sequence attached to them.";
\r
810 if ( ( attached_counter == total_counter ) && ( ext_nodes == ext_nodes_with_seq ) ) {
\r
811 JOptionPane.showMessageDialog( this,
\r
812 "Attached all " + total_counter + " sequences to tree nodes.\n" + s,
\r
813 "All sequences attached",
\r
814 JOptionPane.INFORMATION_MESSAGE );
\r
817 JOptionPane.showMessageDialog( this, "Attached " + attached_counter
\r
818 + " sequences out of a total of " + total_counter + " sequences.\n" + s, attached_counter
\r
819 + " sequences attached", JOptionPane.WARNING_MESSAGE );
\r
823 JOptionPane.showMessageDialog( this, "No maching tree node for any of the " + total_counter
\r
824 + " sequences", "Could not attach any sequences", JOptionPane.ERROR_MESSAGE );
\r
829 void buildAnalysisMenu() {
\r
830 _analysis_menu = MainFrame.createMenu( "Analysis", getConfiguration() );
\r
831 _analysis_menu.add( _gsdi_item = new JMenuItem( "GSDI (Generalized Speciation Duplication Inference)" ) );
\r
832 _analysis_menu.add( _gsdir_item = new JMenuItem( "GSDIR (GSDI with re-rooting)" ) );
\r
833 _analysis_menu.add( _load_species_tree_item = new JMenuItem( "Load Species Tree..." ) );
\r
834 customizeJMenuItem( _gsdi_item );
\r
835 customizeJMenuItem( _gsdir_item );
\r
836 customizeJMenuItem( _load_species_tree_item );
\r
837 _analysis_menu.addSeparator();
\r
838 _analysis_menu.add( _lineage_inference = new JMenuItem( INFER_ANCESTOR_TAXONOMIES ) );
\r
839 customizeJMenuItem( _lineage_inference );
\r
840 _lineage_inference.setToolTipText( "Inference of ancestor taxonomies/lineages" );
\r
841 _jmenubar.add( _analysis_menu );
\r
845 void buildFileMenu() {
\r
846 _file_jmenu = MainFrame.createMenu( "File", getConfiguration() );
\r
847 _file_jmenu.add( _open_item = new JMenuItem( "Read Tree from File..." ) );
\r
848 _file_jmenu.addSeparator();
\r
849 _file_jmenu.add( _open_url_item = new JMenuItem( "Read Tree from URL/Webservice..." ) );
\r
850 _file_jmenu.addSeparator();
\r
851 final WebservicesManager webservices_manager = WebservicesManager.getInstance();
\r
852 _load_phylogeny_from_webservice_menu_items = new JMenuItem[ webservices_manager
\r
853 .getAvailablePhylogeniesWebserviceClients().size() ];
\r
854 for( int i = 0; i < webservices_manager.getAvailablePhylogeniesWebserviceClients().size(); ++i ) {
\r
855 final PhylogeniesWebserviceClient client = webservices_manager.getAvailablePhylogeniesWebserviceClient( i );
\r
856 _load_phylogeny_from_webservice_menu_items[ i ] = new JMenuItem( client.getMenuName() );
\r
857 _file_jmenu.add( _load_phylogeny_from_webservice_menu_items[ i ] );
\r
859 if ( getConfiguration().isEditable() ) {
\r
860 _file_jmenu.addSeparator();
\r
861 _file_jmenu.add( _new_item = new JMenuItem( "New" ) );
\r
862 _new_item.setToolTipText( "to create a new tree with one node, as source for manual tree construction" );
\r
864 _file_jmenu.addSeparator();
\r
865 _file_jmenu.add( _save_item = new JMenuItem( "Save Tree As..." ) );
\r
866 _file_jmenu.add( _save_all_item = new JMenuItem( "Save All Trees As..." ) );
\r
867 _save_all_item.setToolTipText( "Write all phylogenies to one file." );
\r
868 _save_all_item.setEnabled( false );
\r
869 _file_jmenu.addSeparator();
\r
870 _file_jmenu.add( _write_to_pdf_item = new JMenuItem( "Export to PDF file ..." ) );
\r
871 if ( AptxUtil.canWriteFormat( "tif" ) || AptxUtil.canWriteFormat( "tiff" ) || AptxUtil.canWriteFormat( "TIF" ) ) {
\r
872 _file_jmenu.add( _write_to_tif_item = new JMenuItem( "Export to TIFF file..." ) );
\r
874 _file_jmenu.add( _write_to_png_item = new JMenuItem( "Export to PNG file..." ) );
\r
875 _file_jmenu.add( _write_to_jpg_item = new JMenuItem( "Export to JPG file..." ) );
\r
876 if ( AptxUtil.canWriteFormat( "gif" ) ) {
\r
877 _file_jmenu.add( _write_to_gif_item = new JMenuItem( "Export to GIF file..." ) );
\r
879 if ( AptxUtil.canWriteFormat( "bmp" ) ) {
\r
880 _file_jmenu.add( _write_to_bmp_item = new JMenuItem( "Export to BMP file..." ) );
\r
882 _file_jmenu.addSeparator();
\r
883 _file_jmenu.add( _print_item = new JMenuItem( "Print..." ) );
\r
884 _file_jmenu.addSeparator();
\r
885 _file_jmenu.add( _close_item = new JMenuItem( "Close Tab" ) );
\r
886 _close_item.setToolTipText( "To close the current pane." );
\r
887 _close_item.setEnabled( true );
\r
888 _file_jmenu.addSeparator();
\r
889 _file_jmenu.add( _exit_item = new JMenuItem( "Exit" ) );
\r
890 customizeJMenuItem( _open_item );
\r
892 .setFont( new Font( _open_item.getFont().getFontName(), Font.BOLD, _open_item.getFont().getSize() + 4 ) );
\r
893 customizeJMenuItem( _open_url_item );
\r
894 for( int i = 0; i < webservices_manager.getAvailablePhylogeniesWebserviceClients().size(); ++i ) {
\r
895 customizeJMenuItem( _load_phylogeny_from_webservice_menu_items[ i ] );
\r
897 customizeJMenuItem( _save_item );
\r
898 if ( getConfiguration().isEditable() ) {
\r
899 customizeJMenuItem( _new_item );
\r
901 customizeJMenuItem( _close_item );
\r
902 customizeJMenuItem( _save_all_item );
\r
903 customizeJMenuItem( _write_to_pdf_item );
\r
904 customizeJMenuItem( _write_to_png_item );
\r
905 customizeJMenuItem( _write_to_jpg_item );
\r
906 customizeJMenuItem( _write_to_gif_item );
\r
907 customizeJMenuItem( _write_to_tif_item );
\r
908 customizeJMenuItem( _write_to_bmp_item );
\r
909 customizeJMenuItem( _print_item );
\r
910 customizeJMenuItem( _exit_item );
\r
911 _jmenubar.add( _file_jmenu );
\r
914 void buildOptionsMenu() {
\r
915 _options_jmenu = MainFrame.createMenu( OPTIONS_HEADER, getConfiguration() );
\r
916 _options_jmenu.addChangeListener( new ChangeListener() {
\r
919 public void stateChanged( final ChangeEvent e ) {
\r
920 MainFrame.setOvPlacementColorChooseMenuItem( _overview_placment_mi, getOptions() );
\r
921 MainFrame.setTextColorChooseMenuItem( _switch_colors_mi, getCurrentTreePanel() );
\r
923 .setTextMinSupportMenuItem( _choose_minimal_confidence_mi, getOptions(), getCurrentTreePanel() );
\r
924 MainFrame.setTextForFontChooserMenuItem( _choose_font_mi, MainFrame
\r
925 .createCurrentFontDesc( getMainPanel().getTreeFontSet() ) );
\r
926 MainFrame.setTextForGraphicsSizeChooserMenuItem( _print_size_mi, getOptions() );
\r
927 MainFrame.setTextForPdfLineWidthChooserMenuItem( _choose_pdf_width_mi, getOptions() );
\r
928 MainFrame.setCycleNodeFillMenuItem( _cycle_node_fill_mi, getOptions() );
\r
929 MainFrame.setCycleNodeShapeMenuItem( _cycle_node_shape_mi, getOptions() );
\r
930 MainFrame.setCycleDataReturnMenuItem( _cycle_data_return, getOptions() );
\r
931 MainFrame.setTextNodeSizeMenuItem( _choose_node_size_mi, getOptions() );
\r
933 getMainPanel().getControlPanel().setVisibilityOfDomainStrucureCB();
\r
934 getMainPanel().getControlPanel().setVisibilityOfX();
\r
936 catch ( final Exception ignore ) {
\r
937 // do nothing, not important.
\r
941 _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( DISPLAY_SUBHEADER ), getConfiguration() ) );
\r
943 .add( _ext_node_dependent_cladogram_rbmi = new JRadioButtonMenuItem( MainFrame.NONUNIFORM_CLADOGRAMS_LABEL ) );
\r
944 _options_jmenu.add( _uniform_cladograms_rbmi = new JRadioButtonMenuItem( MainFrame.UNIFORM_CLADOGRAMS_LABEL ) );
\r
945 _options_jmenu.add( _non_lined_up_cladograms_rbmi = new JRadioButtonMenuItem( NON_LINED_UP_CLADOGRAMS_LABEL ) );
\r
946 _radio_group_1 = new ButtonGroup();
\r
947 _radio_group_1.add( _ext_node_dependent_cladogram_rbmi );
\r
948 _radio_group_1.add( _uniform_cladograms_rbmi );
\r
949 _radio_group_1.add( _non_lined_up_cladograms_rbmi );
\r
950 _options_jmenu.add( _show_overview_cbmi = new JCheckBoxMenuItem( SHOW_OVERVIEW_LABEL ) );
\r
951 _options_jmenu.add( _show_scale_cbmi = new JCheckBoxMenuItem( DISPLAY_SCALE_LABEL ) );
\r
953 .add( _show_default_node_shapes_internal_cbmi = new JCheckBoxMenuItem( DISPLAY_NODE_BOXES_LABEL_INT ) );
\r
955 .add( _show_default_node_shapes_external_cbmi = new JCheckBoxMenuItem( DISPLAY_NODE_BOXES_LABEL_EXT ) );
\r
957 .add( _show_default_node_shapes_for_marked_cbmi = new JCheckBoxMenuItem( MainFrame.DISPLAY_NODE_BOXES_LABEL_MARKED ) );
\r
958 _options_jmenu.add( _line_up_renderable_data_cbmi = new JCheckBoxMenuItem( MainFrame.LINE_UP_RENDERABLE_DATA ) );
\r
959 if ( getConfiguration().doDisplayOption( Configuration.show_domain_architectures ) ) {
\r
960 _options_jmenu.add( _right_line_up_domains_cbmi = new JCheckBoxMenuItem( MainFrame.RIGHT_LINE_UP_DOMAINS ) );
\r
961 _options_jmenu.add( _show_domain_labels = new JCheckBoxMenuItem( MainFrame.SHOW_DOMAIN_LABELS_LABEL ) );
\r
963 _options_jmenu.add( _show_annotation_ref_source = new JCheckBoxMenuItem( SHOW_ANN_REF_SOURCE_LABEL ) );
\r
964 _options_jmenu.add( _show_confidence_stddev_cbmi = new JCheckBoxMenuItem( SHOW_CONF_STDDEV_LABEL ) );
\r
965 _options_jmenu.add( _color_by_taxonomic_group_cbmi = new JCheckBoxMenuItem( COLOR_BY_TAXONOMIC_GROUP ) );
\r
966 _options_jmenu.add( _color_labels_same_as_parent_branch = new JCheckBoxMenuItem( COLOR_LABELS_LABEL ) );
\r
967 _color_labels_same_as_parent_branch.setToolTipText( MainFrame.COLOR_LABELS_TIP );
\r
968 _options_jmenu.add( _abbreviate_scientific_names = new JCheckBoxMenuItem( ABBREV_SN_LABEL ) );
\r
969 _options_jmenu.add( _label_direction_cbmi = new JCheckBoxMenuItem( LABEL_DIRECTION_LABEL ) );
\r
970 _label_direction_cbmi.setToolTipText( LABEL_DIRECTION_TIP );
\r
971 _options_jmenu.add( _screen_antialias_cbmi = new JCheckBoxMenuItem( SCREEN_ANTIALIAS_LABEL ) );
\r
972 _options_jmenu.add( _background_gradient_cbmi = new JCheckBoxMenuItem( BG_GRAD_LABEL ) );
\r
973 _options_jmenu.add( _cycle_node_shape_mi = new JMenuItem( MainFrame.CYCLE_NODE_SHAPE_LABEL ) );
\r
974 _options_jmenu.add( _cycle_node_fill_mi = new JMenuItem( MainFrame.CYCLE_NODE_FILL_LABEL ) );
\r
975 _options_jmenu.add( _choose_node_size_mi = new JMenuItem( MainFrame.CHOOSE_NODE_SIZE_LABEL ) );
\r
976 _options_jmenu.add( _choose_minimal_confidence_mi = new JMenuItem( "" ) );
\r
977 _options_jmenu.add( _overview_placment_mi = new JMenuItem( "" ) );
\r
978 _options_jmenu.add( _switch_colors_mi = new JMenuItem( "" ) );
\r
979 _options_jmenu.add( _choose_font_mi = new JMenuItem( "" ) );
\r
980 _options_jmenu.addSeparator();
\r
981 _options_jmenu.add( _cycle_data_return = new JMenuItem( "Cycle Data Return" ) );
\r
982 _options_jmenu.addSeparator();
\r
983 _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( SEARCH_SUBHEADER ), getConfiguration() ) );
\r
984 _options_jmenu.add( _search_case_senstive_cbmi = new JCheckBoxMenuItem( SEARCH_CASE_SENSITIVE_LABEL ) );
\r
985 _options_jmenu.add( _search_whole_words_only_cbmi = new JCheckBoxMenuItem( SEARCH_TERMS_ONLY_LABEL ) );
\r
986 _options_jmenu.add( _search_with_regex_cbmi = new JCheckBoxMenuItem( MainFrame.SEARCH_REGEX_LABEL ) );
\r
987 _search_with_regex_cbmi.setToolTipText( MainFrame.SEARCH_WITH_REGEX_TIP );
\r
988 _options_jmenu.add( _inverse_search_result_cbmi = new JCheckBoxMenuItem( INVERSE_SEARCH_RESULT_LABEL ) );
\r
989 _options_jmenu.addSeparator();
\r
990 _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( "Graphics Export & Printing:" ),
\r
991 getConfiguration() ) );
\r
992 _options_jmenu.add( _antialias_print_cbmi = new JCheckBoxMenuItem( "Antialias" ) );
\r
993 _options_jmenu.add( _print_black_and_white_cbmi = new JCheckBoxMenuItem( "Export in Black and White" ) );
\r
995 .add( _print_using_actual_size_cbmi = new JCheckBoxMenuItem( "Use Current Image Size for PDF export and Printing" ) );
\r
997 .add( _graphics_export_using_actual_size_cbmi = new JCheckBoxMenuItem( "Use Current Image Size for PNG, JPG, and GIF export" ) );
\r
999 .add( _graphics_export_visible_only_cbmi = new JCheckBoxMenuItem( "Limit to Visible ('Screenshot') for PNG, JPG, and GIF export" ) );
\r
1000 _options_jmenu.add( _print_size_mi = new JMenuItem( "" ) );
\r
1001 _options_jmenu.add( _choose_pdf_width_mi = new JMenuItem( "" ) );
\r
1002 _options_jmenu.addSeparator();
\r
1003 _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( "Newick/NHX/Nexus Input:" ), getConfiguration() ) );
\r
1005 .add( _internal_number_are_confidence_for_nh_parsing_cbmi = new JCheckBoxMenuItem( "Internal Node Names are Confidence Values" ) );
\r
1006 _options_jmenu.add( _replace_underscores_cbmi = new JCheckBoxMenuItem( "Replace Underscores with Spaces" ) );
\r
1008 .add( _allow_errors_in_distance_to_parent_cbmi = new JCheckBoxMenuItem( "Ignore Distance Values Format Errors" ) );
\r
1009 _options_jmenu.add( _extract_taxonomy_no_rbmi = new JRadioButtonMenuItem( "No Taxonomy Extraction" ) );
\r
1011 .add( _extract_taxonomy_pfam_strict_rbmi = new JRadioButtonMenuItem( "Extract Taxonomy Codes/Ids from Pfam-style Node Names" ) );
\r
1013 .add( _extract_taxonomy_pfam_relaxed_rbmi = new JRadioButtonMenuItem( "Extract Taxonomy Codes/Ids from Pfam-style like Node Names" ) );
\r
1015 .add( _extract_taxonomy_agressive_rbmi = new JRadioButtonMenuItem( "Extract Taxonomy Codes/Ids/Scientific Names from Node Names" ) );
\r
1016 _extract_taxonomy_pfam_strict_rbmi
\r
1017 .setToolTipText( "To extract taxonomy codes/ids from node names in the form of e.g. \"BCL2_MOUSE/123-304\" or \"BCL2_10090/123-304\"" );
\r
1018 _extract_taxonomy_pfam_relaxed_rbmi
\r
1019 .setToolTipText( "To extract taxonomy codes/ids from node names in the form of e.g. \"bax_MOUSE\" or \"bax_10090\"" );
\r
1020 _extract_taxonomy_agressive_rbmi
\r
1021 .setToolTipText( "To extract taxonomy codes/ids or scientific names from node names in the form of e.g. \"MOUSE\" or \"10090\" or \"xyz_Nematostella_vectensis\"" );
\r
1022 _radio_group_2 = new ButtonGroup();
\r
1023 _radio_group_2.add( _extract_taxonomy_no_rbmi );
\r
1024 _radio_group_2.add( _extract_taxonomy_pfam_strict_rbmi );
\r
1025 _radio_group_2.add( _extract_taxonomy_pfam_relaxed_rbmi );
\r
1026 _radio_group_2.add( _extract_taxonomy_agressive_rbmi );
\r
1027 _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( "Newick/Nexus Output:" ), getConfiguration() ) );
\r
1029 .add( _use_brackets_for_conf_in_nh_export_cbmi = new JCheckBoxMenuItem( USE_BRACKETS_FOR_CONF_IN_NH_LABEL ) );
\r
1030 _use_brackets_for_conf_in_nh_export_cbmi
\r
1031 .setToolTipText( "e.g. \"0.1[90]\" for a branch with support 90 and a length of 0.1" );
\r
1033 .add( _use_internal_names_for_conf_in_nh_export_cbmi = new JCheckBoxMenuItem( USE_INTERNAL_NAMES_FOR_CONF_IN_NH_LABEL ) );
\r
1034 customizeJMenuItem( _choose_font_mi );
\r
1035 customizeJMenuItem( _choose_minimal_confidence_mi );
\r
1036 customizeJMenuItem( _switch_colors_mi );
\r
1037 customizeJMenuItem( _print_size_mi );
\r
1038 customizeJMenuItem( _choose_pdf_width_mi );
\r
1039 customizeJMenuItem( _overview_placment_mi );
\r
1040 customizeCheckBoxMenuItem( _show_default_node_shapes_external_cbmi, getOptions()
\r
1041 .isShowDefaultNodeShapesExternal() );
\r
1042 customizeCheckBoxMenuItem( _show_default_node_shapes_internal_cbmi, getOptions()
\r
1043 .isShowDefaultNodeShapesInternal() );
\r
1044 customizeCheckBoxMenuItem( _show_default_node_shapes_for_marked_cbmi, getOptions()
\r
1045 .isShowDefaultNodeShapesForMarkedNodes() );
\r
1046 customizeJMenuItem( _cycle_node_shape_mi );
\r
1047 customizeJMenuItem( _cycle_node_fill_mi );
\r
1048 customizeJMenuItem( _choose_node_size_mi );
\r
1049 customizeJMenuItem( _cycle_data_return );
\r
1050 customizeCheckBoxMenuItem( _color_labels_same_as_parent_branch, getOptions().isColorLabelsSameAsParentBranch() );
\r
1051 customizeCheckBoxMenuItem( _color_by_taxonomic_group_cbmi, getOptions().isColorByTaxonomicGroup() );
\r
1052 customizeCheckBoxMenuItem( _screen_antialias_cbmi, getOptions().isAntialiasScreen() );
\r
1053 customizeCheckBoxMenuItem( _background_gradient_cbmi, getOptions().isBackgroundColorGradient() );
\r
1054 customizeCheckBoxMenuItem( _show_domain_labels, getOptions().isShowDomainLabels() );
\r
1055 customizeCheckBoxMenuItem( _show_annotation_ref_source, getOptions().isShowAnnotationRefSource() );
\r
1056 customizeCheckBoxMenuItem( _abbreviate_scientific_names, getOptions().isAbbreviateScientificTaxonNames() );
\r
1057 customizeCheckBoxMenuItem( _search_case_senstive_cbmi, getOptions().isSearchCaseSensitive() );
\r
1058 customizeCheckBoxMenuItem( _show_scale_cbmi, getOptions().isShowScale() );
\r
1059 customizeRadioButtonMenuItem( _non_lined_up_cladograms_rbmi,
\r
1060 getOptions().getCladogramType() == CLADOGRAM_TYPE.NON_LINED_UP );
\r
1061 customizeRadioButtonMenuItem( _uniform_cladograms_rbmi,
\r
1062 getOptions().getCladogramType() == CLADOGRAM_TYPE.TOTAL_NODE_SUM_DEP );
\r
1063 customizeRadioButtonMenuItem( _ext_node_dependent_cladogram_rbmi,
\r
1064 getOptions().getCladogramType() == CLADOGRAM_TYPE.EXT_NODE_SUM_DEP );
\r
1065 customizeCheckBoxMenuItem( _show_overview_cbmi, getOptions().isShowOverview() );
\r
1066 customizeCheckBoxMenuItem( _label_direction_cbmi,
\r
1067 getOptions().getNodeLabelDirection() == NODE_LABEL_DIRECTION.RADIAL );
\r
1068 customizeCheckBoxMenuItem( _antialias_print_cbmi, getOptions().isAntialiasPrint() );
\r
1069 customizeCheckBoxMenuItem( _print_black_and_white_cbmi, getOptions().isPrintBlackAndWhite() );
\r
1070 customizeCheckBoxMenuItem( _internal_number_are_confidence_for_nh_parsing_cbmi, getOptions()
\r
1071 .isInternalNumberAreConfidenceForNhParsing() );
\r
1072 customizeRadioButtonMenuItem( _extract_taxonomy_no_rbmi,
\r
1073 getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.NO );
\r
1074 customizeRadioButtonMenuItem( _extract_taxonomy_pfam_strict_rbmi,
\r
1075 getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.PFAM_STYLE_STRICT );
\r
1076 customizeRadioButtonMenuItem( _extract_taxonomy_pfam_relaxed_rbmi,
\r
1077 getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.PFAM_STYLE_RELAXED );
\r
1078 customizeRadioButtonMenuItem( _extract_taxonomy_agressive_rbmi,
\r
1079 getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.AGGRESSIVE );
\r
1080 customizeCheckBoxMenuItem( _replace_underscores_cbmi, getOptions().isReplaceUnderscoresInNhParsing() );
\r
1081 customizeCheckBoxMenuItem( _allow_errors_in_distance_to_parent_cbmi, getOptions()
\r
1082 .isReplaceUnderscoresInNhParsing() );
\r
1083 customizeCheckBoxMenuItem( _search_with_regex_cbmi, getOptions().isSearchWithRegex() );
\r
1084 customizeCheckBoxMenuItem( _search_whole_words_only_cbmi, getOptions().isMatchWholeTermsOnly() );
\r
1085 customizeCheckBoxMenuItem( _inverse_search_result_cbmi, getOptions().isInverseSearchResult() );
\r
1086 customizeCheckBoxMenuItem( _graphics_export_visible_only_cbmi, getOptions().isGraphicsExportVisibleOnly() );
\r
1087 customizeCheckBoxMenuItem( _print_using_actual_size_cbmi, getOptions().isPrintUsingActualSize() );
\r
1088 customizeCheckBoxMenuItem( _graphics_export_using_actual_size_cbmi, getOptions()
\r
1089 .isGraphicsExportUsingActualSize() );
\r
1090 customizeCheckBoxMenuItem( _show_confidence_stddev_cbmi, getOptions().isShowConfidenceStddev() );
\r
1091 customizeCheckBoxMenuItem( _use_brackets_for_conf_in_nh_export_cbmi, getOptions()
\r
1092 .getNhConversionSupportValueStyle() == NH_CONVERSION_SUPPORT_VALUE_STYLE.IN_SQUARE_BRACKETS );
\r
1093 customizeCheckBoxMenuItem( _use_internal_names_for_conf_in_nh_export_cbmi, getOptions()
\r
1094 .getNhConversionSupportValueStyle() == NH_CONVERSION_SUPPORT_VALUE_STYLE.AS_INTERNAL_NODE_NAMES );
\r
1095 customizeCheckBoxMenuItem( _line_up_renderable_data_cbmi, getOptions().isLineUpRendarableNodeData() );
\r
1096 customizeCheckBoxMenuItem( _right_line_up_domains_cbmi, getOptions().isRightLineUpDomains() );
\r
1097 _jmenubar.add( _options_jmenu );
\r
1100 void buildPhylogeneticInferenceMenu() {
\r
1101 final InferenceManager im = getInferenceManager();
\r
1102 _inference_menu = MainFrame.createMenu( "Inference", getConfiguration() );
\r
1103 _inference_menu.add( _inference_from_msa_item = new JMenuItem( "From Multiple Sequence Alignment..." ) );
\r
1104 customizeJMenuItem( _inference_from_msa_item );
\r
1105 _inference_from_msa_item.setToolTipText( "Basic phylogenetic inference from MSA" );
\r
1106 if ( im.canDoMsa() ) {
\r
1107 _inference_menu.add( _inference_from_seqs_item = new JMenuItem( "From Unaligned Sequences..." ) );
\r
1108 customizeJMenuItem( _inference_from_seqs_item );
\r
1109 _inference_from_seqs_item
\r
1110 .setToolTipText( "Basic phylogenetic inference including multiple sequence alignment" );
\r
1114 .add( _inference_from_seqs_item = new JMenuItem( "From Unaligned Sequences (no program found)" ) );
\r
1115 customizeJMenuItem( _inference_from_seqs_item );
\r
1116 _inference_from_seqs_item.setEnabled( false );
\r
1118 _jmenubar.add( _inference_menu );
\r
1121 void buildToolsMenu() {
\r
1122 _tools_menu = createMenu( "Tools", getConfiguration() );
\r
1123 _tools_menu.add( _confcolor_item = new JMenuItem( "Colorize Branches Depending on Confidence" ) );
\r
1124 customizeJMenuItem( _confcolor_item );
\r
1125 _tools_menu.add( _color_rank_jmi = new JMenuItem( "Colorize Subtrees via Taxonomic Rank" ) );
\r
1126 customizeJMenuItem( _color_rank_jmi );
\r
1127 _color_rank_jmi.setToolTipText( "for example, at \"Class\" level, colorize mammal specific subtree red" );
\r
1128 _tools_menu.add( _taxcolor_item = new JMenuItem( "Taxonomy Colorize Branches" ) );
\r
1129 customizeJMenuItem( _taxcolor_item );
\r
1130 _tools_menu.addSeparator();
\r
1131 _tools_menu.add( _remove_visual_styles_item = new JMenuItem( "Delete All Visual Styles From Nodes" ) );
\r
1132 _remove_visual_styles_item
\r
1133 .setToolTipText( "To remove all node visual styles (fonts, colors) from the current phylogeny" );
\r
1134 customizeJMenuItem( _remove_visual_styles_item );
\r
1135 _tools_menu.add( _remove_branch_color_item = new JMenuItem( "Delete All Colors From Branches" ) );
\r
1136 _remove_branch_color_item.setToolTipText( "To remove all branch color values from the current phylogeny" );
\r
1137 customizeJMenuItem( _remove_branch_color_item );
\r
1138 _tools_menu.addSeparator();
\r
1139 _tools_menu.add( _annotate_item = new JMenuItem( "Annotate Sequences of Selected Nodes" ) );
\r
1140 customizeJMenuItem( _annotate_item );
\r
1141 _tools_menu.addSeparator();
\r
1142 _tools_menu.add( _midpoint_root_item = new JMenuItem( "Midpoint-Root" ) );
\r
1143 customizeJMenuItem( _midpoint_root_item );
\r
1144 _tools_menu.addSeparator();
\r
1145 _tools_menu.add( _delete_selected_nodes_item = new JMenuItem( "Delete Selected Nodes" ) );
\r
1146 _delete_selected_nodes_item.setToolTipText( "To delete all selected external nodes" );
\r
1147 customizeJMenuItem( _delete_selected_nodes_item );
\r
1148 _tools_menu.add( _delete_not_selected_nodes_item = new JMenuItem( "Retain Selected Nodes" ) );
\r
1149 _delete_not_selected_nodes_item.setToolTipText( "To delete all not selected external nodes" );
\r
1150 customizeJMenuItem( _delete_not_selected_nodes_item );
\r
1151 _tools_menu.addSeparator();
\r
1152 _tools_menu.add( _collapse_species_specific_subtrees = new JMenuItem( "Collapse Species-Specific Subtrees" ) );
\r
1153 customizeJMenuItem( _collapse_species_specific_subtrees );
\r
1154 _collapse_species_specific_subtrees.setToolTipText( "To (reversibly) collapse species-specific subtrees" );
\r
1156 .add( _collapse_below_threshold = new JMenuItem( "Collapse Branches with Confidence Below Threshold into Multifurcations" ) );
\r
1157 customizeJMenuItem( _collapse_below_threshold );
\r
1158 _collapse_below_threshold
\r
1159 .setToolTipText( "To (permanently) collapse branches with confidence values below a threshold into multifurcations (in the case of multiple confidences per branch: without at least one confidence value above a threshold)" );
\r
1162 .add( _collapse_below_branch_length = new JMenuItem( "Collapse Branches with Branch Lengths Below Threshold into Multifurcations" ) );
\r
1163 customizeJMenuItem( _collapse_below_branch_length );
\r
1164 _collapse_below_branch_length
\r
1165 .setToolTipText( "To (permanently) collapse branches with branches with branch lengths below a threshold into multifurcations" );
\r
1167 _tools_menu.addSeparator();
\r
1169 .add( _extract_tax_code_from_node_names_jmi = new JMenuItem( "Extract Taxonomic Data from Node Names" ) );
\r
1170 customizeJMenuItem( _extract_tax_code_from_node_names_jmi );
\r
1171 _extract_tax_code_from_node_names_jmi
\r
1172 .setToolTipText( "To extract SwissProt/Uniprot taxonomic codes (mnemonics) from nodes names in the form of 'xyz_CAEEL', Uniprot/NCBI identifiers form of 'xyz_6239', or scientific names form of 'xyz_Caenorhabditis_elegans'" );
\r
1174 .add( _move_node_names_to_tax_sn_jmi = new JMenuItem( "Transfer Node Names to Taxonomic Scientific Names" ) );
\r
1175 customizeJMenuItem( _move_node_names_to_tax_sn_jmi );
\r
1176 _move_node_names_to_tax_sn_jmi.setToolTipText( "To interpret node names as taxonomic scientific names" );
\r
1177 _tools_menu.add( _move_node_names_to_seq_names_jmi = new JMenuItem( "Transfer Node Names to Sequence Names" ) );
\r
1178 customizeJMenuItem( _move_node_names_to_seq_names_jmi );
\r
1179 _move_node_names_to_seq_names_jmi.setToolTipText( "To interpret node names as sequence (protein, gene) names" );
\r
1180 _tools_menu.addSeparator();
\r
1181 _tools_menu.add( _obtain_seq_information_jmi = new JMenuItem( "Obtain Sequence Information" ) );
\r
1182 customizeJMenuItem( _obtain_seq_information_jmi );
\r
1183 _obtain_seq_information_jmi.setToolTipText( "To add additional sequence information" );
\r
1185 .add( _obtain_detailed_taxonomic_information_jmi = new JMenuItem( OBTAIN_DETAILED_TAXONOMIC_INFORMATION ) );
\r
1186 customizeJMenuItem( _obtain_detailed_taxonomic_information_jmi );
\r
1187 _obtain_detailed_taxonomic_information_jmi
\r
1188 .setToolTipText( "To add additional taxonomic information (from UniProt Taxonomy)" );
\r
1190 .add( _obtain_detailed_taxonomic_information_deleting_jmi = new JMenuItem( "Obtain Detailed Taxonomic Information (deletes nodes!)" ) );
\r
1191 customizeJMenuItem( _obtain_detailed_taxonomic_information_deleting_jmi );
\r
1192 _obtain_detailed_taxonomic_information_deleting_jmi
\r
1193 .setToolTipText( "To add additional taxonomic information, deletes nodes for which taxonomy cannot found (from UniProt Taxonomy)" );
\r
1194 _tools_menu.addSeparator();
\r
1195 _tools_menu.add( _read_values_jmi = new JMenuItem( "Attach Vector/Expression Values" ) );
\r
1196 customizeJMenuItem( _read_values_jmi );
\r
1197 _read_values_jmi.setToolTipText( "To attach vector (e.g. gene expression) values to tree nodes (beta)" );
\r
1198 _jmenubar.add( _tools_menu );
\r
1199 _tools_menu.add( _read_seqs_jmi = new JMenuItem( "Attach Molecular Sequences" ) );
\r
1200 customizeJMenuItem( _read_seqs_jmi );
\r
1202 .setToolTipText( "To attach molecular sequences to tree nodes (from Fasta-formatted file) (beta)" );
\r
1203 _jmenubar.add( _tools_menu );
\r
1208 if ( isUnsavedDataPresent() ) {
\r
1209 final int r = JOptionPane.showConfirmDialog( this,
\r
1210 "Exit despite potentially unsaved changes?",
\r
1212 JOptionPane.YES_NO_OPTION );
\r
1213 if ( r != JOptionPane.YES_OPTION ) {
\r
1220 private void closeCurrentPane() {
\r
1221 if ( getMainPanel().getCurrentTreePanel() != null ) {
\r
1222 if ( getMainPanel().getCurrentTreePanel().isEdited() ) {
\r
1223 final int r = JOptionPane.showConfirmDialog( this,
\r
1224 "Close tab despite potentially unsaved changes?",
\r
1226 JOptionPane.YES_NO_OPTION );
\r
1227 if ( r != JOptionPane.YES_OPTION ) {
\r
1231 getMainPanel().closeCurrentPane();
\r
1232 activateSaveAllIfNeeded();
\r
1236 private void collapse( final Phylogeny phy ) {
\r
1237 final PhylogenyNodeIterator it = phy.iteratorPostorder();
\r
1238 final List<PhylogenyNode> to_be_removed = new ArrayList<PhylogenyNode>();
\r
1239 double min_support = Double.MAX_VALUE;
\r
1240 boolean conf_present = false;
\r
1241 while ( it.hasNext() ) {
\r
1242 final PhylogenyNode n = it.next();
\r
1243 if ( !n.isExternal() && !n.isRoot() ) {
\r
1244 final List<Confidence> c = n.getBranchData().getConfidences();
\r
1245 if ( ( c != null ) && ( c.size() > 0 ) ) {
\r
1246 conf_present = true;
\r
1248 for( final Confidence confidence : c ) {
\r
1249 if ( confidence.getValue() > max ) {
\r
1250 max = confidence.getValue();
\r
1253 if ( max < getMinNotCollapseConfidenceValue() ) {
\r
1254 to_be_removed.add( n );
\r
1256 if ( max < min_support ) {
\r
1257 min_support = max;
\r
1262 if ( conf_present ) {
\r
1263 for( final PhylogenyNode node : to_be_removed ) {
\r
1264 PhylogenyMethods.removeNode( node, phy );
\r
1266 if ( to_be_removed.size() > 0 ) {
\r
1267 phy.externalNodesHaveChanged();
\r
1268 phy.clearHashIdToNodeMap();
\r
1269 phy.recalculateNumberOfExternalDescendants( true );
\r
1270 getCurrentTreePanel().resetNodeIdToDistToLeafMap();
\r
1271 getCurrentTreePanel().updateSetOfCollapsedExternalNodes();
\r
1272 getCurrentTreePanel().calculateLongestExtNodeInfo();
\r
1273 getCurrentTreePanel().setNodeInPreorderToNull();
\r
1274 getCurrentTreePanel().recalculateMaxDistanceToRoot();
\r
1275 getCurrentTreePanel().resetPreferredSize();
\r
1276 getCurrentTreePanel().setEdited( true );
\r
1277 getCurrentTreePanel().repaint();
\r
1280 if ( to_be_removed.size() > 0 ) {
\r
1281 JOptionPane.showMessageDialog( this, "Collapsed " + to_be_removed.size()
\r
1282 + " branches with\nconfidence values below " + getMinNotCollapseConfidenceValue(), "Collapsed "
\r
1283 + to_be_removed.size() + " branches", JOptionPane.INFORMATION_MESSAGE );
\r
1286 JOptionPane.showMessageDialog( this, "No branch collapsed,\nminimum confidence value per branch is "
\r
1287 + min_support, "No branch collapsed", JOptionPane.INFORMATION_MESSAGE );
\r
1291 JOptionPane.showMessageDialog( this,
\r
1292 "No branch collapsed because no confidence values present",
\r
1293 "No confidence values present",
\r
1294 JOptionPane.INFORMATION_MESSAGE );
\r
1298 private void collapseBelowBranchLengthThreshold() {
\r
1299 if ( getCurrentTreePanel() != null ) {
\r
1300 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1301 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1302 final String s = ( String ) JOptionPane
\r
1303 .showInputDialog( this,
\r
1304 "Please enter the minimum branch length value\n",
\r
1305 "Minimal Branch Length Value",
\r
1306 JOptionPane.QUESTION_MESSAGE,
\r
1309 getMinNotCollapseBlValue() );
\r
1310 if ( !ForesterUtil.isEmpty( s ) ) {
\r
1311 boolean success = true;
\r
1313 final String m_str = s.trim();
\r
1314 if ( !ForesterUtil.isEmpty( m_str ) ) {
\r
1316 m = Double.parseDouble( m_str );
\r
1318 catch ( final Exception ex ) {
\r
1325 if ( success && ( m >= 0.0 ) ) {
\r
1326 setMinNotCollapseBlValue( m );
\r
1327 collapseBl( phy );
\r
1334 private void collapseBelowThreshold() {
\r
1335 if ( getCurrentTreePanel() != null ) {
\r
1336 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1337 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1338 final String s = ( String ) JOptionPane.showInputDialog( this,
\r
1339 "Please enter the minimum confidence value\n",
\r
1340 "Minimal Confidence Value",
\r
1341 JOptionPane.QUESTION_MESSAGE,
\r
1344 getMinNotCollapseConfidenceValue() );
\r
1345 if ( !ForesterUtil.isEmpty( s ) ) {
\r
1346 boolean success = true;
\r
1348 final String m_str = s.trim();
\r
1349 if ( !ForesterUtil.isEmpty( m_str ) ) {
\r
1351 m = Double.parseDouble( m_str );
\r
1353 catch ( final Exception ex ) {
\r
1360 if ( success && ( m >= 0.0 ) ) {
\r
1361 setMinNotCollapseConfidenceValue( m );
\r
1369 private void collapseBl( final Phylogeny phy ) {
\r
1370 final PhylogenyNodeIterator it = phy.iteratorPostorder();
\r
1371 final List<PhylogenyNode> to_be_removed = new ArrayList<PhylogenyNode>();
\r
1372 double min_bl = Double.MAX_VALUE;
\r
1373 boolean bl_present = false;
\r
1374 while ( it.hasNext() ) {
\r
1375 final PhylogenyNode n = it.next();
\r
1376 if ( !n.isExternal() && !n.isRoot() ) {
\r
1377 final double bl = n.getDistanceToParent();
\r
1378 if ( bl != PhylogenyDataUtil.BRANCH_LENGTH_DEFAULT ) {
\r
1379 bl_present = true;
\r
1380 if ( bl < getMinNotCollapseBlValue() ) {
\r
1381 to_be_removed.add( n );
\r
1383 if ( bl < min_bl ) {
\r
1389 if ( bl_present ) {
\r
1390 for( final PhylogenyNode node : to_be_removed ) {
\r
1391 PhylogenyMethods.removeNode( node, phy );
\r
1393 if ( to_be_removed.size() > 0 ) {
\r
1394 phy.externalNodesHaveChanged();
\r
1395 phy.clearHashIdToNodeMap();
\r
1396 phy.recalculateNumberOfExternalDescendants( true );
\r
1397 getCurrentTreePanel().resetNodeIdToDistToLeafMap();
\r
1398 getCurrentTreePanel().updateSetOfCollapsedExternalNodes();
\r
1399 getCurrentTreePanel().calculateLongestExtNodeInfo();
\r
1400 getCurrentTreePanel().setNodeInPreorderToNull();
\r
1401 getCurrentTreePanel().recalculateMaxDistanceToRoot();
\r
1402 getCurrentTreePanel().resetPreferredSize();
\r
1403 getCurrentTreePanel().setEdited( true );
\r
1404 getCurrentTreePanel().repaint();
\r
1407 if ( to_be_removed.size() > 0 ) {
\r
1408 JOptionPane.showMessageDialog( this, "Collapsed " + to_be_removed.size()
\r
1409 + " branches with\nbranch length values below " + getMinNotCollapseBlValue(), "Collapsed "
\r
1410 + to_be_removed.size() + " branches", JOptionPane.INFORMATION_MESSAGE );
\r
1413 JOptionPane.showMessageDialog( this,
\r
1414 "No branch collapsed,\nminimum branch length is " + min_bl,
\r
1415 "No branch collapsed",
\r
1416 JOptionPane.INFORMATION_MESSAGE );
\r
1420 JOptionPane.showMessageDialog( this,
\r
1421 "No branch collapsed because no branch length values present",
\r
1422 "No branch length values present",
\r
1423 JOptionPane.INFORMATION_MESSAGE );
\r
1427 private PhyloXmlParser createPhyloXmlParser() {
\r
1428 PhyloXmlParser xml_parser = null;
\r
1429 if ( getConfiguration().isValidatePhyloXmlAgainstSchema() ) {
\r
1431 xml_parser = PhyloXmlParser.createPhyloXmlParserXsdValidating();
\r
1433 catch ( final Exception e ) {
\r
1434 JOptionPane.showMessageDialog( this,
\r
1435 e.getLocalizedMessage(),
\r
1436 "failed to create validating XML parser",
\r
1437 JOptionPane.WARNING_MESSAGE );
\r
1440 if ( xml_parser == null ) {
\r
1441 xml_parser = PhyloXmlParser.createPhyloXmlParser();
\r
1443 return xml_parser;
\r
1446 public void end() {
\r
1447 _mainpanel.terminate();
\r
1448 _contentpane.removeAll();
\r
1449 setVisible( false );
\r
1453 void executeLineageInference() {
\r
1454 if ( ( _mainpanel.getCurrentPhylogeny() == null ) || ( _mainpanel.getCurrentPhylogeny().isEmpty() ) ) {
\r
1457 if ( !_mainpanel.getCurrentPhylogeny().isRooted() ) {
\r
1458 JOptionPane.showMessageDialog( this,
\r
1459 "Phylogeny is not rooted.",
\r
1460 "Cannot infer ancestral taxonomies",
\r
1461 JOptionPane.ERROR_MESSAGE );
\r
1464 final AncestralTaxonomyInferrer inferrer = new AncestralTaxonomyInferrer( this,
\r
1465 _mainpanel.getCurrentTreePanel(),
\r
1466 _mainpanel.getCurrentPhylogeny()
\r
1468 new Thread( inferrer ).start();
\r
1471 private void executePhyleneticInference( final boolean from_unaligned_seqs ) {
\r
1472 final PhyloInferenceDialog dialog = new PhyloInferenceDialog( this,
\r
1473 getPhylogeneticInferenceOptions(),
\r
1474 from_unaligned_seqs );
\r
1475 dialog.activate();
\r
1476 if ( dialog.getValue() == JOptionPane.OK_OPTION ) {
\r
1477 if ( !from_unaligned_seqs ) {
\r
1478 if ( getMsa() != null ) {
\r
1479 final PhylogeneticInferrer inferrer = new PhylogeneticInferrer( getMsa(),
\r
1480 getPhylogeneticInferenceOptions()
\r
1482 new Thread( inferrer ).start();
\r
1485 JOptionPane.showMessageDialog( this,
\r
1486 "No multiple sequence alignment selected",
\r
1487 "Phylogenetic Inference Not Launched",
\r
1488 JOptionPane.WARNING_MESSAGE );
\r
1492 if ( getSeqs() != null ) {
\r
1493 final PhylogeneticInferrer inferrer = new PhylogeneticInferrer( getSeqs(),
\r
1494 getPhylogeneticInferenceOptions()
\r
1496 new Thread( inferrer ).start();
\r
1499 JOptionPane.showMessageDialog( this,
\r
1500 "No input sequences selected",
\r
1501 "Phylogenetic Inference Not Launched",
\r
1502 JOptionPane.WARNING_MESSAGE );
\r
1509 removeAllTextFrames();
\r
1510 _mainpanel.terminate();
\r
1511 _contentpane.removeAll();
\r
1512 setVisible( false );
\r
1514 // System.exit( 0 ); //TODO reconfirm that this is OK, then remove.
\r
1517 private void extractTaxDataFromNodeNames() throws PhyloXmlDataFormatException {
\r
1518 final StringBuilder sb = new StringBuilder();
\r
1519 final StringBuilder sb_failed = new StringBuilder();
\r
1521 int counter_failed = 0;
\r
1522 if ( getCurrentTreePanel() != null ) {
\r
1523 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1524 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1525 final PhylogenyNodeIterator it = phy.iteratorExternalForward();
\r
1526 while ( it.hasNext() ) {
\r
1527 final PhylogenyNode n = it.next();
\r
1528 final String name = n.getName().trim();
\r
1529 if ( !ForesterUtil.isEmpty( name ) ) {
\r
1530 final String nt = ParserUtils.extractTaxonomyDataFromNodeName( n,
\r
1531 TAXONOMY_EXTRACTION.AGGRESSIVE );
\r
1532 if ( !ForesterUtil.isEmpty( nt ) ) {
\r
1533 if ( counter < 15 ) {
\r
1534 sb.append( name + ": " + nt + "\n" );
\r
1536 else if ( counter == 15 ) {
\r
1537 sb.append( "...\n" );
\r
1542 if ( counter_failed < 15 ) {
\r
1543 sb_failed.append( name + "\n" );
\r
1545 else if ( counter_failed == 15 ) {
\r
1546 sb_failed.append( "...\n" );
\r
1552 if ( counter > 0 ) {
\r
1553 String failed = "";
\r
1554 String all = "all ";
\r
1555 if ( counter_failed > 0 ) {
\r
1557 failed = "\nCould not extract taxonomic data for " + counter_failed
\r
1558 + " named external nodes:\n" + sb_failed;
\r
1560 JOptionPane.showMessageDialog( this,
\r
1561 "Extracted taxonomic data from " + all + counter
\r
1562 + " named external nodes:\n" + sb.toString() + failed,
\r
1563 "Taxonomic Data Extraction Completed",
\r
1564 counter_failed > 0 ? JOptionPane.WARNING_MESSAGE
\r
1565 : JOptionPane.INFORMATION_MESSAGE );
\r
1569 .showMessageDialog( this,
\r
1570 "Could not extract any taxonomic data.\nMaybe node names are empty\n"
\r
1571 + "or not in the forms \"XYZ_CAEEL\", \"XYZ_6239\", or \"XYZ_Caenorhabditis_elegans\"\n"
\r
1572 + "or nodes already have taxonomic data?\n",
\r
1573 "No Taxonomic Data Extracted",
\r
1574 JOptionPane.ERROR_MESSAGE );
\r
1581 public MainPanel getMainPanel() {
\r
1582 return _mainpanel;
\r
1585 private double getMinNotCollapseBlValue() {
\r
1586 return _min_not_collapse_bl;
\r
1589 private double getMinNotCollapseConfidenceValue() {
\r
1590 return _min_not_collapse;
\r
1593 public Msa getMsa() {
\r
1597 public File getMsaFile() {
\r
1601 private PhylogeneticInferenceOptions getPhylogeneticInferenceOptions() {
\r
1602 if ( _phylogenetic_inference_options == null ) {
\r
1603 _phylogenetic_inference_options = new PhylogeneticInferenceOptions();
\r
1605 return _phylogenetic_inference_options;
\r
1608 public List<MolecularSequence> getSeqs() {
\r
1612 public File getSeqsFile() {
\r
1613 return _seqs_file;
\r
1616 private boolean isUnsavedDataPresent() {
\r
1617 final List<TreePanel> tps = getMainPanel().getTreePanels();
\r
1618 for( final TreePanel tp : tps ) {
\r
1619 if ( tp.isEdited() ) {
\r
1626 private void moveNodeNamesToSeqNames() throws PhyloXmlDataFormatException {
\r
1627 if ( getCurrentTreePanel() != null ) {
\r
1628 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1629 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1631 .transferNodeNameToField( phy, PhylogenyMethods.PhylogenyNodeField.SEQUENCE_NAME, false );
\r
1636 private void moveNodeNamesToTaxSn() throws PhyloXmlDataFormatException {
\r
1637 if ( getCurrentTreePanel() != null ) {
\r
1638 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1639 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1640 PhylogenyMethods.transferNodeNameToField( phy,
\r
1641 PhylogenyMethods.PhylogenyNodeField.TAXONOMY_SCIENTIFIC_NAME,
\r
1647 private void newTree() {
\r
1648 final Phylogeny[] phys = new Phylogeny[ 1 ];
\r
1649 final Phylogeny phy = new Phylogeny();
\r
1650 final PhylogenyNode node = new PhylogenyNode();
\r
1651 phy.setRoot( node );
\r
1652 phy.setRooted( true );
\r
1654 AptxUtil.addPhylogeniesToTabs( phys, "", "", getConfiguration(), getMainPanel() );
\r
1655 _mainpanel.getControlPanel().showWhole();
\r
1656 _mainpanel.getCurrentTreePanel().setPhylogenyGraphicsType( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );
\r
1657 _mainpanel.getOptions().setPhylogenyGraphicsType( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );
\r
1658 if ( getMainPanel().getMainFrame() == null ) {
\r
1659 // Must be "E" applet version.
\r
1660 ( ( ArchaeopteryxE ) ( ( MainPanelApplets ) getMainPanel() ).getApplet() )
\r
1661 .setSelectedTypeInTypeMenu( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );
\r
1664 getMainPanel().getMainFrame().setSelectedTypeInTypeMenu( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );
\r
1666 activateSaveAllIfNeeded();
\r
1670 private void obtainDetailedTaxonomicInformation() {
\r
1671 if ( getCurrentTreePanel() != null ) {
\r
1672 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1673 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1674 final TaxonomyDataManager t = new TaxonomyDataManager( this,
\r
1675 _mainpanel.getCurrentTreePanel(),
\r
1679 new Thread( t ).start();
\r
1684 private void obtainDetailedTaxonomicInformationDelete() {
\r
1685 if ( getCurrentTreePanel() != null ) {
\r
1686 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1687 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1688 final TaxonomyDataManager t = new TaxonomyDataManager( this,
\r
1689 _mainpanel.getCurrentTreePanel(),
\r
1693 new Thread( t ).start();
\r
1698 private void obtainSequenceInformation() {
\r
1699 if ( getCurrentTreePanel() != null ) {
\r
1700 final Phylogeny phy = getCurrentTreePanel().getPhylogeny();
\r
1701 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1702 final SequenceDataRetriver u = new SequenceDataRetriver( this,
\r
1703 _mainpanel.getCurrentTreePanel(),
\r
1705 new Thread( u ).start();
\r
1710 private void preProcessTreesUponReading( final Phylogeny[] phys ) {
\r
1711 for( final Phylogeny phy : phys ) {
\r
1712 if ( ( phy != null ) && !phy.isEmpty() ) {
\r
1713 for( final PhylogenyNodeIterator it = phy.iteratorPreorder(); it.hasNext(); ) {
\r
1714 final PhylogenyNode n = it.next();
\r
1715 if ( n.isExternal() ) {
\r
1716 if ( n.getNodeData().isHasSequence() ) {
\r
1717 final Sequence s = n.getNodeData().getSequence();
\r
1718 if ( ForesterUtil.isEmpty( s.getGeneName() ) || s.getGeneName().startsWith( "LOC" ) ) {
\r
1719 if ( ( s.getAccession() != null )
\r
1720 && !ForesterUtil.isEmpty( s.getAccession().getValue() ) ) {
\r
1721 s.setGeneName( s.getAccession().getValue() );
\r
1723 else if ( !ForesterUtil.isEmpty( n.getName() ) ) {
\r
1724 s.setGeneName( n.getName() );
\r
1734 private void print() {
\r
1735 if ( ( getCurrentTreePanel() == null ) || ( getCurrentTreePanel().getPhylogeny() == null )
\r
1736 || getCurrentTreePanel().getPhylogeny().isEmpty() ) {
\r
1739 if ( !getOptions().isPrintUsingActualSize() ) {
\r
1740 getCurrentTreePanel().calcParametersForPainting( getOptions().getPrintSizeX() - 80,
\r
1741 getOptions().getPrintSizeY() - 140 );
\r
1742 getCurrentTreePanel().resetPreferredSize();
\r
1743 getCurrentTreePanel().repaint();
\r
1745 final String job_name = Constants.PRG_NAME;
\r
1746 boolean error = false;
\r
1747 String printer_name = null;
\r
1749 printer_name = Printer.print( getCurrentTreePanel(), job_name );
\r
1751 catch ( final Exception e ) {
\r
1753 JOptionPane.showMessageDialog( this, e.getMessage(), "Printing Error", JOptionPane.ERROR_MESSAGE );
\r
1755 if ( !error && ( printer_name != null ) ) {
\r
1756 String msg = "Printing data sent to printer";
\r
1757 if ( printer_name.length() > 1 ) {
\r
1758 msg += " [" + printer_name + "]";
\r
1760 JOptionPane.showMessageDialog( this, msg, "Printing...", JOptionPane.INFORMATION_MESSAGE );
\r
1762 if ( !getOptions().isPrintUsingActualSize() ) {
\r
1763 getControlPanel().showWhole();
\r
1767 public void readMsaFromFile() {
\r
1768 // Set an initial directory if none set yet
\r
1769 final File my_dir = getCurrentDir();
\r
1770 _msa_filechooser.setMultiSelectionEnabled( false );
\r
1771 // Open file-open dialog and set current directory
\r
1772 if ( my_dir != null ) {
\r
1773 _msa_filechooser.setCurrentDirectory( my_dir );
\r
1775 final int result = _msa_filechooser.showOpenDialog( _contentpane );
\r
1776 // All done: get the msa
\r
1777 final File file = _msa_filechooser.getSelectedFile();
\r
1778 setCurrentDir( _msa_filechooser.getCurrentDirectory() );
\r
1779 if ( ( file != null ) && !file.isDirectory() && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
1780 setMsaFile( null );
\r
1784 final InputStream is = new FileInputStream( file );
\r
1785 if ( FastaParser.isLikelyFasta( file ) ) {
\r
1786 msa = FastaParser.parseMsa( is );
\r
1789 msa = GeneralMsaParser.parse( is );
\r
1792 catch ( final MsaFormatException e ) {
\r
1794 JOptionPane.showMessageDialog( this,
\r
1795 e.getLocalizedMessage(),
\r
1796 "Multiple sequence alignment format error",
\r
1797 JOptionPane.ERROR_MESSAGE );
\r
1800 catch ( final IOException e ) {
\r
1802 JOptionPane.showMessageDialog( this,
\r
1803 e.getLocalizedMessage(),
\r
1804 "Failed to read multiple sequence alignment",
\r
1805 JOptionPane.ERROR_MESSAGE );
\r
1808 catch ( final IllegalArgumentException e ) {
\r
1810 JOptionPane.showMessageDialog( this,
\r
1811 e.getLocalizedMessage(),
\r
1812 "Unexpected error during reading of multiple sequence alignment",
\r
1813 JOptionPane.ERROR_MESSAGE );
\r
1816 catch ( final Exception e ) {
\r
1818 e.printStackTrace();
\r
1819 JOptionPane.showMessageDialog( this,
\r
1820 e.getLocalizedMessage(),
\r
1821 "Unexpected error during reading of multiple sequence alignment",
\r
1822 JOptionPane.ERROR_MESSAGE );
\r
1825 if ( ( msa == null ) || ( msa.getNumberOfSequences() < 1 ) ) {
\r
1826 JOptionPane.showMessageDialog( this,
\r
1827 "Multiple sequence alignment is empty",
\r
1828 "Illegal Multiple Sequence Alignment",
\r
1829 JOptionPane.ERROR_MESSAGE );
\r
1832 if ( msa.getNumberOfSequences() < 4 ) {
\r
1833 JOptionPane.showMessageDialog( this,
\r
1834 "Multiple sequence alignment needs to contain at least 3 sequences",
\r
1835 "Illegal multiple sequence alignment",
\r
1836 JOptionPane.ERROR_MESSAGE );
\r
1839 if ( msa.getLength() < 2 ) {
\r
1840 JOptionPane.showMessageDialog( this,
\r
1841 "Multiple sequence alignment needs to contain at least 2 residues",
\r
1842 "Illegal multiple sequence alignment",
\r
1843 JOptionPane.ERROR_MESSAGE );
\r
1847 setMsaFile( _msa_filechooser.getSelectedFile() );
\r
1852 private void readPhylogeniesFromFile() {
\r
1853 boolean exception = false;
\r
1854 Phylogeny[] phys = null;
\r
1855 // Set an initial directory if none set yet
\r
1856 final File my_dir = getCurrentDir();
\r
1857 _open_filechooser.setMultiSelectionEnabled( true );
\r
1858 // Open file-open dialog and set current directory
\r
1859 if ( my_dir != null ) {
\r
1860 _open_filechooser.setCurrentDirectory( my_dir );
\r
1862 final int result = _open_filechooser.showOpenDialog( _contentpane );
\r
1863 // All done: get the file
\r
1864 final File[] files = _open_filechooser.getSelectedFiles();
\r
1865 setCurrentDir( _open_filechooser.getCurrentDirectory() );
\r
1866 boolean nhx_or_nexus = false;
\r
1867 if ( ( files != null ) && ( files.length > 0 ) && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
1868 for( final File file : files ) {
\r
1869 if ( ( file != null ) && !file.isDirectory() ) {
\r
1870 if ( _mainpanel.getCurrentTreePanel() != null ) {
\r
1871 _mainpanel.getCurrentTreePanel().setWaitCursor();
\r
1874 _mainpanel.setWaitCursor();
\r
1876 if ( ( _open_filechooser.getFileFilter() == MainFrame.nhfilter )
\r
1877 || ( _open_filechooser.getFileFilter() == MainFrame.nhxfilter ) ) {
\r
1879 final NHXParser nhx = new NHXParser();
\r
1880 setSpecialOptionsForNhxParser( nhx );
\r
1881 phys = PhylogenyMethods.readPhylogenies( nhx, file );
\r
1882 nhx_or_nexus = true;
\r
1884 catch ( final Exception e ) {
\r
1886 exceptionOccuredDuringOpenFile( e );
\r
1889 else if ( _open_filechooser.getFileFilter() == MainFrame.xmlfilter ) {
\r
1890 warnIfNotPhyloXmlValidation( getConfiguration() );
\r
1892 final PhyloXmlParser xml_parser = createPhyloXmlParser();
\r
1893 phys = PhylogenyMethods.readPhylogenies( xml_parser, file );
\r
1895 catch ( final Exception e ) {
\r
1897 exceptionOccuredDuringOpenFile( e );
\r
1900 else if ( _open_filechooser.getFileFilter() == MainFrame.tolfilter ) {
\r
1902 phys = PhylogenyMethods.readPhylogenies( new TolParser(), file );
\r
1904 catch ( final Exception e ) {
\r
1906 exceptionOccuredDuringOpenFile( e );
\r
1909 else if ( _open_filechooser.getFileFilter() == MainFrame.nexusfilter ) {
\r
1911 final NexusPhylogeniesParser nex = new NexusPhylogeniesParser();
\r
1912 setSpecialOptionsForNexParser( nex );
\r
1913 phys = PhylogenyMethods.readPhylogenies( nex, file );
\r
1914 nhx_or_nexus = true;
\r
1916 catch ( final Exception e ) {
\r
1918 exceptionOccuredDuringOpenFile( e );
\r
1924 final PhylogenyParser parser = ParserUtils
\r
1925 .createParserDependingOnFileType( file, getConfiguration()
\r
1926 .isValidatePhyloXmlAgainstSchema() );
\r
1927 if ( parser instanceof NexusPhylogeniesParser ) {
\r
1928 final NexusPhylogeniesParser nex = ( NexusPhylogeniesParser ) parser;
\r
1929 setSpecialOptionsForNexParser( nex );
\r
1930 nhx_or_nexus = true;
\r
1932 else if ( parser instanceof NHXParser ) {
\r
1933 final NHXParser nhx = ( NHXParser ) parser;
\r
1934 setSpecialOptionsForNhxParser( nhx );
\r
1935 nhx_or_nexus = true;
\r
1937 else if ( parser instanceof PhyloXmlParser ) {
\r
1938 warnIfNotPhyloXmlValidation( getConfiguration() );
\r
1940 phys = PhylogenyMethods.readPhylogenies( parser, file );
\r
1942 catch ( final Exception e ) {
\r
1944 exceptionOccuredDuringOpenFile( e );
\r
1947 if ( _mainpanel.getCurrentTreePanel() != null ) {
\r
1948 _mainpanel.getCurrentTreePanel().setArrowCursor();
\r
1951 _mainpanel.setArrowCursor();
\r
1953 if ( !exception && ( phys != null ) && ( phys.length > 0 ) ) {
\r
1954 boolean one_desc = false;
\r
1955 if ( nhx_or_nexus ) {
\r
1956 for( final Phylogeny phy : phys ) {
\r
1957 if ( getOptions().isInternalNumberAreConfidenceForNhParsing() ) {
\r
1958 PhylogenyMethods.transferInternalNodeNamesToConfidence( phy, "" );
\r
1960 if ( PhylogenyMethods.getMinimumDescendentsPerInternalNodes( phy ) == 1 ) {
\r
1966 if ( PREPROCESS_TREES ) {
\r
1967 preProcessTreesUponReading( phys );
\r
1969 AptxUtil.addPhylogeniesToTabs( phys,
\r
1971 file.getAbsolutePath(),
\r
1972 getConfiguration(),
\r
1974 _mainpanel.getControlPanel().showWhole();
\r
1975 if ( nhx_or_nexus && one_desc ) {
\r
1977 .showMessageDialog( this,
\r
1978 "One or more trees contain (a) node(s) with one descendant, "
\r
1979 + ForesterUtil.LINE_SEPARATOR
\r
1980 + "possibly indicating illegal parentheses within node names.",
\r
1981 "Warning: Possible Error in New Hampshire Formatted Data",
\r
1982 JOptionPane.WARNING_MESSAGE );
\r
1988 activateSaveAllIfNeeded();
\r
1992 void readPhylogeniesFromURL() {
\r
1994 Phylogeny[] phys = null;
\r
1995 final String message = "Please enter a complete URL, for example \"http://purl.org/phylo/treebase/phylows/study/TB2:S15480?format=nexus\"";
\r
1996 final String url_string = JOptionPane.showInputDialog( this,
\r
1998 "Use URL/webservice to obtain a phylogeny",
\r
1999 JOptionPane.QUESTION_MESSAGE );
\r
2000 boolean nhx_or_nexus = false;
\r
2001 if ( ( url_string != null ) && ( url_string.length() > 0 ) ) {
\r
2003 url = new URL( url_string );
\r
2004 PhylogenyParser parser = null;
\r
2005 if ( url.getHost().toLowerCase().indexOf( "tolweb" ) >= 0 ) {
\r
2006 parser = new TolParser();
\r
2009 parser = ParserUtils.createParserDependingOnUrlContents( url, getConfiguration()
\r
2010 .isValidatePhyloXmlAgainstSchema() );
\r
2012 if ( parser instanceof NexusPhylogeniesParser ) {
\r
2013 nhx_or_nexus = true;
\r
2015 else if ( parser instanceof NHXParser ) {
\r
2016 nhx_or_nexus = true;
\r
2018 if ( _mainpanel.getCurrentTreePanel() != null ) {
\r
2019 _mainpanel.getCurrentTreePanel().setWaitCursor();
\r
2022 _mainpanel.setWaitCursor();
\r
2024 final PhylogenyFactory factory = ParserBasedPhylogenyFactory.getInstance();
\r
2025 phys = factory.create( url.openStream(), parser );
\r
2027 catch ( final MalformedURLException e ) {
\r
2028 JOptionPane.showMessageDialog( this,
\r
2029 "Malformed URL: " + url + "\n" + e.getLocalizedMessage(),
\r
2031 JOptionPane.ERROR_MESSAGE );
\r
2033 catch ( final IOException e ) {
\r
2034 JOptionPane.showMessageDialog( this,
\r
2035 "Could not read from " + url + "\n"
\r
2036 + ForesterUtil.wordWrap( e.getLocalizedMessage(), 80 ),
\r
2037 "Failed to read URL",
\r
2038 JOptionPane.ERROR_MESSAGE );
\r
2040 catch ( final Exception e ) {
\r
2041 JOptionPane.showMessageDialog( this,
\r
2042 ForesterUtil.wordWrap( e.getLocalizedMessage(), 80 ),
\r
2043 "Unexpected Exception",
\r
2044 JOptionPane.ERROR_MESSAGE );
\r
2047 if ( _mainpanel.getCurrentTreePanel() != null ) {
\r
2048 _mainpanel.getCurrentTreePanel().setArrowCursor();
\r
2051 _mainpanel.setArrowCursor();
\r
2054 if ( ( phys != null ) && ( phys.length > 0 ) ) {
\r
2055 if ( nhx_or_nexus && getOptions().isInternalNumberAreConfidenceForNhParsing() ) {
\r
2056 for( final Phylogeny phy : phys ) {
\r
2057 PhylogenyMethods.transferInternalNodeNamesToConfidence( phy, "" );
\r
2060 AptxUtil.addPhylogeniesToTabs( phys,
\r
2061 new File( url.getFile() ).getName(),
\r
2062 new File( url.getFile() ).toString(),
\r
2063 getConfiguration(),
\r
2065 _mainpanel.getControlPanel().showWhole();
\r
2068 activateSaveAllIfNeeded();
\r
2072 public void readSeqsFromFileforPI() {
\r
2073 // Set an initial directory if none set yet
\r
2074 final File my_dir = getCurrentDir();
\r
2075 _seqs_pi_filechooser.setMultiSelectionEnabled( false );
\r
2076 // Open file-open dialog and set current directory
\r
2077 if ( my_dir != null ) {
\r
2078 _seqs_pi_filechooser.setCurrentDirectory( my_dir );
\r
2080 final int result = _seqs_pi_filechooser.showOpenDialog( _contentpane );
\r
2081 // All done: get the seqs
\r
2082 final File file = _seqs_pi_filechooser.getSelectedFile();
\r
2083 setCurrentDir( _seqs_pi_filechooser.getCurrentDirectory() );
\r
2084 if ( ( file != null ) && !file.isDirectory() && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
2085 setSeqsFile( null );
\r
2087 List<MolecularSequence> seqs = null;
\r
2089 if ( FastaParser.isLikelyFasta( new FileInputStream( file ) ) ) {
\r
2090 seqs = FastaParser.parse( new FileInputStream( file ) );
\r
2091 for( final MolecularSequence seq : seqs ) {
\r
2092 System.out.println( SequenceWriter.toFasta( seq, 60 ) );
\r
2099 catch ( final MsaFormatException e ) {
\r
2101 JOptionPane.showMessageDialog( this,
\r
2102 e.getLocalizedMessage(),
\r
2103 "Multiple sequence file format error",
\r
2104 JOptionPane.ERROR_MESSAGE );
\r
2107 catch ( final IOException e ) {
\r
2109 JOptionPane.showMessageDialog( this,
\r
2110 e.getLocalizedMessage(),
\r
2111 "Failed to read multiple sequence file",
\r
2112 JOptionPane.ERROR_MESSAGE );
\r
2115 catch ( final IllegalArgumentException e ) {
\r
2117 JOptionPane.showMessageDialog( this,
\r
2118 e.getLocalizedMessage(),
\r
2119 "Unexpected error during reading of multiple sequence file",
\r
2120 JOptionPane.ERROR_MESSAGE );
\r
2123 catch ( final Exception e ) {
\r
2125 e.printStackTrace();
\r
2126 JOptionPane.showMessageDialog( this,
\r
2127 e.getLocalizedMessage(),
\r
2128 "Unexpected error during reading of multiple sequence file",
\r
2129 JOptionPane.ERROR_MESSAGE );
\r
2132 if ( ( seqs == null ) || ( seqs.size() < 1 ) ) {
\r
2133 JOptionPane.showMessageDialog( this,
\r
2134 "Multiple sequence file is empty",
\r
2135 "Illegal multiple sequence file",
\r
2136 JOptionPane.ERROR_MESSAGE );
\r
2139 if ( seqs.size() < 4 ) {
\r
2140 JOptionPane.showMessageDialog( this,
\r
2141 "Multiple sequence file needs to contain at least 3 sequences",
\r
2142 "Illegal multiple sequence file",
\r
2143 JOptionPane.ERROR_MESSAGE );
\r
2146 // if ( msa.getLength() < 2 ) {
\r
2147 // JOptionPane.showMessageDialog( this,
\r
2148 // "Multiple sequence alignment needs to contain at least 2 residues",
\r
2149 // "Illegal multiple sequence file",
\r
2150 // JOptionPane.ERROR_MESSAGE );
\r
2154 setSeqsFile( _seqs_pi_filechooser.getSelectedFile() );
\r
2159 private void readSpeciesTreeFromFile() {
\r
2160 Phylogeny t = null;
\r
2161 boolean exception = false;
\r
2162 final File my_dir = getCurrentDir();
\r
2163 _open_filechooser_for_species_tree.setSelectedFile( new File( "" ) );
\r
2164 if ( my_dir != null ) {
\r
2165 _open_filechooser_for_species_tree.setCurrentDirectory( my_dir );
\r
2167 final int result = _open_filechooser_for_species_tree.showOpenDialog( _contentpane );
\r
2168 final File file = _open_filechooser_for_species_tree.getSelectedFile();
\r
2169 if ( ( file != null ) && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
2170 if ( _open_filechooser_for_species_tree.getFileFilter() == MainFrame.xmlfilter ) {
\r
2172 final Phylogeny[] trees = PhylogenyMethods.readPhylogenies( PhyloXmlParser
\r
2173 .createPhyloXmlParserXsdValidating(), file );
\r
2176 catch ( final Exception e ) {
\r
2178 exceptionOccuredDuringOpenFile( e );
\r
2181 else if ( _open_filechooser_for_species_tree.getFileFilter() == MainFrame.tolfilter ) {
\r
2183 final Phylogeny[] trees = PhylogenyMethods.readPhylogenies( new TolParser(), file );
\r
2186 catch ( final Exception e ) {
\r
2188 exceptionOccuredDuringOpenFile( e );
\r
2194 final Phylogeny[] trees = PhylogenyMethods.readPhylogenies( PhyloXmlParser
\r
2195 .createPhyloXmlParserXsdValidating(), file );
\r
2198 catch ( final Exception e ) {
\r
2200 exceptionOccuredDuringOpenFile( e );
\r
2203 if ( !exception && ( t != null ) && !t.isRooted() ) {
\r
2206 JOptionPane.showMessageDialog( this,
\r
2207 "Species tree is not rooted",
\r
2208 "Species tree not loaded",
\r
2209 JOptionPane.ERROR_MESSAGE );
\r
2211 if ( !exception && ( t != null ) ) {
\r
2212 final Set<Taxonomy> tax_set = new HashSet<Taxonomy>();
\r
2213 for( final PhylogenyNodeIterator it = t.iteratorExternalForward(); it.hasNext(); ) {
\r
2214 final PhylogenyNode node = it.next();
\r
2215 if ( !node.getNodeData().isHasTaxonomy() ) {
\r
2219 .showMessageDialog( this,
\r
2220 "Species tree contains external node(s) without taxonomy information",
\r
2221 "Species tree not loaded",
\r
2222 JOptionPane.ERROR_MESSAGE );
\r
2226 if ( tax_set.contains( node.getNodeData().getTaxonomy() ) ) {
\r
2229 JOptionPane.showMessageDialog( this,
\r
2231 + node.getNodeData().getTaxonomy().asSimpleText()
\r
2232 + "] is not unique in species tree",
\r
2233 "Species tree not loaded",
\r
2234 JOptionPane.ERROR_MESSAGE );
\r
2238 tax_set.add( node.getNodeData().getTaxonomy() );
\r
2243 if ( !exception && ( t != null ) ) {
\r
2244 setSpeciesTree( t );
\r
2245 JOptionPane.showMessageDialog( this,
\r
2246 "Species tree successfully loaded",
\r
2247 "Species tree loaded",
\r
2248 JOptionPane.INFORMATION_MESSAGE );
\r
2250 _contentpane.repaint();
\r
2255 private void setArrowCursor() {
\r
2257 _mainpanel.getCurrentTreePanel().setArrowCursor();
\r
2259 catch ( final Exception ex ) {
\r
2264 private void setMinNotCollapseBlValue( final double min_not_collapse_bl ) {
\r
2265 _min_not_collapse_bl = min_not_collapse_bl;
\r
2268 private void setMinNotCollapseConfidenceValue( final double min_not_collapse ) {
\r
2269 _min_not_collapse = min_not_collapse;
\r
2272 void setMsa( final Msa msa ) {
\r
2276 void setMsaFile( final File msa_file ) {
\r
2277 _msa_file = msa_file;
\r
2280 private void setPhylogeneticInferenceOptions( final PhylogeneticInferenceOptions phylogenetic_inference_options ) {
\r
2281 _phylogenetic_inference_options = phylogenetic_inference_options;
\r
2284 void setSeqs( final List<MolecularSequence> seqs ) {
\r
2288 void setSeqsFile( final File seqs_file ) {
\r
2289 _seqs_file = seqs_file;
\r
2292 private void setSpecialOptionsForNexParser( final NexusPhylogeniesParser nex ) {
\r
2293 nex.setReplaceUnderscores( getOptions().isReplaceUnderscoresInNhParsing() );
\r
2294 nex.setTaxonomyExtraction( getOptions().getTaxonomyExtraction() );
\r
2297 private void setSpecialOptionsForNhxParser( final NHXParser nhx ) {
\r
2298 nhx.setReplaceUnderscores( getOptions().isReplaceUnderscoresInNhParsing() );
\r
2299 nhx.setTaxonomyExtraction( getOptions().getTaxonomyExtraction() );
\r
2300 nhx.setAllowErrorsInDistanceToParent( getOptions().isAllowErrorsInDistanceToParent() );
\r
2303 private void writeAllToFile() {
\r
2304 if ( ( getMainPanel().getTabbedPane() == null ) || ( getMainPanel().getTabbedPane().getTabCount() < 1 ) ) {
\r
2307 final File my_dir = getCurrentDir();
\r
2308 if ( my_dir != null ) {
\r
2309 _save_filechooser.setCurrentDirectory( my_dir );
\r
2311 _save_filechooser.setSelectedFile( new File( "" ) );
\r
2312 final int result = _save_filechooser.showSaveDialog( _contentpane );
\r
2313 final File file = _save_filechooser.getSelectedFile();
\r
2314 setCurrentDir( _save_filechooser.getCurrentDirectory() );
\r
2315 if ( ( file != null ) && ( result == JFileChooser.APPROVE_OPTION ) ) {
\r
2316 if ( file.exists() ) {
\r
2317 final int i = JOptionPane.showConfirmDialog( this,
\r
2318 file + " already exists. Overwrite?",
\r
2320 JOptionPane.OK_CANCEL_OPTION,
\r
2321 JOptionPane.WARNING_MESSAGE );
\r
2322 if ( i != JOptionPane.OK_OPTION ) {
\r
2329 catch ( final Exception e ) {
\r
2330 JOptionPane.showMessageDialog( this,
\r
2331 "Failed to delete: " + file,
\r
2333 JOptionPane.WARNING_MESSAGE );
\r
2337 final int count = getMainPanel().getTabbedPane().getTabCount();
\r
2338 final List<Phylogeny> trees = new ArrayList<Phylogeny>();
\r
2339 for( int i = 0; i < count; ++i ) {
\r
2340 final Phylogeny phy = getMainPanel().getPhylogeny( i );
\r
2341 if ( ForesterUtil.isEmpty( phy.getName() )
\r
2342 && !ForesterUtil.isEmpty( getMainPanel().getTabbedPane().getTitleAt( i ) ) ) {
\r
2343 phy.setName( getMainPanel().getTabbedPane().getTitleAt( i ) );
\r
2346 getMainPanel().getTreePanels().get( i ).setEdited( false );
\r
2348 final PhylogenyWriter writer = new PhylogenyWriter();
\r
2350 writer.toPhyloXML( file, trees, 0, ForesterUtil.LINE_SEPARATOR );
\r
2352 catch ( final IOException e ) {
\r
2353 JOptionPane.showMessageDialog( this,
\r
2354 "Failed to write to: " + file,
\r
2356 JOptionPane.WARNING_MESSAGE );
\r
2360 } // MainFrameApplication.
\r