X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Farchaeopteryx%2Ftools%2FPhylogeneticInferrer.java;h=6cf66148e526bc0f97ad1ff2f434e196c3b5eb7d;hb=bd30826b30945fec52ffd33f65f9456b6ef57e77;hp=981e3c80e22a11607c65418878ecdd6b7a0974a7;hpb=505270dd6bd8892fe00658607de06e5f030b11db;p=jalview.git diff --git a/forester/java/src/org/forester/archaeopteryx/tools/PhylogeneticInferrer.java b/forester/java/src/org/forester/archaeopteryx/tools/PhylogeneticInferrer.java index 981e3c8..6cf6614 100644 --- a/forester/java/src/org/forester/archaeopteryx/tools/PhylogeneticInferrer.java +++ b/forester/java/src/org/forester/archaeopteryx/tools/PhylogeneticInferrer.java @@ -58,12 +58,13 @@ import org.forester.sequence.Sequence; import org.forester.tools.ConfidenceAssessor; import org.forester.util.ForesterUtil; -public class PhylogeneticInferrer implements Runnable { +public class PhylogeneticInferrer extends RunnableProcess { private Msa _msa; private final MainFrameApplication _mf; private final PhylogeneticInferenceOptions _options; private final List _seqs; + private final boolean DEBUG = true; public final static String MSA_FILE_SUFFIX = ".aln"; public final static String PWD_FILE_SUFFIX = ".pwd"; @@ -87,9 +88,11 @@ public class PhylogeneticInferrer implements Runnable { private Msa inferMsa() throws IOException, InterruptedException { final File temp_seqs_file = File.createTempFile( "__msa__temp__", ".fasta" ); - System.out.println(); - System.out.println( "temp file: " + temp_seqs_file ); - System.out.println(); + if ( DEBUG ) { + System.out.println(); + System.out.println( "temp file: " + temp_seqs_file ); + System.out.println(); + } //final File temp_seqs_file = new File( _options.getTempDir() + ForesterUtil.FILE_SEPARATOR + "s.fasta" ); final BufferedWriter writer = new BufferedWriter( new FileWriter( temp_seqs_file ) ); SequenceWriter.writeSeqs( _seqs, writer, SEQ_FORMAT.FASTA, 100 ); @@ -142,7 +145,7 @@ public class PhylogeneticInferrer implements Runnable { e.printStackTrace(); } } - final NeighborJoining nj = new NeighborJoining(); + final NeighborJoining nj = NeighborJoining.createInstance(); final Phylogeny phy = nj.execute( m ); PhylogeneticInferrer.extractFastaInformation( phy ); return phy; @@ -153,22 +156,42 @@ public class PhylogeneticInferrer implements Runnable { if ( ( _msa == null ) && ( _seqs == null ) ) { throw new IllegalArgumentException( "cannot run phylogenetic analysis with null msa and seq array" ); } + start( _mf, "phylogenetic inference" ); if ( _msa == null ) { Msa msa = null; try { msa = inferMsa(); } catch ( final IOException e ) { + end( _mf ); JOptionPane.showMessageDialog( _mf, - "Could not create multiple sequence alignment with " - + _options.getMsaPrg() + "\nand the following parameters:\n\"" - + _options.getMsaPrgParameters() + "\"\nError:" + "Could not create multiple sequence alignment with \"" + + _options.getMsaPrg() + "\" and the following parameters:\n\"" + + _options.getMsaPrgParameters() + "\"\nError: " + e.getLocalizedMessage(), "Failed to Calculate MSA", JOptionPane.ERROR_MESSAGE ); + if ( DEBUG ) { + e.printStackTrace(); + } + return; + } + catch ( final Exception e ) { + end( _mf ); + JOptionPane.showMessageDialog( _mf, + "Could not create multiple sequence alignment with \"" + + _options.getMsaPrg() + "\" and the following parameters:\n\"" + + _options.getMsaPrgParameters() + "\"\nError: " + + e.getLocalizedMessage(), + "Unexpected Exception During MSA Calculation", + JOptionPane.ERROR_MESSAGE ); + if ( DEBUG ) { + e.printStackTrace(); + } return; } if ( msa == null ) { + end( _mf ); JOptionPane.showMessageDialog( _mf, "Could not create multiple sequence alignment with " + _options.getMsaPrg() + "\nand the following parameters:\n\"" @@ -177,14 +200,17 @@ public class PhylogeneticInferrer implements Runnable { JOptionPane.ERROR_MESSAGE ); return; } - System.out.println( msa.toString() ); - System.out.println( MsaMethods.calcBasicGapinessStatistics( msa ).toString() ); + if ( DEBUG ) { + System.out.println( msa.toString() ); + System.out.println( MsaMethods.calcBasicGapinessStatistics( msa ).toString() ); + } final MsaMethods msa_tools = MsaMethods.createInstance(); if ( _options.isExecuteMsaProcessing() ) { msa = msa_tools.removeGapColumns( _options.getMsaProcessingMaxAllowedGapRatio(), _options.getMsaProcessingMinAllowedLength(), msa ); if ( msa == null ) { + end( _mf ); JOptionPane.showMessageDialog( _mf, "Less than two sequences longer than " + _options.getMsaProcessingMinAllowedLength() @@ -194,9 +220,11 @@ public class PhylogeneticInferrer implements Runnable { return; } } - System.out.println( msa_tools.getIgnoredSequenceIds() ); - System.out.println( msa.toString() ); - System.out.println( MsaMethods.calcBasicGapinessStatistics( msa ).toString() ); + if ( DEBUG ) { + System.out.println( msa_tools.getIgnoredSequenceIds() ); + System.out.println( msa.toString() ); + System.out.println( MsaMethods.calcBasicGapinessStatistics( msa ).toString() ); + } _msa = msa; } final int n = _options.getBootstrapSamples(); @@ -214,7 +242,8 @@ public class PhylogeneticInferrer implements Runnable { ConfidenceAssessor.evaluate( "bootstrap", eval_phys, master_phy, true, 1 ); } _mf.getMainPanel().addPhylogenyInNewTab( master_phy, _mf.getConfiguration(), "nj", "njpath" ); - _mf.getMainPanel().getCurrentTreePanel().setArrowCursor(); + // _mf.getMainPanel().getCurrentTreePanel().setArrowCursor(); + end( _mf ); JOptionPane.showMessageDialog( _mf, "Inference successfully completed", "Inference Completed",