fixed issue with UTF8 encoding.
[jalview.git] / forester / java / src / org / forester / archaeopteryx / MainFrameApplication.java
1 // $Id:\r
2 // FORESTER -- software libraries and applications\r
3 // for evolutionary biology research and applications.\r
4 //\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
9 //\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
14 //\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
19 //\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
23 //\r
24 // Contact: phylosoft @ gmail . com\r
25 // WWW: https://sites.google.com/site/cmzmasek/home/software/forester\r
26 \r
27 package org.forester.archaeopteryx;\r
28 \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
46 \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
61 \r
62 import org.forester.analysis.TaxonomyDataManager;\r
63 import org.forester.archaeopteryx.Options.CLADOGRAM_TYPE;\r
64 import org.forester.archaeopteryx.Options.NODE_LABEL_DIRECTION;\r
65 import org.forester.archaeopteryx.Options.PHYLOGENY_GRAPHICS_TYPE;\r
66 import org.forester.archaeopteryx.tools.InferenceManager;\r
67 import org.forester.archaeopteryx.tools.PhyloInferenceDialog;\r
68 import org.forester.archaeopteryx.tools.PhylogeneticInferenceOptions;\r
69 import org.forester.archaeopteryx.tools.PhylogeneticInferrer;\r
70 import org.forester.archaeopteryx.tools.SequenceDataRetriver;\r
71 import org.forester.archaeopteryx.webservices.PhylogeniesWebserviceClient;\r
72 import org.forester.archaeopteryx.webservices.WebservicesManager;\r
73 import org.forester.io.parsers.FastaParser;\r
74 import org.forester.io.parsers.GeneralMsaParser;\r
75 import org.forester.io.parsers.PhylogenyParser;\r
76 import org.forester.io.parsers.nexus.NexusPhylogeniesParser;\r
77 import org.forester.io.parsers.nhx.NHXParser;\r
78 import org.forester.io.parsers.nhx.NHXParser.TAXONOMY_EXTRACTION;\r
79 import org.forester.io.parsers.phyloxml.PhyloXmlDataFormatException;\r
80 import org.forester.io.parsers.phyloxml.PhyloXmlParser;\r
81 import org.forester.io.parsers.phyloxml.PhyloXmlUtil;\r
82 import org.forester.io.parsers.tol.TolParser;\r
83 import org.forester.io.parsers.util.ParserUtils;\r
84 import org.forester.io.writers.SequenceWriter;\r
85 import org.forester.msa.Msa;\r
86 import org.forester.msa.MsaFormatException;\r
87 import org.forester.phylogeny.Phylogeny;\r
88 import org.forester.phylogeny.PhylogenyMethods;\r
89 import org.forester.phylogeny.PhylogenyNode;\r
90 import org.forester.phylogeny.PhylogenyNode.NH_CONVERSION_SUPPORT_VALUE_STYLE;\r
91 import org.forester.phylogeny.data.Confidence;\r
92 import org.forester.phylogeny.data.PhylogenyDataUtil;\r
93 import org.forester.phylogeny.data.Sequence;\r
94 import org.forester.phylogeny.data.Taxonomy;\r
95 import org.forester.phylogeny.factories.ParserBasedPhylogenyFactory;\r
96 import org.forester.phylogeny.factories.PhylogenyFactory;\r
97 import org.forester.phylogeny.iterators.PhylogenyNodeIterator;\r
98 import org.forester.sequence.MolecularSequence;\r
99 import org.forester.util.BasicDescriptiveStatistics;\r
100 import org.forester.util.BasicTable;\r
101 import org.forester.util.BasicTableParser;\r
102 import org.forester.util.DescriptiveStatistics;\r
103 import org.forester.util.ForesterUtil;\r
104 \r
105 public final class MainFrameApplication extends MainFrame {\r
106 \r
107     private final static int             FRAME_X_SIZE                    = 800;\r
108     private final static int             FRAME_Y_SIZE                    = 800;\r
109     // Filters for the file-open dialog (classes defined in this file)\r
110     private static final long            serialVersionUID                = -799735726778865234L;\r
111     private static final boolean         PREPROCESS_TREES                = false;\r
112     private final JFileChooser           _values_filechooser;\r
113     private final JFileChooser           _sequences_filechooser;\r
114     private final JFileChooser           _open_filechooser;\r
115     private final JFileChooser           _msa_filechooser;\r
116     private final JFileChooser           _seqs_pi_filechooser;\r
117     private final JFileChooser           _open_filechooser_for_species_tree;\r
118     // Application-only print menu items\r
119     private JMenuItem                    _collapse_below_threshold;\r
120     private JMenuItem                    _collapse_below_branch_length;\r
121     private ButtonGroup                  _radio_group_1;\r
122     private ButtonGroup                  _radio_group_2;\r
123     // Others:\r
124     double                               _min_not_collapse               = AptxConstants.MIN_NOT_COLLAPSE_DEFAULT;\r
125     double                               _min_not_collapse_bl            = 0.001;\r
126     // Phylogeny Inference menu\r
127     private JMenu                        _inference_menu;\r
128     private JMenuItem                    _inference_from_msa_item;\r
129     private JMenuItem                    _inference_from_seqs_item;\r
130     // Phylogeny Inference\r
131     private PhylogeneticInferenceOptions _phylogenetic_inference_options = null;\r
132     private Msa                          _msa                            = null;\r
133     private File                         _msa_file                       = null;\r
134     private List<MolecularSequence>      _seqs                           = null;\r
135     private File                         _seqs_file                      = null;\r
136     JMenuItem                            _read_values_jmi;\r
137     JMenuItem                            _read_seqs_jmi;\r
138 \r
139     private MainFrameApplication( final Phylogeny[] phys, final Configuration config ) {\r
140         _configuration = config;\r
141         if ( _configuration == null ) {\r
142             throw new IllegalArgumentException( "configuration is null" );\r
143         }\r
144         setVisible( false );\r
145         setOptions( Options.createInstance( _configuration ) );\r
146         _mainpanel = new MainPanel( _configuration, this );\r
147         _open_filechooser = null;\r
148         _open_filechooser_for_species_tree = null;\r
149         _save_filechooser = null;\r
150         _writetopdf_filechooser = null;\r
151         _writetographics_filechooser = null;\r
152         _msa_filechooser = null;\r
153         _seqs_pi_filechooser = null;\r
154         _values_filechooser = null;\r
155         _sequences_filechooser = null;\r
156         _jmenubar = new JMenuBar();\r
157         buildFileMenu();\r
158         buildTypeMenu();\r
159         _contentpane = getContentPane();\r
160         _contentpane.setLayout( new BorderLayout() );\r
161         _contentpane.add( _mainpanel, BorderLayout.CENTER );\r
162         // App is this big\r
163         setSize( MainFrameApplication.FRAME_X_SIZE, MainFrameApplication.FRAME_Y_SIZE );\r
164         // The window listener\r
165         setDefaultCloseOperation( WindowConstants.DO_NOTHING_ON_CLOSE );\r
166         addWindowListener( new WindowAdapter() {\r
167 \r
168             @Override\r
169             public void windowClosing( final WindowEvent e ) {\r
170                 exit();\r
171             }\r
172         } );\r
173         //   setVisible( true );\r
174         if ( ( phys != null ) && ( phys.length > 0 ) ) {\r
175             AptxUtil.addPhylogeniesToTabs( phys, "", null, _configuration, _mainpanel );\r
176             validate();\r
177             getMainPanel().getControlPanel().showWholeAll();\r
178             getMainPanel().getControlPanel().showWhole();\r
179         }\r
180         //activateSaveAllIfNeeded();\r
181         // ...and its children\r
182         _contentpane.repaint();\r
183     }\r
184 \r
185     private MainFrameApplication( final Phylogeny[] phys, final Configuration config, final String title ) {\r
186         this( phys, config, title, null );\r
187     }\r
188 \r
189     private MainFrameApplication( final Phylogeny[] phys,\r
190                                   final Configuration config,\r
191                                   final String title,\r
192                                   final File current_dir ) {\r
193         super();\r
194         _configuration = config;\r
195         if ( _configuration == null ) {\r
196             throw new IllegalArgumentException( "configuration is null" );\r
197         }\r
198         try {\r
199             boolean synth_exception = false;\r
200             if ( AptxConstants.__SYNTH_LF ) {\r
201                 try {\r
202                     final SynthLookAndFeel synth = new SynthLookAndFeel();\r
203                     synth.load( MainFrameApplication.class.getResourceAsStream( "/resources/synth_look_and_feel_1.xml" ),\r
204                                 MainFrameApplication.class );\r
205                     UIManager.setLookAndFeel( synth );\r
206                 }\r
207                 catch ( final Exception ex ) {\r
208                     synth_exception = true;\r
209                     ForesterUtil.printWarningMessage( AptxConstants.PRG_NAME,\r
210                                                       "could not create synth look and feel: "\r
211                                                               + ex.getLocalizedMessage() );\r
212                 }\r
213             }\r
214             if ( !AptxConstants.__SYNTH_LF || synth_exception ) {\r
215                 if ( _configuration.isUseNativeUI() ) {\r
216                     UIManager.setLookAndFeel( UIManager.getSystemLookAndFeelClassName() );\r
217                 }\r
218                 else {\r
219                     UIManager.setLookAndFeel( UIManager.getCrossPlatformLookAndFeelClassName() );\r
220                 }\r
221             }\r
222             //UIManager.setLookAndFeel( "com.sun.java.swing.plaf.nimbus.NimbusLookAndFeel" );\r
223         }\r
224         catch ( final UnsupportedLookAndFeelException e ) {\r
225             AptxUtil.dieWithSystemError( "unsupported look and feel: " + e.toString() );\r
226         }\r
227         catch ( final ClassNotFoundException e ) {\r
228             AptxUtil.dieWithSystemError( "class not found exception: " + e.toString() );\r
229         }\r
230         catch ( final InstantiationException e ) {\r
231             AptxUtil.dieWithSystemError( "instantiation exception: " + e.toString() );\r
232         }\r
233         catch ( final IllegalAccessException e ) {\r
234             AptxUtil.dieWithSystemError( "illegal access exception: " + e.toString() );\r
235         }\r
236         if ( ( current_dir != null ) && current_dir.canRead() && current_dir.isDirectory() ) {\r
237             setCurrentDir( current_dir );\r
238         }\r
239         // hide until everything is ready\r
240         setVisible( false );\r
241         setOptions( Options.createInstance( _configuration ) );\r
242         setInferenceManager( InferenceManager.createInstance( _configuration ) );\r
243         setPhylogeneticInferenceOptions( PhylogeneticInferenceOptions.createInstance( _configuration ) );\r
244         // set title\r
245         setTitle( AptxConstants.PRG_NAME + " " + AptxConstants.VERSION + " (" + AptxConstants.PRG_DATE + ")" );\r
246         _mainpanel = new MainPanel( _configuration, this );\r
247         // The file dialogs\r
248         _open_filechooser = new JFileChooser();\r
249         _open_filechooser.setMultiSelectionEnabled( true );\r
250         _open_filechooser.addChoosableFileFilter( MainFrame.xmlfilter );\r
251         _open_filechooser.addChoosableFileFilter( MainFrame.nhxfilter );\r
252         _open_filechooser.addChoosableFileFilter( MainFrame.nhfilter );\r
253         _open_filechooser.addChoosableFileFilter( MainFrame.nexusfilter );\r
254         _open_filechooser.addChoosableFileFilter( MainFrame.tolfilter );\r
255         _open_filechooser.addChoosableFileFilter( _open_filechooser.getAcceptAllFileFilter() );\r
256         _open_filechooser.setFileFilter( MainFrame.defaultfilter );\r
257         _open_filechooser_for_species_tree = new JFileChooser();\r
258         _open_filechooser_for_species_tree.setMultiSelectionEnabled( false );\r
259         _open_filechooser_for_species_tree.addChoosableFileFilter( MainFrame.xmlfilter );\r
260         _open_filechooser_for_species_tree.addChoosableFileFilter( MainFrame.tolfilter );\r
261         _open_filechooser_for_species_tree.setFileFilter( MainFrame.xmlfilter );\r
262         // Msa:\r
263         _msa_filechooser = new JFileChooser();\r
264         _msa_filechooser.setName( "Read Multiple Sequence Alignment File" );\r
265         _msa_filechooser.setMultiSelectionEnabled( false );\r
266         _msa_filechooser.addChoosableFileFilter( _msa_filechooser.getAcceptAllFileFilter() );\r
267         _msa_filechooser.addChoosableFileFilter( MainFrame.msafilter );\r
268         // Seqs:\r
269         _seqs_pi_filechooser = new JFileChooser();\r
270         _seqs_pi_filechooser.setName( "Read Sequences File" );\r
271         _seqs_pi_filechooser.setMultiSelectionEnabled( false );\r
272         _seqs_pi_filechooser.addChoosableFileFilter( _seqs_pi_filechooser.getAcceptAllFileFilter() );\r
273         _seqs_pi_filechooser.addChoosableFileFilter( MainFrame.seqsfilter );\r
274         // Expression\r
275         _values_filechooser = new JFileChooser();\r
276         _values_filechooser.setMultiSelectionEnabled( false );\r
277         // Sequences\r
278         _sequences_filechooser = new JFileChooser();\r
279         _sequences_filechooser.setMultiSelectionEnabled( false );\r
280         try {\r
281             final String home_dir = System.getProperty( "user.home" );\r
282             _open_filechooser.setCurrentDirectory( new File( home_dir ) );\r
283             _open_filechooser_for_species_tree.setCurrentDirectory( new File( home_dir ) );\r
284             _msa_filechooser.setCurrentDirectory( new File( home_dir ) );\r
285             _seqs_pi_filechooser.setCurrentDirectory( new File( home_dir ) );\r
286             _values_filechooser.setCurrentDirectory( new File( home_dir ) );\r
287             _sequences_filechooser.setCurrentDirectory( new File( home_dir ) );\r
288         }\r
289         catch ( final Exception e ) {\r
290             e.printStackTrace();\r
291             // Do nothing. Not important.\r
292         }\r
293         // build the menu bar\r
294         _jmenubar = new JMenuBar();\r
295         if ( !_configuration.isUseNativeUI() ) {\r
296             _jmenubar.setBackground( getConfiguration().getGuiMenuBackgroundColor() );\r
297         }\r
298         buildFileMenu();\r
299         if ( AptxConstants.__ALLOW_PHYLOGENETIC_INFERENCE ) {\r
300             buildPhylogeneticInferenceMenu();\r
301         }\r
302         buildAnalysisMenu();\r
303         buildToolsMenu();\r
304         buildViewMenu();\r
305         buildFontSizeMenu();\r
306         buildOptionsMenu();\r
307         buildTypeMenu();\r
308         buildHelpMenu();\r
309         setJMenuBar( _jmenubar );\r
310         _jmenubar.add( _help_jmenu );\r
311         _contentpane = getContentPane();\r
312         _contentpane.setLayout( new BorderLayout() );\r
313         _contentpane.add( _mainpanel, BorderLayout.CENTER );\r
314         // App is this big\r
315         setSize( MainFrameApplication.FRAME_X_SIZE, MainFrameApplication.FRAME_Y_SIZE );\r
316         //        addWindowFocusListener( new WindowAdapter() {\r
317         //\r
318         //            @Override\r
319         //            public void windowGainedFocus( WindowEvent e ) {\r
320         //                requestFocusInWindow();\r
321         //            }\r
322         //        } );\r
323         // The window listener\r
324         setDefaultCloseOperation( WindowConstants.DO_NOTHING_ON_CLOSE );\r
325         addWindowListener( new WindowAdapter() {\r
326 \r
327             @Override\r
328             public void windowClosing( final WindowEvent e ) {\r
329                 if ( isUnsavedDataPresent() ) {\r
330                     final int r = JOptionPane.showConfirmDialog( null,\r
331                                                                  "Exit despite potentially unsaved changes?",\r
332                                                                  "Exit?",\r
333                                                                  JOptionPane.YES_NO_OPTION );\r
334                     if ( r != JOptionPane.YES_OPTION ) {\r
335                         return;\r
336                     }\r
337                 }\r
338                 else {\r
339                     final int r = JOptionPane.showConfirmDialog( null,\r
340                                                                  "Exit Archaeopteryx?",\r
341                                                                  "Exit?",\r
342                                                                  JOptionPane.YES_NO_OPTION );\r
343                     if ( r != JOptionPane.YES_OPTION ) {\r
344                         return;\r
345                     }\r
346                 }\r
347                 exit();\r
348             }\r
349         } );\r
350         // The component listener\r
351         addComponentListener( new ComponentAdapter() {\r
352 \r
353             @Override\r
354             public void componentResized( final ComponentEvent e ) {\r
355                 if ( _mainpanel.getCurrentTreePanel() != null ) {\r
356                     _mainpanel.getCurrentTreePanel().calcParametersForPainting( _mainpanel.getCurrentTreePanel()\r
357                                                                                 .getWidth(),\r
358                                                                                 _mainpanel.getCurrentTreePanel()\r
359                                                                                 .getHeight() );\r
360                 }\r
361             }\r
362         } );\r
363         requestFocusInWindow();\r
364         // addKeyListener( this );\r
365         setVisible( true );\r
366         if ( ( phys != null ) && ( phys.length > 0 ) ) {\r
367             AptxUtil.addPhylogeniesToTabs( phys, title, null, _configuration, _mainpanel );\r
368             validate();\r
369             getMainPanel().getControlPanel().showWholeAll();\r
370             getMainPanel().getControlPanel().showWhole();\r
371         }\r
372         activateSaveAllIfNeeded();\r
373         // ...and its children\r
374         _contentpane.repaint();\r
375         System.gc();\r
376     }\r
377 \r
378     private MainFrameApplication( final Phylogeny[] phys, final String config_file, final String title ) {\r
379         // Reads the config file (false, false => not url, not applet):\r
380         this( phys, new Configuration( config_file, false, false, true ), title );\r
381     }\r
382 \r
383     @Override\r
384     public void actionPerformed( final ActionEvent e ) {\r
385         try {\r
386             super.actionPerformed( e );\r
387             final Object o = e.getSource();\r
388             // Handle app-specific actions here:\r
389             if ( o == _open_item ) {\r
390                 readPhylogeniesFromFile();\r
391             }\r
392             if ( o == _open_url_item ) {\r
393                 readPhylogeniesFromURL();\r
394             }\r
395             else if ( o == _new_item ) {\r
396                 newTree();\r
397             }\r
398             else if ( o == _close_item ) {\r
399                 closeCurrentPane();\r
400             }\r
401             else if ( o == _load_species_tree_item ) {\r
402                 readSpeciesTreeFromFile();\r
403             }\r
404             else if ( o == _obtain_detailed_taxonomic_information_jmi ) {\r
405                 if ( isSubtreeDisplayed() ) {\r
406                     return;\r
407                 }\r
408                 obtainDetailedTaxonomicInformation();\r
409             }\r
410             else if ( o == _obtain_detailed_taxonomic_information_deleting_jmi ) {\r
411                 if ( isSubtreeDisplayed() ) {\r
412                     return;\r
413                 }\r
414                 obtainDetailedTaxonomicInformationDelete();\r
415             }\r
416             else if ( o == _obtain_seq_information_jmi ) {\r
417                 obtainSequenceInformation();\r
418             }\r
419             else if ( o == _read_values_jmi ) {\r
420                 if ( isSubtreeDisplayed() ) {\r
421                     return;\r
422                 }\r
423                 addExpressionValuesFromFile();\r
424             }\r
425             else if ( o == _read_seqs_jmi ) {\r
426                 if ( isSubtreeDisplayed() ) {\r
427                     return;\r
428                 }\r
429                 addSequencesFromFile();\r
430             }\r
431             else if ( o == _move_node_names_to_tax_sn_jmi ) {\r
432                 moveNodeNamesToTaxSn();\r
433             }\r
434             else if ( o == _move_node_names_to_seq_names_jmi ) {\r
435                 moveNodeNamesToSeqNames();\r
436             }\r
437             else if ( o == _extract_tax_code_from_node_names_jmi ) {\r
438                 extractTaxDataFromNodeNames();\r
439             }\r
440             else if ( o == _internal_number_are_confidence_for_nh_parsing_cbmi ) {\r
441                 updateOptions( getOptions() );\r
442             }\r
443             else if ( o == _replace_underscores_cbmi ) {\r
444                 if ( ( _extract_taxonomy_no_rbmi != null ) && !_extract_taxonomy_no_rbmi.isSelected() ) {\r
445                     _extract_taxonomy_no_rbmi.setSelected( true );\r
446                 }\r
447                 updateOptions( getOptions() );\r
448             }\r
449             else if ( o == _allow_errors_in_distance_to_parent_cbmi ) {\r
450                 updateOptions( getOptions() );\r
451             }\r
452             else if ( o == _collapse_below_threshold ) {\r
453                 if ( isSubtreeDisplayed() ) {\r
454                     return;\r
455                 }\r
456                 collapseBelowThreshold();\r
457             }\r
458             else if ( o == _collapse_below_branch_length ) {\r
459                 if ( isSubtreeDisplayed() ) {\r
460                     return;\r
461                 }\r
462                 collapseBelowBranchLengthThreshold();\r
463             }\r
464             else if ( ( o == _extract_taxonomy_pfam_strict_rbmi ) || ( o == _extract_taxonomy_pfam_relaxed_rbmi )\r
465                     || ( o == _extract_taxonomy_agressive_rbmi ) ) {\r
466                 if ( _replace_underscores_cbmi != null ) {\r
467                     _replace_underscores_cbmi.setSelected( false );\r
468                 }\r
469                 updateOptions( getOptions() );\r
470             }\r
471             else if ( o == _extract_taxonomy_no_rbmi ) {\r
472                 updateOptions( getOptions() );\r
473             }\r
474             else if ( o == _inference_from_msa_item ) {\r
475                 executePhyleneticInference( false );\r
476             }\r
477             else if ( o == _inference_from_seqs_item ) {\r
478                 executePhyleneticInference( true );\r
479             }\r
480             _contentpane.repaint();\r
481         }\r
482         catch ( final Exception ex ) {\r
483             AptxUtil.unexpectedException( ex );\r
484         }\r
485         catch ( final Error err ) {\r
486             AptxUtil.unexpectedError( err );\r
487         }\r
488     }\r
489 \r
490     public void end() {\r
491         _mainpanel.terminate();\r
492         _contentpane.removeAll();\r
493         setVisible( false );\r
494         dispose();\r
495     }\r
496 \r
497     @Override\r
498     public MainPanel getMainPanel() {\r
499         return _mainpanel;\r
500     }\r
501 \r
502     public Msa getMsa() {\r
503         return _msa;\r
504     }\r
505 \r
506     public File getMsaFile() {\r
507         return _msa_file;\r
508     }\r
509 \r
510     public List<MolecularSequence> getSeqs() {\r
511         return _seqs;\r
512     }\r
513 \r
514     public File getSeqsFile() {\r
515         return _seqs_file;\r
516     }\r
517 \r
518     public void readMsaFromFile() {\r
519         // Set an initial directory if none set yet\r
520         final File my_dir = getCurrentDir();\r
521         _msa_filechooser.setMultiSelectionEnabled( false );\r
522         // Open file-open dialog and set current directory\r
523         if ( my_dir != null ) {\r
524             _msa_filechooser.setCurrentDirectory( my_dir );\r
525         }\r
526         final int result = _msa_filechooser.showOpenDialog( _contentpane );\r
527         // All done: get the msa\r
528         final File file = _msa_filechooser.getSelectedFile();\r
529         setCurrentDir( _msa_filechooser.getCurrentDirectory() );\r
530         if ( ( file != null ) && !file.isDirectory() && ( result == JFileChooser.APPROVE_OPTION ) ) {\r
531             setMsaFile( null );\r
532             setMsa( null );\r
533             Msa msa = null;\r
534             try {\r
535                 final InputStream is = new FileInputStream( file );\r
536                 if ( FastaParser.isLikelyFasta( file ) ) {\r
537                     msa = FastaParser.parseMsa( is );\r
538                 }\r
539                 else {\r
540                     msa = GeneralMsaParser.parse( is );\r
541                 }\r
542             }\r
543             catch ( final MsaFormatException e ) {\r
544                 setArrowCursor();\r
545                 JOptionPane.showMessageDialog( this,\r
546                                                e.getLocalizedMessage(),\r
547                                                "Multiple sequence alignment format error",\r
548                                                JOptionPane.ERROR_MESSAGE );\r
549                 return;\r
550             }\r
551             catch ( final IOException e ) {\r
552                 setArrowCursor();\r
553                 JOptionPane.showMessageDialog( this,\r
554                                                e.getLocalizedMessage(),\r
555                                                "Failed to read multiple sequence alignment",\r
556                                                JOptionPane.ERROR_MESSAGE );\r
557                 return;\r
558             }\r
559             catch ( final IllegalArgumentException e ) {\r
560                 setArrowCursor();\r
561                 JOptionPane.showMessageDialog( this,\r
562                                                e.getLocalizedMessage(),\r
563                                                "Unexpected error during reading of multiple sequence alignment",\r
564                                                JOptionPane.ERROR_MESSAGE );\r
565                 return;\r
566             }\r
567             catch ( final Exception e ) {\r
568                 setArrowCursor();\r
569                 e.printStackTrace();\r
570                 JOptionPane.showMessageDialog( this,\r
571                                                e.getLocalizedMessage(),\r
572                                                "Unexpected error during reading of multiple sequence alignment",\r
573                                                JOptionPane.ERROR_MESSAGE );\r
574                 return;\r
575             }\r
576             if ( ( msa == null ) || ( msa.getNumberOfSequences() < 1 ) ) {\r
577                 JOptionPane.showMessageDialog( this,\r
578                                                "Multiple sequence alignment is empty",\r
579                                                "Illegal Multiple Sequence Alignment",\r
580                                                JOptionPane.ERROR_MESSAGE );\r
581                 return;\r
582             }\r
583             if ( msa.getNumberOfSequences() < 4 ) {\r
584                 JOptionPane.showMessageDialog( this,\r
585                                                "Multiple sequence alignment needs to contain at least 3 sequences",\r
586                                                "Illegal multiple sequence alignment",\r
587                                                JOptionPane.ERROR_MESSAGE );\r
588                 return;\r
589             }\r
590             if ( msa.getLength() < 2 ) {\r
591                 JOptionPane.showMessageDialog( this,\r
592                                                "Multiple sequence alignment needs to contain at least 2 residues",\r
593                                                "Illegal multiple sequence alignment",\r
594                                                JOptionPane.ERROR_MESSAGE );\r
595                 return;\r
596             }\r
597             System.gc();\r
598             setMsaFile( _msa_filechooser.getSelectedFile() );\r
599             setMsa( msa );\r
600         }\r
601     }\r
602 \r
603     public void readSeqsFromFileforPI() {\r
604         // Set an initial directory if none set yet\r
605         final File my_dir = getCurrentDir();\r
606         _seqs_pi_filechooser.setMultiSelectionEnabled( false );\r
607         // Open file-open dialog and set current directory\r
608         if ( my_dir != null ) {\r
609             _seqs_pi_filechooser.setCurrentDirectory( my_dir );\r
610         }\r
611         final int result = _seqs_pi_filechooser.showOpenDialog( _contentpane );\r
612         // All done: get the seqs\r
613         final File file = _seqs_pi_filechooser.getSelectedFile();\r
614         setCurrentDir( _seqs_pi_filechooser.getCurrentDirectory() );\r
615         if ( ( file != null ) && !file.isDirectory() && ( result == JFileChooser.APPROVE_OPTION ) ) {\r
616             setSeqsFile( null );\r
617             setSeqs( null );\r
618             List<MolecularSequence> seqs = null;\r
619             try {\r
620                 if ( FastaParser.isLikelyFasta( new FileInputStream( file ) ) ) {\r
621                     seqs = FastaParser.parse( new FileInputStream( file ) );\r
622                     for( final MolecularSequence seq : seqs ) {\r
623                         System.out.println( SequenceWriter.toFasta( seq, 60 ) );\r
624                     }\r
625                 }\r
626                 else {\r
627                     //TODO error\r
628                 }\r
629             }\r
630             catch ( final MsaFormatException e ) {\r
631                 setArrowCursor();\r
632                 JOptionPane.showMessageDialog( this,\r
633                                                e.getLocalizedMessage(),\r
634                                                "Multiple sequence file format error",\r
635                                                JOptionPane.ERROR_MESSAGE );\r
636                 return;\r
637             }\r
638             catch ( final IOException e ) {\r
639                 setArrowCursor();\r
640                 JOptionPane.showMessageDialog( this,\r
641                                                e.getLocalizedMessage(),\r
642                                                "Failed to read multiple sequence file",\r
643                                                JOptionPane.ERROR_MESSAGE );\r
644                 return;\r
645             }\r
646             catch ( final IllegalArgumentException e ) {\r
647                 setArrowCursor();\r
648                 JOptionPane.showMessageDialog( this,\r
649                                                e.getLocalizedMessage(),\r
650                                                "Unexpected error during reading of multiple sequence file",\r
651                                                JOptionPane.ERROR_MESSAGE );\r
652                 return;\r
653             }\r
654             catch ( final Exception e ) {\r
655                 setArrowCursor();\r
656                 e.printStackTrace();\r
657                 JOptionPane.showMessageDialog( this,\r
658                                                e.getLocalizedMessage(),\r
659                                                "Unexpected error during reading of multiple sequence file",\r
660                                                JOptionPane.ERROR_MESSAGE );\r
661                 return;\r
662             }\r
663             if ( ( seqs == null ) || ( seqs.size() < 1 ) ) {\r
664                 JOptionPane.showMessageDialog( this,\r
665                                                "Multiple sequence file is empty",\r
666                                                "Illegal multiple sequence file",\r
667                                                JOptionPane.ERROR_MESSAGE );\r
668                 return;\r
669             }\r
670             if ( seqs.size() < 4 ) {\r
671                 JOptionPane.showMessageDialog( this,\r
672                                                "Multiple sequence file needs to contain at least 3 sequences",\r
673                                                "Illegal multiple sequence file",\r
674                                                JOptionPane.ERROR_MESSAGE );\r
675                 return;\r
676             }\r
677             //  if ( msa.getLength() < 2 ) {\r
678             //       JOptionPane.showMessageDialog( this,\r
679             //                                      "Multiple sequence alignment needs to contain at least 2 residues",\r
680             //                                      "Illegal multiple sequence file",\r
681             //                                      JOptionPane.ERROR_MESSAGE );\r
682             //       return;\r
683             //   }\r
684             System.gc();\r
685             setSeqsFile( _seqs_pi_filechooser.getSelectedFile() );\r
686             setSeqs( seqs );\r
687         }\r
688     }\r
689 \r
690     private void addExpressionValuesFromFile() {\r
691         if ( ( getCurrentTreePanel() == null ) || ( getCurrentTreePanel().getPhylogeny() == null ) ) {\r
692             JOptionPane.showMessageDialog( this,\r
693                                            "Need to load evolutionary tree first",\r
694                                            "Can Not Read Expression Values",\r
695                                            JOptionPane.WARNING_MESSAGE );\r
696             return;\r
697         }\r
698         final File my_dir = getCurrentDir();\r
699         if ( my_dir != null ) {\r
700             _values_filechooser.setCurrentDirectory( my_dir );\r
701         }\r
702         final int result = _values_filechooser.showOpenDialog( _contentpane );\r
703         final File file = _values_filechooser.getSelectedFile();\r
704         if ( ( file != null ) && ( file.length() > 0 ) && ( result == JFileChooser.APPROVE_OPTION ) ) {\r
705             BasicTable<String> t = null;\r
706             try {\r
707                 t = BasicTableParser.parse( file, '\t' );\r
708                 if ( t.getNumberOfColumns() < 2 ) {\r
709                     t = BasicTableParser.parse( file, ',' );\r
710                 }\r
711                 if ( t.getNumberOfColumns() < 2 ) {\r
712                     t = BasicTableParser.parse( file, ' ' );\r
713                 }\r
714             }\r
715             catch ( final IOException e ) {\r
716                 JOptionPane.showMessageDialog( this,\r
717                                                e.getMessage(),\r
718                                                "Could Not Read Expression Value Table",\r
719                                                JOptionPane.ERROR_MESSAGE );\r
720                 return;\r
721             }\r
722             if ( t.getNumberOfColumns() < 2 ) {\r
723                 JOptionPane.showMessageDialog( this,\r
724                                                "Table contains " + t.getNumberOfColumns() + " column(s)",\r
725                                                "Problem with Expression Value Table",\r
726                                                JOptionPane.ERROR_MESSAGE );\r
727                 return;\r
728             }\r
729             if ( t.getNumberOfRows() < 1 ) {\r
730                 JOptionPane.showMessageDialog( this,\r
731                                                "Table contains zero rows",\r
732                                                "Problem with Expression Value Table",\r
733                                                JOptionPane.ERROR_MESSAGE );\r
734                 return;\r
735             }\r
736             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
737             if ( t.getNumberOfRows() != phy.getNumberOfExternalNodes() ) {\r
738                 JOptionPane.showMessageDialog( this,\r
739                                                "Table contains " + t.getNumberOfRows() + " rows, but tree contains "\r
740                                                        + phy.getNumberOfExternalNodes() + " external nodes",\r
741                                                        "Warning",\r
742                                                        JOptionPane.WARNING_MESSAGE );\r
743             }\r
744             final DescriptiveStatistics stats = new BasicDescriptiveStatistics();\r
745             int not_found = 0;\r
746             for( final PhylogenyNodeIterator iter = phy.iteratorPreorder(); iter.hasNext(); ) {\r
747                 final PhylogenyNode node = iter.next();\r
748                 final String node_name = node.getName();\r
749                 if ( !ForesterUtil.isEmpty( node_name ) ) {\r
750                     int row = -1;\r
751                     try {\r
752                         row = t.findRow( node_name );\r
753                     }\r
754                     catch ( final IllegalArgumentException e ) {\r
755                         JOptionPane\r
756                         .showMessageDialog( this,\r
757                                             e.getMessage(),\r
758                                             "Error Mapping Node Identifiers to Expression Value Identifiers",\r
759                                             JOptionPane.ERROR_MESSAGE );\r
760                         return;\r
761                     }\r
762                     if ( row < 0 ) {\r
763                         if ( node.isExternal() ) {\r
764                             not_found++;\r
765                         }\r
766                         continue;\r
767                     }\r
768                     final List<Double> l = new ArrayList<Double>();\r
769                     for( int col = 1; col < t.getNumberOfColumns(); ++col ) {\r
770                         double d = -100;\r
771                         try {\r
772                             d = Double.parseDouble( t.getValueAsString( col, row ) );\r
773                         }\r
774                         catch ( final NumberFormatException e ) {\r
775                             JOptionPane.showMessageDialog( this,\r
776                                                            "Could not parse \"" + t.getValueAsString( col, row )\r
777                                                            + "\" into a decimal value",\r
778                                                            "Issue with Expression Value Table",\r
779                                                            JOptionPane.ERROR_MESSAGE );\r
780                             return;\r
781                         }\r
782                         stats.addValue( d );\r
783                         l.add( d );\r
784                     }\r
785                     if ( !l.isEmpty() ) {\r
786                         if ( node.getNodeData().getProperties() != null ) {\r
787                             node.getNodeData().getProperties()\r
788                             .removePropertiesWithGivenReferencePrefix( PhyloXmlUtil.VECTOR_PROPERTY_REF );\r
789                         }\r
790                         node.getNodeData().setVector( l );\r
791                     }\r
792                 }\r
793             }\r
794             if ( not_found > 0 ) {\r
795                 JOptionPane.showMessageDialog( this, "Could not fine expression values for " + not_found\r
796                                                + " external node(s)", "Warning", JOptionPane.WARNING_MESSAGE );\r
797             }\r
798             getCurrentTreePanel().setStatisticsForExpressionValues( stats );\r
799         }\r
800     }\r
801 \r
802     private void addSequencesFromFile() {\r
803         if ( ( getCurrentTreePanel() == null ) || ( getCurrentTreePanel().getPhylogeny() == null ) ) {\r
804             JOptionPane.showMessageDialog( this,\r
805                                            "Need to load evolutionary tree first",\r
806                                            "Can Not Read Sequences",\r
807                                            JOptionPane.WARNING_MESSAGE );\r
808             return;\r
809         }\r
810         final File my_dir = getCurrentDir();\r
811         if ( my_dir != null ) {\r
812             _sequences_filechooser.setCurrentDirectory( my_dir );\r
813         }\r
814         final int result = _sequences_filechooser.showOpenDialog( _contentpane );\r
815         final File file = _sequences_filechooser.getSelectedFile();\r
816         List<MolecularSequence> seqs = null;\r
817         if ( ( file != null ) && !file.isDirectory() && ( result == JFileChooser.APPROVE_OPTION ) ) {\r
818             try {\r
819                 if ( FastaParser.isLikelyFasta( new FileInputStream( file ) ) ) {\r
820                     seqs = FastaParser.parse( new FileInputStream( file ) );\r
821                 }\r
822                 else {\r
823                     JOptionPane.showMessageDialog( this,\r
824                                                    "Format does not appear to be Fasta",\r
825                                                    "Multiple sequence file format error",\r
826                                                    JOptionPane.ERROR_MESSAGE );\r
827                     return;\r
828                 }\r
829             }\r
830             catch ( final MsaFormatException e ) {\r
831                 setArrowCursor();\r
832                 JOptionPane.showMessageDialog( this,\r
833                                                e.getLocalizedMessage(),\r
834                                                "Multiple sequence file format error",\r
835                                                JOptionPane.ERROR_MESSAGE );\r
836                 return;\r
837             }\r
838             catch ( final IOException e ) {\r
839                 setArrowCursor();\r
840                 JOptionPane.showMessageDialog( this,\r
841                                                e.getLocalizedMessage(),\r
842                                                "Failed to read multiple sequence file",\r
843                                                JOptionPane.ERROR_MESSAGE );\r
844                 return;\r
845             }\r
846             catch ( final Exception e ) {\r
847                 setArrowCursor();\r
848                 e.printStackTrace();\r
849                 JOptionPane.showMessageDialog( this,\r
850                                                e.getLocalizedMessage(),\r
851                                                "Unexpected error during reading of multiple sequence file",\r
852                                                JOptionPane.ERROR_MESSAGE );\r
853                 return;\r
854             }\r
855             if ( ( seqs == null ) || ( seqs.size() < 1 ) ) {\r
856                 JOptionPane.showMessageDialog( this,\r
857                                                "Multiple sequence file is empty",\r
858                                                "Empty multiple sequence file",\r
859                                                JOptionPane.ERROR_MESSAGE );\r
860                 setArrowCursor();\r
861                 return;\r
862             }\r
863         }\r
864         if ( seqs != null ) {\r
865             for( final MolecularSequence seq : seqs ) {\r
866                 System.out.println( seq.getIdentifier() );\r
867             }\r
868             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
869             int total_counter = 0;\r
870             int attached_counter = 0;\r
871             for( final MolecularSequence seq : seqs ) {\r
872                 ++total_counter;\r
873                 final String seq_name = seq.getIdentifier();\r
874                 if ( !ForesterUtil.isEmpty( seq_name ) ) {\r
875                     List<PhylogenyNode> nodes = phy.getNodesViaSequenceName( seq_name );\r
876                     if ( nodes.isEmpty() ) {\r
877                         nodes = phy.getNodesViaSequenceSymbol( seq_name );\r
878                     }\r
879                     if ( nodes.isEmpty() ) {\r
880                         nodes = phy.getNodesViaGeneName( seq_name );\r
881                     }\r
882                     if ( nodes.isEmpty() ) {\r
883                         nodes = phy.getNodes( seq_name );\r
884                     }\r
885                     if ( nodes.size() > 1 ) {\r
886                         JOptionPane.showMessageDialog( this,\r
887                                                        "Sequence name \"" + seq_name + "\" is not unique",\r
888                                                        "Sequence name not unique",\r
889                                                        JOptionPane.ERROR_MESSAGE );\r
890                         setArrowCursor();\r
891                         return;\r
892                     }\r
893                     final String[] a = seq_name.split( "\\s" );\r
894                     if ( nodes.isEmpty() && ( a.length > 1 ) ) {\r
895                         final String seq_name_split = a[ 0 ];\r
896                         nodes = phy.getNodesViaSequenceName( seq_name_split );\r
897                         if ( nodes.isEmpty() ) {\r
898                             nodes = phy.getNodesViaSequenceSymbol( seq_name_split );\r
899                         }\r
900                         if ( nodes.isEmpty() ) {\r
901                             nodes = phy.getNodes( seq_name_split );\r
902                         }\r
903                         if ( nodes.size() > 1 ) {\r
904                             JOptionPane.showMessageDialog( this, "Split sequence name \"" + seq_name_split\r
905                                                            + "\" is not unique", "Sequence name not unique", JOptionPane.ERROR_MESSAGE );\r
906                             setArrowCursor();\r
907                             return;\r
908                         }\r
909                     }\r
910                     if ( nodes.size() == 1 ) {\r
911                         ++attached_counter;\r
912                         final PhylogenyNode n = nodes.get( 0 );\r
913                         if ( !n.getNodeData().isHasSequence() ) {\r
914                             n.getNodeData().addSequence( new org.forester.phylogeny.data.Sequence() );\r
915                         }\r
916                         n.getNodeData().getSequence().setMolecularSequence( seq.getMolecularSequenceAsString() );\r
917                         if ( ForesterUtil.isEmpty( n.getNodeData().getSequence().getName() ) ) {\r
918                             n.getNodeData().getSequence().setName( seq_name );\r
919                         }\r
920                     }\r
921                 }\r
922             }\r
923             if ( attached_counter > 0 ) {\r
924                 int ext_nodes = 0;\r
925                 int ext_nodes_with_seq = 0;\r
926                 for( final PhylogenyNodeIterator iter = phy.iteratorExternalForward(); iter.hasNext(); ) {\r
927                     ++ext_nodes;\r
928                     final PhylogenyNode n = iter.next();\r
929                     if ( n.getNodeData().isHasSequence()\r
930                             && !ForesterUtil.isEmpty( n.getNodeData().getSequence().getMolecularSequence() ) ) {\r
931                         ++ext_nodes_with_seq;\r
932                     }\r
933                 }\r
934                 final String s;\r
935                 if ( ext_nodes == ext_nodes_with_seq ) {\r
936                     s = "All " + ext_nodes_with_seq + " external nodes now have a molecular sequence attached to them.";\r
937                 }\r
938                 else {\r
939                     s = ext_nodes_with_seq + " out of " + ext_nodes\r
940                             + " external nodes now have a molecular sequence attached to them.";\r
941                 }\r
942                 if ( ( attached_counter == total_counter ) && ( ext_nodes == ext_nodes_with_seq ) ) {\r
943                     JOptionPane.showMessageDialog( this,\r
944                                                    "Attached all " + total_counter + " sequences to tree nodes.\n" + s,\r
945                                                    "All sequences attached",\r
946                                                    JOptionPane.INFORMATION_MESSAGE );\r
947                 }\r
948                 else {\r
949                     JOptionPane.showMessageDialog( this, "Attached " + attached_counter\r
950                                                    + " sequences out of a total of " + total_counter + " sequences.\n" + s, attached_counter\r
951                                                    + " sequences attached", JOptionPane.WARNING_MESSAGE );\r
952                 }\r
953             }\r
954             else {\r
955                 JOptionPane.showMessageDialog( this, "No maching tree node for any of the " + total_counter\r
956                                                + " sequences", "Could not attach any sequences", JOptionPane.ERROR_MESSAGE );\r
957             }\r
958         }\r
959     }\r
960 \r
961     private void closeCurrentPane() {\r
962         if ( getMainPanel().getCurrentTreePanel() != null ) {\r
963             if ( getMainPanel().getCurrentTreePanel().isEdited() ) {\r
964                 final int r = JOptionPane.showConfirmDialog( this,\r
965                                                              "Close tab despite potentially unsaved changes?",\r
966                                                              "Close Tab?",\r
967                                                              JOptionPane.YES_NO_OPTION );\r
968                 if ( r != JOptionPane.YES_OPTION ) {\r
969                     return;\r
970                 }\r
971             }\r
972             getMainPanel().closeCurrentPane();\r
973             activateSaveAllIfNeeded();\r
974         }\r
975     }\r
976 \r
977     private void collapse( final Phylogeny phy ) {\r
978         final PhylogenyNodeIterator it = phy.iteratorPostorder();\r
979         final List<PhylogenyNode> to_be_removed = new ArrayList<PhylogenyNode>();\r
980         double min_support = Double.MAX_VALUE;\r
981         boolean conf_present = false;\r
982         while ( it.hasNext() ) {\r
983             final PhylogenyNode n = it.next();\r
984             if ( !n.isExternal() && !n.isRoot() ) {\r
985                 final List<Confidence> c = n.getBranchData().getConfidences();\r
986                 if ( ( c != null ) && ( c.size() > 0 ) ) {\r
987                     conf_present = true;\r
988                     double max = 0;\r
989                     for( final Confidence confidence : c ) {\r
990                         if ( confidence.getValue() > max ) {\r
991                             max = confidence.getValue();\r
992                         }\r
993                     }\r
994                     if ( max < getMinNotCollapseConfidenceValue() ) {\r
995                         to_be_removed.add( n );\r
996                     }\r
997                     if ( max < min_support ) {\r
998                         min_support = max;\r
999                     }\r
1000                 }\r
1001             }\r
1002         }\r
1003         if ( conf_present ) {\r
1004             for( final PhylogenyNode node : to_be_removed ) {\r
1005                 PhylogenyMethods.removeNode( node, phy );\r
1006             }\r
1007             if ( to_be_removed.size() > 0 ) {\r
1008                 phy.externalNodesHaveChanged();\r
1009                 phy.clearHashIdToNodeMap();\r
1010                 phy.recalculateNumberOfExternalDescendants( true );\r
1011                 getCurrentTreePanel().resetNodeIdToDistToLeafMap();\r
1012                 getCurrentTreePanel().updateSetOfCollapsedExternalNodes();\r
1013                 getCurrentTreePanel().calculateLongestExtNodeInfo();\r
1014                 getCurrentTreePanel().setNodeInPreorderToNull();\r
1015                 getCurrentTreePanel().recalculateMaxDistanceToRoot();\r
1016                 getCurrentTreePanel().resetPreferredSize();\r
1017                 getCurrentTreePanel().setEdited( true );\r
1018                 getCurrentTreePanel().repaint();\r
1019                 repaint();\r
1020             }\r
1021             if ( to_be_removed.size() > 0 ) {\r
1022                 JOptionPane.showMessageDialog( this, "Collapsed " + to_be_removed.size()\r
1023                                                + " branches with\nconfidence values below " + getMinNotCollapseConfidenceValue(), "Collapsed "\r
1024                                                        + to_be_removed.size() + " branches", JOptionPane.INFORMATION_MESSAGE );\r
1025             }\r
1026             else {\r
1027                 JOptionPane.showMessageDialog( this, "No branch collapsed,\nminimum confidence value per branch is "\r
1028                         + min_support, "No branch collapsed", JOptionPane.INFORMATION_MESSAGE );\r
1029             }\r
1030         }\r
1031         else {\r
1032             JOptionPane.showMessageDialog( this,\r
1033                                            "No branch collapsed because no confidence values present",\r
1034                                            "No confidence values present",\r
1035                                            JOptionPane.INFORMATION_MESSAGE );\r
1036         }\r
1037     }\r
1038 \r
1039     private void collapseBelowBranchLengthThreshold() {\r
1040         if ( getCurrentTreePanel() != null ) {\r
1041             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1042             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1043                 final String s = ( String ) JOptionPane\r
1044                         .showInputDialog( this,\r
1045                                           "Please enter the minimum branch length value\n",\r
1046                                           "Minimal Branch Length Value",\r
1047                                           JOptionPane.QUESTION_MESSAGE,\r
1048                                           null,\r
1049                                           null,\r
1050                                           getMinNotCollapseBlValue() );\r
1051                 if ( !ForesterUtil.isEmpty( s ) ) {\r
1052                     boolean success = true;\r
1053                     double m = 0.0;\r
1054                     final String m_str = s.trim();\r
1055                     if ( !ForesterUtil.isEmpty( m_str ) ) {\r
1056                         try {\r
1057                             m = Double.parseDouble( m_str );\r
1058                         }\r
1059                         catch ( final Exception ex ) {\r
1060                             success = false;\r
1061                         }\r
1062                     }\r
1063                     else {\r
1064                         success = false;\r
1065                     }\r
1066                     if ( success && ( m >= 0.0 ) ) {\r
1067                         setMinNotCollapseBlValue( m );\r
1068                         collapseBl( phy );\r
1069                     }\r
1070                 }\r
1071             }\r
1072         }\r
1073     }\r
1074 \r
1075     private void collapseBelowThreshold() {\r
1076         if ( getCurrentTreePanel() != null ) {\r
1077             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1078             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1079                 final String s = ( String ) JOptionPane.showInputDialog( this,\r
1080                                                                          "Please enter the minimum confidence value\n",\r
1081                                                                          "Minimal Confidence Value",\r
1082                                                                          JOptionPane.QUESTION_MESSAGE,\r
1083                                                                          null,\r
1084                                                                          null,\r
1085                                                                          getMinNotCollapseConfidenceValue() );\r
1086                 if ( !ForesterUtil.isEmpty( s ) ) {\r
1087                     boolean success = true;\r
1088                     double m = 0.0;\r
1089                     final String m_str = s.trim();\r
1090                     if ( !ForesterUtil.isEmpty( m_str ) ) {\r
1091                         try {\r
1092                             m = Double.parseDouble( m_str );\r
1093                         }\r
1094                         catch ( final Exception ex ) {\r
1095                             success = false;\r
1096                         }\r
1097                     }\r
1098                     else {\r
1099                         success = false;\r
1100                     }\r
1101                     if ( success && ( m >= 0.0 ) ) {\r
1102                         setMinNotCollapseConfidenceValue( m );\r
1103                         collapse( phy );\r
1104                     }\r
1105                 }\r
1106             }\r
1107         }\r
1108     }\r
1109 \r
1110     private void collapseBl( final Phylogeny phy ) {\r
1111         final PhylogenyNodeIterator it = phy.iteratorPostorder();\r
1112         final List<PhylogenyNode> to_be_removed = new ArrayList<PhylogenyNode>();\r
1113         double min_bl = Double.MAX_VALUE;\r
1114         boolean bl_present = false;\r
1115         while ( it.hasNext() ) {\r
1116             final PhylogenyNode n = it.next();\r
1117             if ( !n.isExternal() && !n.isRoot() ) {\r
1118                 final double bl = n.getDistanceToParent();\r
1119                 if ( bl != PhylogenyDataUtil.BRANCH_LENGTH_DEFAULT ) {\r
1120                     bl_present = true;\r
1121                     if ( bl < getMinNotCollapseBlValue() ) {\r
1122                         to_be_removed.add( n );\r
1123                     }\r
1124                     if ( bl < min_bl ) {\r
1125                         min_bl = bl;\r
1126                     }\r
1127                 }\r
1128             }\r
1129         }\r
1130         if ( bl_present ) {\r
1131             for( final PhylogenyNode node : to_be_removed ) {\r
1132                 PhylogenyMethods.removeNode( node, phy );\r
1133             }\r
1134             if ( to_be_removed.size() > 0 ) {\r
1135                 phy.externalNodesHaveChanged();\r
1136                 phy.clearHashIdToNodeMap();\r
1137                 phy.recalculateNumberOfExternalDescendants( true );\r
1138                 getCurrentTreePanel().resetNodeIdToDistToLeafMap();\r
1139                 getCurrentTreePanel().updateSetOfCollapsedExternalNodes();\r
1140                 getCurrentTreePanel().calculateLongestExtNodeInfo();\r
1141                 getCurrentTreePanel().setNodeInPreorderToNull();\r
1142                 getCurrentTreePanel().recalculateMaxDistanceToRoot();\r
1143                 getCurrentTreePanel().resetPreferredSize();\r
1144                 getCurrentTreePanel().setEdited( true );\r
1145                 getCurrentTreePanel().repaint();\r
1146                 repaint();\r
1147             }\r
1148             if ( to_be_removed.size() > 0 ) {\r
1149                 JOptionPane.showMessageDialog( this, "Collapsed " + to_be_removed.size()\r
1150                                                + " branches with\nbranch length values below " + getMinNotCollapseBlValue(), "Collapsed "\r
1151                                                        + to_be_removed.size() + " branches", JOptionPane.INFORMATION_MESSAGE );\r
1152             }\r
1153             else {\r
1154                 JOptionPane.showMessageDialog( this,\r
1155                                                "No branch collapsed,\nminimum branch length is " + min_bl,\r
1156                                                "No branch collapsed",\r
1157                                                JOptionPane.INFORMATION_MESSAGE );\r
1158             }\r
1159         }\r
1160         else {\r
1161             JOptionPane.showMessageDialog( this,\r
1162                                            "No branch collapsed because no branch length values present",\r
1163                                            "No branch length values present",\r
1164                                            JOptionPane.INFORMATION_MESSAGE );\r
1165         }\r
1166     }\r
1167 \r
1168     private PhyloXmlParser createPhyloXmlParser() {\r
1169         PhyloXmlParser xml_parser = null;\r
1170         if ( getConfiguration().isValidatePhyloXmlAgainstSchema() ) {\r
1171             try {\r
1172                 xml_parser = PhyloXmlParser.createPhyloXmlParserXsdValidating();\r
1173             }\r
1174             catch ( final Exception e ) {\r
1175                 JOptionPane.showMessageDialog( this,\r
1176                                                e.getLocalizedMessage(),\r
1177                                                "failed to create validating XML parser",\r
1178                                                JOptionPane.WARNING_MESSAGE );\r
1179             }\r
1180         }\r
1181         if ( xml_parser == null ) {\r
1182             xml_parser = PhyloXmlParser.createPhyloXmlParser();\r
1183         }\r
1184         return xml_parser;\r
1185     }\r
1186 \r
1187     private void executePhyleneticInference( final boolean from_unaligned_seqs ) {\r
1188         final PhyloInferenceDialog dialog = new PhyloInferenceDialog( this,\r
1189                                                                       getPhylogeneticInferenceOptions(),\r
1190                                                                       from_unaligned_seqs );\r
1191         dialog.activate();\r
1192         if ( dialog.getValue() == JOptionPane.OK_OPTION ) {\r
1193             if ( !from_unaligned_seqs ) {\r
1194                 if ( getMsa() != null ) {\r
1195                     final PhylogeneticInferrer inferrer = new PhylogeneticInferrer( getMsa(),\r
1196                                                                                     getPhylogeneticInferenceOptions()\r
1197                                                                                     .copy(), this );\r
1198                     new Thread( inferrer ).start();\r
1199                 }\r
1200                 else {\r
1201                     JOptionPane.showMessageDialog( this,\r
1202                                                    "No multiple sequence alignment selected",\r
1203                                                    "Phylogenetic Inference Not Launched",\r
1204                                                    JOptionPane.WARNING_MESSAGE );\r
1205                 }\r
1206             }\r
1207             else {\r
1208                 if ( getSeqs() != null ) {\r
1209                     final PhylogeneticInferrer inferrer = new PhylogeneticInferrer( getSeqs(),\r
1210                                                                                     getPhylogeneticInferenceOptions()\r
1211                                                                                     .copy(), this );\r
1212                     new Thread( inferrer ).start();\r
1213                 }\r
1214                 else {\r
1215                     JOptionPane.showMessageDialog( this,\r
1216                                                    "No input sequences selected",\r
1217                                                    "Phylogenetic Inference Not Launched",\r
1218                                                    JOptionPane.WARNING_MESSAGE );\r
1219                 }\r
1220             }\r
1221         }\r
1222     }\r
1223 \r
1224     private void extractTaxDataFromNodeNames() throws PhyloXmlDataFormatException {\r
1225         final StringBuilder sb = new StringBuilder();\r
1226         final StringBuilder sb_failed = new StringBuilder();\r
1227         int counter = 0;\r
1228         int counter_failed = 0;\r
1229         if ( getCurrentTreePanel() != null ) {\r
1230             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1231             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1232                 final PhylogenyNodeIterator it = phy.iteratorExternalForward();\r
1233                 while ( it.hasNext() ) {\r
1234                     final PhylogenyNode n = it.next();\r
1235                     final String name = n.getName().trim();\r
1236                     if ( !ForesterUtil.isEmpty( name ) ) {\r
1237                         final String nt = ParserUtils.extractTaxonomyDataFromNodeName( n,\r
1238                                                                                        TAXONOMY_EXTRACTION.AGGRESSIVE );\r
1239                         if ( !ForesterUtil.isEmpty( nt ) ) {\r
1240                             if ( counter < 15 ) {\r
1241                                 sb.append( name + ": " + nt + "\n" );\r
1242                             }\r
1243                             else if ( counter == 15 ) {\r
1244                                 sb.append( "...\n" );\r
1245                             }\r
1246                             counter++;\r
1247                         }\r
1248                         else {\r
1249                             if ( counter_failed < 15 ) {\r
1250                                 sb_failed.append( name + "\n" );\r
1251                             }\r
1252                             else if ( counter_failed == 15 ) {\r
1253                                 sb_failed.append( "...\n" );\r
1254                             }\r
1255                             counter_failed++;\r
1256                         }\r
1257                     }\r
1258                 }\r
1259                 if ( counter > 0 ) {\r
1260                     String failed = "";\r
1261                     String all = "all ";\r
1262                     if ( counter_failed > 0 ) {\r
1263                         all = "";\r
1264                         failed = "\nCould not extract taxonomic data for " + counter_failed\r
1265                                 + " named external nodes:\n" + sb_failed;\r
1266                     }\r
1267                     JOptionPane.showMessageDialog( this,\r
1268                                                    "Extracted taxonomic data from " + all + counter\r
1269                                                    + " named external nodes:\n" + sb.toString() + failed,\r
1270                                                    "Taxonomic Data Extraction Completed",\r
1271                                                    counter_failed > 0 ? JOptionPane.WARNING_MESSAGE\r
1272                                                            : JOptionPane.INFORMATION_MESSAGE );\r
1273                 }\r
1274                 else {\r
1275                     JOptionPane\r
1276                     .showMessageDialog( this,\r
1277                                         "Could not extract any taxonomic data.\nMaybe node names are empty\n"\r
1278                                                 + "or not in the forms \"XYZ_CAEEL\", \"XYZ_6239\", or \"XYZ_Caenorhabditis_elegans\"\n"\r
1279                                                 + "or nodes already have taxonomic data?\n",\r
1280                                                 "No Taxonomic Data Extracted",\r
1281                                                 JOptionPane.ERROR_MESSAGE );\r
1282                 }\r
1283             }\r
1284         }\r
1285     }\r
1286 \r
1287     private double getMinNotCollapseBlValue() {\r
1288         return _min_not_collapse_bl;\r
1289     }\r
1290 \r
1291     private double getMinNotCollapseConfidenceValue() {\r
1292         return _min_not_collapse;\r
1293     }\r
1294 \r
1295     private PhylogeneticInferenceOptions getPhylogeneticInferenceOptions() {\r
1296         if ( _phylogenetic_inference_options == null ) {\r
1297             _phylogenetic_inference_options = new PhylogeneticInferenceOptions();\r
1298         }\r
1299         return _phylogenetic_inference_options;\r
1300     }\r
1301 \r
1302     private boolean isUnsavedDataPresent() {\r
1303         final List<TreePanel> tps = getMainPanel().getTreePanels();\r
1304         for( final TreePanel tp : tps ) {\r
1305             if ( tp.isEdited() ) {\r
1306                 return true;\r
1307             }\r
1308         }\r
1309         return false;\r
1310     }\r
1311 \r
1312     private void moveNodeNamesToSeqNames() throws PhyloXmlDataFormatException {\r
1313         if ( getCurrentTreePanel() != null ) {\r
1314             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1315             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1316                 PhylogenyMethods\r
1317                 .transferNodeNameToField( phy, PhylogenyMethods.PhylogenyNodeField.SEQUENCE_NAME, false );\r
1318             }\r
1319         }\r
1320     }\r
1321 \r
1322     private void moveNodeNamesToTaxSn() throws PhyloXmlDataFormatException {\r
1323         if ( getCurrentTreePanel() != null ) {\r
1324             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1325             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1326                 PhylogenyMethods.transferNodeNameToField( phy,\r
1327                                                           PhylogenyMethods.PhylogenyNodeField.TAXONOMY_SCIENTIFIC_NAME,\r
1328                                                           false );\r
1329             }\r
1330         }\r
1331     }\r
1332 \r
1333     private void newTree() {\r
1334         final Phylogeny[] phys = new Phylogeny[ 1 ];\r
1335         final Phylogeny phy = new Phylogeny();\r
1336         final PhylogenyNode node = new PhylogenyNode();\r
1337         phy.setRoot( node );\r
1338         phy.setRooted( true );\r
1339         phys[ 0 ] = phy;\r
1340         AptxUtil.addPhylogeniesToTabs( phys, "", "", getConfiguration(), getMainPanel() );\r
1341         _mainpanel.getControlPanel().showWhole();\r
1342         _mainpanel.getCurrentTreePanel().setPhylogenyGraphicsType( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );\r
1343         _mainpanel.getOptions().setPhylogenyGraphicsType( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );\r
1344         if ( getMainPanel().getMainFrame() == null ) {\r
1345             // Must be "E" applet version.\r
1346             ( ( ArchaeopteryxE ) ( ( MainPanelApplets ) getMainPanel() ).getApplet() )\r
1347             .setSelectedTypeInTypeMenu( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );\r
1348         }\r
1349         else {\r
1350             getMainPanel().getMainFrame().setSelectedTypeInTypeMenu( PHYLOGENY_GRAPHICS_TYPE.RECTANGULAR );\r
1351         }\r
1352         activateSaveAllIfNeeded();\r
1353         System.gc();\r
1354     }\r
1355 \r
1356     private void obtainDetailedTaxonomicInformation() {\r
1357         if ( getCurrentTreePanel() != null ) {\r
1358             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1359             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1360                 final TaxonomyDataManager t = new TaxonomyDataManager( this,\r
1361                                                                        _mainpanel.getCurrentTreePanel(),\r
1362                                                                        phy.copy(),\r
1363                                                                        false,\r
1364                                                                        true );\r
1365                 new Thread( t ).start();\r
1366             }\r
1367         }\r
1368     }\r
1369 \r
1370     private void obtainDetailedTaxonomicInformationDelete() {\r
1371         if ( getCurrentTreePanel() != null ) {\r
1372             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1373             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1374                 final TaxonomyDataManager t = new TaxonomyDataManager( this,\r
1375                                                                        _mainpanel.getCurrentTreePanel(),\r
1376                                                                        phy.copy(),\r
1377                                                                        true,\r
1378                                                                        true );\r
1379                 new Thread( t ).start();\r
1380             }\r
1381         }\r
1382     }\r
1383 \r
1384     private void obtainSequenceInformation() {\r
1385         if ( getCurrentTreePanel() != null ) {\r
1386             final Phylogeny phy = getCurrentTreePanel().getPhylogeny();\r
1387             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1388                 final SequenceDataRetriver u = new SequenceDataRetriver( this,\r
1389                                                                          _mainpanel.getCurrentTreePanel(),\r
1390                                                                          phy.copy() );\r
1391                 new Thread( u ).start();\r
1392             }\r
1393         }\r
1394     }\r
1395 \r
1396     private void preProcessTreesUponReading( final Phylogeny[] phys ) {\r
1397         for( final Phylogeny phy : phys ) {\r
1398             if ( ( phy != null ) && !phy.isEmpty() ) {\r
1399                 for( final PhylogenyNodeIterator it = phy.iteratorPreorder(); it.hasNext(); ) {\r
1400                     final PhylogenyNode n = it.next();\r
1401                     if ( n.isExternal() ) {\r
1402                         if ( n.getNodeData().isHasSequence() ) {\r
1403                             final Sequence s = n.getNodeData().getSequence();\r
1404                             if ( ForesterUtil.isEmpty( s.getGeneName() ) || s.getGeneName().startsWith( "LOC" ) ) {\r
1405                                 if ( ( s.getAccession() != null )\r
1406                                         && !ForesterUtil.isEmpty( s.getAccession().getValue() ) ) {\r
1407                                     s.setGeneName( s.getAccession().getValue() );\r
1408                                 }\r
1409                                 else if ( !ForesterUtil.isEmpty( n.getName() ) ) {\r
1410                                     s.setGeneName( n.getName() );\r
1411                                 }\r
1412                             }\r
1413                         }\r
1414                     }\r
1415                 }\r
1416             }\r
1417         }\r
1418     }\r
1419 \r
1420     private void readPhylogeniesFromFile() {\r
1421         boolean exception = false;\r
1422         Phylogeny[] phys = null;\r
1423         // Set an initial directory if none set yet\r
1424         final File my_dir = getCurrentDir();\r
1425         // Open file-open dialog and set current directory\r
1426         if ( my_dir != null ) {\r
1427             _open_filechooser.setCurrentDirectory( my_dir );\r
1428         }\r
1429         final int result = _open_filechooser.showOpenDialog( _contentpane );\r
1430         // All done: get the file\r
1431         final File[] files = _open_filechooser.getSelectedFiles();\r
1432         setCurrentDir( _open_filechooser.getCurrentDirectory() );\r
1433         boolean nhx_or_nexus = false;\r
1434         if ( ( files != null ) && ( files.length > 0 ) && ( result == JFileChooser.APPROVE_OPTION ) ) {\r
1435             for( final File file : files ) {\r
1436                 if ( ( file != null ) && !file.isDirectory() ) {\r
1437                     if ( _mainpanel.getCurrentTreePanel() != null ) {\r
1438                         _mainpanel.getCurrentTreePanel().setWaitCursor();\r
1439                     }\r
1440                     else {\r
1441                         _mainpanel.setWaitCursor();\r
1442                     }\r
1443                     if ( ( _open_filechooser.getFileFilter() == MainFrame.nhfilter )\r
1444                             || ( _open_filechooser.getFileFilter() == MainFrame.nhxfilter ) ) {\r
1445                         try {\r
1446                             final NHXParser nhx = new NHXParser();\r
1447                             setSpecialOptionsForNhxParser( nhx );\r
1448                             phys = PhylogenyMethods.readPhylogenies( nhx, file );\r
1449                             nhx_or_nexus = true;\r
1450                         }\r
1451                         catch ( final Exception e ) {\r
1452                             exception = true;\r
1453                             exceptionOccuredDuringOpenFile( e );\r
1454                         }\r
1455                     }\r
1456                     else if ( _open_filechooser.getFileFilter() == MainFrame.xmlfilter ) {\r
1457                         warnIfNotPhyloXmlValidation( getConfiguration() );\r
1458                         try {\r
1459                             final PhyloXmlParser xml_parser = createPhyloXmlParser();\r
1460                             phys = PhylogenyMethods.readPhylogenies( xml_parser, file );\r
1461                         }\r
1462                         catch ( final Exception e ) {\r
1463                             exception = true;\r
1464                             exceptionOccuredDuringOpenFile( e );\r
1465                         }\r
1466                     }\r
1467                     else if ( _open_filechooser.getFileFilter() == MainFrame.tolfilter ) {\r
1468                         try {\r
1469                             phys = PhylogenyMethods.readPhylogenies( new TolParser(), file );\r
1470                         }\r
1471                         catch ( final Exception e ) {\r
1472                             exception = true;\r
1473                             exceptionOccuredDuringOpenFile( e );\r
1474                         }\r
1475                     }\r
1476                     else if ( _open_filechooser.getFileFilter() == MainFrame.nexusfilter ) {\r
1477                         try {\r
1478                             final NexusPhylogeniesParser nex = new NexusPhylogeniesParser();\r
1479                             setSpecialOptionsForNexParser( nex );\r
1480                             phys = PhylogenyMethods.readPhylogenies( nex, file );\r
1481                             nhx_or_nexus = true;\r
1482                         }\r
1483                         catch ( final Exception e ) {\r
1484                             exception = true;\r
1485                             exceptionOccuredDuringOpenFile( e );\r
1486                         }\r
1487                     }\r
1488                     // "*.*":\r
1489                     else {\r
1490                         try {\r
1491                             final PhylogenyParser parser = ParserUtils\r
1492                                     .createParserDependingOnFileType( file, getConfiguration()\r
1493                                                                       .isValidatePhyloXmlAgainstSchema() );\r
1494                             if ( parser instanceof NexusPhylogeniesParser ) {\r
1495                                 final NexusPhylogeniesParser nex = ( NexusPhylogeniesParser ) parser;\r
1496                                 setSpecialOptionsForNexParser( nex );\r
1497                                 nhx_or_nexus = true;\r
1498                             }\r
1499                             else if ( parser instanceof NHXParser ) {\r
1500                                 final NHXParser nhx = ( NHXParser ) parser;\r
1501                                 setSpecialOptionsForNhxParser( nhx );\r
1502                                 nhx_or_nexus = true;\r
1503                             }\r
1504                             else if ( parser instanceof PhyloXmlParser ) {\r
1505                                 warnIfNotPhyloXmlValidation( getConfiguration() );\r
1506                             }\r
1507                             phys = PhylogenyMethods.readPhylogenies( parser, file );\r
1508                         }\r
1509                         catch ( final Exception e ) {\r
1510                             exception = true;\r
1511                             exceptionOccuredDuringOpenFile( e );\r
1512                         }\r
1513                     }\r
1514                     if ( _mainpanel.getCurrentTreePanel() != null ) {\r
1515                         _mainpanel.getCurrentTreePanel().setArrowCursor();\r
1516                     }\r
1517                     else {\r
1518                         _mainpanel.setArrowCursor();\r
1519                     }\r
1520                     if ( !exception && ( phys != null ) && ( phys.length > 0 ) ) {\r
1521                         boolean one_desc = false;\r
1522                         if ( nhx_or_nexus ) {\r
1523                             for( final Phylogeny phy : phys ) {\r
1524                                 if ( getOptions().isInternalNumberAreConfidenceForNhParsing() ) {\r
1525                                     PhylogenyMethods.transferInternalNodeNamesToConfidence( phy, "" );\r
1526                                 }\r
1527                                 if ( PhylogenyMethods.getMinimumDescendentsPerInternalNodes( phy ) == 1 ) {\r
1528                                     one_desc = true;\r
1529                                     break;\r
1530                                 }\r
1531                             }\r
1532                         }\r
1533                         if ( PREPROCESS_TREES ) {\r
1534                             preProcessTreesUponReading( phys );\r
1535                         }\r
1536                         AptxUtil.addPhylogeniesToTabs( phys,\r
1537                                                        file.getName(),\r
1538                                                        file.getAbsolutePath(),\r
1539                                                        getConfiguration(),\r
1540                                                        getMainPanel() );\r
1541                         _mainpanel.getControlPanel().showWhole();\r
1542                         if ( nhx_or_nexus && one_desc ) {\r
1543                             JOptionPane\r
1544                             .showMessageDialog( this,\r
1545                                                 "One or more trees contain (a) node(s) with one descendant, "\r
1546                                                         + ForesterUtil.LINE_SEPARATOR\r
1547                                                         + "possibly indicating illegal parentheses within node names.",\r
1548                                                         "Warning: Possible Error in New Hampshire Formatted Data",\r
1549                                                         JOptionPane.WARNING_MESSAGE );\r
1550                         }\r
1551                     }\r
1552                 }\r
1553             }\r
1554         }\r
1555         activateSaveAllIfNeeded();\r
1556         System.gc();\r
1557     }\r
1558 \r
1559     private void readSpeciesTreeFromFile() {\r
1560         Phylogeny t = null;\r
1561         boolean exception = false;\r
1562         final File my_dir = getCurrentDir();\r
1563         _open_filechooser_for_species_tree.setSelectedFile( new File( "" ) );\r
1564         if ( my_dir != null ) {\r
1565             _open_filechooser_for_species_tree.setCurrentDirectory( my_dir );\r
1566         }\r
1567         final int result = _open_filechooser_for_species_tree.showOpenDialog( _contentpane );\r
1568         final File file = _open_filechooser_for_species_tree.getSelectedFile();\r
1569         if ( ( file != null ) && ( result == JFileChooser.APPROVE_OPTION ) ) {\r
1570             if ( _open_filechooser_for_species_tree.getFileFilter() == MainFrame.xmlfilter ) {\r
1571                 try {\r
1572                     final Phylogeny[] trees = PhylogenyMethods.readPhylogenies( PhyloXmlParser\r
1573                                                                                 .createPhyloXmlParserXsdValidating(), file );\r
1574                     t = trees[ 0 ];\r
1575                 }\r
1576                 catch ( final Exception e ) {\r
1577                     exception = true;\r
1578                     exceptionOccuredDuringOpenFile( e );\r
1579                 }\r
1580             }\r
1581             else if ( _open_filechooser_for_species_tree.getFileFilter() == MainFrame.tolfilter ) {\r
1582                 try {\r
1583                     final Phylogeny[] trees = PhylogenyMethods.readPhylogenies( new TolParser(), file );\r
1584                     t = trees[ 0 ];\r
1585                 }\r
1586                 catch ( final Exception e ) {\r
1587                     exception = true;\r
1588                     exceptionOccuredDuringOpenFile( e );\r
1589                 }\r
1590             }\r
1591             // "*.*":\r
1592             else {\r
1593                 try {\r
1594                     final Phylogeny[] trees = PhylogenyMethods.readPhylogenies( PhyloXmlParser\r
1595                                                                                 .createPhyloXmlParserXsdValidating(), file );\r
1596                     t = trees[ 0 ];\r
1597                 }\r
1598                 catch ( final Exception e ) {\r
1599                     exception = true;\r
1600                     exceptionOccuredDuringOpenFile( e );\r
1601                 }\r
1602             }\r
1603             if ( !exception && ( t != null ) && !t.isRooted() ) {\r
1604                 exception = true;\r
1605                 t = null;\r
1606                 JOptionPane.showMessageDialog( this,\r
1607                                                "Species tree is not rooted",\r
1608                                                "Species tree not loaded",\r
1609                                                JOptionPane.ERROR_MESSAGE );\r
1610             }\r
1611             if ( !exception && ( t != null ) ) {\r
1612                 final Set<Taxonomy> tax_set = new HashSet<Taxonomy>();\r
1613                 for( final PhylogenyNodeIterator it = t.iteratorExternalForward(); it.hasNext(); ) {\r
1614                     final PhylogenyNode node = it.next();\r
1615                     if ( !node.getNodeData().isHasTaxonomy() ) {\r
1616                         exception = true;\r
1617                         t = null;\r
1618                         JOptionPane\r
1619                         .showMessageDialog( this,\r
1620                                             "Species tree contains external node(s) without taxonomy information",\r
1621                                             "Species tree not loaded",\r
1622                                             JOptionPane.ERROR_MESSAGE );\r
1623                         break;\r
1624                     }\r
1625                     else {\r
1626                         if ( tax_set.contains( node.getNodeData().getTaxonomy() ) ) {\r
1627                             exception = true;\r
1628                             t = null;\r
1629                             JOptionPane.showMessageDialog( this,\r
1630                                                            "Taxonomy ["\r
1631                                                                    + node.getNodeData().getTaxonomy().asSimpleText()\r
1632                                                                    + "] is not unique in species tree",\r
1633                                                                    "Species tree not loaded",\r
1634                                                                    JOptionPane.ERROR_MESSAGE );\r
1635                             break;\r
1636                         }\r
1637                         else {\r
1638                             tax_set.add( node.getNodeData().getTaxonomy() );\r
1639                         }\r
1640                     }\r
1641                 }\r
1642             }\r
1643             if ( !exception && ( t != null ) ) {\r
1644                 setSpeciesTree( t );\r
1645                 JOptionPane.showMessageDialog( this,\r
1646                                                "Species tree successfully loaded",\r
1647                                                "Species tree loaded",\r
1648                                                JOptionPane.INFORMATION_MESSAGE );\r
1649             }\r
1650             _contentpane.repaint();\r
1651             System.gc();\r
1652         }\r
1653     }\r
1654 \r
1655     private void setArrowCursor() {\r
1656         try {\r
1657             _mainpanel.getCurrentTreePanel().setArrowCursor();\r
1658         }\r
1659         catch ( final Exception ex ) {\r
1660             // Do nothing.\r
1661         }\r
1662     }\r
1663 \r
1664     private void setMinNotCollapseBlValue( final double min_not_collapse_bl ) {\r
1665         _min_not_collapse_bl = min_not_collapse_bl;\r
1666     }\r
1667 \r
1668     private void setMinNotCollapseConfidenceValue( final double min_not_collapse ) {\r
1669         _min_not_collapse = min_not_collapse;\r
1670     }\r
1671 \r
1672     private void setPhylogeneticInferenceOptions( final PhylogeneticInferenceOptions phylogenetic_inference_options ) {\r
1673         _phylogenetic_inference_options = phylogenetic_inference_options;\r
1674     }\r
1675 \r
1676     private void setSpecialOptionsForNexParser( final NexusPhylogeniesParser nex ) {\r
1677         nex.setReplaceUnderscores( getOptions().isReplaceUnderscoresInNhParsing() );\r
1678         nex.setTaxonomyExtraction( getOptions().getTaxonomyExtraction() );\r
1679     }\r
1680 \r
1681     private void setSpecialOptionsForNhxParser( final NHXParser nhx ) {\r
1682         nhx.setReplaceUnderscores( getOptions().isReplaceUnderscoresInNhParsing() );\r
1683         nhx.setTaxonomyExtraction( getOptions().getTaxonomyExtraction() );\r
1684         nhx.setAllowErrorsInDistanceToParent( getOptions().isAllowErrorsInDistanceToParent() );\r
1685     }\r
1686 \r
1687     void buildAnalysisMenu() {\r
1688         _analysis_menu = MainFrame.createMenu( "Analysis", getConfiguration() );\r
1689         _analysis_menu.add( _gsdi_item = new JMenuItem( "GSDI (Generalized Speciation Duplication Inference)" ) );\r
1690         _analysis_menu.add( _gsdir_item = new JMenuItem( "GSDIR (GSDI with re-rooting)" ) );\r
1691         _analysis_menu.add( _load_species_tree_item = new JMenuItem( "Load Species Tree..." ) );\r
1692         customizeJMenuItem( _gsdi_item );\r
1693         customizeJMenuItem( _gsdir_item );\r
1694         customizeJMenuItem( _load_species_tree_item );\r
1695         _analysis_menu.addSeparator();\r
1696         _analysis_menu.add( _lineage_inference = new JMenuItem( INFER_ANCESTOR_TAXONOMIES ) );\r
1697         customizeJMenuItem( _lineage_inference );\r
1698         _lineage_inference.setToolTipText( "Inference of ancestor taxonomies/lineages" );\r
1699         _jmenubar.add( _analysis_menu );\r
1700     }\r
1701 \r
1702     @Override\r
1703     void buildFileMenu() {\r
1704         _file_jmenu = MainFrame.createMenu( "File", getConfiguration() );\r
1705         _file_jmenu.add( _open_item = new JMenuItem( "Read Tree from File..." ) );\r
1706         _file_jmenu.addSeparator();\r
1707         _file_jmenu.add( _open_url_item = new JMenuItem( "Read Tree from URL/Webservice..." ) );\r
1708         _file_jmenu.addSeparator();\r
1709         final WebservicesManager webservices_manager = WebservicesManager.getInstance();\r
1710         _load_phylogeny_from_webservice_menu_items = new JMenuItem[ webservices_manager\r
1711                                                                     .getAvailablePhylogeniesWebserviceClients().size() ];\r
1712         for( int i = 0; i < webservices_manager.getAvailablePhylogeniesWebserviceClients().size(); ++i ) {\r
1713             final PhylogeniesWebserviceClient client = webservices_manager.getAvailablePhylogeniesWebserviceClient( i );\r
1714             _load_phylogeny_from_webservice_menu_items[ i ] = new JMenuItem( client.getMenuName() );\r
1715             _file_jmenu.add( _load_phylogeny_from_webservice_menu_items[ i ] );\r
1716         }\r
1717         if ( getConfiguration().isEditable() ) {\r
1718             _file_jmenu.addSeparator();\r
1719             _file_jmenu.add( _new_item = new JMenuItem( "New" ) );\r
1720             _new_item.setToolTipText( "to create a new tree with one node, as source for manual tree construction" );\r
1721         }\r
1722         _file_jmenu.addSeparator();\r
1723         _file_jmenu.add( _save_item = new JMenuItem( "Save Tree As..." ) );\r
1724         _file_jmenu.add( _save_all_item = new JMenuItem( "Save All Trees As..." ) );\r
1725         _save_all_item.setToolTipText( "Write all phylogenies to one file." );\r
1726         _save_all_item.setEnabled( false );\r
1727         _file_jmenu.addSeparator();\r
1728         _file_jmenu.add( _write_to_pdf_item = new JMenuItem( "Export to PDF file ..." ) );\r
1729         if ( AptxUtil.canWriteFormat( "tif" ) || AptxUtil.canWriteFormat( "tiff" ) || AptxUtil.canWriteFormat( "TIF" ) ) {\r
1730             _file_jmenu.add( _write_to_tif_item = new JMenuItem( "Export to TIFF file..." ) );\r
1731         }\r
1732         _file_jmenu.add( _write_to_png_item = new JMenuItem( "Export to PNG file..." ) );\r
1733         _file_jmenu.add( _write_to_jpg_item = new JMenuItem( "Export to JPG file..." ) );\r
1734         if ( AptxUtil.canWriteFormat( "gif" ) ) {\r
1735             _file_jmenu.add( _write_to_gif_item = new JMenuItem( "Export to GIF file..." ) );\r
1736         }\r
1737         if ( AptxUtil.canWriteFormat( "bmp" ) ) {\r
1738             _file_jmenu.add( _write_to_bmp_item = new JMenuItem( "Export to BMP file..." ) );\r
1739         }\r
1740         _file_jmenu.addSeparator();\r
1741         _file_jmenu.add( _print_item = new JMenuItem( "Print..." ) );\r
1742         _file_jmenu.addSeparator();\r
1743         _file_jmenu.add( _close_item = new JMenuItem( "Close Tab" ) );\r
1744         _close_item.setToolTipText( "To close the current pane." );\r
1745         _close_item.setEnabled( true );\r
1746         _file_jmenu.addSeparator();\r
1747         _file_jmenu.add( _exit_item = new JMenuItem( "Exit" ) );\r
1748         customizeJMenuItem( _open_item );\r
1749         _open_item\r
1750         .setFont( new Font( _open_item.getFont().getFontName(), Font.BOLD, _open_item.getFont().getSize() + 4 ) );\r
1751         customizeJMenuItem( _open_url_item );\r
1752         for( int i = 0; i < webservices_manager.getAvailablePhylogeniesWebserviceClients().size(); ++i ) {\r
1753             customizeJMenuItem( _load_phylogeny_from_webservice_menu_items[ i ] );\r
1754         }\r
1755         customizeJMenuItem( _save_item );\r
1756         if ( getConfiguration().isEditable() ) {\r
1757             customizeJMenuItem( _new_item );\r
1758         }\r
1759         customizeJMenuItem( _close_item );\r
1760         customizeJMenuItem( _save_all_item );\r
1761         customizeJMenuItem( _write_to_pdf_item );\r
1762         customizeJMenuItem( _write_to_png_item );\r
1763         customizeJMenuItem( _write_to_jpg_item );\r
1764         customizeJMenuItem( _write_to_gif_item );\r
1765         customizeJMenuItem( _write_to_tif_item );\r
1766         customizeJMenuItem( _write_to_bmp_item );\r
1767         customizeJMenuItem( _print_item );\r
1768         customizeJMenuItem( _exit_item );\r
1769         _jmenubar.add( _file_jmenu );\r
1770     }\r
1771 \r
1772     void buildOptionsMenu() {\r
1773         _options_jmenu = MainFrame.createMenu( OPTIONS_HEADER, getConfiguration() );\r
1774         _options_jmenu.addChangeListener( new ChangeListener() {\r
1775 \r
1776             @Override\r
1777             public void stateChanged( final ChangeEvent e ) {\r
1778                 MainFrame.setOvPlacementColorChooseMenuItem( _overview_placment_mi, getOptions() );\r
1779                 MainFrame.setTextColorChooseMenuItem( _switch_colors_mi, getCurrentTreePanel() );\r
1780                 MainFrame\r
1781                 .setTextMinSupportMenuItem( _choose_minimal_confidence_mi, getOptions(), getCurrentTreePanel() );\r
1782                 MainFrame.setTextForFontChooserMenuItem( _choose_font_mi, MainFrame\r
1783                                                          .createCurrentFontDesc( getMainPanel().getTreeFontSet() ) );\r
1784                 MainFrame.setTextForGraphicsSizeChooserMenuItem( _print_size_mi, getOptions() );\r
1785                 MainFrame.setTextForPdfLineWidthChooserMenuItem( _choose_pdf_width_mi, getOptions() );\r
1786                 MainFrame.setCycleNodeFillMenuItem( _cycle_node_fill_mi, getOptions() );\r
1787                 MainFrame.setCycleNodeShapeMenuItem( _cycle_node_shape_mi, getOptions() );\r
1788                 MainFrame.setCycleDataReturnMenuItem( _cycle_data_return, getOptions() );\r
1789                 MainFrame.setTextNodeSizeMenuItem( _choose_node_size_mi, getOptions() );\r
1790                 try {\r
1791                     getMainPanel().getControlPanel().setVisibilityOfDomainStrucureCB();\r
1792                     getMainPanel().getControlPanel().setVisibilityOfX();\r
1793                 }\r
1794                 catch ( final Exception ignore ) {\r
1795                     // do nothing, not important.\r
1796                 }\r
1797             }\r
1798         } );\r
1799         _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( DISPLAY_SUBHEADER ), getConfiguration() ) );\r
1800         _options_jmenu\r
1801         .add( _ext_node_dependent_cladogram_rbmi = new JRadioButtonMenuItem( MainFrame.NONUNIFORM_CLADOGRAMS_LABEL ) );\r
1802         _options_jmenu.add( _uniform_cladograms_rbmi = new JRadioButtonMenuItem( MainFrame.UNIFORM_CLADOGRAMS_LABEL ) );\r
1803         _options_jmenu.add( _non_lined_up_cladograms_rbmi = new JRadioButtonMenuItem( NON_LINED_UP_CLADOGRAMS_LABEL ) );\r
1804         _radio_group_1 = new ButtonGroup();\r
1805         _radio_group_1.add( _ext_node_dependent_cladogram_rbmi );\r
1806         _radio_group_1.add( _uniform_cladograms_rbmi );\r
1807         _radio_group_1.add( _non_lined_up_cladograms_rbmi );\r
1808         _options_jmenu.add( _show_overview_cbmi = new JCheckBoxMenuItem( SHOW_OVERVIEW_LABEL ) );\r
1809         _options_jmenu.add( _show_scale_cbmi = new JCheckBoxMenuItem( DISPLAY_SCALE_LABEL ) );\r
1810         _options_jmenu\r
1811         .add( _show_default_node_shapes_internal_cbmi = new JCheckBoxMenuItem( DISPLAY_NODE_BOXES_LABEL_INT ) );\r
1812         _options_jmenu\r
1813         .add( _show_default_node_shapes_external_cbmi = new JCheckBoxMenuItem( DISPLAY_NODE_BOXES_LABEL_EXT ) );\r
1814         _options_jmenu\r
1815         .add( _show_default_node_shapes_for_marked_cbmi = new JCheckBoxMenuItem( MainFrame.DISPLAY_NODE_BOXES_LABEL_MARKED ) );\r
1816         _options_jmenu.add( _line_up_renderable_data_cbmi = new JCheckBoxMenuItem( MainFrame.LINE_UP_RENDERABLE_DATA ) );\r
1817         if ( getConfiguration().doDisplayOption( Configuration.show_domain_architectures ) ) {\r
1818             _options_jmenu.add( _right_line_up_domains_cbmi = new JCheckBoxMenuItem( MainFrame.RIGHT_LINE_UP_DOMAINS ) );\r
1819             _options_jmenu.add( _show_domain_labels = new JCheckBoxMenuItem( MainFrame.SHOW_DOMAIN_LABELS_LABEL ) );\r
1820         }\r
1821         _options_jmenu.add( _show_annotation_ref_source = new JCheckBoxMenuItem( SHOW_ANN_REF_SOURCE_LABEL ) );\r
1822         _options_jmenu.add( _show_confidence_stddev_cbmi = new JCheckBoxMenuItem( SHOW_CONF_STDDEV_LABEL ) );\r
1823         _options_jmenu.add( _color_by_taxonomic_group_cbmi = new JCheckBoxMenuItem( COLOR_BY_TAXONOMIC_GROUP ) );\r
1824         _options_jmenu.add( _color_labels_same_as_parent_branch = new JCheckBoxMenuItem( COLOR_LABELS_LABEL ) );\r
1825         _color_labels_same_as_parent_branch.setToolTipText( MainFrame.COLOR_LABELS_TIP );\r
1826         _options_jmenu.add( _abbreviate_scientific_names = new JCheckBoxMenuItem( ABBREV_SN_LABEL ) );\r
1827         _options_jmenu.add( _label_direction_cbmi = new JCheckBoxMenuItem( LABEL_DIRECTION_LABEL ) );\r
1828         _label_direction_cbmi.setToolTipText( LABEL_DIRECTION_TIP );\r
1829         _options_jmenu.add( _screen_antialias_cbmi = new JCheckBoxMenuItem( SCREEN_ANTIALIAS_LABEL ) );\r
1830         _options_jmenu.add( _background_gradient_cbmi = new JCheckBoxMenuItem( BG_GRAD_LABEL ) );\r
1831         _options_jmenu.add( _cycle_node_shape_mi = new JMenuItem( MainFrame.CYCLE_NODE_SHAPE_LABEL ) );\r
1832         _options_jmenu.add( _cycle_node_fill_mi = new JMenuItem( MainFrame.CYCLE_NODE_FILL_LABEL ) );\r
1833         _options_jmenu.add( _choose_node_size_mi = new JMenuItem( MainFrame.CHOOSE_NODE_SIZE_LABEL ) );\r
1834         _options_jmenu.add( _choose_minimal_confidence_mi = new JMenuItem( "" ) );\r
1835         _options_jmenu.add( _overview_placment_mi = new JMenuItem( "" ) );\r
1836         _options_jmenu.add( _switch_colors_mi = new JMenuItem( "" ) );\r
1837         _options_jmenu.add( _choose_font_mi = new JMenuItem( "" ) );\r
1838         _options_jmenu.addSeparator();\r
1839         _options_jmenu.add( _cycle_data_return = new JMenuItem( "Cycle Data Return" ) );\r
1840         _options_jmenu.addSeparator();\r
1841         _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( SEARCH_SUBHEADER ), getConfiguration() ) );\r
1842         _options_jmenu.add( _search_case_senstive_cbmi = new JCheckBoxMenuItem( SEARCH_CASE_SENSITIVE_LABEL ) );\r
1843         _options_jmenu.add( _search_whole_words_only_cbmi = new JCheckBoxMenuItem( SEARCH_TERMS_ONLY_LABEL ) );\r
1844         _options_jmenu.add( _search_with_regex_cbmi = new JCheckBoxMenuItem( MainFrame.SEARCH_REGEX_LABEL ) );\r
1845         _search_with_regex_cbmi.setToolTipText( MainFrame.SEARCH_WITH_REGEX_TIP );\r
1846         _options_jmenu.add( _inverse_search_result_cbmi = new JCheckBoxMenuItem( INVERSE_SEARCH_RESULT_LABEL ) );\r
1847         _options_jmenu.addSeparator();\r
1848         _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( "Graphics Export & Printing:" ),\r
1849                                                       getConfiguration() ) );\r
1850         _options_jmenu.add( _antialias_print_cbmi = new JCheckBoxMenuItem( "Antialias" ) );\r
1851         _options_jmenu.add( _print_black_and_white_cbmi = new JCheckBoxMenuItem( "Export in Black and White" ) );\r
1852         _options_jmenu\r
1853         .add( _print_using_actual_size_cbmi = new JCheckBoxMenuItem( "Use Current Image Size for PDF export and Printing" ) );\r
1854         _options_jmenu\r
1855         .add( _graphics_export_using_actual_size_cbmi = new JCheckBoxMenuItem( "Use Current Image Size for PNG, JPG, and GIF export" ) );\r
1856         _options_jmenu\r
1857         .add( _graphics_export_visible_only_cbmi = new JCheckBoxMenuItem( "Limit to Visible ('Screenshot') for PNG, JPG, and GIF export" ) );\r
1858         _options_jmenu.add( _print_size_mi = new JMenuItem( "" ) );\r
1859         _options_jmenu.add( _choose_pdf_width_mi = new JMenuItem( "" ) );\r
1860         _options_jmenu.addSeparator();\r
1861         _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( "Newick/NHX/Nexus Read:" ), getConfiguration() ) );\r
1862         _options_jmenu\r
1863         .add( _internal_number_are_confidence_for_nh_parsing_cbmi = new JCheckBoxMenuItem( "Internal Node Names are Confidence Values" ) );\r
1864         _options_jmenu.add( _replace_underscores_cbmi = new JCheckBoxMenuItem( "Replace Underscores with Spaces" ) );\r
1865         _options_jmenu\r
1866         .add( _allow_errors_in_distance_to_parent_cbmi = new JCheckBoxMenuItem( "Ignore Distance Values Format Errors" ) );\r
1867         _options_jmenu.add( _extract_taxonomy_no_rbmi = new JRadioButtonMenuItem( "No Taxonomy Extraction" ) );\r
1868         _options_jmenu\r
1869         .add( _extract_taxonomy_pfam_strict_rbmi = new JRadioButtonMenuItem( "Extract Taxonomy Codes/Ids from Pfam-style Node Names" ) );\r
1870         _options_jmenu\r
1871         .add( _extract_taxonomy_pfam_relaxed_rbmi = new JRadioButtonMenuItem( "Extract Taxonomy Codes/Ids from Pfam-style like Node Names" ) );\r
1872         _options_jmenu\r
1873         .add( _extract_taxonomy_agressive_rbmi = new JRadioButtonMenuItem( "Extract Taxonomy Codes/Ids/Scientific Names from Node Names" ) );\r
1874         _extract_taxonomy_pfam_strict_rbmi\r
1875         .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
1876         _extract_taxonomy_pfam_relaxed_rbmi\r
1877         .setToolTipText( "To extract taxonomy codes/ids from node names in the form of e.g. \"bax_MOUSE\" or \"bax_10090\"" );\r
1878         _extract_taxonomy_agressive_rbmi\r
1879         .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
1880         _radio_group_2 = new ButtonGroup();\r
1881         _radio_group_2.add( _extract_taxonomy_no_rbmi );\r
1882         _radio_group_2.add( _extract_taxonomy_pfam_strict_rbmi );\r
1883         _radio_group_2.add( _extract_taxonomy_pfam_relaxed_rbmi );\r
1884         _radio_group_2.add( _extract_taxonomy_agressive_rbmi );\r
1885         _options_jmenu.add( customizeMenuItemAsLabel( new JMenuItem( "Newick/Nexus Save:" ), getConfiguration() ) );\r
1886         _options_jmenu\r
1887         .add( _use_brackets_for_conf_in_nh_export_cbmi = new JCheckBoxMenuItem( USE_BRACKETS_FOR_CONF_IN_NH_LABEL ) );\r
1888         _use_brackets_for_conf_in_nh_export_cbmi\r
1889         .setToolTipText( "e.g. \"0.1[90]\" for a branch with support 90 and a length of 0.1" );\r
1890         _options_jmenu\r
1891         .add( _use_internal_names_for_conf_in_nh_export_cbmi = new JCheckBoxMenuItem( USE_INTERNAL_NAMES_FOR_CONF_IN_NH_LABEL ) );\r
1892         customizeJMenuItem( _choose_font_mi );\r
1893         customizeJMenuItem( _choose_minimal_confidence_mi );\r
1894         customizeJMenuItem( _switch_colors_mi );\r
1895         customizeJMenuItem( _print_size_mi );\r
1896         customizeJMenuItem( _choose_pdf_width_mi );\r
1897         customizeJMenuItem( _overview_placment_mi );\r
1898         customizeCheckBoxMenuItem( _show_default_node_shapes_external_cbmi, getOptions()\r
1899                                    .isShowDefaultNodeShapesExternal() );\r
1900         customizeCheckBoxMenuItem( _show_default_node_shapes_internal_cbmi, getOptions()\r
1901                                    .isShowDefaultNodeShapesInternal() );\r
1902         customizeCheckBoxMenuItem( _show_default_node_shapes_for_marked_cbmi, getOptions()\r
1903                                    .isShowDefaultNodeShapesForMarkedNodes() );\r
1904         customizeJMenuItem( _cycle_node_shape_mi );\r
1905         customizeJMenuItem( _cycle_node_fill_mi );\r
1906         customizeJMenuItem( _choose_node_size_mi );\r
1907         customizeJMenuItem( _cycle_data_return );\r
1908         customizeCheckBoxMenuItem( _color_labels_same_as_parent_branch, getOptions().isColorLabelsSameAsParentBranch() );\r
1909         customizeCheckBoxMenuItem( _color_by_taxonomic_group_cbmi, getOptions().isColorByTaxonomicGroup() );\r
1910         customizeCheckBoxMenuItem( _screen_antialias_cbmi, getOptions().isAntialiasScreen() );\r
1911         customizeCheckBoxMenuItem( _background_gradient_cbmi, getOptions().isBackgroundColorGradient() );\r
1912         customizeCheckBoxMenuItem( _show_domain_labels, getOptions().isShowDomainLabels() );\r
1913         customizeCheckBoxMenuItem( _show_annotation_ref_source, getOptions().isShowAnnotationRefSource() );\r
1914         customizeCheckBoxMenuItem( _abbreviate_scientific_names, getOptions().isAbbreviateScientificTaxonNames() );\r
1915         customizeCheckBoxMenuItem( _search_case_senstive_cbmi, getOptions().isSearchCaseSensitive() );\r
1916         customizeCheckBoxMenuItem( _show_scale_cbmi, getOptions().isShowScale() );\r
1917         customizeRadioButtonMenuItem( _non_lined_up_cladograms_rbmi,\r
1918                                       getOptions().getCladogramType() == CLADOGRAM_TYPE.NON_LINED_UP );\r
1919         customizeRadioButtonMenuItem( _uniform_cladograms_rbmi,\r
1920                                       getOptions().getCladogramType() == CLADOGRAM_TYPE.TOTAL_NODE_SUM_DEP );\r
1921         customizeRadioButtonMenuItem( _ext_node_dependent_cladogram_rbmi,\r
1922                                       getOptions().getCladogramType() == CLADOGRAM_TYPE.EXT_NODE_SUM_DEP );\r
1923         customizeCheckBoxMenuItem( _show_overview_cbmi, getOptions().isShowOverview() );\r
1924         customizeCheckBoxMenuItem( _label_direction_cbmi,\r
1925                                    getOptions().getNodeLabelDirection() == NODE_LABEL_DIRECTION.RADIAL );\r
1926         customizeCheckBoxMenuItem( _antialias_print_cbmi, getOptions().isAntialiasPrint() );\r
1927         customizeCheckBoxMenuItem( _print_black_and_white_cbmi, getOptions().isPrintBlackAndWhite() );\r
1928         customizeCheckBoxMenuItem( _internal_number_are_confidence_for_nh_parsing_cbmi, getOptions()\r
1929                                    .isInternalNumberAreConfidenceForNhParsing() );\r
1930         customizeRadioButtonMenuItem( _extract_taxonomy_no_rbmi,\r
1931                                       getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.NO );\r
1932         customizeRadioButtonMenuItem( _extract_taxonomy_pfam_strict_rbmi,\r
1933                                       getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.PFAM_STYLE_STRICT );\r
1934         customizeRadioButtonMenuItem( _extract_taxonomy_pfam_relaxed_rbmi,\r
1935                                       getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.PFAM_STYLE_RELAXED );\r
1936         customizeRadioButtonMenuItem( _extract_taxonomy_agressive_rbmi,\r
1937                                       getOptions().getTaxonomyExtraction() == TAXONOMY_EXTRACTION.AGGRESSIVE );\r
1938         customizeCheckBoxMenuItem( _replace_underscores_cbmi, getOptions().isReplaceUnderscoresInNhParsing() );\r
1939         customizeCheckBoxMenuItem( _allow_errors_in_distance_to_parent_cbmi, getOptions()\r
1940                                    .isReplaceUnderscoresInNhParsing() );\r
1941         customizeCheckBoxMenuItem( _search_with_regex_cbmi, getOptions().isSearchWithRegex() );\r
1942         customizeCheckBoxMenuItem( _search_whole_words_only_cbmi, getOptions().isMatchWholeTermsOnly() );\r
1943         customizeCheckBoxMenuItem( _inverse_search_result_cbmi, getOptions().isInverseSearchResult() );\r
1944         customizeCheckBoxMenuItem( _graphics_export_visible_only_cbmi, getOptions().isGraphicsExportVisibleOnly() );\r
1945         customizeCheckBoxMenuItem( _print_using_actual_size_cbmi, getOptions().isPrintUsingActualSize() );\r
1946         customizeCheckBoxMenuItem( _graphics_export_using_actual_size_cbmi, getOptions()\r
1947                                    .isGraphicsExportUsingActualSize() );\r
1948         customizeCheckBoxMenuItem( _show_confidence_stddev_cbmi, getOptions().isShowConfidenceStddev() );\r
1949         customizeCheckBoxMenuItem( _use_brackets_for_conf_in_nh_export_cbmi, getOptions()\r
1950                                    .getNhConversionSupportValueStyle() == NH_CONVERSION_SUPPORT_VALUE_STYLE.IN_SQUARE_BRACKETS );\r
1951         customizeCheckBoxMenuItem( _use_internal_names_for_conf_in_nh_export_cbmi, getOptions()\r
1952                                    .getNhConversionSupportValueStyle() == NH_CONVERSION_SUPPORT_VALUE_STYLE.AS_INTERNAL_NODE_NAMES );\r
1953         customizeCheckBoxMenuItem( _line_up_renderable_data_cbmi, getOptions().isLineUpRendarableNodeData() );\r
1954         customizeCheckBoxMenuItem( _right_line_up_domains_cbmi, getOptions().isRightLineUpDomains() );\r
1955         _jmenubar.add( _options_jmenu );\r
1956     }\r
1957 \r
1958     void buildPhylogeneticInferenceMenu() {\r
1959         final InferenceManager im = getInferenceManager();\r
1960         _inference_menu = MainFrame.createMenu( "Inference", getConfiguration() );\r
1961         _inference_menu.add( _inference_from_msa_item = new JMenuItem( "From Multiple Sequence Alignment..." ) );\r
1962         customizeJMenuItem( _inference_from_msa_item );\r
1963         _inference_from_msa_item.setToolTipText( "Basic phylogenetic inference from MSA" );\r
1964         if ( im.canDoMsa() ) {\r
1965             _inference_menu.add( _inference_from_seqs_item = new JMenuItem( "From Unaligned Sequences..." ) );\r
1966             customizeJMenuItem( _inference_from_seqs_item );\r
1967             _inference_from_seqs_item\r
1968             .setToolTipText( "Basic phylogenetic inference including multiple sequence alignment" );\r
1969         }\r
1970         else {\r
1971             _inference_menu\r
1972             .add( _inference_from_seqs_item = new JMenuItem( "From Unaligned Sequences (no program found)" ) );\r
1973             customizeJMenuItem( _inference_from_seqs_item );\r
1974             _inference_from_seqs_item.setEnabled( false );\r
1975         }\r
1976         _jmenubar.add( _inference_menu );\r
1977     }\r
1978 \r
1979     void buildToolsMenu() {\r
1980         _tools_menu = createMenu( "Tools", getConfiguration() );\r
1981         _tools_menu.add( _confcolor_item = new JMenuItem( "Colorize Branches Depending on Confidence" ) );\r
1982         customizeJMenuItem( _confcolor_item );\r
1983         _tools_menu.add( _color_rank_jmi = new JMenuItem( "Colorize Subtrees via Taxonomic Rank" ) );\r
1984         customizeJMenuItem( _color_rank_jmi );\r
1985         _color_rank_jmi.setToolTipText( "for example, at \"Class\" level, colorize mammal specific subtree red" );\r
1986         _tools_menu.add( _taxcolor_item = new JMenuItem( "Taxonomy Colorize Branches" ) );\r
1987         customizeJMenuItem( _taxcolor_item );\r
1988         _tools_menu.addSeparator();\r
1989         _tools_menu.add( _remove_visual_styles_item = new JMenuItem( "Delete All Visual Styles From Nodes" ) );\r
1990         _remove_visual_styles_item\r
1991         .setToolTipText( "To remove all node visual styles (fonts, colors) from the current phylogeny" );\r
1992         customizeJMenuItem( _remove_visual_styles_item );\r
1993         _tools_menu.add( _remove_branch_color_item = new JMenuItem( "Delete All Colors From Branches" ) );\r
1994         _remove_branch_color_item.setToolTipText( "To remove all branch color values from the current phylogeny" );\r
1995         customizeJMenuItem( _remove_branch_color_item );\r
1996         _tools_menu.addSeparator();\r
1997         _tools_menu.add( _annotate_item = new JMenuItem( "Annotate Sequences of Selected Nodes" ) );\r
1998         customizeJMenuItem( _annotate_item );\r
1999         _tools_menu.addSeparator();\r
2000         _tools_menu.add( _midpoint_root_item = new JMenuItem( "Midpoint-Root" ) );\r
2001         customizeJMenuItem( _midpoint_root_item );\r
2002         _tools_menu.addSeparator();\r
2003         _tools_menu.add( _delete_selected_nodes_item = new JMenuItem( "Delete Selected Nodes" ) );\r
2004         _delete_selected_nodes_item.setToolTipText( "To delete all selected external nodes" );\r
2005         customizeJMenuItem( _delete_selected_nodes_item );\r
2006         _tools_menu.add( _delete_not_selected_nodes_item = new JMenuItem( "Retain Selected Nodes" ) );\r
2007         _delete_not_selected_nodes_item.setToolTipText( "To delete all not selected external nodes" );\r
2008         customizeJMenuItem( _delete_not_selected_nodes_item );\r
2009         _tools_menu.addSeparator();\r
2010         _tools_menu.add( _collapse_species_specific_subtrees = new JMenuItem( "Collapse Species-Specific Subtrees" ) );\r
2011         customizeJMenuItem( _collapse_species_specific_subtrees );\r
2012         _collapse_species_specific_subtrees.setToolTipText( "To (reversibly) collapse species-specific subtrees" );\r
2013         _tools_menu\r
2014         .add( _collapse_below_threshold = new JMenuItem( "Collapse Branches with Confidence Below Threshold into Multifurcations" ) );\r
2015         customizeJMenuItem( _collapse_below_threshold );\r
2016         _collapse_below_threshold\r
2017         .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
2018         //\r
2019         _tools_menu\r
2020         .add( _collapse_below_branch_length = new JMenuItem( "Collapse Branches with Branch Lengths Below Threshold into Multifurcations" ) );\r
2021         customizeJMenuItem( _collapse_below_branch_length );\r
2022         _collapse_below_branch_length\r
2023         .setToolTipText( "To (permanently) collapse branches with branches with branch lengths below a threshold into multifurcations" );\r
2024         //\r
2025         _tools_menu.addSeparator();\r
2026         _tools_menu\r
2027         .add( _extract_tax_code_from_node_names_jmi = new JMenuItem( "Extract Taxonomic Data from Node Names" ) );\r
2028         customizeJMenuItem( _extract_tax_code_from_node_names_jmi );\r
2029         _extract_tax_code_from_node_names_jmi\r
2030         .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
2031         _tools_menu\r
2032         .add( _move_node_names_to_tax_sn_jmi = new JMenuItem( "Transfer Node Names to Taxonomic Scientific Names" ) );\r
2033         customizeJMenuItem( _move_node_names_to_tax_sn_jmi );\r
2034         _move_node_names_to_tax_sn_jmi.setToolTipText( "To interpret node names as taxonomic scientific names" );\r
2035         _tools_menu.add( _move_node_names_to_seq_names_jmi = new JMenuItem( "Transfer Node Names to Sequence Names" ) );\r
2036         customizeJMenuItem( _move_node_names_to_seq_names_jmi );\r
2037         _move_node_names_to_seq_names_jmi.setToolTipText( "To interpret node names as sequence (protein, gene) names" );\r
2038         _tools_menu.addSeparator();\r
2039         _tools_menu.add( _obtain_seq_information_jmi = new JMenuItem( "Obtain Sequence Information" ) );\r
2040         customizeJMenuItem( _obtain_seq_information_jmi );\r
2041         _obtain_seq_information_jmi.setToolTipText( "To add additional sequence information" );\r
2042         _tools_menu\r
2043         .add( _obtain_detailed_taxonomic_information_jmi = new JMenuItem( OBTAIN_DETAILED_TAXONOMIC_INFORMATION ) );\r
2044         customizeJMenuItem( _obtain_detailed_taxonomic_information_jmi );\r
2045         _obtain_detailed_taxonomic_information_jmi\r
2046         .setToolTipText( "To add additional taxonomic information (from UniProt Taxonomy)" );\r
2047         _tools_menu\r
2048         .add( _obtain_detailed_taxonomic_information_deleting_jmi = new JMenuItem( "Obtain Detailed Taxonomic Information (deletes nodes!)" ) );\r
2049         customizeJMenuItem( _obtain_detailed_taxonomic_information_deleting_jmi );\r
2050         _obtain_detailed_taxonomic_information_deleting_jmi\r
2051         .setToolTipText( "To add additional taxonomic information, deletes nodes for which taxonomy cannot found (from UniProt Taxonomy)" );\r
2052         _tools_menu.addSeparator();\r
2053         _tools_menu.add( _read_values_jmi = new JMenuItem( "Attach Vector/Expression Values" ) );\r
2054         customizeJMenuItem( _read_values_jmi );\r
2055         _read_values_jmi.setToolTipText( "To attach vector (e.g. gene expression) values to tree nodes (beta)" );\r
2056         _jmenubar.add( _tools_menu );\r
2057         _tools_menu.add( _read_seqs_jmi = new JMenuItem( "Attach Molecular Sequences" ) );\r
2058         customizeJMenuItem( _read_seqs_jmi );\r
2059         _read_seqs_jmi\r
2060         .setToolTipText( "To attach molecular sequences to tree nodes (from Fasta-formatted file) (beta)" );\r
2061         _jmenubar.add( _tools_menu );\r
2062     }\r
2063 \r
2064     @Override\r
2065     void close() {\r
2066         if ( isUnsavedDataPresent() ) {\r
2067             final int r = JOptionPane.showConfirmDialog( this,\r
2068                                                          "Exit despite potentially unsaved changes?",\r
2069                                                          "Exit?",\r
2070                                                          JOptionPane.YES_NO_OPTION );\r
2071             if ( r != JOptionPane.YES_OPTION ) {\r
2072                 return;\r
2073             }\r
2074         }\r
2075         exit();\r
2076     }\r
2077 \r
2078     void exit() {\r
2079         removeAllTextFrames();\r
2080         _mainpanel.terminate();\r
2081         _contentpane.removeAll();\r
2082         setVisible( false );\r
2083         dispose();\r
2084         // System.exit( 0 ); //TODO reconfirm that this is OK, then remove.\r
2085     }\r
2086 \r
2087     void readPhylogeniesFromURL() {\r
2088         URL url = null;\r
2089         Phylogeny[] phys = null;\r
2090         final String message = "Please enter a complete URL, for example \"http://purl.org/phylo/treebase/phylows/study/TB2:S15480?format=nexus\"";\r
2091         final String url_string = JOptionPane.showInputDialog( this,\r
2092                                                                message,\r
2093                                                                "Use URL/webservice to obtain a phylogeny",\r
2094                                                                JOptionPane.QUESTION_MESSAGE );\r
2095         boolean nhx_or_nexus = false;\r
2096         if ( ( url_string != null ) && ( url_string.length() > 0 ) ) {\r
2097             try {\r
2098                 url = new URL( url_string );\r
2099                 PhylogenyParser parser = null;\r
2100                 if ( url.getHost().toLowerCase().indexOf( "tolweb" ) >= 0 ) {\r
2101                     parser = new TolParser();\r
2102                 }\r
2103                 else {\r
2104                     parser = ParserUtils.createParserDependingOnUrlContents( url, getConfiguration()\r
2105                                                                              .isValidatePhyloXmlAgainstSchema() );\r
2106                 }\r
2107                 if ( parser instanceof NexusPhylogeniesParser ) {\r
2108                     nhx_or_nexus = true;\r
2109                 }\r
2110                 else if ( parser instanceof NHXParser ) {\r
2111                     nhx_or_nexus = true;\r
2112                 }\r
2113                 if ( _mainpanel.getCurrentTreePanel() != null ) {\r
2114                     _mainpanel.getCurrentTreePanel().setWaitCursor();\r
2115                 }\r
2116                 else {\r
2117                     _mainpanel.setWaitCursor();\r
2118                 }\r
2119                 final PhylogenyFactory factory = ParserBasedPhylogenyFactory.getInstance();\r
2120                 phys = factory.create( url.openStream(), parser );\r
2121             }\r
2122             catch ( final MalformedURLException e ) {\r
2123                 JOptionPane.showMessageDialog( this,\r
2124                                                "Malformed URL: " + url + "\n" + e.getLocalizedMessage(),\r
2125                                                "Malformed URL",\r
2126                                                JOptionPane.ERROR_MESSAGE );\r
2127             }\r
2128             catch ( final IOException e ) {\r
2129                 JOptionPane.showMessageDialog( this,\r
2130                                                "Could not read from " + url + "\n"\r
2131                                                        + ForesterUtil.wordWrap( e.getLocalizedMessage(), 80 ),\r
2132                                                        "Failed to read URL",\r
2133                                                        JOptionPane.ERROR_MESSAGE );\r
2134             }\r
2135             catch ( final Exception e ) {\r
2136                 JOptionPane.showMessageDialog( this,\r
2137                                                ForesterUtil.wordWrap( e.getLocalizedMessage(), 80 ),\r
2138                                                "Unexpected Exception",\r
2139                                                JOptionPane.ERROR_MESSAGE );\r
2140             }\r
2141             finally {\r
2142                 if ( _mainpanel.getCurrentTreePanel() != null ) {\r
2143                     _mainpanel.getCurrentTreePanel().setArrowCursor();\r
2144                 }\r
2145                 else {\r
2146                     _mainpanel.setArrowCursor();\r
2147                 }\r
2148             }\r
2149             if ( ( phys != null ) && ( phys.length > 0 ) ) {\r
2150                 if ( nhx_or_nexus && getOptions().isInternalNumberAreConfidenceForNhParsing() ) {\r
2151                     for( final Phylogeny phy : phys ) {\r
2152                         PhylogenyMethods.transferInternalNodeNamesToConfidence( phy, "" );\r
2153                     }\r
2154                 }\r
2155                 AptxUtil.addPhylogeniesToTabs( phys,\r
2156                                                new File( url.getFile() ).getName(),\r
2157                                                new File( url.getFile() ).toString(),\r
2158                                                getConfiguration(),\r
2159                                                getMainPanel() );\r
2160                 _mainpanel.getControlPanel().showWhole();\r
2161             }\r
2162         }\r
2163         activateSaveAllIfNeeded();\r
2164         System.gc();\r
2165     }\r
2166 \r
2167     void setMsa( final Msa msa ) {\r
2168         _msa = msa;\r
2169     }\r
2170 \r
2171     void setMsaFile( final File msa_file ) {\r
2172         _msa_file = msa_file;\r
2173     }\r
2174 \r
2175     void setSeqs( final List<MolecularSequence> seqs ) {\r
2176         _seqs = seqs;\r
2177     }\r
2178 \r
2179     void setSeqsFile( final File seqs_file ) {\r
2180         _seqs_file = seqs_file;\r
2181     }\r
2182 \r
2183     public static MainFrameApplication createInstance( final Phylogeny[] phys, final Configuration config ) {\r
2184         return new MainFrameApplication( phys, config );\r
2185     }\r
2186 \r
2187     public static MainFrame createInstance( final Phylogeny[] phys,\r
2188                                             final Configuration config,\r
2189                                             final String title,\r
2190                                             final File current_dir ) {\r
2191         return new MainFrameApplication( phys, config, title, current_dir );\r
2192     }\r
2193 \r
2194     static MainFrame createInstance( final Phylogeny[] phys, final Configuration config, final String title ) {\r
2195         return new MainFrameApplication( phys, config, title );\r
2196     }\r
2197 \r
2198     static MainFrame createInstance( final Phylogeny[] phys, final String config_file_name, final String title ) {\r
2199         return new MainFrameApplication( phys, config_file_name, title );\r
2200     }\r
2201 \r
2202     static void warnIfNotPhyloXmlValidation( final Configuration c ) {\r
2203         if ( !c.isValidatePhyloXmlAgainstSchema() ) {\r
2204             JOptionPane\r
2205             .showMessageDialog( null,\r
2206                                 ForesterUtil\r
2207                                 .wordWrap( "phyloXML XSD-based validation is turned off [enable with line 'validate_against_phyloxml_xsd_schem: true' in configuration file]",\r
2208                                            80 ),\r
2209                                            "Warning",\r
2210                                            JOptionPane.WARNING_MESSAGE );\r
2211         }\r
2212     }\r
2213 } // MainFrameApplication.\r